[ncl-talk] Problem with panel plot (size difference)
Adam Phillips
asphilli at ucar.edu
Mon Dec 11 10:14:57 MST 2017
Hi Tomoko,
NCL uses the space needed to draw the first panel to allocate the required
space for each of the following panels. Thus, if your last column has
bigger plots (and labelbars are considered part of the plot in this case),
NCL will shrink the plot to fit into the allocated space. You have two
choices:
1) Tell NCL to allocate space by looking at a different panel by setting
gsnPanelScalePlotIndex = 1 (or 2). See example #24 here:
https://www.ncl.ucar.edu/Applications/panel.shtml#ex24
2) You can explicitly set the sizes and locations of each panel by setting
gsnPanelXF and gsnPanelYF, using information returned from gsnPanelDebug.
See example 19 here:
https://www.ncl.ucar.edu/Applications/panel.shtml#ex19
Hope that helps! If not, or if you have further questions please respond to
the ncl-talk email list.
Adam
On Fri, Dec 8, 2017 at 10:54 PM, Tomoko Koyama <Tomoko.Koyama at colorado.edu>
wrote:
> When I run the following script to make a panel, figures of the rightmost
> appear smaller.
> I guess that’s because I only add legends to the rightmost, but is there
> any way to have all the figures in the same size?
>
> Thank you in advance,
> Tomoko
>
> ;-------------------------------------------------------------------------
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
>
> ;;; Constant information
>
> model = "CCSM4"
>
> weeks = (/2,4,6,8,10,12/)
> nwk = dimsizes(weeks)
>
> ; SOM size
> xdim = my_xdim
> ydim = my_ydim
>
> ; Which node to process
> xcomp = my_comp_x
> ycomp = my_comp_y
> som_node = xcomp+(ycomp*xdim)
>
> ; Directories
> indir = "/work/koyama_sice/composites/"
> outdir = "/work/koyama_sice/img_panel/"
>
> ; Arrays for plotting
> ext_var1 = new((/nwk,180,180/), "float")
> ext_var2 = new((/nwk,180,180/), "float")
> ext_var3 = new((/nwk,180,180/), "float")
> ext_var4 = new((/nwk,180,180/), "float")
> ext_var5 = new((/nwk,180,180/), "float")
> ext_var6 = new((/nwk,180,180/), "float")
> ext_var7 = new((/nwk,180,180/), "float")
> ext_var8 = new((/nwk,180,180/), "float")
>
> ss_var1 = new((/nwk,180,180/), "float")
> ss_var2 = new((/nwk,180,180/), "float")
> ss_var3 = new((/nwk,180,180/), "float")
> ss_var4 = new((/nwk,180,180/), "float")
> ss_var5 = new((/nwk,180,180/), "float")
> ss_var6 = new((/nwk,180,180/), "float")
> ss_var7 = new((/nwk,180,180/), "float")
> ss_var8 = new((/nwk,180,180/), "float")
>
> ;-------------------------------------------------------------------------
> ; Get information for what was masked out
> dirlsm="/work/koyama_sice/Land_Mask/"
> fnlsm=addfile(dirlsm+"ERAI.LandMask.180x180.org.nc","r")
> lsm=fnlsm->Land_Mask
>
> ;-------------------------------------------------------------------------
> ; Grab near-surface temperature data
> fili = systemfunc("cd " + indir + " ; ls CCSM4.tas_composite_my_
> xdimxmy_ydim_wk*.nc")
> fils = systemfunc("cd " + indir + " ; ls CCSM4.tas_statsig_my_xdimxmy_
> ydim_wk*.nc")
>
> do iwk=0,nwk-1
> fn = addfile(indir+fili(weeks(iwk)),"r")
> var = fn->tas
> lat = fn->lat2d
> lon = fn->lon2d
> ext_var1(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> do iwk=0,nwk-1
> fn = addfile(indir+fils(weeks(iwk)),"r")
> var = fn->alpha
> ss_var1(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> ;-------------------------------------------------------------------------
> ; Grab sea level pressure data
> fili = systemfunc("cd " + indir + " ; ls CCSM4.psl_composite_my_
> xdimxmy_ydim_wk*.nc")
> fils = systemfunc("cd " + indir + " ; ls CCSM4.psl_statsig_my_xdimxmy_
> ydim_wk*.nc")
>
> do iwk=0,nwk-1
> fn = addfile(indir+fili(weeks(iwk)),"r")
> var = fn->psl
> ext_var2(iwk,:,:) = var(som_node,:,:)/100.
> delete(var)
> end do
>
> do iwk=0,nwk-1
> fn = addfile(indir+fils(weeks(iwk)),"r")
> var = fn->alpha
> ss_var2(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> ;-------------------------------------------------------------------------
> ; Grab 1000 - 500 hPa thickness data
> fili = systemfunc("cd " + indir + " ; ls CCSM4.thickness_composite_my_
> xdimxmy_ydim_wk*.nc")
> fils = systemfunc("cd " + indir + " ; ls CCSM4.thickness_statsig_my_
> xdimxmy_ydim_wk*.nc")
>
> do iwk=0,nwk-1
> fn = addfile(indir+fili(weeks(iwk)),"r")
> var = fn->thickness
> ext_var3(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> do iwk=0,nwk-1
> fn = addfile(indir+fils(weeks(iwk)),"r")
> var = fn->alpha
> ss_var3(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> ;-------------------------------------------------------------------------
> ; Grab sea ice concentration data
> fili = systemfunc("cd " + indir + " ; ls CCSM4.sic_composite_my_
> xdimxmy_ydim_wk*.nc")
> fils = systemfunc("cd " + indir + " ; ls CCSM4.sic_statsig_my_xdimxmy_
> ydim_wk*.nc")
>
> do iwk=0,nwk-1
> fn = addfile(indir+fili(weeks(iwk)),"r")
> var = fn->sic
> ext_var4(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> do iwk=0,nwk-1
> fn = addfile(indir+fils(weeks(iwk)),"r")
> var = fn->alpha
> ss_var4(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> ;-------------------------------------------------------------------------
> ; Grab sea ice thickness data
> fili = systemfunc("cd " + indir + " ; ls CCSM4.sit_composite_my_
> xdimxmy_ydim_wk*.nc")
> fils = systemfunc("cd " + indir + " ; ls CCSM4.sit_statsig_my_xdimxmy_
> ydim_wk*.nc")
>
> do iwk=0,nwk-1
> fn = addfile(indir+fili(weeks(iwk)),"r")
> var = fn->sit
> ext_var5(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> do iwk=0,nwk-1
> fn = addfile(indir+fils(weeks(iwk)),"r")
> var = fn->alpha
> ss_var5(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> ;-------------------------------------------------------------------------
> ; Grab total surface energy flux
> fili = systemfunc("cd " + indir + " ; ls CCSM4.sfcnf_composite_my_
> xdimxmy_ydim_wk*.nc")
> fils = systemfunc("cd " + indir + " ; ls CCSM4.sfcnf_statsig_my_
> xdimxmy_ydim_wk*.nc")
>
> do iwk=0,nwk-1
> fn = addfile(indir+fili(weeks(iwk)),"r")
> var = fn->sfcnf
> ext_var6(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> do iwk=0,nwk-1
> fn = addfile(indir+fils(weeks(iwk)),"r")
> var = fn->alpha
> ss_var6(iwk,:,:) = var(som_node,:,:)
> delete(var)
> end do
>
> ;-------------------------------------------------------------------------
> ; Grab geopotential height
> fili = systemfunc("cd " + indir + " ; ls CCSM4.zg_composite_my_xdimxmy_
> ydim_wk*.nc")
> fils = systemfunc("cd " + indir + " ; ls CCSM4.zg_statsig_my_xdimxmy_
> ydim_wk*.nc")
>
> do iwk=0,nwk-1
> fn = addfile(indir+fili(weeks(iwk)),"r")
> var = fn->zg
> ext_var7(iwk,:,:) = var(som_node,3,:,:) ; 500 hPa
> ext_var8(iwk,:,:) = var(som_node,4,:,:) ; 250 hPa
> delete(var)
> end do
>
> do iwk=0,nwk-1
> fn = addfile(indir+fils(weeks(iwk)),"r")
> var = fn->alpha
> ss_var7(iwk,:,:) = var(som_node,3,:,:)
> ss_var8(iwk,:,:) = var(som_node,4,:,:)
> delete(var)
> end do
>
> ;-------------------------------------------------------------------------
> ;;; Make plot
>
> fnout="CCSM4_seq_my_xdimxmy_ydim_my_comp_x_my_comp_y.pdf"
> wks = gsn_open_wks("pdf", outdir+fnout)
> gsn_define_colormap(wks,"BlRe")
> additional_colors = (/"gray85"/)
> gsn_merge_colormaps(wks,"BlRe",additional_colors)
>
> ; Map resources
> res = True
> res at mpLandFillColor = "gray85"
> res at mpFillOn = True
> res at mpGeophysicalLineColor = "black" ; color of continents
> res at mpGeophysicalLineThicknessF = 2.
> res at gsnPolar = "NH" ; northern hemisphere
> plot
> res at mpMinLatF = 30 ; to focus in on NP
>
> ; Contour info
> res at cnFillOn = True
> res at cnFillMode = "CellFill"
> res at cnLineLabelsOn = False
> res at gsnSpreadColors = True ; use total colormap
> res at cnLevelSelectionMode = "ManualLevels"
> res at gsnSpreadColorStart = 35
> res at gsnSpreadColorEnd = -35
> res at gsnAddCyclic = False
> res at cnInfoLabelOn = False ; no contour info label
>
> res at lbLabelBarOn = False ; turn off individual cb's
> res at lbOrientation = "Vertical"
> res at lbTitleAngleF = 90.
> res at lbTitlePosition = "Right"
> res at lbTitleDirection = "Across"
> ;res at lbLabelFontHeightF = 0.025
> res at gsnRightString = " "
> res at gsnLeftString = " "
>
> res at trGridType = "TriangularMesh"
> res at sfXArray = lon
> res at sfYArray = lat
> res at gsnMaximize = True ; max plot area
> res at gsnFrame = False ; do not advance frame
> res at gsnDraw = False
> res at tiXAxisString = ""
> res at tiYAxisString = ""
> res at gsnRightStringFontHeightF = 0.025
> res at gsnLeftStringFontHeightF = 0.025
>
> ; Overplot
> res_ss = True
> res_ss at gsnDraw = False
> res_ss at gsnFrame = False
> res_ss at cnFillOn = True
> res_ss at cnFillMode = "CellFill"
> res_ss at cnFillOpacityF = 0.2
> res_ss at cnLinesOn = False
> res_ss at gsnTickMarksOn = False
> res_ss at gsnAddCyclic = False
> res_ss at cnLineThicknessF = 0.
> res_ss at cnLineLabelsOn = False
> res_ss at cnLineColor = "black"
> res_ss at cnInfoLabelOn = False
> res_ss at lbLabelBarOn = False
> res_ss at cnLevelSelectionMode = "ExplicitLevels"
> res_ss at cnLevels = 89.9
> res_ss at cnFillColors = (/"white","black"/)
> res_ss at cnMonoLineDashPattern = False
> res_ss at cnLineDashPattern = 0
> res_ss at gsnMaximize = True
> res_ss at trGridType = "TriangularMesh"
> res_ss at sfXArray = lon
> res_ss at sfYArray = lat
>
> opt = True
> opt at gsnShadeFillType = "pattern"
> opt at gsnShadeHigh = 17
>
> maxLev = 18
>
> plot = new(48,graphic)
>
> ;-------------------------------------------------------------------------
> ; 1. tas
> res1 = res
> maximum=2.4
> minimum=-1.*maximum
> mnmxint = nice_mnmxintvl(minimum,maximum,maxLev,True)
> res1 at cnLinesOn = False
> res1 at cnMinLevelValF = mnmxint(0)
> res1 at cnMaxLevelValF = mnmxint(1)
> res1 at cnLevelSpacingF = mnmxint(2)
> res1 at lbTitleString = "Near-surface temperature"
>
> ; 2. psl
> res2 = res
> maximum=4.
> minimum=-1.*maximum
> mnmxint = nice_mnmxintvl(minimum,maximum,maxLev,True)
> res2 at cnLinesOn = True
> res2 at cnMinLevelValF = mnmxint(0)
> res2 at cnMaxLevelValF = mnmxint(1)
> res2 at cnLevelSpacingF = mnmxint(2)
> res2 at lbTitleString = "Sea level pressure"
>
> ; 3. 1000 - 500 hPa thickness
> res3 = res
> maximum=28.
> minimum=-1.*maximum
> mnmxint = nice_mnmxintvl(minimum,maximum,maxLev,True)
> res3 at cnLinesOn = False
> res3 at cnMinLevelValF = mnmxint(0)
> res3 at cnMaxLevelValF = mnmxint(1)
> res3 at cnLevelSpacingF = mnmxint(2)
> res3 at lbTitleString = "1000-500hPa thickness"
>
> ; 4. sic
> ; Mask land areas
> ext_var4 = mask(ext_var4,lsm.ge.0.5,False)
>
> res4 = res
> maximum=12.
> minimum=-1.*maximum
> mnmxint = nice_mnmxintvl(minimum,maximum,maxLev,True)
> cmapt=read_colormap_file("BlRe")
> cmap=cmapt(::-1,:)
> res4 at gsnSpreadColors = False ; use total colormap
> res4 at cnSpanFillPalette = True
> res4 at cnFillPalette = cmap(33:63,:)
> res4 at cnLinesOn = False
> res4 at cnMinLevelValF = mnmxint(0)
> res4 at cnMaxLevelValF = mnmxint(1)
> res4 at cnLevelSpacingF = mnmxint(2)
> res4 at lbTitleString = "Sea ice concentration"
>
> ; 5. sit
> ; Mask land areas
> ext_var5 = mask(ext_var5,lsm.ge.0.5,False)
>
> res5 = res
> maximum=0.24
> minimum=-1.*maximum
> mnmxint = nice_mnmxintvl(minimum,maximum,maxLev,True)
> res5 at gsnSpreadColors = False ; use total colormap
> res5 at cnSpanFillPalette = True
> res5 at cnFillPalette = cmap(33:63,:)
> res5 at cnLinesOn = False
> res5 at cnMinLevelValF = mnmxint(0)
> res5 at cnMaxLevelValF = mnmxint(1)
> res5 at cnLevelSpacingF = mnmxint(2)
> res5 at lbTitleString = "Sea ice thickness"
>
> ; 6. sfcnf
> res6 = res
> maximum=28.
> minimum=-1.*maximum
> mnmxint = nice_mnmxintvl(minimum,maximum,maxLev,True)
> res6 at cnLinesOn = False
> res6 at cnMinLevelValF = mnmxint(0)
> res6 at cnMaxLevelValF = mnmxint(1)
> res6 at cnLevelSpacingF = mnmxint(2)
> res6 at lbTitleString = "Total surface energy flux"
>
> ; 7. zg500
> res7 = res
> maximum=40.
> minimum=-1.*maximum
> mnmxint = nice_mnmxintvl(minimum,maximum,maxLev,True)
> res7 at cnLinesOn = False
> res7 at cnMinLevelValF = mnmxint(0)
> res7 at cnMaxLevelValF = mnmxint(1)
> res7 at cnLevelSpacingF = mnmxint(2)
> res7 at lbTitleString = "Geopotential height @ 500 hPa"
>
> ; 8. zg250
> res8 = res
> maximum=60.
> minimum=-1.*maximum
> mnmxint = nice_mnmxintvl(minimum,maximum,maxLev,True)
> res8 at cnLinesOn = False
> res8 at cnMinLevelValF = mnmxint(0)
> res8 at cnMaxLevelValF = mnmxint(1)
> res8 at cnLevelSpacingF = mnmxint(2)
> res8 at lbTitleString = "Geopotential height @ 250 hPa"
>
> res11 = res1
> res12 = res1
> res13 = res1
> res14 = res1
> res15 = res1
> res16 = res1
> res21 = res2
> res22 = res2
> res23 = res2
> res24 = res2
> res25 = res2
> res26 = res2
> res31 = res3
> res32 = res3
> res33 = res3
> res34 = res3
> res35 = res3
> res36 = res3
> res41 = res4
> res42 = res4
> res43 = res4
> res44 = res4
> res45 = res4
> res46 = res4
> res51 = res5
> res52 = res5
> res53 = res5
> res54 = res5
> res55 = res5
> res56 = res5
> res61 = res6
> res62 = res6
> res63 = res6
> res64 = res6
> res65 = res6
> res66 = res6
> res71 = res7
> res72 = res7
> res73 = res7
> res74 = res7
> res75 = res7
> res76 = res7
> res81 = res8
> res82 = res8
> res83 = res8
> res84 = res8
> res85 = res8
> res86 = res8
>
> ;---1st row plots
> res11 at gsnLeftString = "week 2"
> res12 at gsnLeftString = "week 4"
> res13 at gsnLeftString = "week 6"
> res14 at gsnLeftString = "week 8"
> res15 at gsnLeftString = "week 10"
> res16 at gsnLeftString = "week 12"
>
> ;---Rightmost plots
> res16 at gsnRightString = "K"
> res26 at gsnRightString = "hPa"
> res36 at gsnRightString = "m"
> res46 at gsnRightString = "%"
> res56 at gsnRightString = "m"
> res66 at gsnRightString = "W m~S~-2~N~"
> res76 at gsnRightString = "m"
> res86 at gsnRightString = "m"
>
> res16 at lbLabelBarOn = True
> res26 at lbLabelBarOn = True
> res36 at lbLabelBarOn = True
> res46 at lbLabelBarOn = True
> res56 at lbLabelBarOn = True
> res66 at lbLabelBarOn = True
> res76 at lbLabelBarOn = True
> res86 at lbLabelBarOn = True
>
> ; 1st row
> ; 1. tas
> plot(0) = gsn_csm_contour_map_polar(wks,ext_var1(0,:,:),res11)
> overlay_plot = gsn_csm_contour(wks,ss_var1(0,:,:),res_ss)
> overlay(plot(0),overlay_plot)
> plot(1) = gsn_csm_contour_map_polar(wks,ext_var1(1,:,:),res12)
> overlay_plot = gsn_csm_contour(wks,ss_var1(1,:,:),res_ss)
> overlay(plot(1),overlay_plot)
> plot(2) = gsn_csm_contour_map_polar(wks,ext_var1(2,:,:),res13)
> overlay_plot = gsn_csm_contour(wks,ss_var1(2,:,:),res_ss)
> overlay(plot(2),overlay_plot)
> plot(3) = gsn_csm_contour_map_polar(wks,ext_var1(3,:,:),res14)
> overlay_plot = gsn_csm_contour(wks,ss_var1(3,:,:),res_ss)
> overlay(plot(3),overlay_plot)
> plot(4) = gsn_csm_contour_map_polar(wks,ext_var1(4,:,:),res15)
> overlay_plot = gsn_csm_contour(wks,ss_var1(4,:,:),res_ss)
> overlay(plot(4),overlay_plot)
> plot(5) = gsn_csm_contour_map_polar(wks,ext_var1(5,:,:),res16)
> overlay_plot = gsn_csm_contour(wks,ss_var1(5,:,:),res_ss)
> overlay(plot(5),overlay_plot)
>
> ; 2nd row
> plot(6) = gsn_csm_contour_map_polar(wks,ext_var2(0,:,:),res21)
> overlay_plot = gsn_csm_contour(wks,ss_var2(0,:,:),res_ss)
> overlay(plot(6),overlay_plot)
> plot(7) = gsn_csm_contour_map_polar(wks,ext_var2(1,:,:),res22)
> overlay_plot = gsn_csm_contour(wks,ss_var2(1,:,:),res_ss)
> overlay(plot(7),overlay_plot)
> plot(8) = gsn_csm_contour_map_polar(wks,ext_var2(2,:,:),res23)
> overlay_plot = gsn_csm_contour(wks,ss_var2(2,:,:),res_ss)
> overlay(plot(8),overlay_plot)
> plot(9) = gsn_csm_contour_map_polar(wks,ext_var2(3,:,:),res24)
> overlay_plot = gsn_csm_contour(wks,ss_var2(3,:,:),res_ss)
> overlay(plot(9),overlay_plot)
> plot(10) = gsn_csm_contour_map_polar(wks,ext_var2(4,:,:),res25)
> overlay_plot = gsn_csm_contour(wks,ss_var2(4,:,:),res_ss)
> overlay(plot(10),overlay_plot)
> plot(11) = gsn_csm_contour_map_polar(wks,ext_var2(5,:,:),res26)
> overlay_plot = gsn_csm_contour(wks,ss_var2(5,:,:),res_ss)
> overlay(plot(11),overlay_plot)
>
> ; 3rd row
> plot(12) = gsn_csm_contour_map_polar(wks,ext_var3(0,:,:),res31)
> overlay_plot = gsn_csm_contour(wks,ss_var3(0,:,:),res_ss)
> overlay(plot(12),overlay_plot)
> plot(13) = gsn_csm_contour_map_polar(wks,ext_var3(1,:,:),res32)
> overlay_plot = gsn_csm_contour(wks,ss_var3(1,:,:),res_ss)
> overlay(plot(13),overlay_plot)
> plot(14) = gsn_csm_contour_map_polar(wks,ext_var3(2,:,:),res33)
> overlay_plot = gsn_csm_contour(wks,ss_var3(2,:,:),res_ss)
> overlay(plot(14),overlay_plot)
> plot(15) = gsn_csm_contour_map_polar(wks,ext_var3(3,:,:),res34)
> overlay_plot = gsn_csm_contour(wks,ss_var3(3,:,:),res_ss)
> overlay(plot(15),overlay_plot)
> plot(16) = gsn_csm_contour_map_polar(wks,ext_var3(4,:,:),res35)
> overlay_plot = gsn_csm_contour(wks,ss_var3(4,:,:),res_ss)
> overlay(plot(16),overlay_plot)
> plot(17) = gsn_csm_contour_map_polar(wks,ext_var3(5,:,:),res36)
> overlay_plot = gsn_csm_contour(wks,ss_var3(5,:,:),res_ss)
> overlay(plot(17),overlay_plot)
>
> ; 4th row
> plot(18) = gsn_csm_contour_map_polar(wks,ext_var4(0,:,:),res41)
> overlay_plot = gsn_csm_contour(wks,ss_var4(0,:,:),res_ss)
> overlay(plot(18),overlay_plot)
> plot(19) = gsn_csm_contour_map_polar(wks,ext_var4(1,:,:),res42)
> overlay_plot = gsn_csm_contour(wks,ss_var4(1,:,:),res_ss)
> overlay(plot(19),overlay_plot)
> plot(20) = gsn_csm_contour_map_polar(wks,ext_var4(2,:,:),res43)
> overlay_plot = gsn_csm_contour(wks,ss_var4(2,:,:),res_ss)
> overlay(plot(20),overlay_plot)
> plot(21) = gsn_csm_contour_map_polar(wks,ext_var4(3,:,:),res44)
> overlay_plot = gsn_csm_contour(wks,ss_var4(3,:,:),res_ss)
> overlay(plot(21),overlay_plot)
> plot(22) = gsn_csm_contour_map_polar(wks,ext_var4(4,:,:),res45)
> overlay_plot = gsn_csm_contour(wks,ss_var4(4,:,:),res_ss)
> overlay(plot(22),overlay_plot)
> plot(23) = gsn_csm_contour_map_polar(wks,ext_var4(5,:,:),res46)
> overlay_plot = gsn_csm_contour(wks,ss_var4(5,:,:),res_ss)
> overlay(plot(23),overlay_plot)
>
> ; 5th row
> plot(24) = gsn_csm_contour_map_polar(wks,ext_var5(0,:,:),res51)
> overlay_plot = gsn_csm_contour(wks,ss_var5(0,:,:),res_ss)
> overlay(plot(24),overlay_plot)
> plot(25) = gsn_csm_contour_map_polar(wks,ext_var5(1,:,:),res52)
> overlay_plot = gsn_csm_contour(wks,ss_var5(1,:,:),res_ss)
> overlay(plot(25),overlay_plot)
> plot(26) = gsn_csm_contour_map_polar(wks,ext_var5(2,:,:),res53)
> overlay_plot = gsn_csm_contour(wks,ss_var5(2,:,:),res_ss)
> overlay(plot(26),overlay_plot)
> plot(27) = gsn_csm_contour_map_polar(wks,ext_var5(3,:,:),res54)
> overlay_plot = gsn_csm_contour(wks,ss_var5(3,:,:),res_ss)
> overlay(plot(27),overlay_plot)
> plot(28) = gsn_csm_contour_map_polar(wks,ext_var5(4,:,:),res55)
> overlay_plot = gsn_csm_contour(wks,ss_var5(4,:,:),res_ss)
> overlay(plot(28),overlay_plot)
> plot(29) = gsn_csm_contour_map_polar(wks,ext_var5(5,:,:),res56)
> overlay_plot = gsn_csm_contour(wks,ss_var5(5,:,:),res_ss)
> overlay(plot(29),overlay_plot)
>
> ; 6th row
> plot(30) = gsn_csm_contour_map_polar(wks,ext_var6(0,:,:),res61)
> overlay_plot = gsn_csm_contour(wks,ss_var6(0,:,:),res_ss)
> overlay(plot(30),overlay_plot)
> plot(31) = gsn_csm_contour_map_polar(wks,ext_var6(1,:,:),res62)
> overlay_plot = gsn_csm_contour(wks,ss_var4(1,:,:),res_ss)
> overlay(plot(31),overlay_plot)
> plot(32) = gsn_csm_contour_map_polar(wks,ext_var6(2,:,:),res63)
> overlay_plot = gsn_csm_contour(wks,ss_var6(2,:,:),res_ss)
> overlay(plot(32),overlay_plot)
> plot(33) = gsn_csm_contour_map_polar(wks,ext_var6(3,:,:),res64)
> overlay_plot = gsn_csm_contour(wks,ss_var6(3,:,:),res_ss)
> overlay(plot(33),overlay_plot)
> plot(34) = gsn_csm_contour_map_polar(wks,ext_var6(4,:,:),res65)
> overlay_plot = gsn_csm_contour(wks,ss_var6(4,:,:),res_ss)
> overlay(plot(34),overlay_plot)
> plot(35) = gsn_csm_contour_map_polar(wks,ext_var6(5,:,:),res66)
> overlay_plot = gsn_csm_contour(wks,ss_var6(5,:,:),res_ss)
> overlay(plot(35),overlay_plot)
>
> ; 7th row
> plot(36) = gsn_csm_contour_map_polar(wks,ext_var7(0,:,:),res71)
> overlay_plot = gsn_csm_contour(wks,ss_var7(0,:,:),res_ss)
> overlay(plot(36),overlay_plot)
> plot(37) = gsn_csm_contour_map_polar(wks,ext_var7(1,:,:),res72)
> overlay_plot = gsn_csm_contour(wks,ss_var7(1,:,:),res_ss)
> overlay(plot(37),overlay_plot)
> plot(38) = gsn_csm_contour_map_polar(wks,ext_var7(2,:,:),res73)
> overlay_plot = gsn_csm_contour(wks,ss_var7(2,:,:),res_ss)
> overlay(plot(38),overlay_plot)
> plot(39) = gsn_csm_contour_map_polar(wks,ext_var7(3,:,:),res74)
> overlay_plot = gsn_csm_contour(wks,ss_var7(3,:,:),res_ss)
> overlay(plot(39),overlay_plot)
> plot(40) = gsn_csm_contour_map_polar(wks,ext_var7(4,:,:),res75)
> overlay_plot = gsn_csm_contour(wks,ss_var7(4,:,:),res_ss)
> overlay(plot(40),overlay_plot)
> plot(41) = gsn_csm_contour_map_polar(wks,ext_var7(5,:,:),res76)
> overlay_plot = gsn_csm_contour(wks,ss_var7(5,:,:),res_ss)
> overlay(plot(41),overlay_plot)
>
> ; 7th row
> plot(42) = gsn_csm_contour_map_polar(wks,ext_var8(0,:,:),res81)
> overlay_plot = gsn_csm_contour(wks,ss_var7(0,:,:),res_ss)
> overlay(plot(42),overlay_plot)
> plot(43) = gsn_csm_contour_map_polar(wks,ext_var8(1,:,:),res82)
> overlay_plot = gsn_csm_contour(wks,ss_var7(1,:,:),res_ss)
> overlay(plot(43),overlay_plot)
> plot(44) = gsn_csm_contour_map_polar(wks,ext_var8(2,:,:),res83)
> overlay_plot = gsn_csm_contour(wks,ss_var7(2,:,:),res_ss)
> overlay(plot(44),overlay_plot)
> plot(45) = gsn_csm_contour_map_polar(wks,ext_var8(3,:,:),res84)
> overlay_plot = gsn_csm_contour(wks,ss_var7(3,:,:),res_ss)
> overlay(plot(45),overlay_plot)
> plot(46) = gsn_csm_contour_map_polar(wks,ext_var8(4,:,:),res85)
> overlay_plot = gsn_csm_contour(wks,ss_var7(4,:,:),res_ss)
> overlay(plot(46),overlay_plot)
> plot(47) = gsn_csm_contour_map_polar(wks,ext_var8(5,:,:),res86)
> overlay_plot = gsn_csm_contour(wks,ss_var7(5,:,:),res_ss)
> overlay(plot(47),overlay_plot)
>
> ;-------------------------------------------------------------------------
> res_pan = True
> res_pan at gsnPaperOrientation = "portrait"
> res_pan at gsnMaximize = True
> res_pan at gsnPanelLabelBar = False
> res_pan at txString = "CCSM4 Node (my_comp_x,my_comp_y)"
>
> res_pan at gsnPanelYWhiteSpacePercent = 2.
> res_pan at gsnPanelXWhiteSpacePercent = 2.
>
> gsn_panel(wks,plot,(/8,6/),res_pan)
>
> end
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
--
Adam Phillips
Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171211/c3b41788/attachment.html>
More information about the ncl-talk
mailing list