import xarray as xr
import sys
import os

ifile = sys.argv[1]

ds = xr.open_dataset(ifile)
ds = ds.sel(x=ds.x[26:90], y=ds.y[154:204])
ofile = '/home/disk/rocinante/DATA/temp/TNC_stormwater/pnnl/{}'.format(os.path.basename(ifile))
ds.to_netcdf(ofile)