
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">

<HTML>

<HEAD>

  <META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">

  <META NAME="GENERATOR" CONTENT="GtkHTML/3.16.1">

</HEAD>

<BODY>

Hello,<BR>

<BR>

I would like to create Web Coverage Service from one netCDF file in following <A HREF="http://mapserver.gis.umn.edu/docs/howto/WCSServerFormatHowTo">http://mapserver.gis.umn.edu/docs/howto/WCSServerFormatHowTo</A><BR>

<BR>

My netCDF store rainfall matrix every five minutes in the "pluie" variable. This the CDL :<BR>

<BR>

dimensions:

<PRE>

   time = UNLIMITED;   // (496 currently)   // (has coord.var)

   x = 50;   // (has coord.var)

   y = 50;   // (has coord.var)

variables:

   int Lambert_Conformal;

     :grid_mapping_name = "lambert_conformal_conic";

     :standard_parallel = 43.20317, 44.99683; // double

     :longitude_of_central_meridian = 2.337229; // double

     :latitude_of_projection_origin = 44.1; // double

     :false_easting = 600.0; // double

     :false_northing = 3200.0; // double

     :_CoordinateTransformType = "Projection";

     :_CoordinateSystems = "ProjectionCoordinateSystem";

     :_CoordinateAxes = "y x";

     :_CoordinateAxisTypes = "GeoX GeoY";

<B>   double pluie(time=496, y=50, x=50);</B>

<B>     :units = "mm/h";</B>

<B>     :long_name = "Intensité de pluie";</B>

<B>     :coordinates = "lat lon";</B>

<B>     :grid_mapping = "Lambert_Conformal";</B>

<B>     :_CoordinateSystems = "ProjectionCoordinateSystem LatLonCoordinateSystem";</B>

   char LatLonCoordinateSystem;

     :_CoordinateAxes = "time lat lon";

   char ProjectionCoordinateSystem;

     :_CoordinateAxes = "time y x";

     :_CoordinateTransforms = "LambertConformalProjection";

   double time(time=496);

     :long_name = "rainfall time";

     :standard_name = "time";

     :units = "seconds since 1970-01-01 01:00:00";

     :axis = "T";

     :_CoordinateAxisType = "Time";

     :bounds = "time_bnds";

   double x(x=50);

     :units = "km";

     :long_name = "x coordinate of projection";

     :standard_name = "projection_x_coordinate";

     :_CoordinateAxisType = "GeoX";

   double y(y=50);

     :units = "km";

     :long_name = "y coordinate of projection";

     :standard_name = "projection_y_coordinate";

     :_CoordinateAxisType = "GeoY";

   double lat(x=50, y=50);

     :units = "degrees_north";

     :long_name = "latitude coordinate";

     :standard_name = "latitude";

     :_CoordinateAxisType = "Lat";

   double lon(x=50, y=50);

     :units = "degrees_east";

     :long_name = "longitude";

     :standard_name = "longitude coordinate";

     :_CoordinateAxisType = "Lon";


 :Zone_Radar = "Gard-Ceze Amont";

 :history = "Tuesday 19 February 2008, 00:00";

 :Projection_cartographique = "Lambert III";

 :X0 = 700.7; // double

 :Y0 = 3189.1; // double

 :Conventions = "CF-1.0";

}

</PRE>

<BR>

The HowTo of Mapserver Website seems adapt to a set of netcdf files but not for one netCDF file including several bands corresponding to each time step.<BR>

When I use :<BR>

<B>gdal_translate -ot Float64 -of gTiff -b 300 -sds netCDF:"pluie-4-20020908.nc":pluie test.tif</B><BR>

A tif file is created based on the 300th band of "pluie" variable<BR>

<BR>

When I use :<BR>

<B>gdaltindex test.shp NETCDF:pluie-3-20020908.nc:pluie</B><BR>

A shape file is created but without every band, only a record with the netcdf extent and his name<BR>

<BR>

I suppose i have to create a new record for each time's band of netcdf file in the tileindex after creating with gdaltindex ?!<BR>

<BR>

I would like to use the Python Script of Norman Barker <A HREF="http://www.nabble.com/Map-File-Problem--td5541099.html#a5544967">http://www.nabble.com/Map-File-Problem--td5541099.html#a5544967</A><BR>

<BR>

But it is possible or the only solution is to create a tif file for each band ?<BR>

<BR>

Thanks<BR>

<BR>

Vincent<BR>

<BR>

<BR>

<BR>

</BODY>

</HTML>