Re: the longitude problem when read grib format data

Hi Li,

Is it possible for you to provide us with the
"lcld.ensm.apr.cfs.data.grb" file so we can take a look?

If so, please email me separately and I'll give you some information
on where you can put it.

Thanks,

--Mary

On Tue, 5 Aug 2008, lizhou wrote:

> Dear all,
>
> when I read CFS output with grib format, I got wrong longitude. the longitude is "lon_98: [ 0..1.875]" instead of beginning with 0 with 1.875 step. does anyone know what's happen?
>
> Thanks in advance,
> Li
> -------------------------------------
>
>
> The read code is:
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
> in = addfile("data/lcld.ensm.apr.cfs.data.grb","r")
> lcld = in->T_CDC_98_LCY_ave1m
> printVarSummary(lcld)
> end
>
> The output information is:
> Variable: lcld
> Type: float
> Total Size: 14943744 bytes
> 3735936 values
> Number of Dimensions: 4
> Dimensions and sizes: [initial_time0 | 23] x [forecast_time1 | 9] x [lat_98 | 94] x [lon_98 | 192]
> Coordinates:
> initial_time0: [04/01/1981 (00:00)..04/01/2003 (00:00)]
> forecast_time1: [1..9]
> lat_98: [88.54195..-88.54195]
> lon_98: [ 0..1.875]
> Number Of Attributes: 9
> sub_center : Environmental Modeling Center
> center : US National Weather Service - NCEP (WMC)
> long_name : Total cloud cover
> units : %
> _FillValue : -999
> level_indicator : 214
> gds_grid_type : 4
> parameter_table_version : 2
> parameter_number : 71
>
>
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 06 2008 - 07:02:34 MDT