[ncl-talk] Difficulty with RAP coordinate arrays in NCL

[Dennis Shea]

The only 'rap...grb2' file I have has:

      float gridrot_0 ( ygrid_0, xgrid_0 )
         long_name :    vector rotation angle
         GridType :     Lambert Conformal (secant, tangent, conical or
bipolar)
         units :        radians
         formula_u :    Uearth = sin(rot)*Vgrid + cos(rot)*Ugrid
         formula_v :    Vearth = cos(rot)*Vgrid - sin(rot)*Ugrid
         note1 :        u and v components of vector quantities are
resolved relative to grid
         note2 :        apply formulas to derive u and v components
relative to earth

      float gridlat_0 ( ygrid_0, xgrid_0 )
         corners :      ( 16.281, 17.34023, 55.48131, 54.17241 )
         long_name :    latitude
         grid_type :    Lambert Conformal (secant, tangent, conical or
bipolar)
         units :        degrees_north
         Latin2 :       25
         Latin1 :       25
         Dy :   13.545
         Dx :   13.545
         Lov :  265
         Lo1 :  233.862
         La1 :  16.281

      float gridlon_0 ( ygrid_0, xgrid_0 )
[SNIP]

Your coordimates are contained in 2-dimenional variables.
Could you pls ftp the file to

ftp ftp.cgd.ucar.edu
anonymous
ypur_email
cd incoming
put rap_130_20160115_0000_000.grb2
quit

===

Attached is a simple script I have lying around. It does not use shape
files. Try running it with your data.

Note:

  lat2d    = f->gridlat_0   ; these are 2D
  lon2d   = f->gridlon_0

Later,

  res at vfXArray               = lon2d
  res at vfYArray               = lat2d

Also, depending on you usage, you may want to rotate the winds.

===
FYI: I always helps to indicate which version of NCL you are using.

%> ncl -V

===
Also, printVarSummary( ...)    of selected variables.

Good luck




On Fri, Sep 16, 2016 at 12:48 PM, Black, Forrest (LARC-D318)[UNIVERSITIES
SPACE RESEARCH ASSOCIATION] <forrest.black at nasa.gov> wrote:

> Hello,
>
>
>
> I am having difficulty plotting RAP data showing wind and reflectivity
> data overlaying a shapefile of the Northeast US. I looked to a previous
> correspondence you had on the subject
> <http://mailman.ucar.edu/pipermail/ncl-talk/2015-August/003632.html>, but
> have been unable to fix the error. The map plots, but I get the message:
>
>
>
> “check_for_y_lat_coord: Warning: Data either does not contain a valid
> latitude coordinate array or doesn't contain one at all.”
>
>
>
> I uploaded the NCL script along with the relevant files to your ftp. The
> folder I uploaded is “ncl_talk_forrestblack.7z”.  Attached to this email is
> also the script in case that is all you need.
>
>
>
> Thank you for your help,
>
>
>
> -        Forrest
>
>
>
> Forrest Black
>
> NIFS Intern
>
> Crew Systems & Aviation Operations
>
> NASA Langley Research Center
>
> forrest.black at nasa.gov
>
>
>
> _______________________________________________
> 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/20160916/cab65f87/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rap.ncl
Type: text/x-ncl
Size: 3397 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160916/cab65f87/attachment.bin