begin
yrStrt = 1980 yrLast = 2011 season = "JFM" ; choose Dec-Jan-Feb seasonal mean
f = addfile ("/disk_raid5/data/NCEP2-2.5/monthly/pressure/uwnd.mon.mean.nc", "r") TIME = f->time YYYY = cd_calendar(TIME,-1)/100 ; entire file iYYYY = ind(YYYY.ge.yrStrt .and. YYYY.le.yrLast) uwnd = short2flt(f->uwnd(11:394,7,0:36,:)) ;
UWND = month_to_season (uwnd, season) nyrs = dimsizes(UWND&time)
time = TIME(11:394) lat=f->lat(0:36) lon = f->lon level = f->level ; calculate climatology UWND_cli = dim_avg_n_Wrap(UWND, 0)
ff = addfile ("/disk_raid5/data/NCEP2-2.5/monthly/pressure/vwnd.mon.mean.nc", "r") TIME2 = ff->time YYYY2 = cd_calendar(TIME2,-1)/100 ; entire file iYYYY2 = ind(YYYY2.ge.yrStrt .and. YYYY2.le.yrLast) vwnd = short2flt(ff->vwnd(11:394,7,0:36,:)) ;VWND = month_to_season (vwnd, season) VWND = month_to_season (vwnd, season) nyrs = dimsizes(VWND&time) time = TIME(11:394) lat=ff->lat(0:36) lon = ff->lon level = ff->level ; calculate climatology VWND_cli = dim_avg_n_Wrap(VWND, 0)
wks = gsn_open_wks("ps","UV300") gsn_define_colormap(wks,"temp_diff_18lev") ; choose colormap
|