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"""