<div dir="ltr"><div class="gmail_default" style="font-size:small">Kerwyn,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Please respond back to the ncl-talk email list with follow-up questions, and not just to me personally.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It helps if you can provide a clean script, because it&#39;s hard to go through a script that has lots of code commented out.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I took your code and tried to clean it up and use it on one of my own WRF output files.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Are are a few issues:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">[1] You are not providing lat/lon coordinate arrays for all the variables you are plotting. You have several cases of this:<br></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace"><br></font></div><div class="gmail_default" style=""><div class="gmail_default" style=""><font face="monospace, monospace">  froz@lon2d = lon2d</font></div><div class="gmail_default" style=""><font face="monospace, monospace">  froz@lat2d = lat2d</font></div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">which is the right idea, but you need to do this for the variables you&#39;re actually plotting, like &quot;z_plane&quot;:</div><div class="gmail_default" style=""><font face="monospace, monospace"><br></font></div><div class="gmail_default" style=""><div class="gmail_default"><font face="monospace, monospace">  z_plane@lon2d = lon2d</font></div><div class="gmail_default"><font face="monospace, monospace">  z_plane@lat2d = lat2d</font></div><div class="gmail_default"><br></div><div class="gmail_default">[2] You need to always set res@gsnAddCyclic = False. You had this in a few places, but commented out:</div><div class="gmail_default">







<p class=""><font face="monospace, monospace"><span class=""> </span><span class="">; res@</span><span class="">gsnAddCy</span><span class="">clic = False                                                     </span></font></p><p class=""><span class="">You need to uncomment these for every single plot that you are creating that will be overlaid on a map.</span></p></div><div class="gmail_default">[3] The &quot;geo&quot; plot was appearing by itself because you hadn&#39;t set:</div><div class="gmail_default"><br></div><div class="gmail_default"><div class="gmail_default">res5@gsnDraw = False</div><div class="gmail_default">res5@gsnFrame = False</div><div class="gmail_default"><br></div><div class="gmail_default">And, it wasn&#39;t showing up on the &quot;contour&quot; plot because of the issue mentioned in [1].</div><div class="gmail_default"><br></div></div><div class="gmail_default">[4] You have hard-coded your min/max lat/lon resources, which I don&#39;t normally recommend. I suggest using the min/max of your lat2d/lon2d arrays:</div><div class="gmail_default"><br></div><div class="gmail_default"><br></div></div><div class="gmail_default" style=""><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">  res@mpMinLatF     = min(lat2d) ; 37.85                                        </font></div><div class="gmail_default"><font face="monospace, monospace">  res@mpMaxLatF     = max(lat2d) ; 38.50                                        </font></div><div class="gmail_default"><font face="monospace, monospace">  res@mpMinLonF     = min(lon2d) ; -120.0                                       </font></div><div class="gmail_default"><font face="monospace, monospace">  res@mpMaxLonF     = max(lon2d) ; -119.35  </font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif">If these suggestions don&#39;t fix the problem, please post back to ncl-talk.</font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif">Thanks,</font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif"><br></font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif">--Mary</font></div><div class="gmail_default"><font face="arial, helvetica, sans-serif"><br></font></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 4, 2016 at 3:38 PM, Kerwyn Texeira <span dir="ltr">&lt;<a href="mailto:ktish86@gmail.com" target="_blank">ktish86@gmail.com</a>&gt;</span> wrote:<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><div>Hi ncl-talk,<br><br></div>I not getting the errors anymore but I have two problems<br><br></div>1. The geopotential contours are not smooth.  It&#39;s plotting like it&#39;s contours at the surface and not at 600hpa<br><br></div>2. Also the geopotential heights are not overlaying the terrain.  Instead, its plotting separately. <br><br></div>I would like to get smoother geopotential heights at 600hpa and overlay with winds over the terrain zoomed in.<br><br></div>How do I overcome these two issues? Can anyone help with these please?  I have attached the plots. <br><br></div><div><div class="h5">My Script:<br><br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;<br> <br>begin<br>   a = addfile(&quot;./wrfout_d03_2014-01-11_22:40:<a href="http://00.nc" target="_blank">00.nc</a>&quot;,&quot;r&quot;)<br><br>  it = 0<br>  hgt       = wrf_user_getvar(a, &quot;HGT&quot;, it)<br>  hgt@lat2d = wrf_user_getvar(a, &quot;XLAT&quot;, it)<br>  hgt@lon2d = wrf_user_getvar(a, &quot;XLONG&quot;, it)<br>  u         = wrf_user_getvar(a, &quot;ua&quot;, it)<br>  v         = wrf_user_getvar(a, &quot;va&quot;, it)<br>  p         = wrf_user_getvar(a, &quot;pressure&quot;, it)<br>  qc        = wrf_user_getvar(a, &quot;QCLOUD&quot;, it)<br>  qr        = wrf_user_getvar(a, &quot;QRAIN&quot;, it)<br>  qs        = wrf_user_getvar(a, &quot;QSNOW&quot;, it)<br>  qi        = wrf_user_getvar(a, &quot;QICE&quot;, it)<br>  pre       = wrf_user_getvar(a, &quot;RAINNC&quot;, it)<br>   z        = wrf_user_getvar(a, &quot;z&quot;, it)<br><br>  u_wind    = wrf_user_intrp3d(u, p, &quot;h&quot;, 600., 0.0, False)<br>  v_wind    = wrf_user_intrp3d(v, p, &quot;h&quot;, 600., 0.0, False)<br>  qcl       = wrf_user_intrp3d(qc, p, &quot;h&quot;, 600., 0.0, False)<br>  qrn       = wrf_user_intrp3d(qr, p, &quot;h&quot;, 600., 0.0, False)<br>  qsn       = wrf_user_intrp3d(qs, p, &quot;h&quot;, 600., 0.0, False)<br>  qice      = wrf_user_intrp3d(qi, p, &quot;h&quot;, 600., 0.0, False)<br>  z_plane   = wrf_user_intrp3d(z, p, &quot;h&quot;, 600., 0.0, False)<br><br>  qcl = qcl*1000<br>  qrn = qrn*1000<br>  qsn = qsn*1000<br>  qice = qice*1000<br>  precip = pre*0.03937<br>  froz = qice + qsn<br>  melt = qcl + qrn<br>   <br>  froz@lon2d = hgt@lon2d<br>  froz@lat2d = hgt@lat2d<br>  <br>  melt@lon2d = hgt@lon2d<br>  melt@lat2d = hgt@lat2d<br>  <br> ; qcl@lon2d = hgt@lon2d<br> ; qcl@lat2d = hgt@lat2d<br> ; qrn@lon2d = hgt@lon2d<br> ; qrn@lat2d = hgt@lat2d<br> ; qsn@lon2d = hgt@lon2d<br> ; qsn@lat2d = hgt@lat2d<br> ; qice@lon2d = hgt@lon2d<br> ; qice@lat2d = hgt@lat2d<br>  <br><br>   spd = (u_wind*u_wind + v_wind*v_wind)^(0.5)  ;m/s<br>   u_wind = u_wind*1.94384449<br>   v_wind = v_wind*1.94384449<br><br>   u_wind@lon2d =   hgt@lon2d<br>   u_wind@lat2d =   hgt@lat2d<br><br>   v_wind@lon2d =   hgt@lon2d<br>   v_wind@lat2d =   hgt@lat2d<br><br> ; spd = spd*1.94384449<br> ; spd@units = &quot;Wind Speed&quot;<br> ; spd@units = &quot;kts&quot;<br>  <br>    wks_type = &quot;png&quot;<br>    wks_type@wkWidth = 2500<br>    wks_type@wkHeight = 2500<br>    wks = gsn_open_wks(wks_type,&quot;geo&quot;)         ; send graphics to PNG file<br>  <br>  <br>; gsn_define_colormap(wks,&quot;BlAqGrYeOrRe&quot;)<br>  gsn_define_colormap(wks,&quot;matlab_jet&quot;) <br>  <br>  res = True<br>  res@gsnDraw      =  False                   ; do not draw the plot<br>  res@gsnFrame     =  False                   ; do not advance the frame<br>  res@cnLineLabelsOn       = False            ; do not use line labels<br>  res@cnFillOn             = True             ; color fill<br>  res@cnLinesOn            = False            ; do not draw contour lines<br>  <br><br><br>  res@tiMainString = &quot;&quot;<br>  res@pmTickMarkDisplayMode  = &quot;Always&quot;<br>  res@mpProjection  = &quot;CylindricalEquidistant&quot;    ;The default<br>  res@mpDataBaseVersion      = &quot;MediumRes&quot;<br>  res@mpOutlineOn            =True<br>  res@lbOrientation          = &quot;Vertical&quot;<br>  res@tiMainOffsetYF         = -0.03<br>  res@mpFillOn     = False<br>  res@mpOutlineOn  = True                  ; turn the map outline on<br>  res@mpMinLatF     = 37.85<br>  res@mpMaxLatF     = 38.50<br>  res@mpMinLonF     = -120.0<br>  res@mpMaxLonF     = -119.35<br>  res@gsnLeftString = &quot;&quot;<br>  res@gsnCenterString = &quot;Geopotential Height at 600 hpa on Jan 11 at 22:40UTC&quot;<br>  res@gsnStringFontHeightF = 0.020<br>  res@gsnRightString = &quot;&quot;<br>  res@gsnMaximize     = True<br>  res@mpShapeMode     = &quot;FreeAspect&quot;  <br>  res@lbTitleString   = &quot;Terrain (m)&quot;<br>  res@lbTitlePosition = &quot;Right&quot;<br>  res@lbTitleDirection = &quot;Across&quot;<br>  res@lbTitleAngleF    = 90.<br>  res@lbTitleFontHeightF = 0.015<br></div></div>  ;res@trGridType = &quot;TrianglarMesh&quot;<br> ; res@gsnAddCyclic = False <br><div><div class="h5"><br> ;------wind vectors<br>  res2 = True<br>  res2@gsnDraw = False<br>  res2@gsnFrame = False<br>  res2@vcWindBarbLineThicknessF= 3.0<br>  res2@vcRefLengthF= 0.018<br>  res2@vcRefMagnitudeF= 10<br>  res2@vcMinDistanceF = 0.05<br>  res2@vcGlyphStyle = &quot;WindBarb&quot;<br>  res2@gsnLeftString = &quot;&quot;<br>  res2@gsnRightString = &quot;&quot;<br>  res2@vcRefAnnoOn = False<br>  <br><br>;-------Frozen<br>;res3 = True<br>;res3@cnLineColor = &quot;Blue&quot;<br>;res3@gsnContourNegLineDashPattern = 2.0<br>;res3@gsnContourLineThicknessesScale = 2.0<br>;res3@gsnContourPosLineDashPattern = 3.0<br>;res3@gsnContourZeroLineThicknessF = 0.0<br>;solid = gsn_csm_contour(wks,froz,res3)<br><br>;------Melted<br>;res4 = True<br>;res4@cnLineColor = &quot;Black&quot;<br>;res4@gsnContourLineThicknessesScale = 2.0<br>;res4@cnLineLabelBackgroundColor = -1 ; transparent<br>;res4@cnInfoLabelOn = False<br>;res4@cnLevelSelectionMode = &quot;ExplicitLevels&quot;<br>;res4@cnLevels = (/0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0/)<br>;res4@cnLineLabelFontHeightF = 0.01<br>;res4@cnLineLabelDensityF =2 <br>;res4@cnLineLabelsOn = True<br>;liquid = gsn_csm_contour(wks,melt,res4)<br><br>;------Geopotential Height<br>res5 = True<br></div></div>;res5@cnInfoLabelOn = False<br>;res5@cnLineLabelsOn = False<br>res5@cnLinesOn = True<br>res5@gsnLeftString = &quot;&quot;<br>res5@gsnRightString = &quot;&quot;<br>;res5@cnLevelSelectionMode = &quot;ManualLevels&quot;<br>;res5@cnLineColor = &quot;Black&quot;<br>;res5@cnLineLabelBackgroundColor = -1<br>;res5@cnLineLabelPlacementMode = &quot;constant&quot;<br>;res5@cnLineLabelFontColor = &quot;blue&quot;<br>;res5@cnFillOpacityF = 1.0<br>;res5@cnLabelDrawOrder = &quot;PostDraw&quot;<br>;res5@cnLineDrawOrder = &quot;PostDraw&quot;<br>res5@gsnContourLineThicknessesScale = 3.0<br>res5@trGridType = &quot;TriangularMesh&quot;<br>res5@gsnAddCyclic = False<span class=""><br>geo = gsn_csm_contour(wks, z_plane, res5)<br><br><br>contour = gsn_csm_contour_map(wks,hgt,res)<br>vector = gsn_csm_vector(wks,u_wind,v_wind,res2)<br></span>overlay(contour, vector)<br>overlay(contour, geo)<span class=""><br><br>plres = True<br>plres@gsLineColor = &quot;blue&quot;<br>plres@gsLineThicknessF = 3.0<br><br>lat = (/38.05, 38.05/)<br>lon = (/-120.10, -119.35/)<br><br>a = gsn_add_polyline(wks,contour,lon,lat,plres)<br><br>pltres = True<br>pltres@gsLineColor = &quot;Black&quot;<br>pltres@gsLineThicknessF = 3.0<br><br>lat1 = (/38.27, 38.27/)<br>lon1 = (/-120.10, -119.35/)<br><br>b = gsn_add_polyline(wks,contour,lon1,lat1,pltres)<br><br>draw(contour)<br>frame(wks)<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br><br><br>end<br>    <br><br></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, May 4, 2016 at 7:36 AM, Mary Haley <span dir="ltr">&lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">Kerwyn,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If your lat/lon coordinate arrays have missing values, then this causes issues for the contouring.  Usually, by doing what the error message suggests:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">res@trGridType = &quot;TriangularMesh&quot;</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">then what it does (I think) is throw away any data values that have missing lat/lon pairs, and then it uses a triangular mesh algorithm to construct the contours of the remaining data.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">However, this is WRF data, and if any of your lat2d/lon2d arrays contain missing values, then there&#39;s likely an issue with your WRF output file.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I would do this to make sure your lat/lon are okay:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><span><div class="gmail_default"><font face="monospace, monospace">  lat2d = wrf_user_getvar(a, &quot;XLAT&quot;, it)</font></div></span><span><div class="gmail_default"><font face="monospace, monospace">  lon2d = wrf_user_getvar(a, &quot;XLONG&quot;, it)</font></div></span><div class="gmail_default"><font face="monospace, monospace">  printMinMax(lat2d,0)</font></div><div class="gmail_default"><font face="monospace, monospace">  printMinMax(lon2d,0)</font></div><div class="gmail_default"><font face="monospace, monospace">  hgt@lat2d = lat2d</font></div><div class="gmail_default"><font face="monospace, monospace">  hgt@lon2d = lon2d</font></div><div class="gmail_default"><font face="monospace, monospace">  . . .</font></div><div class="gmail_default"><br></div></div><div class="gmail_default" style="font-size:small">If the lat/lon arrays look okay, then the next thing I would look for is to do a printVarSummary on every single variable you are trying to plot, to make sure it has the correct lat2d/lon2d attached. Also, for WRF data, you want to set res@gsnAddCyclic = False for any data being overlaid on a map, so it doesn&#39;t try to add a longitude cyclic point.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Tue, May 3, 2016 at 9:12 PM, Kerwyn Texeira <span dir="ltr">&lt;<a href="mailto:ktish86@gmail.com" target="_blank">ktish86@gmail.com</a>&gt;</span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div><div><div>Hi ncl-talk,<br><br></div>I&#39;m trying to plot potential height at a certain pressure level across a terrain with winds overlayed on a zoomed in plot. I&#39;m seeing stars.  I&quot;m getting a long message that I don&#39;t have any idea what it means.<br><br></div>Error message:<br>warning:ContourPlotDraw: out of range coordinates encountered; standard ORV rendering method may be unreliable;<br> consider setting the resource trGridType to &quot;TriangularMesh&quot; if coordinates contain missing values<br><br></div>My Script:<br><br>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;<br> <br>begin<br>   a = addfile(&quot;./wrfout_d03_2014-01-11_22:40:<a href="http://00.nc" target="_blank">00.nc</a>&quot;,&quot;r&quot;)<br><br>  it = 0<br>  hgt       = wrf_user_getvar(a, &quot;HGT&quot;, it)<br>  hgt@lat2d = wrf_user_getvar(a, &quot;XLAT&quot;, it)<br>  hgt@lon2d = wrf_user_getvar(a, &quot;XLONG&quot;, it)<br>  u         = wrf_user_getvar(a, &quot;ua&quot;, it)<br>  v         = wrf_user_getvar(a, &quot;va&quot;, it)<br>  p         = wrf_user_getvar(a, &quot;pressure&quot;, it)<br>  qc        = wrf_user_getvar(a, &quot;QCLOUD&quot;, it)<br>  qr        = wrf_user_getvar(a, &quot;QRAIN&quot;, it)<br>  qs        = wrf_user_getvar(a, &quot;QSNOW&quot;, it)<br>  qi        = wrf_user_getvar(a, &quot;QICE&quot;, it)<br>  pre       = wrf_user_getvar(a, &quot;RAINNC&quot;, it)<br>   z        = wrf_user_getvar(a, &quot;z&quot;, it)<br><br>  u_wind    = wrf_user_intrp3d(u, p, &quot;h&quot;, 600., 0.0, False)<br>  v_wind    = wrf_user_intrp3d(v, p, &quot;h&quot;, 600., 0.0, False)<br>  qcl       = wrf_user_intrp3d(qc, p, &quot;h&quot;, 600., 0.0, False)<br>  qrn       = wrf_user_intrp3d(qr, p, &quot;h&quot;, 600., 0.0, False)<br>  qsn       = wrf_user_intrp3d(qs, p, &quot;h&quot;, 600., 0.0, False)<br>  qice      = wrf_user_intrp3d(qi, p, &quot;h&quot;, 600., 0.0, False)<br>  z_plane   = wrf_user_intrp3d(z, p, &quot;h&quot;, 600., 0.0, False)<br><br>  qcl = qcl*1000<br>  qrn = qrn*1000<br>  qsn = qsn*1000<br>  qice = qice*1000<br>  precip = pre*0.03937<br>  froz = qice + qsn<br>  melt = qcl + qrn<br>   <br>  froz@lon2d = hgt@lon2d<br>  froz@lat2d = hgt@lat2d<br>  <br>  melt@lon2d = hgt@lon2d<br>  melt@lat2d = hgt@lat2d<br>  <br> ; qcl@lon2d = hgt@lon2d<br> ; qcl@lat2d = hgt@lat2d<br> ; qrn@lon2d = hgt@lon2d<br> ; qrn@lat2d = hgt@lat2d<br> ; qsn@lon2d = hgt@lon2d<br> ; qsn@lat2d = hgt@lat2d<br> ; qice@lon2d = hgt@lon2d<br> ; qice@lat2d = hgt@lat2d<br>  <br><br>   spd = (u_wind*u_wind + v_wind*v_wind)^(0.5)  ;m/s<br>   u_wind = u_wind*1.94384449<br>   v_wind = v_wind*1.94384449<br><br>   u_wind@lon2d =   hgt@lon2d<br>   u_wind@lat2d =   hgt@lat2d<br><br>   v_wind@lon2d =   hgt@lon2d<br>   v_wind@lat2d =   hgt@lat2d<br><br> ; spd = spd*1.94384449<br> ; spd@units = &quot;Wind Speed&quot;<br> ; spd@units = &quot;kts&quot;<br>  <br>    wks_type = &quot;png&quot;<br>    wks_type@wkWidth = 2500<br>    wks_type@wkHeight = 2500<br>    wks = gsn_open_wks(wks_type,&quot;geo&quot;)         ; send graphics to PNG file<br>  <br>  <br>; gsn_define_colormap(wks,&quot;BlAqGrYeOrRe&quot;)<br>  gsn_define_colormap(wks,&quot;matlab_jet&quot;) <br>  <br>  res = True<br>  res@gsnDraw      =  False                   ; do not draw the plot<br>  res@gsnFrame     =  False                   ; do not advance the frame<br>  res@cnLineLabelsOn       = False            ; do not use line labels<br>  res@cnFillOn             = True             ; color fill<br>  res@cnLinesOn            = False            ; do not draw contour lines<br>  <br><br><br>  res@tiMainString = &quot;&quot;<br>  res@pmTickMarkDisplayMode  = &quot;Always&quot;<br>  res@mpProjection  = &quot;CylindricalEquidistant&quot;    ;The default<br>  res@mpDataBaseVersion      = &quot;MediumRes&quot;<br>  res@mpOutlineOn            =True<br>  res@lbOrientation          = &quot;Vertical&quot;<br>  res@tiMainOffsetYF         = -0.03<br>  res@mpFillOn     = False<br>  res@mpOutlineOn  = True                  ; turn the map outline on<br>  res@mpMinLatF     = 37.85<br>  res@mpMaxLatF     = 38.50<br>  res@mpMinLonF     = -120.0<br>  res@mpMaxLonF     = -119.35<br>  res@gsnLeftString = &quot;&quot;<br>  res@gsnCenterString = &quot;Geopotential Height at 600 hpa on Jan 11 at 22:40UTC&quot;<br>  res@gsnStringFontHeightF = 0.020<br>  res@gsnRightString = &quot;&quot;<br>  res@gsnMaximize     = True<br>  res@mpShapeMode     = &quot;FreeAspect&quot;  <br>  res@lbTitleString   = &quot;Terrain (m)&quot;<br>  res@lbTitlePosition = &quot;Right&quot;<br>  res@lbTitleDirection = &quot;Across&quot;<br>  res@lbTitleAngleF    = 90.<br>  res@lbTitleFontHeightF = 0.015<br> <br><br> ;------wind vectors<br>  res2 = True<br>  res2@gsnDraw = False<br>  res2@gsnFrame = False<br>  res2@vcWindBarbLineThicknessF= 3.0<br>  res2@vcRefLengthF= 0.018<br>  res2@vcRefMagnitudeF= 10<br>  res2@vcMinDistanceF = 0.05<br>  res2@vcGlyphStyle = &quot;WindBarb&quot;<br>  res2@gsnLeftString = &quot;&quot;<br>  res2@gsnRightString = &quot;&quot;<br>  res2@vcRefAnnoOn = False<br>  <br><br>;-------Frozen<br>;res3 = True<br>;res3@cnLineColor = &quot;Blue&quot;<br>;res3@gsnContourNegLineDashPattern = 2.0<br>;res3@gsnContourLineThicknessesScale = 2.0<br>;res3@gsnContourPosLineDashPattern = 3.0<br>;res3@gsnContourZeroLineThicknessF = 0.0<br>;solid = gsn_csm_contour(wks,froz,res3)<br><br>;------Melted<br>;res4 = True<br>;res4@cnLineColor = &quot;Black&quot;<br>;res4@gsnContourLineThicknessesScale = 2.0<br>;res4@cnLineLabelBackgroundColor = -1 ; transparent<br>;res4@cnInfoLabelOn = False<br>;res4@cnLevelSelectionMode = &quot;ExplicitLevels&quot;<br>;res4@cnLevels = (/0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0/)<br>;res4@cnLineLabelFontHeightF = 0.01<br>;res4@cnLineLabelDensityF =2 <br>;res4@cnLineLabelsOn = True<br>;liquid = gsn_csm_contour(wks,melt,res4)<br><br>;------Geopotential Height<br>res5 = True<br>res5@cnLineColor = &quot;Black&quot;<br>res5@gsnContourLineThicknessesScale = 3.0<br>geo = gsn_csm_contour(wks, z_plane, res5)<br><br><br>contour = gsn_csm_contour_map(wks,hgt,res)<br>vector = gsn_csm_vector(wks,u_wind,v_wind,res2)<br>overlay(contour, geo)<br>overlay(contour, vector)<br><br>plres = True<br>plres@gsLineColor = &quot;blue&quot;<br>plres@gsLineThicknessF = 3.0<br><br>lat = (/38.05, 38.05/)<br>lon = (/-120.10, -119.35/)<br><br>a = gsn_add_polyline(wks,contour,lon,lat,plres)<br><br>pltres = True<br>pltres@gsLineColor = &quot;Black&quot;<br>pltres@gsLineThicknessF = 3.0<br><br>lat1 = (/38.27, 38.27/)<br>lon1 = (/-120.10, -119.35/)<br><br>b = gsn_add_polyline(wks,contour,lon1,lat1,pltres)<br><br>draw(contour)<br>frame(wks)<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br><br><br>end<br>    <br><br></div>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>