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

[Black, Forrest (LARC-D318)[UNIVERSITIES SPACE RESEARCH ASSOCIATION]]

Hello,

I tried integrating the example script you gave me, but now there is an error on line 24 which is where I made the new additions to the script.
I have uploaded sample data to your ftp and it is named “fblack_sample.grb2”. I am also including two scripts “test6.ncl” is without the 2D variables script that worked (just didn’t display the data correctly) and the updated script “test6b.ncl” that has a problem at line 24.

The NCL version I am using is 6.3.0

The details for the variables are attached in “varDetails.txt”

Thank you,


-        Forrest

Forrest Black
NIFS Intern
Crew Systems & Aviation Operations
NASA Langley Research Center
Building 2168A
Room 2118
Phone: 757-864-2772
Forrest.black at nasa.gov<mailto:Forrest.black at nasa.gov>



From: Dennis Shea [mailto:shea at ucar.edu]
Sent: Friday, September 16, 2016 6:01 PM
To: Black, Forrest (LARC-D318)[UNIVERSITIES SPACE RESEARCH ASSOCIATION] <forrest.black at nasa.gov<mailto:forrest.black at nasa.gov>>
Cc: ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
Subject: Re: [ncl-talk] Difficulty with RAP coordinate arrays in NCL

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<http://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<mailto: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<mailto:forrest.black at nasa.gov>


_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu<mailto: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/20160919/99237d24/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test6.ncl
Type: application/octet-stream
Size: 3587 bytes
Desc: test6.ncl
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160919/99237d24/attachment.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test6b.ncl
Type: application/octet-stream
Size: 3587 bytes
Desc: test6b.ncl
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160919/99237d24/attachment-0001.obj 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: varSummary.txt
Url: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160919/99237d24/attachment.txt