<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'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'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("<a href="http://uvwr.mean.JJAS.nc">uvwr.mean.JJAS.nc</a>", "r") ; open netcdf file<br> input4 = addfile("<a href="http://atanas_topo.nc">atanas_topo.nc</a>", "r")<br><br> u00 = in1->u <br> v00 = in1->v <br> w00 = in1->w<br> rhum00 = in1->r<br> <br> topo1 = input4->HT<br> <br> <br> lat = in1->latitude<br> lon = in1->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 ("eps", "Overlay_Rhum_udiv_topo") ; 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 = ""<br> res@gsnLeftString = ""<br><br> cnres = res<br> cnres@cnFillOn = True ; turn on color fill<br> cnres@cnFillPalette = "MPL_rainbow"<br> cnres@lbOrientation = "Vertical"<br> cnres@pmLabelBarOrthogonalPosF = 0.08<br> cnres@tiYAxisString = "Pressure (hPa)"<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 = "CurlyVector" ; turn on curley vectors<br> vcres@vcMinDistanceF = 0.04 ; thin out vectors<br> vcres@vcMapDirection = False<br> vcres@vcLineArrowThicknessF = 3.0 <br> vcres@vcVectorDrawOrder = "Draw" ; 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> "trYMinF" : ymin<br> "trYMaxF" : ymax<br> end getvalues<br><br> <br>;---Add topo field using a filled polygon.<br> getvalues contour_plot_2<br> "trYMinF" : ymin2<br> "trYMaxF" : ymax2<br> end getvalues<br> <br> <br> <br>;---Create new X,Y arrays that form a closed polygon.<br> nlon = dimsizes(topo&lon)<br> xtopo = new(nlon+3,typeof(topo&lon))<br> ytopo = new(nlon+3,typeof(topo))<br><br> xtopo(0:nlon-1) = topo&lon<br> ytopo(0:nlon-1) = topo<br> xtopo(nlon) = topo&lon(nlon-1)<br> ytopo(nlon) = ymax ; Use actual Y max of contour plot<br> xtopo(nlon+1) = topo&lon(0)<br> ytopo(nlon+1) = ymax<br> xtopo(nlon+2) = topo&lon(0) ; This last point closes<br> ytopo(nlon+2) = topo(0) ; the polygon.<br> <br> <br> nlon2 = dimsizes(topo2&lon)<br> xtopo2 = new(nlon2+3,typeof(topo2&lon))<br> ytopo2 = new(nlon2+3,typeof(topo2))<br><br> xtopo2(0:nlon-1) = topo2&lon<br> ytopo2(0:nlon-1) = topo2<br> xtopo2(nlon) = topo2&lon(nlon-1)<br> ytopo2(nlon) = ymax2 ; Use actual Y max of contour plot<br> xtopo2(nlon+1) = topo2&lon(0)<br> ytopo2(nlon+1) = ymax2<br> xtopo2(nlon+2) = topo2&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 = "gray20"<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"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></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'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"><<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>></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"><<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>></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's look like what I'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"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></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'm sorry I didn'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'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'm not sure what I'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'm using "sellat" 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"><<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>></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"><<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>></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'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"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></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't think I'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's hard to debug this particular script without having the data so I can run it. But, I'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'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 "udiv00" instead.</div><div style="font-size:small"><br></div><div style="font-size:small">Still, it'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'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><use your email address for the password><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"><<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>></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'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>