load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" begin FireFile=asciiread("dlyFireClim.txt",(/214,12/),"float") print("read fire") LiteFile=asciiread("dlyLightClim.txt",(/214,12/),"float") print("read light") PptFile=asciiread("dlyPptClim.txt",(/214,12/),"float") print("read light") fire = new((/9,214/),"float") light = new((/9,214/),"float") ppt = new((/9,214/),"float") do nyr = 1990, 1998 fire(nyr-1990,:) = FireFile(:,nyr-1987) light(nyr-1990,:) = LiteFile(:,nyr-1987) ppt(nyr-1990,:) = PptFile(:,nyr-1987) end do day = PptFile(:,0) ;----------------------------------------------------------------------- wks_type = "x11" wks = gsn_open_wks (wks_type,"climyrly") plot = new(2,graphic) ;------------------------------------------------------------------------ ; ;BOX PLOT RESOURCES resL = True resL@trXMinF = 91 resL@trXMaxF = 304 resL@gsnDraw = False resL@gsnFrame = False resL@tmYROn = False ;Still shows resL@tmYLOn = False resL@tmXBOn = False resL@tmXTOn = False do nyr = 0, 1 dataL = new((/3,214/),float) dataL(0,:) = light(nyr,:)/stddev(light(nyr,:)) dataL(1,:) = fire(nyr,:)/stddev(fire(nyr,:)) dataL(2,:) = ppt(nyr,:) /stddev(ppt(nyr,:)) plot(nyr) = gsn_csm_xy(wks,day,dataL,resL) end do resPan = True resPan@gsnFrame = False resPan@gsnMaximize = True gsn_panel(wks,plot,(/2,1/),resPan) frame(wks) end