case = "wheeler" diro = "./" latN = 15 latS = -latN nDayWin = 96 nDaySkip =1 opt = True opt@debug= True opt@pltType = "png" opt@cnLinesOn = False ; turn off contour lines spd = 1 var = "olr" setfileoption("nc", "FileStructure", "Advanced") ; not need for newer versions of NCL ; The 1675 index used to start the calculations was determined seperately. ; index 1675 corresponds to day 19790101 [ 1 Jan 1979 ] f = addfile("olr.day.mean.nc", "r") olr = short2flt( f->olr(1675:,{latS:latN},:) ) ; NOAA ;olr = f->olr(1675:,{latS:latN},:)) ; Adit Modi printVarSummary(olr) printMinMax(olr,0) print("-----") ymd = cd_calendar(olr&time, -2) ; time -> yyyymmdd dimo = dimsizes(olr) ntim = dimo(0) ; # days on file print("First and Last dates: ymd(0)="+ymd(0)+" ymd(ntim-1)="+ymd(ntim-1)) wkSpaceTime (olr, diro, case, var,latN, spd, nDayWin, nDaySkip, opt )