<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Rabah,<div><br></div><div>So if you don’t put them both both res and pltres you still end up with 2 plots?</div><div><br></div><div>It might also be worth trying to use the gsn_add_polymarker as opposed to using the gsn-polymarker. gsn_polymarker might have the frame internally built. The gsn_add_polymarker does return an identifier so you would do it similar to your polyline call:</div><div><br></div><div>pm_id = gsn_add_polymarker(wks,plot,npres@ref_lon,npres@ref_lat,mpres)</div><div><br></div><div>If this doesn’t work can you attach your geo_em file and script and I will see if I can see what is going on.</div><div><br></div><div>-Alex</div><div><br></div><div><br></div><div><br><div><div>On May 26, 2015, at 7:06 PM, Rabah Hachelaf &lt;<a href="mailto:hachelaf@sca.uqam.ca">hachelaf@sca.uqam.ca</a>&gt; wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr">Hi Alex ,<div><br></div><div>When i add&nbsp;</div><div><div>&nbsp;opts@gsnDraw = False</div><div>&nbsp;opts@gsnFrame = False &nbsp;&nbsp;</div></div><div>&nbsp;or&nbsp;</div><div><div>&nbsp;res@gsnDraw = False</div><div>&nbsp;res@gsnFrame = False &nbsp;&nbsp;</div></div><div><br></div><div>and also&nbsp;</div><div><br></div><div><div>&nbsp;pltres@gsnDraw = False</div><div>&nbsp;pltres@gsnFrame = False&nbsp;</div></div><div><br></div><div>I get only grid with the box and the marker like attached.&nbsp;</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-05-26 20:47 GMT-04:00 Alex Schaefer <span dir="ltr">&lt;<a href="mailto:alexander.schaefer@mines.sdsmt.edu" target="_blank">alexander.schaefer@mines.sdsmt.edu</a>&gt;</span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Rabah,<div><br></div><div>Where in your script did you add the resources suggested by Jared?</div><div><br></div><div>Since you have a bunch of different resource objects you will need to set it in each going into plots controlling plot resources. For you this should include res or opts and likely pltres.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>-Alex</div></font></span><div><div class="h5"><div><br></div><div><br><div><div>On May 26, 2015, at 6:36 PM, Rabah Hachelaf &lt;<a href="mailto:hachelaf@sca.uqam.ca" target="_blank">hachelaf@sca.uqam.ca</a>&gt; wrote:</div><br><blockquote type="cite"><div dir="ltr">Hi Jared,<div><br></div><div>Thanks for your reponse,</div><div><br></div><div>But add &nbsp;</div><div><br></div><div><br></div><div><div style="font-size:12.8000001907349px">res@gsnDraw = False<br></div><span style="font-size:12.8000001907349px">res@gsnFrame = False</span><br></div><div><br></div><div>doesn't change any things.&nbsp;</div><div><br></div><div>Rabah</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-05-26 18:34 GMT-04:00 Jared Lee <span dir="ltr">&lt;<a href="mailto:jaredlee@ucar.edu" target="_blank">jaredlee@ucar.edu</a>&gt;</span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>Hi Rabah,<br><br></div>I believe what you're missing is these two lines where you set other 'res' resources:<br><br></div>res@gsnDraw = False<br></div>res@gsnFrame = False<br><br></div>By default they're both set to True, which would give you two separate images. If you set both to false you should get a single plot with your polymarker overlaid.<br><br></div><div>Hope that helps.<br></div><div><br></div>Jared<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div>On Tue, May 26, 2015 at 4:21 PM, Rabah Hachelaf <span dir="ltr">&lt;<a href="mailto:hachelaf@sca.uqam.ca" target="_blank">hachelaf@sca.uqam.ca</a>&gt;</span> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div dir="ltr"><div><br></div><div><br></div><div>Hello NCL users,</div><div><br></div><div>I am using the script below to plot a box and a marker dot over topography field from WRF outputs,</div><div><br></div><div>My problem is i can not get one figure, i get two separate figures like attached.&nbsp;</div><div><br></div><div>Any suggestion to solve this.</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>; &nbsp; Example script - plot terrain (using some basic WRF scripts)</div><div>; &nbsp; Overwrite basic map background settings</div><div>; &nbsp; November 2008</div><div><br></div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"</div><div><br></div><div>begin</div><div><br></div><div>&nbsp; a = addfile("../<a href="http://geo_em.d02.nc/" target="_blank">geo_em.d02.nc</a>","r") &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Open a file</div><div>&nbsp; filename = "../namelist.wps"</div><div><br></div><div>&nbsp; npres= True</div><div>&nbsp; npres@ref_lat &nbsp; = stringtofloat(systemfunc("grep ref_lat &nbsp; &nbsp;" +filename+ " | cut -f2 -d'=' | cut -f1 -d','" ) )</div><div>&nbsp; npres@ref_lon &nbsp; = stringtofloat(systemfunc("grep ref_lon &nbsp; &nbsp;" +filename+ " | cut -f2 -d'=' | cut -f1 -d','" ) )</div><div><br></div><div>&nbsp;type = "png"</div><div>&nbsp; wks = gsn_open_wks(type,"plt_geo_6") &nbsp; &nbsp; &nbsp; &nbsp; ; Create a plot workstation</div><div><br></div><div><br></div><div>&nbsp; opts = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Set some Basic Plot options</div><div>&nbsp; opts@MainTitle = "GEOGRID FIELDS"</div><div>&nbsp; opts@InitTime = False &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Do not plot time or footers</div><div>&nbsp; opts@Footer = False &nbsp; &nbsp;</div><div><br></div><div>&nbsp; ter = wrf_user_getvar(a,"HGT_M",0) &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Read the variable to memory</div><div>&nbsp; lon = a-&gt;XLONG_U(0,:,:)</div><div>&nbsp; lat = a-&gt;XLAT_U(0,:,:)</div><div>&nbsp;printVarSummary(lat)&nbsp;</div><div>&nbsp;</div><div>&nbsp;print (lat(0,0)+", "+lon(0,0)+",0,")</div><div>&nbsp;print (lat(0,60)+", "+lon(0,60)+",0,")</div><div>&nbsp;print (lat(59,60)+", "+lon(59,60)+",0,")</div><div>&nbsp;print (lat(59,0)+", "+lon(59,0)+",0")</div><div>&nbsp;</div><div>&nbsp;</div><div>&nbsp; res = opts &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Use basic options for this field</div><div>&nbsp; res@cnFillOn = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Create a color fill plot</div><div>&nbsp; res@ContourParameters = (/ 0., 1000., 50. /) ; Set the levels</div><div><br></div><div>&nbsp; contour = wrf_contour(a,wks,ter,res)</div><div><br></div><div>&nbsp; pltres = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; Set plot options</div><div>&nbsp; mpres = True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; Set map options</div><div>&nbsp; mpres@mpGeophysicalLineColor &nbsp; &nbsp; &nbsp;= "Black" &nbsp;; Overwrite basic map settings</div><div>&nbsp; mpres@mpGridLineColor &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = "Black"</div><div>&nbsp; mpres@mpLimbLineColor &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = "Black"</div><div>&nbsp; mpres@mpNationalLineColor &nbsp; &nbsp; &nbsp; &nbsp; = "Black"</div><div>&nbsp; mpres@mpPerimLineColor &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= "Black"</div><div>&nbsp; mpres@mpUSStateLineColor &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= "Black"</div><div>&nbsp; &nbsp; &nbsp;mpres@mpProjection = "LambertConformal"</div><div>&nbsp; &nbsp; &nbsp; &nbsp;mpres@mpLambertParallel1F = 30</div><div>&nbsp; &nbsp; &nbsp; &nbsp;mpres@mpLambertParallel2F = 60&nbsp;</div><div>&nbsp; &nbsp; &nbsp; &nbsp;mpres@mpLambertMeridianF = -98</div><div>&nbsp; &nbsp; &nbsp; &nbsp;mpres@mpGridSpacingF &nbsp; &nbsp; = 0.25</div><div>;************************************************</div><div>; add the box</div><div>;************************************************</div><div>;************************************************</div><div>&nbsp;ypts = (/ lat(10,10), lat(10,50), lat(49,50), &nbsp;lat(49,10), &nbsp;lat(10,10)/)</div><div>&nbsp;xpts = (/ lon(10,10), lon(10,50), lon(49,50), &nbsp;lon(49,10), &nbsp;lon(10,10)/)</div><div>;************************************************</div><div><br></div><div>&nbsp; resp &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= True &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; polyline mods desired</div><div>&nbsp; resp@gsLineColor &nbsp; &nbsp; &nbsp;= "black" &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; color of lines</div><div>&nbsp; resp@gsLineThicknessF = 2.0 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; thickness of lines</div><div>; &nbsp;resp@gsLineLabelString= "test" &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; adds a line label string</div><div>&nbsp; &nbsp;</div><div><br></div><div>&nbsp; &nbsp; &nbsp; &nbsp;</div><div>&nbsp; plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres) ; Plot field over map background</div><div>&nbsp;&nbsp;</div><div>&nbsp;&nbsp;</div><div>&nbsp; dum = new(4,graphic)&nbsp;</div><div>&nbsp; do i = 0 , 3</div><div>&nbsp; &nbsp; dum(i)=gsn_add_polyline(wks,plot,xpts(i:i+1),ypts(i:i+1),resp) &nbsp; &nbsp; &nbsp;</div><div>&nbsp; end do&nbsp;</div><div>;*********************************************************&nbsp;</div><div>&nbsp; mpres@gsMarkerColor = "Black"</div><div>&nbsp; mpres@gsMarkerIndex = 16</div><div>&nbsp; mpres@gsMarkerSizeF = 0.005</div><div>&nbsp; gsn_polymarker(wks,plot,npres@ref_lon,npres@ref_lat,mpres)</div><div>&nbsp; &nbsp;</div><div>&nbsp; draw(plot)</div><div>&nbsp; frame(wks)</div><div>&nbsp;</div><div>end</div><span><font color="#888888"><div><br></div>-- <br><div><div dir="ltr"><div>------------------------------</div>Cordialement,<br>Best regards,<br>Rabah Hachelaf&nbsp;<br></div></div>
</font></span></div>
<br></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><br>-- <br><div><div dir="ltr"><div><div><div><div><div><font face="courier new, monospace">===============================<br>Jared A. Lee, Ph.D.<br></font></div><font face="courier new, monospace">Project Scientist I<br></font></div><font face="courier new, monospace">Research Applications Laboratory<br></font></div><font face="courier new, monospace">National Center for Atmospheric Research<br>Boulder, Colorado, USA<br><br></font></div><font face="courier new, monospace">Email: <a href="mailto:jaredlee@ucar.edu" target="_blank">jaredlee@ucar.edu</a> (w)<br></font></div><font face="courier new, monospace">Phone: 303.497.8485 (w)</font><div><font face="courier new, monospace">Web: <a href="https://staff.ucar.edu/users/jaredlee" target="_blank">https://staff.ucar.edu/users/jaredlee</a><br>===============================</font><br></div></div></div>
</font></span></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr"><div>------------------------------</div>Cordialement,<br>Best regards,<br>Rabah Hachelaf&nbsp;<br></div></div>
</div>
_______________________________________________<br>ncl-talk mailing list<br>List instructions, subscriber options, unsubscribe:<br><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br></blockquote></div><br></div></div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>------------------------------</div>Cordialement,<br>Best regards,<br>Rabah Hachelaf&nbsp;<br></div></div>
</div>
<span>&lt;plt_geo_6.png&gt;</span></blockquote></div><br></div></body></html>