Per #1905, add a call to ensemble_stat using the -ctrl option. I real…

…ize that if a control member is specified, the number of ensemble members should be increased by 1.

test/xml/unit_ensemble_stat.xml

@@ -122,6 +122,34 @@
     </output>
   </test>
 
+  <test name="ensemble_stat_MASK_SID_CTRL">
+    <exec>&MET_BIN;/ensemble_stat</exec>
+    <env>
+      <pair><name>CENSOR_THRESH</name>  <value></value></pair>
+      <pair><name>CENSOR_VAL</name>     <value></value></pair>
+      <pair><name>SKIP_CONST</name>     <value>FALSE</value></pair>
+      <pair><name>OUTPUT_PREFIX</name>  <value>MASK_SID_CTRL</value></pair>
+      <pair><name>CONFIG_DIR</name>     <value>&CONFIG_DIR;</value></pair>
+    </env>
+    <param> \
+      5 \
+      &DATA_DIR_MODEL;/grib1/arw-fer-gep5/arw-fer-gep5_2012040912_F024.grib \
+      &DATA_DIR_MODEL;/grib1/arw-sch-gep2/arw-sch-gep2_2012040912_F024.grib \
+      &DATA_DIR_MODEL;/grib1/arw-sch-gep6/arw-sch-gep6_2012040912_F024.grib \
+      &DATA_DIR_MODEL;/grib1/arw-tom-gep3/arw-tom-gep3_2012040912_F024.grib \
+      &DATA_DIR_MODEL;/grib1/arw-tom-gep7/arw-tom-gep7_2012040912_F024.grib \
+      &CONFIG_DIR;/EnsembleStatConfig_MASK_SID \
+      -ctrl &DATA_DIR_MODEL;/grib1/arw-fer-gep1/arw-fer-gep1_2012040912_F024.grib \
+      -point_obs &OUTPUT_DIR;/ascii2nc/gauge_2012041012_24hr.nc \
+      -outdir &OUTPUT_DIR;/ensemble_stat -v 1
+    </param>
+    <output>
+      <stat>&OUTPUT_DIR;/ensemble_stat/ensemble_stat_MASK_SID_CTRL_20120410_120000V.stat</stat>
+      <stat>&OUTPUT_DIR;/ensemble_stat/ensemble_stat_MASK_SID_CTRL_20120410_120000V_orank.txt</stat>
+      <grid_nc>&OUTPUT_DIR;/ensemble_stat/ensemble_stat_MASK_SID_CTRL_20120410_120000V_ens.nc</grid_nc>
+    </output>
+  </test>
+
   <test name="ensemble_stat_MASK_SID_CENSOR">
     <exec>&MET_BIN;/ensemble_stat</exec>
     <env>