<div dir="ltr"><div><div>Thank you Adam.<br><br></div><br></div>Dipti<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, May 16, 2017 at 3:08 PM, Adam Phillips <span dir="ltr">&lt;<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Dipti,<div>Yes. You would need to set a couple of resources that control the range of the x-axis:</div><div>r@trXMinF = 0.0</div><div>r@trXMaxF = 0.001.  </div><div><br></div><div>If that does not fix the issue please let the ncl-talk email list know.</div><div>Adam </div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Tue, May 16, 2017 at 10:39 AM, Dipti Sharma <span dir="ltr">&lt;<a href="mailto:sdipti596@gmail.com" target="_blank">sdipti596@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">All,<br><br><pre>The below script works fine but I can not change the x-axix range/y-axis range (values). <br>it is 0. to 0.5 now but I want to make it going from 0 to 0.001 for instance. Is it possible to modify it? <br></pre><pre>Thank you in advance for your help on this.<br><br></pre><pre>Dipti <br></pre><pre>;Here is the script<br><br>;spec_4.ncl
;
; Concepts illustrated:
;   - Calculating and plotting spectra
;   - Calculating confidence intervals
;   - Making an axis logarithmic in an XY plot
;*****************************<wbr>*******************
load &quot;$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_code.ncl&quot;
load &quot;$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_csm.ncl&quot;
load &quot;$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/contributed.ncl&quot;
load &quot;$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/shea_util.ncl&quot;
;*****************************<wbr>*******************
begin
;*****************************<wbr>*******************
; variable and file handling
;*****************************<wbr>*******************
   fn  = &quot;SOI_Darwin.nc&quot; ; define filename
   in  = addfile(fn,&quot;r&quot;)                                  ; open netcdf file
   soi  = in-&gt;DSOI                                        ; get data
;*****************************<wbr>*******************
; set function arguments
;*****************************<wbr>*******************
  d   = 0    ; detrending opt: 0=&gt;remove mean 1=&gt;remove mean + detrend
  sm  = 21   ; smooth: should be at least 3 and odd
  pct = 0.10 ; percent taper: (0.0 &lt;= pct &lt;= 1.0) 0.10 common. 
;*****************************<wbr>*******************
; calculate spectrum
;*****************************<wbr>*******************
  sdof = specx_anal(soi,d,sm,pct)
;*****************************<wbr>*******************
; calculate confidence interval [here 5 and 95%]
; return 4 curves to be plotted
;*****************************<wbr>*******************
  splt = specx_ci (sdof, 0.05, 0.95)
;*****************************<wbr>*******************
; plotting
;*****************************<wbr>*******************
   wks  = gsn_open_wks(&quot;ps&quot;,&quot;spec&quot;)              ; Opens a ps file

   r = True                                      ; plot mods desired
   r@tiMainString = &quot;SOI&quot;                        ; title
   r@tiXAxisString = &quot;Frequency (cycles/month)&quot;  ; xaxis
   r@tiYAxisString = &quot;Variance&quot;                  ; yaxis
;*****************************<wbr>******************
; Generate log plot showing &quot;red noise&quot; confidence bounds
; (a) log scaling and (b) the Band Width
;*****************************<wbr>******************
   r@trYLog              = True                 ; log scaling
   r@trYMinF             = 0.10                 ; manually set lower limit
   r@trYMaxF             = 30.0                 ;   &quot;          upper
   r@gsnFrame            = False                ; do not advance frame
   plot  = gsn_csm_xy(wks,sdof@frq, splt,r)

   xf   = (/0.40, 0.40+sdof@bw/)                ; Create band width line
   ys   = 0.75*max(sdof@spcx)                   ; 75% up Y axis
   yv   = (/ys,ys/)               
   rpl  = True                                  ; resources for polyline
   rpl@gsLineThicknessF  = 2                    ; Define line thickness 
   gsn_polyline(wks,plot,xf,yv,rp<wbr>l)             ; Draw BandWidth

   txres= True                                  ; label BW line
   txres@txFontHeightF = 0.015                  ; font height
   txres@txJust        = &quot;CenterLeft&quot;           ; Set lable location
   gsn_text(wks,plot,&quot;BW&quot;,0.41+sd<wbr>of@bw,ys,txres); Label 
   frame (wks)
end

</pre><div class="gmail_extra"><br><div class="gmail_quote">On Sat, May 13, 2017 at 2:00 PM,  <span dir="ltr">&lt;<a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send ncl-talk mailing list submissions to<br>
        <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<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>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:ncl-talk-owner@ucar.edu" target="_blank">ncl-talk-owner@ucar.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of ncl-talk digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: Overlay profile of topography (Appo derbetini)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Fri, 12 May 2017 19:17:56 +0100<br>
From: Appo derbetini &lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt;<br>
Subject: Re: [ncl-talk] Overlay profile of topography<br>
To: Mary Haley &lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;<br>
Cc: NCL &lt;<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>&gt;<br>
Message-ID:<br>
        &lt;<a href="mailto:CAHVpkHj4MrBrGt1Yrcq1RGkVaF1Fy-fG5e2x6zcg9LOvT77S9Q@mail.gmail.com" target="_blank">CAHVpkHj4MrBrGt1Yrcq1RGkVaF1F<wbr>y-fG5e2x6zcg9LOvT77S9Q@mail.gm<wbr>ail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;utf-8&quot;<br>
<br>
Thank you Mary<br>
<br>
<br>
<br>
2017-05-12 16:07 GMT+01:00 Mary Haley &lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;:<br>
<br>
&gt; Hi Appo,<br>
&gt;<br>
&gt; The first thing I tell people when they are trying to overlay plots and<br>
&gt; not seeing both plots is to draw them individually before you call<br>
&gt; overlay.  This way you can make sure the individual plots look okay before<br>
&gt; you get to the overlay call.<br>
&gt;<br>
&gt; All this involves in your case is commenting out the gsnDraw and gsnFrame<br>
&gt; settings.<br>
&gt;<br>
&gt; When you do this, you&#39;ll see that the Y axis on the wind profile ranges<br>
&gt; from 20000 to 100000, while on the humidity plot it ranges from 200 to 1000.<br>
&gt;<br>
&gt; It&#39;s not well-advertised, but under the hood, gsn_csm_pres_hgt converts<br>
&gt; the pressure axis to mb before plotting. (I&#39;ve fixed the documentation to<br>
&gt; make this more clear.)<br>
&gt;<br>
&gt; gsn_csm_vector doesn&#39;t do this, so you will need to scale the pressure<br>
&gt; values yourself before plotting.<br>
&gt;<br>
&gt; I suggest scaling all the arrays, just to be consistent, even for r00:<br>
&gt;<br>
&gt; ;---Convert levels from Pa to hPa for consistent plotting<br>
&gt;<br>
&gt;    r00&amp;plev = r00&amp;plev * 0.01<br>
&gt;    udiv00&amp;plev = udiv00&amp;plev * 0.01<br>
&gt;    omega00&amp;plev = omega00&amp;plev * 0.01<br>
&gt;    r00&amp;plev@units  = &quot;mb&quot;<br>
&gt;    udiv00&amp;plev@units  = &quot;mb&quot;<br>
&gt;    omega00&amp;plev@units = &quot;mb&quot;<br>
&gt;<br>
&gt; --Mary<br>
&gt;<br>
&gt;<br>
&gt; On Wed, May 10, 2017 at 3:02 AM, Appo derbetini &lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt;<br>
&gt; wrote:<br>
&gt;<br>
&gt;&gt; Dear Mary,<br>
&gt;&gt;<br>
&gt;&gt; I am very surprise.<br>
&gt;&gt; I try the script with another dataset.<br>
&gt;&gt; As you can see in the attached plot, only humidity is displayed.<br>
&gt;&gt; Vertical profile of wind is not appearing.<br>
&gt;&gt;<br>
&gt;&gt; What i going wrong?<br>
&gt;&gt;<br>
&gt;&gt; Best regards,<br>
&gt;&gt;<br>
&gt;&gt; Appo<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; 2017-05-08 20:22 GMT+01:00 Appo derbetini &lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt;:<br>
&gt;&gt;<br>
&gt;&gt;&gt; Thank you very much Mary,<br>
&gt;&gt;&gt; Now, I have a perfect plot.<br>
&gt;&gt;&gt; Cheers,<br>
&gt;&gt;&gt; Appo<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; 2017-05-07 16:08 GMT+01:00 Mary Haley &lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;:<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; Hi,<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; If by &quot;optimized for this plot&quot; you mean is the plotting being done<br>
&gt;&gt;&gt;&gt; correctly, then I would say &quot;yes&quot;.<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; However, I did notice that you have two slightly different scales on<br>
&gt;&gt;&gt;&gt; your color bar. The top one goes from 30 to 90 in steps of 5, while the<br>
&gt;&gt;&gt;&gt; bottom one goes from 20 to 90 in steps of 5.<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; If somebody is not paying attention when looking at these images, they<br>
&gt;&gt;&gt;&gt; may not realize that the bottom purple color represents values &lt; 30 in the<br>
&gt;&gt;&gt;&gt; top plot, while it represents values &lt; 20 in the bottom plot.<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; I recommend using the same color scale for both plots, and then using a<br>
&gt;&gt;&gt;&gt; single panel labelbar so that it&#39;s more clear what the colors mean.<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; To make sure you use the same contour levels for both plots, set these<br>
&gt;&gt;&gt;&gt; resources:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;  res@cnLevelSelectionMode = &quot;ManualLevels&quot; ; manually set contour levels<br>
&gt;&gt;&gt;&gt;  res@cnMinLevelValF       = 20.            ; minimum contour level<br>
&gt;&gt;&gt;&gt;  res@cnMaxLevelValF       = 90.            ; maximum contour level<br>
&gt;&gt;&gt;&gt;  res@cnLevelSpacingF        = 5.<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; You then you need to turn off the individual labelbars for both plots, so you can set a common labelbar in the panel:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;  res@lbLabelBarOn = False<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; Now, in the panel resources, set:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; pres@gsnPanelLabelBar = True<br>
&gt;&gt;&gt;&gt; pres2@lbOrientation   = &quot;vertical&quot; ; set this if you want a vertical label bar<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; If you need more customization of your panel plot, then we have some examples on our panel page:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; <a href="http://www.ncl.ucar.edu/Applications/panel.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/panel.shtml</a><br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; --Mary<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt; On Fri, May 5, 2017 at 4:07 AM, Appo derbetini &lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt;<br>
&gt;&gt;&gt;&gt; wrote:<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; Dear Mary,<br>
&gt;&gt;&gt;&gt;&gt; I applied your suggestions.<br>
&gt;&gt;&gt;&gt;&gt; Now, things are going like a charm.<br>
&gt;&gt;&gt;&gt;&gt; Here attached figure obtained and script.<br>
&gt;&gt;&gt;&gt;&gt; I guess that this script is optimized  for this plot?<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; Thank you very much<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; begin<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;     in1  = addfile(&quot;<a href="http://uvwr.mean.JJAS.nc" rel="noreferrer" target="_blank">uvwr.mean.JJAS.nc</a>&quot;, &quot;r&quot;)<br>
&gt;&gt;&gt;&gt;&gt; ; open netcdf file<br>
&gt;&gt;&gt;&gt;&gt;     input4  = addfile(&quot;<a href="http://topo.nc" rel="noreferrer" target="_blank">topo.nc</a>&quot;, &quot;r&quot;)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    u00     =  in1-&gt;u<br>
&gt;&gt;&gt;&gt;&gt;    v00     =  in1-&gt;v<br>
&gt;&gt;&gt;&gt;&gt;    w00     =  in1-&gt;w<br>
&gt;&gt;&gt;&gt;&gt;    rhum00  =  in1-&gt;r<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    topo1 = input4-&gt;HT<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    lat = in1-&gt;latitude<br>
&gt;&gt;&gt;&gt;&gt;    lon = in1-&gt;longitude<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    dv00    = uv2dvF (u00, v00)<br>
&gt;&gt;&gt;&gt;&gt;    uvd00   = dv2uvF (dv00)<br>
&gt;&gt;&gt;&gt;&gt;    ucomp_div00 = uvd00(0, :, : , :, :)<br>
&gt;&gt;&gt;&gt;&gt;    copy_VarCoords(u00, ucomp_div00)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    sellat = 3.0<br>
&gt;&gt;&gt;&gt;&gt;    latmin = -5.0<br>
&gt;&gt;&gt;&gt;&gt;    latmax = 3.0<br>
&gt;&gt;&gt;&gt;&gt;    lonmin = min(lon)<br>
&gt;&gt;&gt;&gt;&gt;    lonmax = 20.0<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    sellat2 = -1.0<br>
&gt;&gt;&gt;&gt;&gt;    latmin2 = -4.0<br>
&gt;&gt;&gt;&gt;&gt;    latmax2 = -1.0<br>
&gt;&gt;&gt;&gt;&gt;    lonmin = min(lon)<br>
&gt;&gt;&gt;&gt;&gt;    lonmax = 20.0<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    udiv00 = dim_avg_n_Wrap(ucomp_div00(:, :, {latmin:latmax},<br>
&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), 2)<br>
&gt;&gt;&gt;&gt;&gt;    omega00 = dim_avg_n_Wrap(w00(:, :, {latmin:latmax},<br>
&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), 2)<br>
&gt;&gt;&gt;&gt;&gt;    r00 = dim_avg_n_Wrap(rhum00(:, :, {latmin:latmax},{lonmin:lonmax<wbr>}),<br>
&gt;&gt;&gt;&gt;&gt; 2)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    udiv002 = dim_avg_n_Wrap(ucomp_div00(:, :, {latmin2:latmax2},<br>
&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), 2)<br>
&gt;&gt;&gt;&gt;&gt;    omega002 = dim_avg_n_Wrap(w00(:, :, {latmin2:latmax2},<br>
&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), 2)<br>
&gt;&gt;&gt;&gt;&gt;    r002 = dim_avg_n_Wrap(rhum00(:, :, {latmin2:latmax2},{lonmin:lonm<wbr>ax}),<br>
&gt;&gt;&gt;&gt;&gt; 2)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    elev = 1013.25*(1 - topo1*0.0065/288.15)^5.25145<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    copy_VarCoords(topo1, elev)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; ;  topo = dim_avg_n_Wrap(elev({latmin:la<wbr>tmax},{lonmin:lonmax}), 0)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    topo = elev({sellat}, {lonmin:lonmax})<br>
&gt;&gt;&gt;&gt;&gt;    topo2 = elev({sellat2}, {lonmin:lonmax})<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; ;---create plot<br>
&gt;&gt;&gt;&gt;&gt;    wks   = gsn_open_wks (&quot;eps&quot;, &quot;Overlay_Rhum_udiv_JJAS_Timmea<wbr>n_5S3N&quot;)<br>
&gt;&gt;&gt;&gt;&gt; ; open workstation<br>
&gt;&gt;&gt;&gt;&gt;    plot = new(2, graphic)<br>
&gt;&gt;&gt;&gt;&gt;    res                     = True                  ; plot mods desired<br>
&gt;&gt;&gt;&gt;&gt;    res@gsnDraw             = False    ; turn off draw<br>
&gt;&gt;&gt;&gt;&gt;    res@gsnFrame            = False    ; turn off frame<br>
&gt;&gt;&gt;&gt;&gt;    res@gsnMaximize         = True<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    res@gsnRightString      = &quot;&quot;<br>
&gt;&gt;&gt;&gt;&gt;    res@gsnLeftString       = &quot;&quot;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    cnres = res<br>
&gt;&gt;&gt;&gt;&gt;    cnres@cnFillOn          = True               ; turn on color fill<br>
&gt;&gt;&gt;&gt;&gt;    cnres@cnFillPalette     = &quot;MPL_rainbow&quot;<br>
&gt;&gt;&gt;&gt;&gt;    cnres@lbOrientation     = &quot;Vertical&quot;<br>
&gt;&gt;&gt;&gt;&gt;    cnres@pmLabelBarOrthogonalPosF = 0.08<br>
&gt;&gt;&gt;&gt;&gt;    cnres@tiYAxisString      = &quot;Pressure (hPa)&quot;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    vcres = res<br>
&gt;&gt;&gt;&gt;&gt;    vcres@vcRefMagnitudeF = 10.0                ; define vector ref mag<br>
&gt;&gt;&gt;&gt;&gt;    vcres@vcRefLengthF    = 0.1              ; define length of vec ref<br>
&gt;&gt;&gt;&gt;&gt;    vcres@vcGlyphStyle    = &quot;CurlyVector&quot;      ; turn on curley vectors<br>
&gt;&gt;&gt;&gt;&gt;    vcres@vcMinDistanceF  = 0.04               ; thin out vectors<br>
&gt;&gt;&gt;&gt;&gt;    vcres@vcMapDirection  = False<br>
&gt;&gt;&gt;&gt;&gt;    vcres@vcLineArrowThicknessF   = 3.0<br>
&gt;&gt;&gt;&gt;&gt;    vcres@vcVectorDrawOrder  = &quot;Draw&quot;        ; draw vectors last<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    contour_plota1  =  gsn_csm_pres_hgt(wks, r00(0, :,<br>
&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), cnres)<br>
&gt;&gt;&gt;&gt;&gt;    vector_plota1   = gsn_csm_vector(wks, udiv00(0, :,<br>
&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    contour_plota2  =  gsn_csm_pres_hgt(wks, r002(0, :,<br>
&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), cnres)<br>
&gt;&gt;&gt;&gt;&gt;    vector_plota2   = gsn_csm_vector(wks, udiv002(0, :,<br>
&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), -200.0*omega002(0, :, {lonmin:lonmax}), vcres )<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; ;---Add topo field using a filled polygon.<br>
&gt;&gt;&gt;&gt;&gt;    getvalues contour_plota1<br>
&gt;&gt;&gt;&gt;&gt;       &quot;trYMinF&quot; : ymin<br>
&gt;&gt;&gt;&gt;&gt;       &quot;trYMaxF&quot; : ymax<br>
&gt;&gt;&gt;&gt;&gt;    end getvalues<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; ;---Create new X,Y arrays that form a closed polygon.<br>
&gt;&gt;&gt;&gt;&gt;    nlon  = dimsizes(topo&amp;lon)<br>
&gt;&gt;&gt;&gt;&gt;    xtopo = new(nlon+3,typeof(topo&amp;lon))<br>
&gt;&gt;&gt;&gt;&gt;    ytopo = new(nlon+3,typeof(topo))<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    xtopo(0:nlon-1) = topo&amp;lon<br>
&gt;&gt;&gt;&gt;&gt;    ytopo(0:nlon-1) = topo<br>
&gt;&gt;&gt;&gt;&gt;    xtopo(nlon)     = topo&amp;lon(nlon-1)<br>
&gt;&gt;&gt;&gt;&gt;    ytopo(nlon)     = ymax          ; Use actual Y max of contour plot<br>
&gt;&gt;&gt;&gt;&gt;    xtopo(nlon+1)   = topo&amp;lon(0)<br>
&gt;&gt;&gt;&gt;&gt;    ytopo(nlon+1)   = ymax<br>
&gt;&gt;&gt;&gt;&gt;    xtopo(nlon+2)   = topo&amp;lon(0)   ; This last point closes<br>
&gt;&gt;&gt;&gt;&gt;    ytopo(nlon+2)   = topo(0)       ; the polygon.<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; ;---Add topo field using a filled polygon.<br>
&gt;&gt;&gt;&gt;&gt;    getvalues contour_plota2<br>
&gt;&gt;&gt;&gt;&gt;       &quot;trYMinF&quot; : ymin2<br>
&gt;&gt;&gt;&gt;&gt;       &quot;trYMaxF&quot; : ymax2<br>
&gt;&gt;&gt;&gt;&gt;    end getvalues<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; ;---Create new X,Y arrays that form a closed polygon.<br>
&gt;&gt;&gt;&gt;&gt;    nlon  = dimsizes(topo2&amp;lon)<br>
&gt;&gt;&gt;&gt;&gt;    xtopo2 = new(nlon+3,typeof(topo2&amp;lon))<br>
&gt;&gt;&gt;&gt;&gt;    ytopo2 = new(nlon+3,typeof(topo2))<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    xtopo2(0:nlon-1) = topo2&amp;lon<br>
&gt;&gt;&gt;&gt;&gt;    ytopo2(0:nlon-1) = topo2<br>
&gt;&gt;&gt;&gt;&gt;    xtopo2(nlon)     = topo2&amp;lon(nlon-1)<br>
&gt;&gt;&gt;&gt;&gt;    ytopo2(nlon)     = ymax2          ; Use actual Y max of contour plot<br>
&gt;&gt;&gt;&gt;&gt;    xtopo2(nlon+1)   = topo2&amp;lon(0)<br>
&gt;&gt;&gt;&gt;&gt;    ytopo2(nlon+1)   = ymax2<br>
&gt;&gt;&gt;&gt;&gt;    xtopo2(nlon+2)   = topo2&amp;lon(0)   ; This last point closes<br>
&gt;&gt;&gt;&gt;&gt;    ytopo2(nlon+2)   = topo2(0)       ; the polygon.<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; ;---Add the polygon to the contour plot.<br>
&gt;&gt;&gt;&gt;&gt;    gnres = True<br>
&gt;&gt;&gt;&gt;&gt;    gnres@gsFillColor = &quot;gray20&quot;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    id = gsn_add_polygon(wks, contour_plota1, xtopo, ytopo, gnres)<br>
&gt;&gt;&gt;&gt;&gt;    i2 = gsn_add_polygon(wks, contour_plota2, xtopo2, ytopo2, gnres)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;     overlay(contour_plota1, vector_plota1)<br>
&gt;&gt;&gt;&gt;&gt;     overlay(contour_plota2, vector_plota2)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;      plot(0)  = contour_plota1<br>
&gt;&gt;&gt;&gt;&gt;      plot(1)  = contour_plota2<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    ;draw(contour_plot)     ; This draws everything<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;    pres = True<br>
&gt;&gt;&gt;&gt;&gt;    pres@gsnMaximize = True<br>
&gt;&gt;&gt;&gt;&gt;    gsn_panel(wks, plot, (/2, 1/), pres)<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; end<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt; 2017-05-05 0:39 GMT+01:00 Mary Haley &lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;:<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; I&#39;m so sorry I didn&#39;t get back to you sooner on this.  I&#39;ve been<br>
&gt;&gt;&gt;&gt;&gt;&gt; swamped and then was out-of-town for a few days.<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; It&#39;s a little hard to follow what&#39;s going on in this script without<br>
&gt;&gt;&gt;&gt;&gt;&gt; being able to run it.<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; I think the issue is that you are using the &quot;res&quot; variable in your<br>
&gt;&gt;&gt;&gt;&gt;&gt; call to gsn_panel, and &quot;res&quot; has both gsnDraw and gsnFrame set to False.<br>
&gt;&gt;&gt;&gt;&gt;&gt; This will cause *no* plot to appear when you call gsn_panel.  I&#39;m surprised<br>
&gt;&gt;&gt;&gt;&gt;&gt; you are getting any plot at all.<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; I recommend creating a new resource variable for gsn_panel, because<br>
&gt;&gt;&gt;&gt;&gt;&gt; you really don&#39;t want to use the one that you used for the plots themselves:<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; pres = True<br>
&gt;&gt;&gt;&gt;&gt;&gt; pres@gsnMaximize = True<br>
&gt;&gt;&gt;&gt;&gt;&gt; gsn_panel(wks, plot, (/2, 1/), pres)<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; If you<br>
&gt;&gt;&gt;&gt;&gt;&gt; ? continue to have problems, then it would help if you send me the<br>
&gt;&gt;&gt;&gt;&gt;&gt; two data files. If you can&#39;t send the files, then send me the image you&#39;re<br>
&gt;&gt;&gt;&gt;&gt;&gt; getting, and any error messages ?that the script produces.<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; ?Thanks,<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; --Mary<br>
&gt;&gt;&gt;&gt;&gt;&gt; ?<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt; On Tue, Apr 25, 2017 at 12:56 AM, Appo derbetini &lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a><br>
&gt;&gt;&gt;&gt;&gt;&gt; &gt; wrote:<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; Good morning Mary,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; Sorry , i forgot to add ncl-users email to the reply.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; I am attaching the modified version of the script were I&#39;m trying to<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; create a panel plot for two transects at latitude sellat and sellat2.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; Unfortunately, only one plot appears.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; Really, I don&#39;t know what is wrong this script.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; Your help will be appreciated.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; begin<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;     in1  = addfile(&quot;<a href="http://uvwr.mean.JJAS.nc" rel="noreferrer" target="_blank">uvwr.mean.JJAS.nc</a>&quot;,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; &quot;r&quot;)                         ; open netcdf file<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;     input4  = addfile(&quot;<a href="http://atanas_topo.nc" rel="noreferrer" target="_blank">atanas_topo.nc</a>&quot;, &quot;r&quot;)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    u00     =  in1-&gt;u<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    v00     =  in1-&gt;v<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    w00     =  in1-&gt;w<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    rhum00  =  in1-&gt;r<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    topo1 = input4-&gt;HT<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    lat = in1-&gt;latitude<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    lon = in1-&gt;longitude<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    dv00    = uv2dvF (u00, v00)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    uvd00   = dv2uvF (dv00)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ucomp_div00 = uvd00(0, :, : , :, :)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    copy_VarCoords(u00, ucomp_div00)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    sellat = 3.0<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    sellat2 = -2.5<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    latmin = -5.0<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    latmax = 3.0<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    lonmin = min(lon)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    lonmax = 20.0<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    udiv00 = dim_avg_n_Wrap(ucomp_div00(:, :, {latmin:latmax},<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), 2)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    omega00 = dim_avg_n_Wrap(w00(:, :, {latmin:latmax},<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), 2)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    r00 = dim_avg_n_Wrap(rhum00(:, :, {latmin:latmax},{lonmin:lonmax<wbr>}),<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; 2)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ;udiv00  = ucomp_div00(:, :, {sellat}, {lonmin:lonmax})<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ;omega00 = w00(:, :, {sellat}, {lonmin:lonmax})<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ;r00     = rhum00(:, :, {sellat},{lonmin:lonmax})<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    elev = 1013.25*(1 - topo1*0.0065/288.15)^5.25145<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    copy_VarCoords(topo1, elev)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; ;  topo = dim_avg_n_Wrap(elev({latmin:la<wbr>tmax},{lonmin:lonmax}), 0)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    topo  = elev({sellat}, {lonmin:lonmax})<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    topo2 = elev({sellat2}, {lonmin:lonmax})<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; ;---create plot<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    wks   = gsn_open_wks (&quot;eps&quot;, &quot;Overlay_Rhum_udiv_topo&quot;)        ;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; open workstation<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;     plot = new(2, graphic)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    res                     = True                  ; plot mods<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; desired<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    res@gsnDraw             = False    ; turn off draw<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    res@gsnFrame            = False    ; turn off frame<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    res@gsnMaximize         = True<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    res@gsnRightString      = &quot;&quot;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    res@gsnLeftString       = &quot;&quot;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    cnres = res<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    cnres@cnFillOn          = True               ; turn on color fill<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    cnres@cnFillPalette     = &quot;MPL_rainbow&quot;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    cnres@lbOrientation     = &quot;Vertical&quot;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    cnres@pmLabelBarOrthogonalPosF = 0.08<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    cnres@tiYAxisString      = &quot;Pressure (hPa)&quot;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vcres = res<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vcres@vcRefMagnitudeF = 10.0                ; define vector ref<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; mag<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vcres@vcRefLengthF    = 0.1              ; define length of vec<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; ref<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vcres@vcGlyphStyle    = &quot;CurlyVector&quot;      ; turn on curley<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; vectors<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vcres@vcMinDistanceF  = 0.04               ; thin out vectors<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vcres@vcMapDirection  = False<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vcres@vcLineArrowThicknessF   = 3.0<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vcres@vcVectorDrawOrder  = &quot;Draw&quot;        ; draw vectors last<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    contour_plot_1 = gsn_csm_pres_hgt(wks, r00(0, :,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), cnres)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vector_plot_1  = gsn_csm_vector(wks, udiv00(0, :,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    contour_plot_2 = gsn_csm_pres_hgt(wks, r00(0, :,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), cnres)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    vector_plot_2  = gsn_csm_vector(wks, udiv00(0, :,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), vcres )<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; ;---Add topo field using a filled polygon.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    getvalues contour_plot_1<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;       &quot;trYMinF&quot; : ymin<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;       &quot;trYMaxF&quot; : ymax<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    end getvalues<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; ;---Add topo field using a filled polygon.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    getvalues contour_plot_2<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;       &quot;trYMinF&quot; : ymin2<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;       &quot;trYMaxF&quot; : ymax2<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    end getvalues<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; ;---Create new X,Y arrays that form a closed polygon.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    nlon  = dimsizes(topo&amp;lon)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo = new(nlon+3,typeof(topo&amp;lon))<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo = new(nlon+3,typeof(topo))<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo(0:nlon-1) = topo&amp;lon<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo(0:nlon-1) = topo<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo(nlon)     = topo&amp;lon(nlon-1)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo(nlon)     = ymax          ; Use actual Y max of contour plot<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo(nlon+1)   = topo&amp;lon(0)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo(nlon+1)   = ymax<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo(nlon+2)   = topo&amp;lon(0)   ; This last point closes<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo(nlon+2)   = topo(0)       ; the polygon.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    nlon2  = dimsizes(topo2&amp;lon)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo2 = new(nlon2+3,typeof(topo2&amp;lon))<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo2 = new(nlon2+3,typeof(topo2))<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo2(0:nlon-1) = topo2&amp;lon<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo2(0:nlon-1) = topo2<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo2(nlon)     = topo2&amp;lon(nlon-1)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo2(nlon)     = ymax2         ; Use actual Y max of contour<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; plot<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo2(nlon+1)   = topo2&amp;lon(0)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo2(nlon+1)   = ymax2<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    xtopo2(nlon+2)   = topo2&amp;lon(0)   ; This last point closes<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    ytopo2(nlon+2)   = topo2(0)       ; the polygon.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; ;---Add the polygon to the contour plot.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    gnres = True<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    gnres@gsFillColor = &quot;gray20&quot;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    id = gsn_add_polygon(wks, contour_plot_1, xtopo, ytopo, gnres)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    id2 = gsn_add_polygon(wks, contour_plot_2, xtopo2, ytopo2, gnres)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    overlay(contour_plot_1, vector_plot_1)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    overlay(contour_plot_2, vector_plot_2)<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;      draw(contour_plot_1)     ; This draws everything<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;      draw(contour_plot_2)     ; This draws everything<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    plot(0) = contour_plot_1<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;    plot(1) = contour_plot_2<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;     gsn_panel(wks,plot, (/2, 1/), res)               ; now draw as<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; one plot<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; end<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt; 2017-04-24 21:58 GMT+01:00 Mary Haley &lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Hi,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; In the future, please email ncl-talk with follow-up questions so<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; everybody can benefit from the answers.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; What is the exact error message you&#39;re getting?<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; --Mary<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; On Mon, Apr 24, 2017 at 4:08 AM, Appo derbetini &lt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; <a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt; wrote:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Hi Mary,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Since you send me the script that plot masked fields, I am working<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; to do it for many latitude in a panel plot.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Always i got many errors saying that i cannot overlay.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; How to solve it?<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Thank you in advance.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Appo<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; 2017-04-10 8:23 GMT+01:00 Appo derbetini &lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt;:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Thank you very much.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; It&#39;s look like what I&#39;m trying to do.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; regards<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; 2017-04-09 22:11 GMT+01:00 Mary Haley &lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; I&#39;m sorry I didn&#39;t get back to this sooner, but I was having<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; trouble trying to understand what you are doing in the script.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; First, I&#39;m not sure it makes sense to take an average of<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; topographical data across a set of latitudes and plot that as a topo layer.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; But, that aside, I think you just want to plot the topo line as<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; a filled polygon instead of a separate XY plot.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; I&#39;m not sure what I&#39;ve attached is correct, but hopefully it<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; gives you an idea of what to do.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Note that I&#39;m using &quot;sellat&quot; to select a single latitude value,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; rather than averaging across a range of latitudes.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; --Mary<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; On Sun, Apr 9, 2017 at 2:10 AM, Appo derbetini &lt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; <a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt; wrote:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Dear Mary,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Despite applying what you suggest the problem remain the same<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Thank you very much<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; 2017-04-06 8:53 GMT+01:00 Appo derbetini &lt;<a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt;:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Dear Mary,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Datasets are uploaded as indicated.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; What I really trying to do is something like Mask Example 14.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; But I don&#39;t want to interpolate topography to grid of wind.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Here attached figure produced by the script<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Thank you<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; 2017-04-05 15:36 GMT+01:00 Mary Haley &lt;<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>&gt;:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; I don&#39;t think I&#39;ve ever seen this error before.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; It&#39;s hard to debug this particular script without having the<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; data so I can run it. But, I&#39;m wondering if this line might be part of the<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; problem:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;   plot  = gsn_csm_pres_hgt_vector(wks, r00(0, :,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), udiv00(0, :, {lonmin:lonmax}), -200.0*omega00(0, :,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), res )<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; When you do an arithmetic operation on an array, like:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;   -200. * omega00<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; you end up stripping all of the metadata off omega00.  What&#39;s<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; getting passed into gsn_csm_pres_hgt has no metadata. This might be okay,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; however, because I think this routine might just use the metadata from<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; &quot;udiv00&quot; instead.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Still, it&#39;s worth trying this:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;  omega00_scale = omega00(0, :, {lonmin:lonmax})   ; subsets<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; omega00 *and* copies all metadata<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;  omega00_scale = -200.0 * omega00_scale<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;   plot  = gsn_csm_pres_hgt_vector(wks, r00(0, :,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; {lonmin:lonmax}), udiv00(0, :, {lonmin:lonmax}), omega00_scale, res )<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; If you continue to have problems with this script, could you<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; upload your data (if it&#39;s not too large) to our ftp:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; ftp <a href="http://ftp.cgd.ucar.edu" rel="noreferrer" target="_blank">ftp.cgd.ucar.edu</a><br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; anonymous<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; &lt;use your email address for the password&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; cd incoming<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; put <a href="http://uvwr.mean.JJAS.nc" rel="noreferrer" target="_blank">uvwr.mean.JJAS.nc</a><br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; put <a href="http://atanas_topo.nc" rel="noreferrer" target="_blank">atanas_topo.nc</a><br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; quit<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Thanks,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; --Mary<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; On Tue, Apr 4, 2017 at 9:18 AM, Appo derbetini &lt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; <a href="mailto:appopson4@gmail.com" target="_blank">appopson4@gmail.com</a>&gt; wrote:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Dear all,<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; I am trying to overlay topography profil on a profil of wind<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; and relative humidity.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; I don&#39;t want to interpolate topography to show some details.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Atached are script and datasets used for the plot.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; But I having errors.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; warning:LlDataPolygon: point 6.150009,0.000000 outside data<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; domain<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; warning:LlDataPolygon: point 19.983337,0.000000 outside data<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; domain<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; warning:LlDataLineTo: point 6.150009,0.000000 outside data<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; domain<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; warning:LlDataLineTo: point 19.983337,0.000000 outside data<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; domain<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Any help will be appreciated.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; Best regards.<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; ______________________________<wbr>_________________<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; ncl-talk mailing list<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; List instructions, subscriber options, unsubscribe:<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt; <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>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;<br>
&gt;<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170512/89245831/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/piperm<wbr>ail/ncl-talk/attachments/20170<wbr>512/89245831/attachment-0001.<wbr>html</a><br>
<br>
------------------------------<br>
<br>
______________________________<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>
<br>
End of ncl-talk Digest, Vol 162, Issue 26<br>
******************************<wbr>***********<br>
</blockquote></div><br></div></div>
<br>______________________________<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><br clear="all"><div><br></div>-- <br></div></div><span class="HOEnZb"><font color="#888888"><div class="m_4791061683926175533gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist,  </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/<wbr>asphilli/</a>   </font></span><span><font color="#888888"><a href="tel:(303)%20497-1726" value="+13034971726" target="_blank">303-497-1726</a> </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</font></span></div>
</blockquote></div><br></div>