add support of GRIB level type MEANSEA

ChangeLog

 2008-12-??  Uwe Schulzweida  <Uwe.Schulzweida@zmaw.de>
 
+	* add support of GRIB level type MEANSEA
 	* cdf_inq_contents: check zaxis type "depth_blow_sea/land" (bug fix)
 
 2008-11-13  Uwe Schulzweida  <Uwe.Schulzweida@zmaw.de>

src/cdi.h

@@ -124,6 +124,7 @@ extern "C" {
 #define  ZAXIS_TRAJECTORY         9
 #define  ZAXIS_ALTITUDE          10
 #define  ZAXIS_SIGMA             11
+#define  ZAXIS_MEANSEA           12
 
 /* TAXIS types */
 

src/stream_cdf.c

@@ -2033,6 +2033,7 @@ void cdfDefZaxis(int streamID, int zaxisID)
   type   = zaxisInqType(zaxisID);
 
   if ( dimlen == 1 && type == ZAXIS_SURFACE ) return;
+  if ( dimlen == 1 && type == ZAXIS_MEANSEA ) return;
 
   zaxisInqName(zaxisID, axisname);
   /*

src/stream_grb.c

@@ -342,6 +342,11 @@ int gribGetZaxisType(int grb_ltype)
 	zaxistype = ZAXIS_SURFACE;
 	break;
       }
+    case LTYPE_MEANSEA:
+      {
+	zaxistype = ZAXIS_MEANSEA;
+	break;
+      }
     case LTYPE_99:
     case LTYPE_ISOBARIC:
       {
@@ -2029,6 +2034,13 @@ void grbDefLevel(int *isec1, int *isec2, double *fsec2, int zaxisID, int levelID
 	ISEC1_Level2    = 0;
 	break;
       }
+    case ZAXIS_MEANSEA:
+      {
+	ISEC1_LevelType = LTYPE_MEANSEA;
+	ISEC1_Level1    = (int) zaxisInqLevel(zaxisID, levelID);
+	ISEC1_Level2    = 0;
+	break;
+      }
     case ZAXIS_HYBRID:
     case ZAXIS_HYBRID_HALF:
       {

src/zaxis.c

@@ -30,6 +30,7 @@ ZaxistypeEntry[] = {
   {"lev",     "trajectory",               ""},
   {"alt",     "altitude",                 "m"},
   {"lev",     "sigma",                    "level"},
+  {"lev",     "meansea",                  "level"},
 };
 
 static int CDI_MaxZaxistype = sizeof(ZaxistypeEntry) / sizeof(ZaxistypeEntry[0]);