load "/home/cis/nslib/essential_libs.ncl" load "/home/cis/nslib/jeenie.module" begin Data_dir="/apcc01/CIS/data/IDX/OBS" outdir="/apcc01/CIS/figure/indx/" indx="dmi" infil1=Data_dir+"/"+indx+".nc" fin1=addfile(infil1,"r") obs1=fin1->$indx$ obs1=obs1/stddev(obs1) time=obs1&time res=True opt=True wks=gsn_open_wks("x11","timeseries_"+indx) res@gsnYRefLine = 0 res@gsnAboveYRefLineColor="red" res@gsnBelowYRefLineColor="blue" res@xyLineThicknesses = (/0.1/) res@tiYAxisString="Standardized value" time_int=5 time_axis_in_years(res,time,time_int) res@tiMainString = changeCase(indx,"up")+" index " plot1=gsn_csm_xy(wks,time,obs1,res) wks=gsn_open_wks("x11","timeseries_2yr_"+indx) ; Second plot draws for last two years ntim=dimsizes(time) tim2=time(ntim-23:) obs2=obs1(ntim-23:) ntim2=dimsizes(tim2) time_int=3 time_axis_in_years_months(res,tim2,time_int) plot2=gsn_csm_xy(wks,tim2,obs2,res) end