NetCDF question

Don Hooper hoop at COLORADO.EDU
Fri May 12 17:31:56 EDT 2006


Kobi,

Well, like I told you privately, the sdfopen code isn't designed for
station data.  My foolish code assumes there exists more than one
value in the x and y vectors (your "lat" and "lon"), and probably
gets very confused as a result.  You might try writing an xdfopen-style
Data Descriptor File (DDF) (probably only need the DSET, XDEF, and YDEF
lines).  See the COLA GrADS documentation web site for the details.

-Hoop
> From owner-gradsusr at LIST.CINECA.IT  Fri May 12 15:18:22 2006
> From: Kobi Mosquera <kobi at CHAVIN.IGP.GOB.PE>
>
> Hoop,
>
> sorry!!, the "ncdump -c" output is:
>
> ================================================
> netcdf data_igp3 {
> dimensions:
>         time = 72 ;
>         height = 1 ;
>         lat = 1 ;
>         lon = 1 ;
> variables:
>         float time(time) ;
>                 time:units = "hours since 2000-1-1 00:00:0.0 -5:00" ;
>                 time:long_name = "Time" ;
>                 time:delta_t = "0000-00-00 01:00:00" ;
>                 time:axis = "T" ;
>         float height(height) ;
>                 height:actual_range = 0.f, 0.f ;
>                 height:units = "m " ;
>                 height:axis = "Z" ;
>         float lat(lat) ;
>                 lat:actual_range = -5.f, -5.f ;
>                 lat:units = "degrees_north" ;
>                 lat:long_name = "Latitude" ;
>                 lat:axis = "Y" ;
>         float lon(lon) ;
>                 lon:actual_range = -78.f, -78.f ;
>                 lon:units = "degrees_east" ;
>                 lon:long_name = "Longitude" ;
>                 lon:axis = "X" ;
>         float temp_aire(time, height, lat, lon) ;
>                 temp_aire:name = "AT " ;
>                 temp_aire:long_name = "Temperatura del aire " ;
>                 temp_aire:units = "degC  " ;
>                 temp_aire:actual_range = 5.f, 50.f ;
>                 temp_aire:missing_value = -999.f ;
>
> // global attributes:
>                 :Source = "Instituto Geofisico del Peru ( http://www.met.igp.gob.pe )" ;
>                 :File_info = "Contactar a Kobi MOsquera  kobi at chavin.igp.gob.pe " ;
>                 :title = "Temperatura del aire  en la estacion La Calurosa         " ;
>                 :description = "Temperatura del aire  de la estacion La
> Calurosa         " ;
>                 :missing_value = -999.f ;
>                 :Convention = "COARDS" ;
> data:
>
>  lat = -5 ;
>
>  lon = -78 ;
> }
>
> ================================================================
>
> Thanks,
>
> Kobi
>
> --------------------------------------------------------
> Kobi A. Mosquera Vásquez
> Centro de Predicción Numérica del Tiempo y Clima (CPNTC)
> Instituto Geofísico del Perú (IGP)
>
> Calle Badajoz N° 169 – 171
> Urbanización Mayorazgo,IV Etapa
> Lima 3 - Perú
> Telefax (51) 1 3172326
> http://www.met.igp.gob.pe
> --------------------------------------------------------



More information about the gradsusr mailing list