Turn on cycle ice function and bug fix for update GVF. (NOAA-EMC#246)

* Turn on cycle ice function and bug fix for update GVF. * Bug fix and enhance for retro runs.
EricJames-NOAA · Nov 3, 2021 · 808286d · 808286d
1 parent 87a9eb3
commit 808286d
Show file tree

Hide file tree

Showing 4 changed files with 10 additions and 4 deletions.
diff --git a/scripts/exregional_run_prepstart.sh b/scripts/exregional_run_prepstart.sh
@@ -356,9 +356,9 @@ if [ ${SFC_CYC} -eq 1 ] || [ ${SFC_CYC} -eq 2 ] ; then  # cycle surface fields
             ncks --append geolonlat.nc sfc_data.tile7.halo0.nc
             ncrename -v tslb,stc -v smois,smc -v sh2o,slc sfc_data.tile7.halo0.nc
           else
-            mv sfc_data.nc old.sfc_data.nc
+            mv sfc_data.nc gfsice.sfc_data.nc
             mv ${restart_prefix_find}sfc_data.nc sfc_data.nc
-#            ${EXECDIR}/cycle_ice.exe > stdout_cycleICE 2>&1
+            ${EXECDIR}/cycle_ice.exe > stdout_cycleICE 2>&1
           fi
           echo "cycle surface with ${checkfile}" > cycle_surface.done
           if [ ${SAVE_CYCLE_LOG} == "TRUE" ] ; then
@@ -376,8 +376,8 @@ fi
 #
 #-----------------------------------------------------------------------
 if [ ${HH} -eq ${GVF_update_hour} ] && [ ${cycle_type} == "spinup" ]; then
-   latestGVF=$(ls ${GVFOBS_PATH}/GVF-WKL-GLB_v2r3_npp_s*_e${YYYYMMDDm1}_c${YYYYMMDD}*.grib2)
-   latestGVF2=$(ls ${GVFOBS_PATH}/GVF-WKL-GLB_v2r3_npp_s*_e${YYYYMMDDm2}_c${YYYYMMDDm1}*.grib2)
+   latestGVF=$(ls ${GVF_ROOT}/GVF-WKL-GLB_v2r3_npp_s*_e${YYYYMMDDm1}_c${YYYYMMDD}*.grib2)
+   latestGVF2=$(ls ${GVF_ROOT}/GVF-WKL-GLB_v2r3_npp_s*_e${YYYYMMDDm2}_c${YYYYMMDDm1}*.grib2)
    if [ ! -r "${latestGVF}" ]; then
      if [ -r "${latestGVF2}" ]; then
        latestGVF=${latestGVF2}

diff --git a/ush/config.sh.RRFS_CONUS_3km b/ush/config.sh.RRFS_CONUS_3km
@@ -83,6 +83,9 @@ for i in {0..23..3}; do FCST_LEN_HRS_CYCLES[$i]=12; done
 DA_CYCLE_INTERV="1"
 RESTART_INTERVAL="1"
 
+SST_update_hour=01
+GVF_update_hour=04
+
 NCORES_RUN_ANAL=240
 HYBENSMEM_NMIN=66
 HALO_BLEND=20

diff --git a/ush/set_rrfs_config.sh b/ush/set_rrfs_config.sh
@@ -47,6 +47,8 @@ if [[ $DO_RETRO == "TRUE" ]] ; then
     ENKF_FCST=/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/enkf/atm
     AIRCRAFT_REJECT="/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/amdar_reject_lists"
     SFCOBS_USELIST="/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/mesonet_uselists"
+    SST_ROOT="/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/May_2021_data/highres_sst"
+    GVF_ROOT="/mnt/lfs4/BMC/wrfruc/Ruifang.Li/data/May_2021_data/gvf/grib2"
   fi
   if [[ $MACHINE == "hera" ]] ; then
 

diff --git a/ush/templates/FV3LAM_wflow.xml b/ush/templates/FV3LAM_wflow.xml
@@ -615,6 +615,7 @@ MODULES_RUN_TASK_FP script.
 {%- if do_retro %}
            <datadep age="00:00:05:00"><cyclestr>&FG_ROOT;/@Y@m@d@H/ics/gfs_data.tile7.halo0.nc</cyclestr></datadep>
            <datadep age="00:00:05:00"><cyclestr offset="-{{ extrn_mdl_ics_offset_hrs }}:00:00">&FG_ROOT;/@Y@m@d@H/lbcs/gfs_bndy.tile7.001.nc</cyclestr></datadep>
+           <datadep age="00:00:10:00"><cyclestr offset="00:00:00">&FG_ROOT;/@Y@m@d@H/fcst_fv3lam/INPUT/gfs_ctrl.nc</cyclestr></datadep>
 {% else %}
            <and>
              <timedep><cyclestr offset="&START_TIME_SPINUP;">@Y@m@d@H@M00</cyclestr></timedep>