print of average rel err. on demand for some toys

examples/toy_icon_atm.c

@@ -400,6 +400,51 @@ int main (int argc, char *argv[]) {
       cell_in_patch_err_rel[i] = fabs(cell_in_patch_err_abs[i] / cell_data_field[i]);
     }
 
+#define DEBUG
+#ifdef DEBUG
+
+    double f1err_1 = 0.0;
+    double f1err_100 = 0.0;
+    int if1err_1 = 0;
+    int if1err_100 = 0;
+
+    double f2err_1 = 0.0;
+    double f2err_100 = 0.0;
+    int if2err_1 = 0;
+    int if2err_100 = 0;
+
+    for ( int i = 0; i < num_cells; ++i ) {
+
+      if ( cell_in_conserv_err_rel[i] < 100.0 ) {
+    	  f1err_1 += cell_in_conserv_err_rel[i];
+    	  if1err_1++;
+      } else {
+    	  f1err_100 += cell_in_conserv_err_rel[i];
+    	  if1err_100++;
+      }
+
+      if ( cell_in_patch_err_rel[i] < 100.0 ) {
+    	  f2err_1 += cell_in_patch_err_rel[i];
+    	  if2err_1++;
+      } else {
+    	  f2err_100 += cell_in_patch_err_rel[i];
+    	  if2err_100++;
+      }
+    }
+
+    f1err_1 = f1err_1 / (double)if1err_1;
+    f2err_1 = f2err_1 / (double)if2err_1;
+
+    f1err_100 = f1err_100 / (double)if1err_100;
+    f2err_100 = f2err_100 / (double)if2err_100;
+
+    printf ("avg. rel err. for 1st field %f\n", f1err_1);
+    printf ("avg. rel err. for 2nd field %f\n", f2err_1);
+
+    printf ("extreme avg. rel err. for 1st field %f\n", f1err_100);
+    printf ("extreme avg. rel err. for 2nd field %f\n", f2err_100);
+#endif
+
     //----------------------------------------------------------
     // write field to vtk output file
     //----------------------------------------------------------

examples/toy_icon_ocn.c

@@ -402,6 +402,51 @@ int main (int argc, char *argv[]) {
       cell_in_patch_err_rel[i] = fabs(cell_in_patch_err_abs[i] / cell_data_field[i]);
     }
 
+#define DEBUG
+#ifdef DEBUG
+
+    double f1err_1 = 0.0;
+    double f1err_100 = 0.0;
+    int if1err_1 = 0;
+    int if1err_100 = 0;
+
+    double f2err_1 = 0.0;
+    double f2err_100 = 0.0;
+    int if2err_1 = 0;
+    int if2err_100 = 0;
+
+    for ( int i = 0; i < num_cells; ++i ) {
+
+      if ( cell_in_conserv_err_rel[i] < 100.0 ) {
+    	  f1err_1 += cell_in_conserv_err_rel[i];
+    	  if1err_1++;
+      } else {
+    	  f1err_100 += cell_in_conserv_err_rel[i];
+    	  if1err_100++;
+      }
+
+      if ( cell_in_patch_err_rel[i] < 100.0 ) {
+    	  f2err_1 += cell_in_patch_err_rel[i];
+    	  if2err_1++;
+      } else {
+    	  f2err_100 += cell_in_patch_err_rel[i];
+    	  if2err_100++;
+      }
+    }
+
+    f1err_1 = f1err_1 / (double)if1err_1;
+    f2err_1 = f2err_1 / (double)if2err_1;
+
+    f1err_100 = f1err_100 / (double)if1err_100;
+    f2err_100 = f2err_100 / (double)if2err_100;
+
+    printf ("avg. rel err. for 1st field %f\n", f1err_1);
+    printf ("avg. rel err. for 2nd field %f\n", f2err_1);
+
+    printf ("extreme avg. rel err. for 1st field %f\n", f1err_100);
+    printf ("extreme avg. rel err. for 2nd field %f\n", f2err_100);
+#endif
+
     //----------------------------------------------------------
     // write field to vtk output file
     //----------------------------------------------------------