data=old_var[0,:,new_y,new_x]
if var_name in ['W','TABS','QV','QN']:
npz_vars[var_name]=data
print data.shape
out_var[:,:,:]=data[...]
for the_att in old_var.ncattrs():
setattr(out_var,the_att,getattr(old_var,the_att).strip())
for var_name in ['x','y','z']:
old_var=one_d[var_name]
for the_att in old_var.ncattrs():
new_var=out_vars[var_name]
setattr(new_var,the_att,getattr(old_var,the_att))
nc_out.description='subset of BOMEX run created with make_subset.py'
nc_out.in_file=in_file
nc_in.close()
nc_out.close()
var_descrip={'W':{'long_name':'vertical velocity','units':'m/s'},
'TABS':{'long_name':'temperature','units':'K'},
'p':{'long_name':'pressure','units':'hPa'},
'QV':{'long_name':'vapor mixing ratio','units':'g/kg'},
'z':{'long_name':'height','units':'m'},
'x':{'long_name':'xcoord','units':'m'},
'y':{'long_name':'ycoord','units':'m'},
}
var__file="""output from SAM LES run BOMEX_256x256x150_25m_20m_2s_16_0000009240.nc
written by make_sub_set.py on 2012/Jan/30"""
