Merge remote-tracking branch 'remotes/origin/develop' into feature/ha…

…fsv0.3_phase2

parm/hafs_regional_da_C192s1n4.conf

@@ -0,0 +1,93 @@
+# This is a UNIX conf file that contains ALL information relating to
+# the HAFS C192 regional configuration.
+# Note: These are just example settings. Of course, when model resolution
+# changes, one should also modify the model time step and other related
+# dynamics and physics options, so that they work fine for the desired model
+# resolution.
+[config]
+# Specify the domain center Latitude and Longitude
+domlat=22.0
+domlon=-62.0
+# Same as domlat and domlon but vitals formatted
+vitalsdomlat=220N
+vitalsdomlon=620W
+
+# HAFS launcher requires vitals and a storm.
+# This is a default minimum vitals file for a fake storm.
+# The format of the missing value fields was based on the log
+# output returned when running the launcher after the vitals have
+# been "cleaned" up. That is, if  you enter the fields as all -999 values,
+# the "cleaned" output is returned. In essence I'm treating the vitals
+# module as a factory in the sense that it knows the correct format.
+# NHC 00L FAKE 20120710 0000 250N 0800W -99 -99 -999 -999 -099 -9 -99 -999 -999 -999 -999 M
+# fakestormid is a config option set in the launcher and is the
+# last storm id in the list of storms.
+fakestormid=00L
+fakestormname=NATL
+fakestormorg=NHC
+fakestorm_vitals={fakestormorg} {fakestormid} {fakestormname} {YMD} {HH}{min} {vitalsdomlat} {vitalsdomlon} -99 -99 -999 -999 -099 -9 -99 -999 -999 -999 -999 M
+
+[dir]
+PARMforecast={PARMhafs}/forecast/regional
+
+[grid]
+CASE=C192             ;; FV3 resolution
+LEVS=65              ;; Model vertical levels: 65
+gtype=regional       ;; grid type: uniform, stretch, nest, or regional
+# If gridfixdir is provided and the dir exists, then use the pre-generated static grid fix files under gridfixdir
+#gridfixdir={FIXhafs}/fix_fv3/{CASE}
+gridfixdir=/let/hafs_grid/generate/grid
+# Otherwise, will generate the model grid according to the following grid parameters
+# Need for grid types: stretch, nest and regional
+stretch_fac=1.0001   ;; Stretching factor for the grid
+target_lon={domlon}  ;; center longitude of the highest resolution tile
+target_lat={domlat}  ;; center latitude of the highest resolution tile
+# Need for grid types: nest and regional
+refine_ratio=4       ;; Specify the refinement ratio for nest grid
+istart_nest=12
+jstart_nest=42
+iend_nest=371
+jend_nest=341
+
+[grid_ens]
+CASE_ENS=C192            ;; FV3 resolution
+LEVS_ENS=65              ;; Model vertical levels: 65
+gtype_ens=regional       ;; grid type: uniform, stretch, nest, or regional
+gridfixdir_ens=/let/hafs_grid/generate/grid_ens
+stretch_fac_ens=1.0001   ;; Stretching factor for the grid
+target_lon_ens={domlon}  ;; center longitude of the highest resolution tile
+target_lat_ens={domlat}  ;; center latitude of the highest resolution tile
+refine_ratio_ens=4       ;; Specify the refinement ratio for nest grid
+istart_nest_ens=12
+jstart_nest_ens=42
+iend_nest_ens=371
+jend_nest_ens=341
+
+[forecast]
+ccpp_suite_regional=FV3_HAFS_v0_gfdlmp_tedmf
+
+restart_interval="3 6 9"  ;; restart interval in hours for the forecast
+
+layoutx=40
+layouty=30
+npx=721
+npy=601
+npz=64
+
+write_groups=2
+write_tasks_per_group=48
+
+output_grid=regional_latlon
+output_grid_cen_lon=-62.0         ;; central longitude
+output_grid_cen_lat=22.0          ;; central latitude
+output_grid_lon_span=109.8        ;; output domain span for longitude in regular latlon coordinate system (in degrees)
+output_grid_lat_span=72.9         ;; output domain span for latitude in regular latlon coordinate system (in degrees)
+output_grid_dlon=0.1              ;; output grid spacing dlon . . . .
+output_grid_dlat=0.1              ;; output grid spacing dlat . . . .
+
+[forecast_ens]
+restart_interval_ens="3 6 9"  ;; restart interval in hours for the ensda member forecast
+
+[rocotostr]
+FORECAST_RESOURCES=FORECAST_RESOURCES_regional_{forecast/layoutx}x{forecast/layouty}io{forecast/write_groups}x{forecast/write_tasks_per_group}_omp2
+

rocoto/cronjob_hafs_rt.sh

@@ -98,6 +98,16 @@ scrubopt="config.scrub_work=no config.scrub_com=no"
 
 #===============================================================================
 
+# 3DEnVar with GDAS ensembles
+ ${PYTHON3} ./run_hafs.py -t ${dev} 2020082506-2020082512 00L HISTORY \
+     config.EXPT=${EXPT} config.SUBEXPT=${EXPT}_C192s1n4_3denvar \
+     config.run_gsi=yes config.run_envar=yes \
+     gsi.use_bufr_nr=yes \
+     config.NHRS=12 ${scrubopt} \
+     ../parm/hafs_regional_da_C192s1n4.conf
+
+#===============================================================================
+
 date
 
 echo 'cronjob done'

rocoto/hafs_rt_status.sh

@@ -63,15 +63,15 @@ for file in *.xml; do
 
   # Check if rocoto completion task ran successfully or not
 
-  if [[ $if_complete == "1" ]]; then
+  if [ $if_complete == "1" ] || [ $if_complete == "2" ]; then
     echo "ROCOTO SAYS COMPLETION TASK SUCCEEDED"
   else
     echo "ROCOTO SAYS COMPLETION TASK DID NOT SUCCEED"
   fi
 
   # Check the post and product log files
 
-  if [[ $if_complete == "1" ]]; then
+  if [ $if_complete == "1" ] || [ $if_complete == "2" ]; then
     post_log=`cat ${HAFS_out}/${subexpt}/${storm_init}/${sid}/hafs_atm_post.log|grep "post job done"|tail -1`
     prod_log=`cat ${HAFS_out}/${subexpt}/${storm_init}/${sid}/hafs_product.log|grep "product job done"|tail -1`
     if [[ $post_log == "post job done" ]]; then
@@ -96,7 +96,7 @@ for file in *.xml; do
 
   # Check if everything passed
 
-  if [[ $if_complete == "1" ]]; then
+  if [ $if_complete == "1" ] || [ $if_complete == "2" ]; then
     if [[ -e ${storm1_done} && -e ${atcfunix} && -e ${hafsprs_synoptic} ]]; then
     if [[ $post_log == "post job done" ]]; then
     if [[ $prod_log == "product job done" ]]; then