<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><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://topo.nc">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> latmin = -5.0<br> latmax = 3.0<br> lonmin = min(lon)<br> lonmax = 20.0<br> <br> sellat2 = -1.0<br> latmin2 = -4.0<br> latmax2 = -1.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> 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:lonmax}), 2)<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: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_JJAS_Timmean_5S3N") ; 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><br><br> 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 )<br><br><br>;---Add topo field using a filled polygon.<br> getvalues contour_plota1<br> "trYMinF" : ymin<br> "trYMaxF" : ymax<br> end getvalues<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> <br> <br>;---Add topo field using a filled polygon.<br> getvalues contour_plota2<br> "trYMinF" : ymin2<br> "trYMaxF" : ymax2<br> end getvalues<br><br>;---Create new X,Y arrays that form a closed polygon.<br> nlon = dimsizes(topo2&lon)<br> xtopo2 = new(nlon+3,typeof(topo2&lon))<br> ytopo2 = new(nlon+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>;---Add the polygon to the contour plot.<br> gnres = True<br> gnres@gsFillColor = "gray20"<br> <br> 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<br> <br> <br> <br> <br> pres = True<br> pres@gsnMaximize = True<br> gsn_panel(wks, plot, (/2, 1/), pres)<br> <br> <br><br>end <br><div><div><br><br></div></div></div><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 class="gmail_default" 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 class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" 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 class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" 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 class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" 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 class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">pres = True</font></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">pres@gsnMaximize = True</font></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace">gsn_panel(wks, plot, (/2, 1/), pres)</font></div><div class="gmail_default" style="font-size:small"><font face="monospace, monospace"><br></font></div>If you<div class="gmail_default" 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 class="gmail_default" style="font-size:small">Thanks,</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"></div><br><div class="gmail_default" style="font-size:small"><br></div></div></div><div class="HOEnZb"><div class="h5"><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_5105003629348599156HOEnZb"><div class="m_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 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="m_5105003629348599156m_-6385272546512449176HOEnZb"><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_5105003629348599156m_-6385272546512449176HOEnZb"><div class="m_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_5105003629348599156m_-6385272546512449176m_1699758146122632547HOEnZb"><font color="#888888"><br><br></font></span></div><span class="m_5105003629348599156m_-6385272546512449176m_1699758146122632547HOEnZb"><font color="#888888">Appo<br></font></span></div><div class="m_5105003629348599156m_-6385272546512449176m_1699758146122632547HOEnZb"><div class="m_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_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634HOEnZb"><div class="m_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 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_5105003629348599156m_-6385272546512449176m_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_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634m_-1832595038336345784HOEnZb"><div class="m_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_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291HOEnZb"><div class="m_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_5105003629348599156m_-6385272546512449176m_1699758146122632547m_-8427646788359640634m_-1832595038336345784m_2127398401502690291m_-1311893959559757138HOEnZb"><div class="m_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 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_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_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>