netCDF in a mapfile

[Tony Jolibois]

          Hi all,

I try to display  netCDF file with the GDAL driver, but no succeed.

The gdalinfo give a good response, but I don't know how I can specifie
my variable and the depth for the layer in the mapfile (it is a 3
dimension variable).
Can anyone help me ?

Thanks in advance,
Tony

gdal from CVS today
mapserver 4.4.1
on linux AS

gdalinfo :
Driver: netCDF/network Common Data Format
Size is 8, 129
Coordinate System is `'
Origin = (0.000000,0.000000)
Pixel Size = (1.00000000,1.00000000)
Corner Coordinates:
Upper Left  (   0.0000000,   0.0000000)
Lower Left  (       0.000,     129.000)
Upper Right (   8.0000000,   0.0000000)
Lower Right (       8.000,     129.000)
Center      (   4.0000000,  64.5000000)
Band 1 Block=8x1 Type=Float32, ColorInterp=Undefined
...
...

The header of the netCDF is like this :
netcdf mercator_psy2v1_mersea_grid1o8_med_b20050105_f200501189999 {
dimensions:
        depth = 8 ;
        latitude = 129 ;
        longitude = 341 ;
variables:
        float depth(depth) ;
                depth:valid_min = 5.f ;
                depth:valid_max = 2000.f ;
                depth:units = "m" ;
                depth:positive = "down" ;
                depth:unit_long = "Meter" ;
                depth:long_name = "depth" ;
                depth:standard_name = "depth" ;
                depth:axis = "Z" ;
        float latitude(latitude) ;
                latitude:valid_min = 30.f ;
                latitude:valid_max = 46.f ;
                latitude:step = 0.125f ;
                latitude:units = "degrees_north" ;
                latitude:unit_long = "Degrees North" ;
                latitude:long_name = "latitude" ;
                latitude:standard_name = "latitude" ;
                latitude:axis = "Y" ;
        float longitude(longitude) ;
                longitude:valid_min = -6.f ;
                longitude:valid_max = 36.5f ;
                longitude:step = 0.125f ;
                longitude:units = "degrees_east" ;
                longitude:unit_long = "Degrees East" ;
                longitude:long_name = "longitude" ;
                longitude:standard_name = "longitude" ;
                longitude:axis = "X" ;
        float temperature(depth, latitude, longitude) ;
                temperature:long_name = "temperature" ;
                temperature:standard_name =
"sea_water_potential_temperature" ;
                temperature:units = "degC" ;
                temperature:unit_long = "Degrees Celcius" ;
                temperature:missing_value = 1.e+35f ;
                temperature:_FillValue = 1.e+35f ;
                temperature:valid_min = 9.834617f ;
                temperature:valid_max = 20.74316f ;
// global attributes:
                :title = "MERCATOR PSY2V1R1" ;
                :easting = "longitude" ;
                :northing = "latitude" ;
                :history = "2005/01/ 5 19:58:31 MERCATOR Netcdf creation" ;
                :source = "MERCATOR PSY2V1" ;
                :institution = "GIP MERCATOR OCEAN" ;
                :references = "http://www.mercator-ocean.fr" ;
                :comment = " " ;
                :conventions = "CF-1.0" ;
                :domain_name = "MED" ;
                :field_type = "mean" ;
                :field_date = "2005-01-18 00:00:00" ;
                :field_julian_date = 20106 ;
                :julian_day_unit = "days since 1950-01-01 00:00:00" ;
                :forecast_range = "13-day forecast" ;
                :forecast_type = "forecast" ;
                :bulletin_date = "2005-01-05 00:00:00" ;
                :bulletin_type = "operational" ;
                :longitude_min = -6.f ;
                :longitude_max = 36.5f ;
                :latitude_min = 30.f ;
                :latitude_max = 46.f ;
                :z_min = 5.f ;
                :z_max = 2000.f ;
}


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20050120/c2752cd6/attachment.htm>