<div dir="ltr"><div><div><div><div>Good morning Mary,<br></div>Sorry , i forgot to add ncl-users email to the reply. <br><br></div>I am attaching the modified version of the script were I&#39;m trying to create a panel plot for two transects at latitude sellat and sellat2. <br></div>Unfortunately, only one plot appears. <br></div><div>Really, I don&#39;t know what is wrong this script.<br><br></div>Your help will be appreciated.<br>  <br><div><div><div><div><div><br><br>begin<br><br>    in1  = addfile(&quot;<a href="http://uvwr.mean.JJAS.nc">uvwr.mean.JJAS.nc</a>&quot;, &quot;r&quot;)                         ; open netcdf file<br>    input4  = addfile(&quot;<a href="http://atanas_topo.nc">atanas_topo.nc</a>&quot;, &quot;r&quot;)<br><br>   u00     =  in1-&gt;u  <br>   v00     =  in1-&gt;v                              <br>   w00     =  in1-&gt;w<br>   rhum00  =  in1-&gt;r<br>   <br>   topo1 = input4-&gt;HT<br>    <br>   <br>   lat = in1-&gt;latitude<br>   lon = in1-&gt;longitude<br>   <br> <br>   dv00    = uv2dvF (u00, v00)    <br>   uvd00   = dv2uvF (dv00)  <br>   ucomp_div00 = uvd00(0, :, : , :, :)<br>   copy_VarCoords(u00, ucomp_div00)<br><br><br>   sellat = 3.0<br>   sellat2 = -2.5<br>   latmin = -5.0<br>   latmax = 3.0<br>   lonmin = min(lon)<br>   lonmax = 20.0<br>   <br>   udiv00 = dim_avg_n_Wrap(ucomp_div00(:, :, {latmin:latmax}, {lonmin:lonmax}), 2)<br>   omega00 = dim_avg_n_Wrap(w00(:, :, {latmin:latmax}, {lonmin:lonmax}), 2)<br>   r00 = dim_avg_n_Wrap(rhum00(:, :, {latmin:latmax},{lonmin:lonmax}), 2)<br><br>   ;udiv00  = ucomp_div00(:, :, {sellat}, {lonmin:lonmax})<br>   ;omega00 = w00(:, :, {sellat}, {lonmin:lonmax})<br>   ;r00     = rhum00(:, :, {sellat},{lonmin:lonmax})<br>  <br> <br>   elev = 1013.25*(1 - topo1*0.0065/288.15)^5.25145<br>  <br>   copy_VarCoords(topo1, elev)<br><br>;  topo = dim_avg_n_Wrap(elev({latmin:latmax},{lonmin:lonmax}), 0)<br>  <br>   topo  = elev({sellat}, {lonmin:lonmax})<br>   topo2 = elev({sellat2}, {lonmin:lonmax})<br> <br>;---create plot<br>   wks   = gsn_open_wks (&quot;eps&quot;, &quot;Overlay_Rhum_udiv_topo&quot;)        ; open workstation<br>    plot = new(2, graphic)     <br>   res                     = True                  ; plot mods desired<br>   res@gsnDraw             = False    ; turn off draw<br>   res@gsnFrame            = False    ; turn off frame<br>   res@gsnMaximize         = True<br><br>   res@gsnRightString      = &quot;&quot;<br>   res@gsnLeftString       = &quot;&quot;<br><br>   cnres = res<br>   cnres@cnFillOn          = True               ; turn on color fill<br>   cnres@cnFillPalette     = &quot;MPL_rainbow&quot;<br>   cnres@lbOrientation     = &quot;Vertical&quot;<br>   cnres@pmLabelBarOrthogonalPosF = 0.08<br>   cnres@tiYAxisString      = &quot;Pressure (hPa)&quot;<br><br>   vcres = res<br>   vcres@vcRefMagnitudeF = 10.0                ; define vector ref mag<br>   vcres@vcRefLengthF    = 0.1              ; define length of vec ref<br>   vcres@vcGlyphStyle    = &quot;CurlyVector&quot;      ; turn on curley vectors<br>   vcres@vcMinDistanceF  = 0.04               ; thin out vectors<br>   vcres@vcMapDirection  = False<br>   vcres@vcLineArrowThicknessF   = 3.0  <br>   vcres@vcVectorDrawOrder  = &quot;Draw&quot;        ; draw vectors last<br>  <br>   contour_plot_1 = gsn_csm_pres_hgt(wks, r00(0, :, {lonmin:lonmax}), cnres)<br>   vector_plot_1  = gsn_csm_vector(wks, udiv00(0, :, {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )<br>   <br>   contour_plot_2 = gsn_csm_pres_hgt(wks, r00(0, :, {lonmin:lonmax}), cnres)<br>   vector_plot_2  = gsn_csm_vector(wks, udiv00(0, :, {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )<br><br><br>;---Add topo field using a filled polygon.<br>   getvalues contour_plot_1<br>      &quot;trYMinF&quot; : ymin<br>      &quot;trYMaxF&quot; : ymax<br>   end getvalues<br><br> <br>;---Add topo field using a filled polygon.<br>   getvalues contour_plot_2<br>      &quot;trYMinF&quot; : ymin2<br>      &quot;trYMaxF&quot; : ymax2<br>   end getvalues<br>  <br>   <br>   <br>;---Create new X,Y arrays that form a closed polygon.<br>   nlon  = dimsizes(topo&amp;lon)<br>   xtopo = new(nlon+3,typeof(topo&amp;lon))<br>   ytopo = new(nlon+3,typeof(topo))<br><br>   xtopo(0:nlon-1) = topo&amp;lon<br>   ytopo(0:nlon-1) = topo<br>   xtopo(nlon)     = topo&amp;lon(nlon-1)<br>   ytopo(nlon)     = ymax          ; Use actual Y max of contour plot<br>   xtopo(nlon+1)   = topo&amp;lon(0)<br>   ytopo(nlon+1)   = ymax<br>   xtopo(nlon+2)   = topo&amp;lon(0)   ; This last point closes<br>   ytopo(nlon+2)   = topo(0)       ; the polygon.<br>   <br>   <br>   nlon2  = dimsizes(topo2&amp;lon)<br>   xtopo2 = new(nlon2+3,typeof(topo2&amp;lon))<br>   ytopo2 = new(nlon2+3,typeof(topo2))<br><br>   xtopo2(0:nlon-1) = topo2&amp;lon<br>   ytopo2(0:nlon-1) = topo2<br>   xtopo2(nlon)     = topo2&amp;lon(nlon-1)<br>   ytopo2(nlon)     = ymax2         ; Use actual Y max of contour plot<br>   xtopo2(nlon+1)   = topo2&amp;lon(0)<br>   ytopo2(nlon+1)   = ymax2<br>   xtopo2(nlon+2)   = topo2&amp;lon(0)   ; This last point closes<br>   ytopo2(nlon+2)   = topo2(0)       ; the polygon.<br>   <br>  <br>;---Add the polygon to the contour plot.<br>   gnres = True<br>   gnres@gsFillColor = &quot;gray20&quot;<br>  <br>   id = gsn_add_polygon(wks, contour_plot_1, xtopo, ytopo, gnres)<br>   id2 = gsn_add_polygon(wks, contour_plot_2, xtopo2, ytopo2, gnres)<br>   overlay(contour_plot_1, vector_plot_1)<br>   overlay(contour_plot_2, vector_plot_2)<br>   <br><br>     draw(contour_plot_1)     ; This draws everything<br>     draw(contour_plot_2)     ; This draws everything<br>   <br>   plot(0) = contour_plot_1<br>   <br>   plot(1) = contour_plot_2<br>   <br>    gsn_panel(wks,plot, (/2, 1/), res)               ; now draw as one plot<br> <br><br>end <br><br></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-04-24 21:58 GMT+01:00 Mary Haley <span dir="ltr">&lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@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 class="gmail_default" style="font-size:small">Hi,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">In the future, please email ncl-talk with follow-up questions so everybody can benefit from the answers.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">What is the exact error message you&#39;re getting?</div><span class="HOEnZb"><font color="#888888"><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></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 24, 2017 at 4:08 AM, Appo derbetini <span dir="ltr">&lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@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>Hi Mary,<br><br></div>Since you send me the script that plot masked fields, I am working to do it for many latitude in a panel plot. <br></div><div>Always i got many errors saying that i cannot overlay.<br></div><div>How to solve it?<br></div>Thank you in advance.<span class="m_1699758146122632547HOEnZb"><font color="#888888"><br><br></font></span></div><span class="m_1699758146122632547HOEnZb"><font color="#888888">Appo<br></font></span></div><div class="m_1699758146122632547HOEnZb"><div class="m_1699758146122632547h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-04-10 8:23 GMT+01:00 Appo derbetini <span dir="ltr">&lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</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>Thank you very much.<br></div><div>It&#39;s look like what I&#39;m trying to do.<br></div><div>regards<br></div></div><div class="m_1699758146122632547m_-8427646788359640634HOEnZb"><div class="m_1699758146122632547m_-8427646788359640634h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-04-09 22:11 GMT+01:00 Mary Haley <span dir="ltr">&lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@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 class="gmail_default" style="font-size:small">I&#39;m sorry I didn&#39;t get back to this sooner, but I was having trouble trying to understand what you are doing in the script.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">First, I&#39;m not sure it makes sense to take an average of topographical data across a set of latitudes and plot that as a topo layer.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">But, that aside, I think you just want to plot the topo line as a filled polygon instead of a separate XY plot.  <br><br></div><div class="gmail_default" style="font-size:small">I&#39;m not sure what I&#39;ve attached is correct, but hopefully it gives you an idea of what to do.  </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Note that I&#39;m using &quot;sellat&quot; to select a single latitude value, rather than averaging across a range of latitudes.</div><span class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784HOEnZb"><font color="#888888"><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></font></span></div><div class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784HOEnZb"><div class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784h5"><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Apr 9, 2017 at 2:10 AM, Appo derbetini <span dir="ltr">&lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@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>Dear Mary,<br><br></div>Despite applying what you suggest the problem remain the same<br></div>Thank you very much<br></div><div class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291HOEnZb"><div class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-04-06 8:53 GMT+01:00 Appo derbetini <span dir="ltr">&lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</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>Dear Mary,<br><br></div>Datasets are uploaded as indicated. <br><br></div>What I really trying to do is something like Mask Example 14. <br><br></div>But I don&#39;t want to interpolate topography to grid of wind. <br><br></div>Here attached figure produced by the script<br><br></div>Thank you<br><div><div><div><br></div></div></div></div><div class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291m_-1311893959559757138HOEnZb"><div class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291m_-1311893959559757138h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-04-05 15:36 GMT+01:00 Mary Haley <span dir="ltr">&lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@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 class="gmail_default" style="font-size:small">I don&#39;t think I&#39;ve ever seen this error before.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It&#39;s hard to debug this particular script without having the data so I can run it. But, I&#39;m wondering if this line might be part of the problem:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default">  plot  = gsn_csm_pres_hgt_vector(wks, r00(0, :, {lonmin:lonmax}), udiv00(0, :, {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), res )</div><div class="gmail_default"><br><div style="font-size:small">When you do an arithmetic operation on an array, like:</div><div style="font-size:small"><br></div><div style="font-size:small">  -200. * omega00</div><div style="font-size:small"><br></div><div style="font-size:small">you end up stripping all of the metadata off omega00.  What&#39;s getting passed into gsn_csm_pres_hgt has no metadata. This might be okay, however, because I think this routine might just use the metadata from &quot;udiv00&quot; instead.</div><div style="font-size:small"><br></div><div style="font-size:small">Still, it&#39;s worth trying this:</div><div style="font-size:small"><br></div><div style="font-size:small"> omega00_scale = omega00(0, :, {lonmin:lonmax})   ; subsets omega00 *and* copies all metadata</div><div style="font-size:small"> omega00_scale = -200.0 * omega00_scale</div><div style="font-size:small"><br></div><div>  plot  = gsn_csm_pres_hgt_vector(wks, r00(0, :, {lonmin:lonmax}), udiv00(0, :, {lonmin:lonmax}), omega00_scale, res )<br></div><div><br></div><div>If you continue to have problems with this script, could you upload your data (if it&#39;s not too large) to our ftp:</div><div><br></div><div><font face="monospace, monospace">ftp <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a><br>anonymous<br>&lt;use your email address for the password&gt;<br>cd incoming<br>put <a href="http://uvwr.mean.JJAS.nc" target="_blank">uvwr.mean.JJAS.nc</a><br>put <a href="http://atanas_topo.nc" target="_blank">atanas_topo.nc</a><br></font><div><font face="monospace, monospace">quit</font></div></div><div><br></div><div>Thanks,</div><div><br></div><div>--Mary</div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291m_-1311893959559757138m_1425169837060340304h5">On Tue, Apr 4, 2017 at 9:18 AM, Appo derbetini <span dir="ltr">&lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@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 class="m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291m_-1311893959559757138m_1425169837060340304h5"><div dir="ltr"><div><div><div><div><div>Dear all,<br></div>I am trying to overlay topography profil on a profil of wind and relative humidity.<br></div>I don&#39;t want to interpolate topography to show some details. <br><br></div><div>Atached are script and datasets used for the plot.<br></div><div><br></div>But I having errors. <br><br>warning:LlDataPolygon: point 6.150009,0.000000 outside data domain<br>warning:LlDataPolygon: point 19.983337,0.000000 outside data domain<br>warning:LlDataLineTo: point 6.150009,0.000000 outside data domain<br>warning:LlDataLineTo: point 19.983337,0.000000 outside data domain<br><br></div>Any help will be appreciated.<br><br></div>Best regards.</div>
<br></div></div>______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>