[ncl-talk] Overlay profile of topography

[Appo derbetini]

Dear Mary,

Datasets are uploaded as indicated.

What I really trying to do is something like Mask Example 14.

But I don't want to interpolate topography to grid of wind.

Here attached figure produced by the script

Thank you


2017-04-05 15:36 GMT+01:00 Mary Haley <haley at ucar.edu>:

> 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/20170406/fd5520f1/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Overlay.pdf
Type: application/pdf
Size: 89254 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170406/fd5520f1/attachment.pdf