[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