load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl" diri = "/Users/shea/Data/CDC/" fils = systemfunc("cd "+diri+" ; ls air.day.20*.nc") print(fils) pthi = diri+fils f = addfiles(pthi, "r") ; default is 'cat' air = short2flt( f[:]->air ) ; if 'air' is type short air = f[:]->air ; if 'air' is type float or double printVarSummary (air) newlat = fspan(0.00,-3.00,3) newlon = fspan(6.00,53.00,47) newlat@units = "degrees_south" newlon@units = "degrees_east" regrid = linint2_Wrap(air&lon,air&lat,air,True,newlon,newlat,0) printVarSummary(regrid) pnew = (/ 1000., 950., 850., 750., 650., 550., 450., 350., 250., 150., 50./) pnew!0 = "pnew" pnew@units = "mb" linlog = 2 air03 = int2p_n_Wrap(regrid&level,regrid,pnew,linlog,1) printVarSummary(air03) system("/bin/rm -f Air03.nc") fout=addfile("Air03.nc","c") fout@title="Air03" fout->AIR03 = air03