<div dir="ltr"><div>Hi,<br>I am trying to make a panel plot of two data fields of the same dimensions: time x lat x lon. I want to show one data field overlay above the other. Here is the part of script for plotting. <br></div>
<div>My questions are<br></div><div>1. why does the base map not show?<br></div><div>2. how do I remove the text string "CONTOUR FROM ....."?<br></div><div>3. how do I draw the second labelbar for the contour?<br>
</div>Thanks!<br><br>;-----------------------------------<br><br> cmap = (/ (/255, 255, 255/) ,\<br> (/ 0, 0, 0/) ,\<br> (/143, 97, 84 /) ,\ ;0<br> (/219, 205, 202/) ,\ ;1<br> (/246, 229, 149/) ,\ ;2<br>
(/162, 235, 130/) ,\ ;3<br> (/94, 224, 116/) ,\ ;4<br> (/0, 201, 50 /) ,\ ;5<br> (/255, 255, 255/) ,\ ;6<br> (/215, 239, 249/) ,\ ;7<br> (/145, 209, 242/) ,\ ;8<br>
(/55, 117, 185/) ,\ ;9<br> (/30, 61, 124/) /)/255. ;10<br><br>;---------------------------------------<br> print("create plots")<br> res1 = True<br> res1@gsnDraw = False <br>
res1@gsnFrame = False <br> res1@gsnMaximize = False<br><br> res1@cnFillOn = True <br> res1@cnLinesOn = False <br> res1@cnLineLabelsOn = False <br>
res1@cnFillMode = "RasterFill"<br> res1@cnLevelSelectionMode = "ExplicitLevels"<br><br> res1@gsnAddCyclic = False <br> res1@lbLabelBarOn = False<br><br> res1@tmXTOn = False<br>
res1@tmYROn = False<br> res1@tmXBOn = True<br> res1@tmYLOn = True<br> res1@tmXBValues = ispan(toint(lonL),toint(lonR),5)<br> res1@tmXBLabels = res1@tmXBValues+"E"<br>
res1@tmXBMinorValues = ispan(toint(lonL),toint(lonR),1)<br> res1@tmXBLabelFontHeightF = 0.025<br> res1@tmYLLabelFontHeightF = res1@tmXBLabelFontHeightF<br><br> res2 = res1<br><br> res1@mpLimitMode = "LatLon" <br>
res1@mpMinLatF = latS<br> res1@mpMinLonF = lonL<br> res1@mpMaxLatF = latN<br> res1@mpMaxLonF = lonR<br><br> res1@mpFillOn = False <br> res1@mpOutlineOn = True <br>
res1@mpDataBaseVersion = "MediumRes"<br> res1@mpOutlineBoundarySets = "AllBoundaries" <br> res1@mpNationalLineThicknessF = 1.5<br> res1@mpGeophysicalLineThicknessF = 1.5<br> res1@mpGridAndLimbOn = True<br>
res1@mpGridSpacingF = 5.0<br> res1@mpGridLineThicknessF = 0.5<br> res1@mpGridLineDashPattern = 2<br><br> res1@cnLevels = (/.1, .15, .2, .3, .4/)<br> res1@cnFillColors = (/2, 3, 4, 5, 6, 7/) <br>
<br> type = "newpng"<br> outfile = "Figure_1"<br> wks = gsn_open_wks(type,outfile)<br> gsn_define_colormap(wks,cmap)<br><br><br> plts1 = new(nmon,graphic)<br> do i = 0, nmon-1<br><br> if (i.le.nmon-4) then<br>
res1@tmXBOn = False<br> else<br> res1@tmXBOn = True<br> end if<br><br> if (i%3.eq.0) then<br> res1@tmYLOn = True<br> else<br> res1@tmYLOn = False<br> end if<br><br> plts1(i)=gsn_csm_contour_map_ce(wks,ndvi_dom(i,:,:),res1)<br>
end do<br><br> res2@tmXBOn = False<br> res2@tmYLOn = False<br> res2@cnLevels = (/.2, .4, .6, .8/)<br> res2@cnFillColors = (/-1, 9, 10, 11, 12/)<br><br> rtxt = True<br>
rtxt@txJust = "CenterLeft"<br> rtxt@txFontHeightF = 0.035<br> rtxt@PerimOn = False<br> rtxt@txBackgroundFillColor= "white"<br><br> plts2 = new(nmon,graphic)<br>
txts = new(nmon,graphic)<br> do i = 0, nmon-1<br><br> plts2(i)= gsn_csm_contour(wks,snow_dom(i,:,:),res2) <br> txts(i) = gsn_add_text(wks,plts2(i),monstr(i),50.5,53.5,rtxt)<br><br> overlay(plts1(i),plts2(i))<br>
end do<br><br> resP = True<br> resP@gsnMaximize = True<br> resP@gsnPanelLabelBar = True<br> resP@gsnPanelCenter = False <br> resP@gsnFrame = False<br>
<br> resP@pmLabelBarOrthogonalPosF = -0.01<br><br> resP@gsnPanelBottom = 0.05<br> resP@gsnPanelTop = 0.98<br> resP@gsnPanelLeft = 0.05<br> resP@gsnPanelRight = 0.98 <br><br>
resP@gsnPanelScalePlotIndex = 4<br> resP@gsnPanelXWhiteSpacePercent = 5<br> resP@gsnPanelYWhiteSpacePercent = 2<br><br> gsn_panel(wks,plts1,(/4,3/),resP)<br> frame(wks)<br><div><br clear="all"><br>-- <br><font><i style="font-family:arial narrow,sans-serif">Xin Xi<br>
School of Earth and Atmospheric Sciences<br>Georgia Institute of Technology<br>311 Ferst Drive, Atlanta, GA 30332-0340</i></font><br>
</div></div>