<div dir="ltr">Thank you Mary<br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-05-12 16:07 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 Appo,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The first thing I tell people when they are trying to overlay plots and not seeing both plots is to draw them individually before you call overlay. This way you can make sure the individual plots look okay before you get to the overlay call.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">All this involves in your case is commenting out the gsnDraw and gsnFrame settings.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">When you do this, you'll see that the Y axis on the wind profile ranges from 20000 to 100000, while on the humidity plot it ranges from 200 to 1000.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It's not well-advertised, but under the hood, gsn_csm_pres_hgt converts the pressure axis to mb before plotting. (I've fixed the documentation to make this more clear.)</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">gsn_csm_vector doesn't do this, so you will need to scale the pressure values yourself before plotting.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I suggest scaling all the arrays, just to be consistent, even for r00:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default"><font face="monospace, monospace">;---Convert levels from Pa to hPa for consistent plotting </font></div><div class="gmail_default"><font face="monospace, monospace"> r00&plev = r00&plev * 0.01</font></div><div class="gmail_default"><font face="monospace, monospace"> udiv00&plev = udiv00&plev * 0.01</font></div><div class="gmail_default"><font face="monospace, monospace"> omega00&plev = omega00&plev * 0.01</font></div><div class="gmail_default"><font face="monospace, monospace"> r00&plev@units = "mb"</font></div><div class="gmail_default"><font face="monospace, monospace"> udiv00&plev@units = "mb"</font></div><div class="gmail_default"><font face="monospace, monospace"> omega00&plev@units = "mb"</font></div></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 Wed, May 10, 2017 at 3:02 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><div><div><div><div>Dear Mary,<br><br></div>I am very surprise. <br></div>I try the script with another dataset. <br></div>As you can see in the attached plot, only humidity is displayed.<br></div>Vertical profile of wind is not appearing. <br><br></div>What i going wrong?<br><br></div>Best regards,<br><br></div>Appo<div><div class="m_-5904197720154203962h5"><br><div><div><div><div><div><div><div><div><div><br><br><br><div class="gmail_extra"><br><div class="gmail_quote">2017-05-08 20:22 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">Thank you very much Mary,<div>Now, I have a perfect plot.</div><div>Cheers,</div><div>Appo</div></div><div class="m_-5904197720154203962m_5167464244441849060HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-05-07 16:08 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 style="font-size:small">Hi,</div><div style="font-size:small"><br></div><div style="font-size:small">If by "optimized for this plot" you mean is the plotting being done correctly, then I would say "yes".</div><div style="font-size:small"><br></div><div style="font-size:small">However, I did notice that you have two slightly different scales on your color bar. The top one goes from 30 to 90 in steps of 5, while the bottom one goes from 20 to 90 in steps of 5.</div><div style="font-size:small"><br></div><div style="font-size:small">If somebody is not paying attention when looking at these images, they may not realize that the bottom purple color represents values < 30 in the top plot, while it represents values < 20 in the bottom plot.</div><div style="font-size:small"><br></div><div style="font-size:small">I recommend using the same color scale for both plots, and then using a single panel labelbar so that it's more clear what the colors mean.</div><div style="font-size:small"><br></div><div style="font-size:small">To make sure you use the same contour levels for both plots, set these resources:</div><div><pre style="font-size:small;color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap"><font face="monospace, monospace"> res@cnLevelSelectionMode = "ManualLevels" ; manually set contour levels
res@cnMinLevelValF = 20. ; minimum contour level
res@cnMaxLevelValF = 90. ; maximum contour level
</font>res@cnLevelSpacingF = 5. </pre><pre style="font-size:small;color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap"><font face="arial, helvetica, sans-serif">You then you need to turn off the individual labelbars for both plots, so you can set a common labelbar in the panel:</font></pre><pre style="font-size:small;color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap"> res@lbLabelBarOn = False</pre><pre style="font-size:small;color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap"><font face="arial, helvetica, sans-serif">Now, in the panel resources, set:</font></pre><pre style="font-size:small;color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap">pres@gsnPanelLabelBar = True
pres2@lbOrientation = "vertical" ; set this if you want a vertical label bar<div><br></div></pre><pre style="font-size:small;color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap"><font face="arial, helvetica, sans-serif">If you need more customization of your panel plot, then we have some examples on our panel page:</font></pre><pre style="word-wrap:break-word"><font face="arial, helvetica, sans-serif"><font color="#000000"><span style="white-space:pre-wrap"><a href="http://www.ncl.ucar.edu/Applications/panel.shtml" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/panel.shtml</a><span class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627HOEnZb"><font color="#888888"><br></font></span></span></font></font></pre><span class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627HOEnZb"><font color="#888888"><pre style="word-wrap:break-word"><font face="arial, helvetica, sans-serif"><font color="#000000"><span style="white-space:pre-wrap">--Mary</span></font></font></pre><pre style="word-wrap:break-word"><font face="arial, helvetica, sans-serif"><font color="#000000"><span style="white-space:pre-wrap"><br></span></font></font></pre><pre style="font-size:small;color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap"><font face="arial, helvetica, sans-serif"><br></font></pre></font></span></div></div><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627h5"><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 5, 2017 at 4:07 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><div><div><span style="color:rgb(39,78,19)"><span style="background-color:rgb(255,255,255)"><font size="4">Dear Mary,<br></font></span></span></div><span style="color:rgb(39,78,19)"><span style="background-color:rgb(255,255,255)"><font size="4">I applied your suggestions. <br></font></span></span></div><span style="color:rgb(39,78,19)"><span style="background-color:rgb(255,255,255)"><font size="4">Now, things are going like a charm. <br></font></span></span></div><span style="color:rgb(39,78,19)"><span style="background-color:rgb(255,255,255)"><font size="4">Here attached figure obtained and script.<br></font></span></span></div><span style="color:rgb(39,78,19)"><span style="background-color:rgb(255,255,255)"><font size="4">I guess that this script is optimized for this plot?<br><br></font></span></span></div><span style="color:rgb(7,55,99)"><span style="color:rgb(39,78,19)"><span style="background-color:rgb(255,255,255)"><font size="4">Thank you very much</font></span></span><br></span><span><br>begin<br><br> in1 = addfile("<a href="http://uvwr.mean.JJAS.nc" target="_blank">uvwr.mean.JJAS.nc</a>", "r") ; open netcdf file<br></span> input4 = addfile("<a href="http://topo.nc" target="_blank">topo.nc</a>", "r")<span><br><br> u00 = in1->u <br> v00 = in1->v <wbr> <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></span><span> latmin = -5.0<br> latmax = 3.0<br> lonmin = min(lon)<br> lonmax = 20.0<br> <br></span> sellat2 = -1.0<br> latmin2 = -4.0<br> latmax2 = -1.0<span><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<wbr>}), 2)<br><br></span> udiv002 = dim_avg_n_Wrap(ucomp_div00(:, :, {latmin2:latmax2}, {lonmin:lonmax}), 2)<br> omega002 = dim_avg_n_Wrap(w00(:, :, {latmin2:latmax2}, {lonmin:lonmax}), 2)<br> r002 = dim_avg_n_Wrap(rhum00(:, :, {latmin2:latmax2},{lonmin:lonm<wbr>ax}), 2)<span><br><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:la<wbr>tmax},{lonmin:lonmax}), 0)<br> <br> topo = elev({sellat}, {lonmin:lonmax})<br> topo2 = elev({sellat2}, {lonmin:lonmax})<br> <br>;---create plot<br></span> wks = gsn_open_wks ("eps", "Overlay_Rhum_udiv_JJAS_Timmea<wbr>n_5S3N") ; open workstation<span><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><br><br></span> contour_plota1 = gsn_csm_pres_hgt(wks, r00(0, :, {lonmin:lonmax}), cnres)<br> vector_plota1 = gsn_csm_vector(wks, udiv00(0, :, {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )<br> <br> contour_plota2 = gsn_csm_pres_hgt(wks, r002(0, :, {lonmin:lonmax}), cnres)<br> vector_plota2 = gsn_csm_vector(wks, udiv002(0, :, {lonmin:lonmax}), -200.0*omega002(0, :, {lonmin:lonmax}), vcres )<span><br><br><br>;---Add topo field using a filled polygon.<br></span> getvalues contour_plota1<span><br> "trYMinF" : ymin<br> "trYMaxF" : ymax<br> end getvalues<br><br></span><span>;---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> <br> <br></span><span>;---Add topo field using a filled polygon.<br></span> getvalues contour_plota2<span><br> "trYMinF" : ymin2<br> "trYMaxF" : ymax2<br> end getvalues<br><br>;---Create new X,Y arrays that form a closed polygon.<br></span> nlon = dimsizes(topo2&lon)<br> xtopo2 = new(nlon+3,typeof(topo2&lon))<br> ytopo2 = new(nlon+3,typeof(topo2))<span><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>;---Add the polygon to the contour plot.<br> gnres = True<br> gnres@gsFillColor = "gray20"<br> <br></span> id = gsn_add_polygon(wks, contour_plota1, xtopo, ytopo, gnres)<br> i2 = gsn_add_polygon(wks, contour_plota2, xtopo2, ytopo2, gnres)<br> <br> overlay(contour_plota1, vector_plota1)<br> overlay(contour_plota2, vector_plota2)<br> <br> plot(0) = contour_plota1 <br> plot(1) = contour_plota2 <br><br> ;draw(contour_plot) ; This draws everything<span><br> <br> <br> <br> <br> pres = True<br> pres@gsnMaximize = True<br> gsn_panel(wks, plot, (/2, 1/), pres)<br> <br> <br><br></span>end <br><div><div><br><br></div></div></div><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-05-05 0:39 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 style="font-size:small">I'm so sorry I didn't get back to you sooner on this. I've been swamped and then was out-of-town for a few days.</div><div style="font-size:small"><br></div><div style="font-size:small">It's a little hard to follow what's going on in this script without being able to run it. </div><div style="font-size:small"><br></div><div style="font-size:small">I think the issue is that you are using the "res" variable in your call to gsn_panel, and "res" has both gsnDraw and gsnFrame set to False. This will cause *no* plot to appear when you call gsn_panel. I'm surprised you are getting any plot at all.</div><div style="font-size:small"><br></div><div style="font-size:small">I recommend creating a new resource variable for gsn_panel, because you really don't want to use the one that you used for the plots themselves:</div><div style="font-size:small"><br></div><div style="font-size:small"><font face="monospace, monospace">pres = True</font></div><div style="font-size:small"><font face="monospace, monospace">pres@gsnMaximize = True</font></div><div style="font-size:small"><font face="monospace, monospace">gsn_panel(wks, plot, (/2, 1/), pres)</font></div><div style="font-size:small"><font face="monospace, monospace"><br></font></div>If you<div style="font-size:small;display:inline"> continue to have problems, then it would help if you send me the two data files. If you can't send the files, then send me the image you're getting, and any error messages that the script produces.</div><div><br></div><div><div style="font-size:small">Thanks,</div><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div><div style="font-size:small"></div><br><div style="font-size:small"><br></div></div></div><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 25, 2017 at 12:56 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><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" target="_blank">uvwr.mean.JJAS.nc</a>", "r") ; open netcdf file<br> input4 = addfile("<a href="http://atanas_topo.nc" target="_blank">atanas_topo.nc</a>", "r")<br><br> u00 = in1->u <br> v00 = in1->v <wbr> <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<wbr>}), 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:la<wbr>tmax},{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") <wbr> ; 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="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156h5"><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 style="font-size:small">Hi,</div><div style="font-size:small"><br></div><div style="font-size:small">In the future, please email ncl-talk with follow-up questions so everybody can benefit from the answers.</div><div style="font-size:small"><br></div><div style="font-size:small">What is the exact error message you're getting?</div><span class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176HOEnZb"><font color="#888888"><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div><div style="font-size:small"><br></div></font></span></div><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176h5"><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_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_1699758146122632547HOEnZb"><font color="#888888"><br><br></font></span></div><span class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_1699758146122632547HOEnZb"><font color="#888888">Appo<br></font></span></div><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_1699758146122632547HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_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_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_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 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 style="font-size:small"><br></div><div 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 style="font-size:small"><br></div><div 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 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 style="font-size:small"><br></div><div 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_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634m_-1832595038336345784HOEnZb"><font color="#888888"><div style="font-size:small"><br></div><div style="font-size:small">--Mary</div><div style="font-size:small"><br></div></font></span></div><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634m_-1832595038336345784HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_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_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_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_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291m_-1311893959559757138HOEnZb"><div class="m_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_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 style="font-size:small">I don't think I've ever seen this error before.</div><div style="font-size:small"><br></div><div 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 style="font-size:small"><br></div><div> plot = gsn_csm_pres_hgt_vector(wks, r00(0, :, {lonmin:lonmax}), udiv00(0, :, {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), res )</div><div><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_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_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_-5904197720154203962m_5167464244441849060m_-2035199396315559627m_-6467908691799575580m_-5509234554777509720m_5105003629348599156m_-6385272546512449176m_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>
</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></div></div></div></div></div></div></div></div></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>