[netcdfgroup] NetCdf on Panoply

• To: netcdfgroup@xxxxxxxxxxxxxxxx
• Subject: [netcdfgroup] NetCdf on Panoply
• From: Vincent THIERION <vincent.thierion@xxxxxx>
• Date: Tue, 22 Jan 2008 16:08:05 +0100

Hello,

I created netcdf file to store and vizualize radar rainfall data (time serie). 
Radar data is store in two-dimensional matrix with X and Y axis corresponding 
to Lambert Conformal Conic projection. 

I have created all data axis (Lambert Conformal Conic coordinates) and stored 
them in one-dimensional matrix.

I can vizualize it thanks to ncBrowse software. 
However, in using Panoply software, I can't vizualize it. I can read on 
http://www.giss.nasa.gov/tools/panoply/help/ref_datafiles.html

"Lambert conformal conic: Panoply recognizes a Lambert conformal conic 
data..."

Is it necessary to create data axis in longitude/latitude based on the LCC 
coordinates ? Is correspondance performed automatically ?

This is the corresponding CDL:

netCDF file "pluie2.nc"
--------
netcdf file:/home/vincent/pluie2.nc {
 dimensions:
   Date = UNLIMITED;   // (3 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 = 600000.0; // double
     :false_northing = 3200000.0; // double
     :_CoordinateTransformType = "Projection";
     :_CoordinateSystems = "ProjectionCoordinateSystem";
     :_CoordinateAxes = "y x";
   double pluie(Date=3, x=50, y=50);
     :units = "mm";
     :long_name = "Cumul de pluie sur 5 minutes";
     :coordinates = "lat lon";
     :grid_mapping = "Lambert_Conformal";
     :_CoordinateSystems = "ProjectionCoordinateSystem 
LatLonCoordinateSystem";
   char LatLonCoordinateSystem;
     :_CoordinateAxes = "Date lat lon";
   char ProjectionCoordinateSystem;
     :_CoordinateAxes = "Date x y";
     :_CoordinateTransforms = "LambertConformalProjection";
   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";
     :standard_name = "latitude";
     :_CoordinateAxisType = "Lat";
   double lon(x=50, y=50);
     :units = "degrees_east";
     :long_name = "longitude";
     :standard_name = "longitude";
     :_CoordinateAxisType = "Lon";
   double Date(Date=3);
     :long_name = "Date et heure de l'image Radar";
     :units = "seconds since 1970-01-01 01:00:00";
     :axis = "T";
     :_CoordinateAxisType = "Time";

 :X0 = 742.7; // double
 :Y0 = 3232.1; // double
}

Thanks a lot for your help.
-- 
THIERION Vincent                
Eleve-chercheur
LGEI - Ecole des Mines d'Alès
Tel : 04-66-78-27-62
Mail : vincent.thierion@xxxxxx