load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" begin FireFile=asciiread("/fire4/bhall/DryTstorm/dlyFireClim.txt",(/214,12/),"float") print("read fire") LiteFile=asciiread("/fire4/bhall/DryTstorm/dlyLightClim.txt",(/214,12/),"float") print("read light") PptFile=asciiread("/fire4/bhall/DryTstorm/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@trYMinF = 0 resL@trYMaxF = 50000 resL@gsnDraw = False resL@gsnFrame = False resL@gsnMaximize = True resL@tmYROn = False ;Still shows resL@tmYLOn = False resL@tmXBOn = False resL@tmXTOn = False resR = True resR@tmYROn = False ;Still shows resL@tmYROn = False ;Still shows do nyr = 0, 1 dataL = new(214,float) dataL = light(nyr,:) dataR = new((/2,214/),float) dataR(0,:) = fire(nyr,:) dataR(1,:) = ppt(nyr,:) / 5 plot(nyr) = gsn_csm_xy2(wks,day,dataL,dataR,resL,resR) end do resPan = True resPan@gsnFrame = False gsn_panel(wks,plot,(/2,1/),resPan) frame(wks) end