[ncl-talk] Overlay profile of topography

[Mary Haley]

I don't think I've ever seen this error before.

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:

  plot  = gsn_csm_pres_hgt_vector(wks, r00(0, :, {lonmin:lonmax}),
udiv00(0, :, {lonmin:lonmax}), -200.0*omega00(0, :, {lonmin:lonmax}), res )

When you do an arithmetic operation on an array, like:

  -200. * omega00

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.

Still, it's worth trying this:

 omega00_scale = omega00(0, :, {lonmin:lonmax})   ; subsets omega00 *and*
copies all metadata
 omega00_scale = -200.0 * omega00_scale

  plot  = gsn_csm_pres_hgt_vector(wks, r00(0, :, {lonmin:lonmax}),
udiv00(0, :, {lonmin:lonmax}), omega00_scale, res )

If you continue to have problems with this script, could you upload your
data (if it's not too large) to our ftp:

ftp ftp.cgd.ucar.edu
anonymous
<use your email address for the password>
cd incoming
put uvwr.mean.JJAS.nc
put atanas_topo.nc
quit

Thanks,

--Mary


On Tue, Apr 4, 2017 at 9:18 AM, Appo derbetini <appopson4 at gmail.com> wrote:

> Dear all,
> I am trying to overlay topography profil on a profil of wind and relative
> humidity.
> I don't want to interpolate topography to show some details.
>
> Atached are script and datasets used for the plot.
>
> But I having errors.
>
> warning:LlDataPolygon: point 6.150009,0.000000 outside data domain
> warning:LlDataPolygon: point 19.983337,0.000000 outside data domain
> warning:LlDataLineTo: point 6.150009,0.000000 outside data domain
> warning:LlDataLineTo: point 19.983337,0.000000 outside data domain
>
> Any help will be appreciated.
>
> Best regards.
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170405/f268ef05/attachment-0001.html