<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="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="h5">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="h5"><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>