import xarray as xr
import sys

gcm = sys.argv[1]
wrf = xr.open_dataset('{}_PREC_interp.nc'.format(gcm))
#wrf = wrf.sel(time=wrf.time[skips:])
zz = wrf.sel(time=wrf.time[0:1])
zz['time'] = zz.time - 3600000000000
zz['PREC'] = zz.PREC * 0
ds = xr.concat([zz, wrf], dim='time')
ds.to_netcdf('merge/{}_PREC_interp.nc'.format(gcm))
wrfd = wrf.resample(time='D').sum()
wrfd.to_netcdf('merge/{}_PREC_interp_daily.nc'.format(gcm))