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