check_for_y_lat_coord: Warning: Data either does not contain ...

From: <antti.pessi_at_nyahnyahspammersnyahnyah>
Date: Thu Sep 16 2010 - 11:50:32 MDT

Hi,

 

I'm trying to plot regular WRF grib2 output that comes out from WRF
Post-Processor. I get these warnings about lat/lon coordinate array and
units. The units should be correct in the file (degrees_north), but the
lat/lon variable name is gridlat_0 / gridlon_0. How do I make NCL to
find the correct coordinate array? I have seen this problem before in
this list, but I couldn't make it work using those tips.

 

NCL code simplified:

-----------------------------

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

 

begin

  a = addfile("1009141400_f0100-co-d01.grb2","r")

 

  temp = a->TMP_P0_L100_GLC0({50000},:,:)

  uwnd = a->UGRD_P0_L100_GLC0({50000},:,:)

 

  wks = gsn_open_wks("x11","overlay")

  plot = gsn_csm_contour_map_overlay(wks,temp,uwnd,res,sres)

end

 

NCL error message:

----------------------------

(0) check_for_y_lat_coord: Warning: Data either does not contain a
valid latitude coordinate array or doesn't contain one at all.

(0) A valid latitude coordinate array should have a 'units'
attribute equal to one of the following values:

(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg
north'

(0) check_for_lon_coord: Warning: Data either does not contain a
valid longitude coordinate array or doesn't contain one at all.

(0) A valid longitude coordinate array should have a 'units'
attribute equal to one of the following values:

(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'

 

Some output from ncl_filedump grib2file:

--------------------------------------------------------

   file global attributes:

   dimensions:

      ygrid_0 = 156

      xgrid_0 = 162

      lv_ISBL0 = 39

.................

      float TMP_P0_L100_GLC0 ( lv_ISBL0, ygrid_0, xgrid_0 )

         center : US National Weather Service - NCEP (WMC)

         production_status : Operational products

         long_name : Temperature

         units : K

         _FillValue : 1e+20

         coordinates : gridlat_0 gridlon_0

         grid_type : Lambert Conformal can be secant or tangent,
conical or bipolar

         parameter_discipline_and_category : Meteorological products,
Temperature

         parameter_template_discipline_category_number : ( 0, 0,
0, 0 )

         level_type : Isobaric surface (Pa)

..........

      float gridlat_0 ( ygrid_0, xgrid_0 )

         corners : ( 37.203, 37.20306, 42.77142, 42.77136 )

         long_name : latitude

         grid_type : Lambert Conformal (secant, tangent, conical or
bipolar)

         units : degrees_north

...............

      float gridlon_0 ( ygrid_0, xgrid_0 )

         corners : ( -108.441, -101.1767, -100.8674, -108.7505 )

         long_name : longitude

         grid_type : Lambert Conformal (secant, tangent, conical or
bipolar)

         units : degrees_east

 

 

Any help appreciated.

 

Thanks,

Antti

 

----------------------------------------------

Antti T. Pessi, PhD

Meteorologist

Vaisala Inc.

194 South Taylor Avenue

Louisville, Colorado 80027

Email: antti.pessi@vaisala.com

Phone: +1 303-402-4758

Fax: +1 303-499-1767

http://www.vaisala.com

----------------------------------------------

 

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 16 11:50:49 2010

This archive was generated by hypermail 2.1.8 : Mon Sep 20 2010 - 15:45:07 MDT