Per #2395, write the new VCNT columns to the output and document the …

…additions to the VL1L2, VAL1L2, and VCNT columns.
dtcenter · Feb 13, 2024 · 1d387ad · 1d387ad
1 parent e58899e
commit 1d387ad
Show file tree

Hide file tree

Showing 2 changed files with 55 additions and 2 deletions.
diff --git a/docs/Users_Guide/point-stat.rst b/docs/Users_Guide/point-stat.rst
@@ -1285,6 +1285,15 @@ The first set of header columns are common to all of the output files generated
   * - 34
     - O_SPEED_BAR
     - Mean observed wind speed
+  * - 35
+    - DIR_ME
+    - Mean wind direction difference, from -180 to 180 degrees
+  * - 36
+    - DIR_MAE
+    - Mean absolute wind direction difference
+  * - 37
+    - DIR_MSE
+    - Mean squared wind direction difference
 
 .. _table_PS_format_info_VAL1L2:
 
@@ -1331,6 +1340,15 @@ The first set of header columns are common to all of the output files generated
   * - 34
     - OA_SPEED_BAR
     - Mean observed wind speed anomaly
+  * - 35
+    - DIRA_ME
+    - Mean wind direction anomaly difference, from -180 to 180 degrees
+  * - 36
+    - DIRA_MAE
+    - Mean absolute wind direction anomaly difference
+  * - 37
+    - DIRA_MSE
+    - Mean squared wind direction anomaly difference
 
 .. _table_PS_format_info_VCNT:
 
@@ -1410,6 +1428,18 @@ The first set of header columns are common to all of the output files generated
   * - 85-87
     - ANOM_CORR_UNCNTR, :raw-html:`<br />` ANOM_CORR_UNCNTR_BCL, :raw-html:`<br />` ANOM_CORR_UNCNTR_BCU
     - Uncentered vector Anomaly Correlation excluding mean error including bootstrap upper and lower confidence limits
+  * - 88-90
+    - DIR_ME, :raw-html:`<br />` DIR_ME_BCL, :raw-html:`<br />` DIR_ME_BCU
+    - Mean direction difference, from -180 to 180 degrees, including bootstrap upper and lower confidence limits
+  * - 91-93
+    - DIR_MAE, :raw-html:`<br />` DIR_MAE_BCL, :raw-html:`<br />` DIR_MAE_BCU
+    - Mean absolute direction difference including bootstrap upper and lower confidence limits
+  * - 94-96
+    - DIR_MSE, :raw-html:`<br />` DIR_MSE_BCL, :raw-html:`<br />` DIR_MSE_BCU
+    - Mean squared direction difference including bootstrap upper and lower confidence limits
+  * - 97-99
+    - DIR_RMSE, :raw-html:`<br />` DIR_RMSE_BCL, :raw-html:`<br />` DIR_RMSE_BCU
+    - Root mean squared direction difference including bootstrap upper and lower confidence limits
 
 .. _table_PS_format_info_MPR:
 

diff --git a/src/libcode/vx_stat_out/stat_columns.cc b/src/libcode/vx_stat_out/stat_columns.cc
@@ -3105,10 +3105,17 @@ void write_vcnt_cols(const VL1L2Info &vcnt_info, int i,
    //    SPD_ABSERR,       SPD_ABSERR_BCL,       SPD_ABSERR_BCU,
    //    DIR_ERR,          DIR_ERR_BCL,          DIR_ERR_BCU,
    //    DIR_ABSERR,       DIR_ABSERR_BCL,       DIR_ABSERR_BCU,
-   //    ANOM_CORR,        ANOM_CORR_NCL,        ANOM_CORR_NCU,       ANOM_CORR_BCL,   ANOM_CORR_BCU
-   //    ANOM_CORR_UNCNTR, ANOM_CORR_UNCNTR_BCL, ANOM_CORR_UNCNTR_BCU
+   //    ANOM_CORR,        ANOM_CORR_NCL,        ANOM_CORR_NCU,
+   //                      ANOM_CORR_BCL,        ANOM_CORR_BCU
+   //    ANOM_CORR_UNCNTR, ANOM_CORR_UNCNTR_BCL, ANOM_CORR_UNCNTR_BCU,
+   //    DIR_ME,           DIR_ME_BCL,           DIR_ME_BCU,
+   //    DIR_MAE,          DIR_MAE_BCL,          DIR_MAE_BCU,
+   //    DIR_MSE,          DIR_MSE_BCL,          DIR_MSE_BCU,
+   //    DIR_RMSE,         DIR_RMSE_BCL,         DIR_RMSE_BCU
    //
 
+   // TODO: MET #963 Compute the VCNT bootstrap CI's
+
    at.set_entry(r, c++, max(vcnt_info.vcount,          // TOTAL
                             vcnt_info.vacount));
 
@@ -3194,6 +3201,22 @@ void write_vcnt_cols(const VL1L2Info &vcnt_info, int i,
    at.set_entry(r, c++, (string)na_str);               // ANOM_CORR_UNCNTR_BCL
    at.set_entry(r, c++, (string)na_str);               // ANOM_CORR_UNCNTR_BCU
 
+   at.set_entry(r, c++, vcnt_info.DIR_ME.v);           // DIR_ME
+   at.set_entry(r, c++, (string)na_str);               // DIR_ME_BCL
+   at.set_entry(r, c++, (string)na_str);               // DIR_ME_BCU
+
+   at.set_entry(r, c++, vcnt_info.DIR_MAE.v);          // DIR_MAE
+   at.set_entry(r, c++, (string)na_str);               // DIR_MAE_BCL
+   at.set_entry(r, c++, (string)na_str);               // DIR_MAE_BCU
+
+   at.set_entry(r, c++, vcnt_info.DIR_MSE.v);          // DIR_MSE
+   at.set_entry(r, c++, (string)na_str);               // DIR_MSE_BCL
+   at.set_entry(r, c++, (string)na_str);               // DIR_MSE_BCU
+
+   at.set_entry(r, c++, vcnt_info.DIR_RMSE.v);         // DIR_RMSE
+   at.set_entry(r, c++, (string)na_str);               // DIR_RMSE_BCL
+   at.set_entry(r, c++, (string)na_str);               // DIR_RMSE_BCU
+
    return;
 }