[ncl-talk] Problem with panel plot (size difference)

Tomoko Koyama Tomoko.Koyama at Colorado.EDU
Fri Dec 8 22:54:28 MST 2017


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



More information about the ncl-talk mailing list