Bugfix #2979 develop MTD Grid

internal/test_unit/python/unit.py

@@ -234,7 +234,7 @@ def unit(test_xml, file_log=None, cmd_only=False, noexit=False, memchk=False, ca
 
             #   # on failure, print the problematic test and exit, if requested
             if not (ret_ok and out_ok):
-                logger.info("\n".join(set_envs) + cmd + cmd_outs + "\n".join(unset_envs) + "\n")
+                logger.info("\n".join(set_envs) + "\n" + cmd + "\n" + cmd_outs + "\n".join(unset_envs) + "\n")
                 if not noexit:
                     sys.exit(1)
 

internal/test_util/tools/other/mode_time_domain/Makefile.am

@@ -38,8 +38,6 @@ test_velocity_LDADD = \
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_partition.o \
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_read_data.o \
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_txt_output.o \
-	${top_builddir}/src/tools/other/mode_time_domain/mtd-nc_grid.o \
-	${top_builddir}/src/tools/other/mode_time_domain/mtd-nc_utils_local.o \
 	-lvx_pxm \
 	-lvx_plot_util \
 	-lvx_nav \

internal/test_util/tools/other/mode_time_domain/Makefile.in

@@ -125,8 +125,6 @@ test_velocity_DEPENDENCIES = ${top_builddir}/src/tools/other/mode_time_domain/mt
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_partition.o \
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_read_data.o \
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_txt_output.o \
-	${top_builddir}/src/tools/other/mode_time_domain/mtd-nc_grid.o \
-	${top_builddir}/src/tools/other/mode_time_domain/mtd-nc_utils_local.o \
 	$(am__DEPENDENCIES_1) $(am__DEPENDENCIES_1) \
 	$(am__DEPENDENCIES_1) $(am__DEPENDENCIES_1) \
 	$(am__DEPENDENCIES_1) $(am__DEPENDENCIES_1)
@@ -378,8 +376,6 @@ test_velocity_LDADD = \
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_partition.o \
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_read_data.o \
 	${top_builddir}/src/tools/other/mode_time_domain/mtd-mtd_txt_output.o \
-	${top_builddir}/src/tools/other/mode_time_domain/mtd-nc_grid.o \
-	${top_builddir}/src/tools/other/mode_time_domain/mtd-nc_utils_local.o \
 	-lvx_pxm \
 	-lvx_plot_util \
 	-lvx_nav \

src/basic/vx_math/is_bad_data.h

@@ -80,6 +80,14 @@ inline int is_eq(float a, float b) {
    return is_eq((double)a, (double)b);
 }
 
+inline int is_eq(double a, float b) {
+   return is_eq(a, (double)b);
+}
+
+inline int is_eq(float a, double b) {
+   return is_eq((double)a, b);
+}
+
 template <typename T>
 inline int is_eq(T a, T b) {
    return (a == b);

src/libcode/vx_statistics/compute_stats.cc

@@ -1593,7 +1593,7 @@ void compute_aggregated_seeps_grid(const DataPlane &fcst_dp, const DataPlane &ob
          float obs_value = obs_dp.get(ix, iy);
          fcst_cat = obs_cat = bad_data_int;
          seeps_score = bad_data_double;
-         if (!is_eq(fcst_value, -9999.0) && !is_eq(obs_value, -9999.0)) {
+         if (!is_bad_data(fcst_value) && !is_bad_data(obs_value)) {
             SeepsScore *seeps_mpr = seeps_climo->get_record(ix, iy, fcst_value, obs_value);
             if (seeps_mpr != nullptr) {
                fcst_cat = seeps_mpr->fcst_cat;

src/tools/other/mode_time_domain/2d_att.cc

@@ -7,7 +7,6 @@
 // *=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*
 
 
-
 ////////////////////////////////////////////////////////////////////////
 
 
@@ -399,15 +398,7 @@ a.ObjectNumber = obj_number;
 
 moments = mask_2d.calc_2d_moments();
 
-if ( moments.N == 0 )  {
-
-   // mlog << Error << "\n\n  calc_2d_single_atts() -> empty object!\n\n";
-
-   // exit ( 1 );
-
-   return a;
-
-}
+if ( moments.N == 0 )  return a;
 
 a.Xbar = (moments.Sx)/(moments.N);
 a.Ybar = (moments.Sy)/(moments.N);
@@ -430,7 +421,8 @@ values = new float [a.Area];
 
 if ( !values )  {
 
-   mlog << Error << "\n\n  calc_2d_single_atts() -> memory allocation error\n\n";
+   mlog << Error << "\ncalc_2d_single_atts() -> "
+        << "memory allocation error\n\n";
 
    exit ( 1 );
 
@@ -626,6 +618,3 @@ return;
 ////////////////////////////////////////////////////////////////////////
 
 
-
-
-

src/tools/other/mode_time_domain/2d_att_array.cc

@@ -7,20 +7,6 @@
 // *=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*
 
 
-
-////////////////////////////////////////////////////////////////////////
-
-
-   //
-   //  Warning:  This file is machine generated
-   //
-   //            Do not edit by hand
-   //
-   //
-   //  Created by arraygen on September 17, 2015   10:11 am   MDT
-   //
-
-
 ////////////////////////////////////////////////////////////////////////
 
 
@@ -127,13 +113,10 @@ void SingleAtt2DArray::clear()
 if ( e )  { delete [] e;  e = (SingleAtt2D *) nullptr; }
 
 
-
 Nelements = 0;
 
 Nalloc = 0;
 
-// AllocInc = 50;   //  don't reset AllocInc
-
 
 return;
 
@@ -174,7 +157,8 @@ SingleAtt2D * u = new SingleAtt2D [N];
 
 if ( !u )  {
 
-   mlog << Error << "SingleAtt2DArray::extend(int) -> memory allocation error\n\n";
+   mlog << Error << "\nSingleAtt2DArray::extend(int) -> "
+        << "memory allocation error\n\n";
 
    exit ( 1 );
 
@@ -238,7 +222,8 @@ void SingleAtt2DArray::set_alloc_inc(int N)
 
 if ( N < 0 )  {
 
-   mlog << Error << "SingleAtt2DArray::set_alloc_int(int) -> bad value ... " << N << "\n\n";
+   mlog << Error << "\nSingleAtt2DArray::set_alloc_int(int) -> "
+        << "bad value ... " << N << "\n\n";
 
    exit ( 1 );
 
@@ -299,7 +284,8 @@ SingleAtt2D & SingleAtt2DArray::operator[](int N) const
 
 if ( (N < 0) || (N >= Nelements) )  {
 
-   mlog << Error << "\n\n  SingleAtt2DArray::operator[](int) -> range check error ... " << N << "\n\n";
+   mlog << Error << "\nSingleAtt2DArray::operator[](int) -> "
+        << "range check error ... " << N << "\n\n";
 
    exit ( 1 );
 }
@@ -344,8 +330,8 @@ unixtime SingleAtt2DArray::valid_time(int index) const
 
 if ( (index < 0) || (index >= Nelements) )  {
 
-   mlog << Error
-        << "SingleAtt2DArray::valid_time(int) const -> range check error\n\n";
+   mlog << Error << "\nSingleAtt2DArray::valid_time(int) const -> "
+        << "range check error\n\n";
 
    exit ( 1 );
 
@@ -365,8 +351,8 @@ int SingleAtt2DArray::lead_time(int index) const
 
 if ( (index < 0) || (index >= Nelements) )  {
 
-   mlog << Error
-        << "SingleAtt2DArray::lead_time(int) const -> range check error\n\n";
+   mlog << Error << "\nSingleAtt2DArray::lead_time(int) const -> "
+        << "range check error\n\n";
 
    exit ( 1 );
 
@@ -386,8 +372,8 @@ int SingleAtt2DArray::time_index(int index) const
 
 if ( (index < 0) || (index >= Nelements) )  {
 
-   mlog << Error
-        << "SingleAtt2DArray::time_index(int) const -> range check error\n\n";
+   mlog << Error << "\nSingleAtt2DArray::time_index(int) const -> "
+        << "range check error\n\n";
 
    exit ( 1 );
 
@@ -401,5 +387,3 @@ return e[index].time_index();
 ////////////////////////////////////////////////////////////////////////
 
 
-
-

src/tools/other/mode_time_domain/2d_moments.cc

@@ -7,7 +7,6 @@
 // *=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*
 
 
-
 ////////////////////////////////////////////////////////////////////////
 
 
@@ -25,12 +24,10 @@ using namespace std;
 
 ////////////////////////////////////////////////////////////////////////
 
-
    //
    //  Code for class Mtd_2D_Moments
    //
 
-
 ////////////////////////////////////////////////////////////////////////
 
 
@@ -177,7 +174,8 @@ void Mtd_2D_Moments::centralize()
 
 if ( N == 0 )  {
 
-   mlog << Error << "\n\n  Mtd_2D_Moments::centralize() -> no data!\n\n";
+   mlog << Error << "\nMtd_2D_Moments::centralize() -> "
+        << "no data!\n\n";
 
    exit ( 1 );
 
@@ -220,7 +218,8 @@ double Mtd_2D_Moments::calc_2D_axis_plane_angle() const
 
 if ( ! IsCentralized )  {
 
-   mlog << Error << "\n\n  Mtd_2D_Moments::calc_2D_axis_plane_angle() const -> moments must be centralized first!\n\n";
+   mlog << Error << "\nMtd_2D_Moments::calc_2D_axis_plane_angle() const -> "
+        << "moments must be centralized first!\n\n";
 
    exit ( 1 );
 
@@ -242,6 +241,3 @@ return angle;
 ////////////////////////////////////////////////////////////////////////
 
 
-
-
-

src/tools/other/mode_time_domain/2d_moments.h

@@ -57,7 +57,6 @@ class Mtd_2D_Moments {
 
       int area() const;
 
-
          //
          //  do stuff
          //
@@ -66,11 +65,8 @@ class Mtd_2D_Moments {
 
       void centralize();
 
-
       double calc_2D_axis_plane_angle() const;
 
-
-
 };