Feature 1815 level value

docs/Users_Guide/appendixA.rst

@@ -147,9 +147,10 @@ Q. How do I choose a time slice in a NetCDF file?
 
 A.
 When processing NetCDF files, the level information needs to be
-specified to tell MET which 2D slice of data to use. There is
-currently no way to explicitly define which time slice to use
-other than selecting the time index.
+specified to tell MET which 2D slice of data to use. The index is selected from
+a value when it starts with "@" for vertical level (pressure or height)
+and time. The actual time, @YYYYMMDD_HHMM, is allowed instead of selecting
+the time index.
 
 Let's use plot_data_plane as an example:
 
@@ -160,6 +161,11 @@ Let's use plot_data_plane as an example:
 		obs.ps \ 
 		'name="APCP"; level="(5,*,*)";'
 		
+		plot_data_plane \
+		gtg_obs_forecast.20130730.i00.f00.nc \
+		altitude_20000.ps \
+		'name = "edr"; level = "(@20130730_0000,@20000,*,*)";'
+		
 Assuming that the first array is the time, this will select the 6-th
 time slice of the APCP data and plot it since these indices are 0-based.
 

docs/Users_Guide/config_options.rst

@@ -991,7 +991,8 @@ File-format specific settings for the "field" entry:
 
       * (i,...,j,*,*) for a single field, where i,...,j specifies fixed
         dimension values and *,* specifies the two dimensions for the
-        gridded field. For example:
+        gridded field. @ specifies the vertical level value or time value
+        instead of offset, (i,...,@NNN,*,*). For example:
 
       .. code-block:: none
 
@@ -1006,6 +1007,17 @@ File-format specific settings for the "field" entry:
              }
            ];
 
+        field = [
+             {
+               name       = "QVAPOR";
+               level      = "(@20220601_1200,@850,*,*)";
+             },
+             {
+               name       = "TMP_P850_ENS_MEAN";
+               level      = [ "(*,*)" ];
+             }
+           ];
+
   * Python (using PYTHON_NUMPY or PYTHON_XARRAY):
 
     * The Python interface for MET is described in Appendix F of the MET

internal/test_unit/xml/unit_plot_data_plane.xml

@@ -273,6 +273,20 @@
     </output>
   </test>
 
+  <test name="plot_data_plane_NCCF_latlon_by_value">
+    <exec>&MET_BIN;/plot_data_plane</exec>
+    <param> \
+      &DATA_DIR_MODEL;/nccf/gtg/latlon/gtg_obs_forecast.20130730.i00.f00.nc \
+      &OUTPUT_DIR;/plot_data_plane/gtg_obs_forecast.20130730.i00.f00.NCCF_latlon_20000.ps \
+      'name = "edr"; level = "(0,@20000,*,*)";' \
+      -title "NCCF Latitude/Longitude Level 0" \
+      -v 1
+    </param>
+    <output>
+      <ps>&OUTPUT_DIR;/plot_data_plane/gtg_obs_forecast.20130730.i00.f00.NCCF_latlon_20000.ps</ps>
+    </output>
+  </test>
+
   <test name="plot_data_plane_NCCF_north_to_south">
     <exec>&MET_BIN;/plot_data_plane</exec>
     <param> \
@@ -292,7 +306,7 @@
     <param> \
       &DATA_DIR_MODEL;/nccf/GloTEC_TEC_2015_03_17.nc \
       &OUTPUT_DIR;/plot_data_plane/plot_data_plane_NCCF_time.ps \
-      'name="TEC"; level="(20150317_000500,*,*)"; file_type=NETCDF_NCCF;' \
+      'name="TEC"; level="(@20150317_000500,*,*)"; file_type=NETCDF_NCCF;' \
       -title "NCCF GloTEC - Total Electron Content" \
       -v 1
     </param>
@@ -403,7 +417,7 @@
     <param> \
       &DATA_DIR_MODEL;/easm/pr_day_MPI-ESM-MR_rcp85_r1i1p1_20060101-20091231.nc \
       &OUTPUT_DIR;/plot_data_plane/EaSM_CMIP5_pr_day_MPI-ESM-MR_rcp85_r1i1p1_20060101_12_time.ps \
-      'name="pr"; level="(20060102_000000,*,*)";' \
+      'name="pr"; level="(@20060102_000000,*,*)";' \
       -v 4
     </param>
     <output>

src/basic/vx_util/util_constants.h

@@ -119,6 +119,7 @@ static const double const_gop   = 9.80616; // from The Ceaseless Wind
 static const double const_rd    = 287.0;   // kg/k dry gas constant
 
 static const int vx_data2d_star = -12345;
+static const int vx_data2d_dim_by_value = -123456; // apply the value instead of offset for slicing
 
 ////////////////////////////////////////////////////////////////////////
 

src/libcode/vx_data2d/level_info.cc

@@ -70,7 +70,6 @@ void LevelInfo::init_from_scratch() {
 
    clear();
 
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
@@ -88,9 +87,8 @@ void LevelInfo::assign(const LevelInfo &l) {
    Upper   = l.upper();
    Lower   = l.lower();
    Increment = l.increment();
-   time_as_offset = l.is_time_as_offset();
+   Is_offset = l.is_offset();
 
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
@@ -106,9 +104,8 @@ void LevelInfo::clear() {
    Upper  = 0.0;
    Lower  = 0.0;
    Increment = 0.0;
-   time_as_offset = true;
+   Is_offset = true;
 
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
@@ -125,78 +122,67 @@ void LevelInfo::dump(ostream &out) const {
        << "  Upper     = " << Upper << "\n"
        << "  Increment = " << Increment << "\n";
 
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_type(LevelType lt) {
    Type = lt;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_type_num(int i) {
    TypeNum = i;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_req_name(const char *str) {
    ReqName = str;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_name(const char *str) {
    Name = str;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_units(const char *str) {
    Units = str;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_upper(double u) {
    Upper = u;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_lower(double l) {
    Lower = l;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_range(double l, double u) {
    Lower = l;
    Upper = u;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
 void LevelInfo::set_increment(double i) {
    Increment = i;
-   return;
 }
 
 ///////////////////////////////////////////////////////////////////////////////
 
-void LevelInfo::set_time_as_offset(bool b) {
-   time_as_offset = b;
-   return;
+void LevelInfo::set_is_offset(bool b) {
+   Is_offset = b;
 }
 
 ///////////////////////////////////////////////////////////////////////////////

src/libcode/vx_data2d/level_info.h

@@ -48,7 +48,7 @@ class LevelInfo
       double Upper;         // Upper level limit
       double Lower;         // Lower level limit
       double Increment;     // Increment (time: seconds, 0 for no increment)
-      bool   time_as_offset;// default: true, false: the (time) value instead
+      bool   Is_offset;     // default: true, false: the value instead
                             // of the offset at Lower and Upper
 
       void init_from_scratch();
@@ -77,7 +77,7 @@ class LevelInfo
       double       upper()       const;
       double       lower()       const;
       double       increment()   const;
-      bool         is_time_as_offset()const;
+      bool         is_offset()    const;
 
       //
       // set stuff
@@ -92,7 +92,7 @@ class LevelInfo
       void set_lower(double l);
       void set_range(double l, double u);
       void set_increment(double i);
-      void set_time_as_offset(bool b);
+      void set_is_offset(bool b);
 
 };
 
@@ -106,7 +106,7 @@ inline ConcatString LevelInfo::units()    const { return(Units);   }
 inline double       LevelInfo::upper()    const { return(Upper);   }
 inline double       LevelInfo::lower()    const { return(Lower);   }
 inline double       LevelInfo::increment()const { return(Increment);}
-inline bool         LevelInfo::is_time_as_offset()const { return(time_as_offset);}
+inline bool         LevelInfo::is_offset()const { return(Is_offset);}
 
 ///////////////////////////////////////////////////////////////////////////////