<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 12 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:"Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"\@Malgun Gothic";
panose-1:2 11 5 3 2 0 0 2 0 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.apple-tab-span
{mso-style-name:apple-tab-span;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Hello,
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Try this, please.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Id = new(3, graphic)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Id(0) = gsn_add_shapefile_polylines(wks,plot(0),shp_name1,lnres)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Id(1) = gsn_add_shapefile_polylines(wks,plot(1),shp_name1,lnres)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Id(2) = gsn_add_shapefile_polylines(wks,plot(2),shp_name1,lnres)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif"">Cheers.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> ncl-talk-bounces@ucar.edu [mailto:ncl-talk-bounces@ucar.edu]
<b>On Behalf Of </b>grace<br>
<b>Sent:</b> 28 July 2016 03:47<br>
<b>To:</b> ncl-talk<br>
<b>Subject:</b> [ncl-talk] The shp only appear on the last plot of a (/3,1/)panel<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Hi:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> All,I have write a script to plot a (/3,1/)panel,but the shp I added only appeared on the last plot.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">The description only worked on the first plot.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><img width="344" height="809" id="_x0000_i1025" src="cid:image001.jpg@01D1E8BA.E6FD8330"><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">This is my script:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">begin<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; The WRF ARW input file. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; This needs to have a ".nc" appended, so just do it.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> a = addfile("/public/home/huanglei/data/20160724/oldwrf/wrfout_d03_2016-07-23_12:00:00"+".nc","r")<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">; We generate plots, but what kind do we prefer?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; type = "x11"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> type = "pdf"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; type = "ps"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; type = "ncgm"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> wks = gsn_open_wks(type,"plt_Precip_oldwrf3km")<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> plot = new(3,graphic)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">; Set some basic resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> res = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres = True ; Map resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres@mpOutlineOn = False ; Turn off map outlines<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres@mpFillOn = False ; Turn off map fill<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> mpres@mpGridAndLimbOn = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ;res@mpProjection = "Lambert"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> pltres = True ; Plot resources<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> pltres@PanelPlot = True ; Tells wrf_map_overlays not to remove overlays<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">; What times and how many time steps are in the data set?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> FirstTime = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> times = wrf_user_getvar(a,"times",-1) ; get all times in the file<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ntimes = dimsizes(times) ; number of times in the file<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; print(times)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; exit<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> it_end = 29<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; First get the variables we will need <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> ; Get non-convective, convective and total precipitation of 5km <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_exp = wrf_user_getvar(a,"RAINNC",it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_con = wrf_user_getvar(a,"RAINC",it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot = rain_exp + rain_con<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot@description = "Total Precipitation"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><span class="apple-tab-span"> </span>;calculate the precipitation <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> it_start = 22<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_exp_save = wrf_user_getvar(a,"RAINNC",it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_con_save = wrf_user_getvar(a,"RAINC",it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot_save = rain_exp_save + rain_con_save<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> times_sav = times(it_start)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot_tend = rain_tot - rain_tot_save<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rainc_tend = rain_con - rain_con_save ; CUMULUS PRECIPITATION<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> rainnc_tend= rain_exp - rain_exp_save ; SCALE PRECIPITATION<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> rain_tot_tend@description = "Precipitation of 3km(old wrf)"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> ; Plotting options for Precipitation<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r = res <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@UnitLabel = "mm"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnLevelSelectionMode = "ExplicitLevels"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnLevels = (/ .1, .2, .4, .8, 1.6, 3.2, 6.4, \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 12.8, 25.6, 51.2/)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnFillColors = (/"White","DarkOliveGreen1", \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> "DarkOliveGreen3","Chartreuse", \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> "Chartreuse3","Green","ForestGreen", \<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> "Yellow","Orange","Red","Violet"/)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnInfoLabelOn = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnConstFLabelOn = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@cnFillOn = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@vpHeightF = 0.1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@vpWidthF = 0.9 <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@lgLegendOn = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@gsnDraw = False <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@gsnFrame = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@lbLabelBarOn = False <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@Footer = False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> opts_r@NoHeaderFooter =True <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; Precipitation Tendencies<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; opts_r@SubFieldTitle = "from " + times(it_start) + " to " + times(it_end)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> contour_tend = wrf_contour(a,wks, rain_tot_tend,opts_r) ; total (color)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> contour_rainc_tend = wrf_contour(a,wks, rainc_tend,opts_r) ; total cumulus precipitation (color)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> contour_rainnc_tend = wrf_contour(a,wks, rainnc_tend,opts_r) ; total scale precipitation(color)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> delete(opts_r)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> ; MAKE PLOTS <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> plot(0) = wrf_map_overlays(a,wks,contour_tend,pltres,mpres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> plot(1) = wrf_map_overlays(a,wks,contour_rainc_tend,pltres,mpres) <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> plot(2) = wrf_map_overlays(a,wks,contour_rainnc_tend,pltres,mpres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;>============================================================<<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">; add China map<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">;>------------------------------------------------------------<<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> shp_name1 = "/public/home/huanglei/map/xian.shp"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> lnres = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> lnres@gsLineColor = "gray25"<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> lnres@gsLineThicknessF = 0.5 <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> id = gsn_add_shapefile_polylines(wks,plot(0),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id = gsn_add_shapefile_polylines(wks,plot(1),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> id = gsn_add_shapefile_polylines(wks,plot(2),shp_name1,lnres)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">;;;;;;;;;;;;;;;;;;creat panel;;;;;;;;;;;;;;;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnFrame =False<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelLabelBar = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelBottom = 0.05<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelTop = 0.95 ; Make sure not too close to <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@gsnPanelBottom = 0.35 ; edge, so it maximizes better.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> resP@tyPolyDrawList = True<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> gsn_panel(wks,plot,(/3,1/),resP)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">; draw(plot) ; This will draw the map and the shapefile outlines.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> frame(wks)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">end<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif""> How can I slove the problem?</span><o:p></o:p></p>
</div>
</div>
</body>
</html>