preparation for reduce_dim of z axis.

src/cdf_write.c

@@ -986,6 +986,7 @@ void cdf_write_var_data(int fileID, int vlistID, int varID, int ncvarid, int dty
 void cdf_write_var(stream_t *streamptr, int varID, int memtype, const void *data, int nmiss)
 {
   if ( streamptr->accessmode == 0 ) cdfEndDef(streamptr);
+  // if ( streamptr->ncmode == 0 ) { cdfEndDef(streamptr); streamptr->ncmode = 2; }
 
   size_t xsize = 0, ysize = 0;
   size_t size;
@@ -1081,6 +1082,7 @@ void cdf_write_var_chunk(stream_t *streamptr, int varID, int memtype,
                          const int rect[][2], const void *data, int nmiss)
 {
   if ( streamptr->accessmode == 0 ) cdfEndDef(streamptr);
+  // if ( streamptr->ncmode == 0 ) { cdfEndDef(streamptr); streamptr->ncmode = 2; }
 
   int xid = CDI_UNDEFID, yid = CDI_UNDEFID;
   size_t xsize = 0, ysize = 0;
@@ -1182,6 +1184,7 @@ void cdf_write_var_chunk(stream_t *streamptr, int varID, int memtype,
 void cdf_write_var_slice(stream_t *streamptr, int varID, int levelID, int memtype, const void *data, int nmiss)
 {
   if ( streamptr->accessmode == 0 ) cdfEndDef(streamptr);
+  //  if ( streamptr->ncmode == 0 ) { streamptr->ncmode = 2; cdfEndDef(streamptr); }
 
   size_t xsize = 0, ysize = 0;
   size_t start[5];

src/stream_cdf_o.c

@@ -1495,20 +1495,20 @@ void cdfDefZaxis(stream_t *streamptr, int zaxisID)
   size_t dimlen = (size_t)zaxisInqSize(zaxisID);
   int type = zaxisInqType(zaxisID);
 
-  bool is_scalar = false;
+  int ndims = 1;
+
   if ( dimlen == 1 )
     {
-      is_scalar = zaxisInqScalar(zaxisID) > 0;
+      bool is_scalar = zaxisInqScalar(zaxisID) > 0;
       if ( !is_scalar && CDI_cmor_mode )
         {
           is_scalar = true;
           zaxisDefScalar(zaxisID);
         }
-    }
 
-  int ndims = is_scalar ? 0 : 1;
+      if ( is_scalar ) ndims = 0;
+      // if ( CDI_reduce_dim ) return;
 
-  if ( dimlen == 1 )
       switch (type)
         {
         case ZAXIS_SURFACE:
@@ -1526,6 +1526,7 @@ void cdfDefZaxis(stream_t *streamptr, int zaxisID)
         case ZAXIS_MIX_LAYER:
           return;
         }
+    }
 
   zaxisInqName(zaxisID, axisname);