import pandas as pd
import sys
import os
import numpy as np

args = sys.argv
ifile = args[1]
fn = os.path.basename(ifile)
opath = "/home/disk/rocinante/DATA/temp/WRF/scripts/regrid_c2c/longwave/output/chehalis_delta/"
#deltaf = "/home/disk/tsuga2/jswon11/workdir/2021_09_KingCounty-Phase3/scripts/lelv/delta_livneh_15km/low_elev_delta.csv"
deltaf = "/home/disk/tsuga2/jswon11/workdir/2021_09_KingCounty-Phase3/scripts/lelv/delta_che/low_elev_delta.csv"
dtfmt = "%m/%d/%Y-%H:%M:%S"

df = pd.read_csv(ifile, header=None, sep=' ')
df.columns = ['date', 'temp', 'wind', 'rh', 'sw', 'lw', 'prec']
df.index = pd.to_datetime(df.date, format=dtfmt)
df.index = df.index.shift(-8, freq='H')
df.date = df.index.strftime(date_format=dtfmt)
dlt = pd.read_csv(deltaf)
mons = df.index.month.values
tmin = np.array([dlt.Tmin[x-1] for x in mons])
tmax = np.array([dlt.Tmax[x-1] for x in mons])
prec = np.array([dlt.Prec[x-1] for x in mons])
hrs = np.array(df.index.hour.values)
#delta = (tmax+tmin)/2 + (tmax-tmin)/2*np.sin((2*np.pi/24*(hrs+14)))
delta = tmin
df.temp += delta
df.prec *= prec


opath = "{}/{}".format(opath, fn)
df.to_csv(opath, header=False, index=False, sep=' ')