Re: Plotting eta data

Hi Agnes,

Instead of naming your attributes "lat2d_eta" and "lon2d_eta", try
"lat2d" and "lon2d":

GHGT_eta_at_lat2d=eta->gridlat_212
GHGT_eta_at_lon2d=eta->gridlon_212

The the "lat2d", "lon2d" attributes are speicla ones that some of the
gsn_csm_xxxx scripts look for,
if no 1D coordinate variables are available.

--Mary

On Nov 12, 2008, at 10:34 AM, Agnes Lim wrote:

> Hi
>
> I had the following code segment to read the geopotential data from
> an ETA model grib file.
> GHGT_eta=eta->HGT_212_ISBL(it_eta,:,:,:)
> GHGT_eta = GHGT_eta(::-1,:, :)
> GHGT_eta_at_lat2d_eta=eta->gridlat_212
> GHGT_eta_at_lon2d_eta=eta->gridlon_212
> When I tried plotting the data, I get the following complaints
> (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'
>
> (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 ea
>
> Wonder why is that so? I encoutered that for wrf data after some
> processing but doing what i did above remove that problem.
>
> Agnes
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 12 2008 - 10:41:50 MST