Ocean Ensmemble Options

NeilBarton-NOAA · Mar 7, 2024 · 267db5e · 267db5e
1 parent 622683d
commit 267db5e
Show file tree

Hide file tree

Showing 7 changed files with 55 additions and 35 deletions.
diff --git a/parm/config/gefs/config.base b/parm/config/gefs/config.base
@@ -260,10 +260,6 @@ export NMEM_ENS=@NMEM_ENS@
 export ENSMEM=${ENSMEM:-"000"}
 export MEMDIR="mem${ENSMEM}"
 
-# initialize ocean ensemble members with perturbations
-#   if true, only occurs for members greater than zero
-export OCN_ENS_PERTURB_FILES=false
-
 export DOIAU="NO"  # While we are not doing IAU, we may want to warm start w/ IAU in the future
 # Check if cycle is cold starting
 if [[ "${EXP_WARM_START}" = ".false." ]]; then

diff --git a/parm/config/gefs/config.efcs b/parm/config/gefs/config.efcs
@@ -60,6 +60,17 @@ export SPPT_TAU=21600.
 export SPPT_LSCALE=500000.
 export SPPT_LOGIT=".true."
 export SPPT_SFCLIMIT=".true."
+# OCN options
+export DO_OCN_SPPT="YES"
+export DO_OCN_PERT_EPBL="YES"
+export ODA_INCUPD="True"
+export ODA_TEMPINC_VAR='t_pert'
+export ODA_SALTINC_VAR='s_pert'
+export ODA_THK_VAR='h_anl'
+export ODA_UINC_VAR='u_pert'
+export ODA_VINC_VAR='v_pert'
+export ODA_INCUPD_NHOURS=0.0
+
 
 export restart_interval="${restart_interval_gfs}"
 

diff --git a/parm/config/gefs/config.stage_ic b/parm/config/gefs/config.stage_ic
@@ -4,6 +4,10 @@
 
 echo "BEGIN: config.stage_ic"
 
+# initialize ocean ensemble members with perturbations
+#   if true, only occurs for members greater than zero
+export STAGE_OCN_PERTURB_FILES=false
+
 # Get task specific resources
 source "${EXPDIR}/config.resources" stage_ic
 

diff --git a/scripts/exglobal_stage_ic.sh b/scripts/exglobal_stage_ic.sh
@@ -106,7 +106,7 @@ for MEMDIR in "${MEMDIR_ARRAY[@]}"; do
 
     # Ocean Perturbation Files
     # Extra zero on MEMDIR ensure we have a number even if the string is empty
-    if (( 0${MEMDIR:3} > 0 )) && [[ "${OCN_ENS_PERTURB_FILES:-false}" == "true" ]]; then
+    if (( 0${MEMDIR:3} > 0 )) && [[ "${STAGE_OCN_PERTURB_FILES:-false}" == "true" ]]; then
         src="${BASE_CPLIC}/${CPL_OCNIC:-}/${PDY}${cyc}/${MEMDIR}/ocean/${PDY}.${cyc}0000.mom6_increment.nc"
         tgt="${COM_OCEAN_RESTART_PREV}/${PDY}.${cyc}0000.mom6_increment.nc"
         ${NCP} "${src}" "${tgt}"

diff --git a/ush/forecast_postdet.sh b/ush/forecast_postdet.sh
@@ -693,6 +693,11 @@ MOM6_postdet() {
       ${NLN} "${COM_OCEAN_ANALYSIS}/${RUN}.t${cyc}z.ocninc.nc" "${DATA}/INPUT/mom6_increment.nc"
   fi
 
+  # GEFS perturbations
+  if (( 0${MEMDIR:3} > 0 )) && [[ "${DO_INCUPD:-False}" == "True" ]]; then
+     ${NLN} "${COM_OCEAN_RESTART_PREV}/${sPDY}.${scyc}0000.mom6_increment.nc" "${DATA}/INPUT/mom6_increment.nc"
+  fi
+
   # Copy MOM6 fixed files
   ${NCP} "${FIXgfs}/mom6/${OCNRES}/"* "${DATA}/INPUT/"  # TODO: These need to be explicit
 
@@ -707,13 +712,17 @@ MOM6_postdet() {
 
   # If using stochatic parameterizations, create a seed that does not exceed the
   # largest signed integer
-  if [[ "${DO_OCN_SPPT}" = "YES" ]] || [[ "${DO_OCN_PERT_EPBL}" = "YES" ]]; then
-    if [[ ${SET_STP_SEED:-"YES"} = "YES" ]]; then
-      ISEED_OCNSPPT=$(( (current_cycle*1000 + MEMBER*10 + 6) % 2147483647 ))
-      ISEED_EPBL=$(( (current_cycle*1000 + MEMBER*10 + 7) % 2147483647 ))
-    else
-      ISEED=${ISEED:-0}
-    fi
+  if [[ ${DO_OCN_SPPT} = "YES" ]]; then
+    OCNSPPT=${OCNSPPT:-0.8,0.4,0.2,0.08,0.04}
+    OCNSPPT_TAU=${OCNSPPT_TAU:-2.16E4,2.592E5,2.592E6,7.776E6,3.1536E7}
+    OCNSPPT_LSCALE=${OCNSPPT_LSCALE:-500.E3,1000.E3,2000.E3,2000.E3,2000.E3}
+    ISEED_OCNSPPT=$((CDATE*10000 + ${MEMBER#0}*100 + 8)),$((CDATE*10000 + ${MEMBER#0}*100 + 9)),$((CDATE*10000 + ${MEMBER#0}*100 + 10)),$((CDATE*10000 + ${MEMBER#0}*100 + 11)),$((CDATE*10000 + ${MEMBER#0}*100 + 12))
+  fi
+  if [[ ${DO_OCN_PERT_EPBL} = "YES" ]]; then
+    EPBL=${EPBL:-0.8,0.4,0.2,0.08,0.04}
+    EPBL_TAU=${EPBL_TAU:-2.16E4,2.592E5,2.592E6,7.776E6,3.1536E7}
+    EPBL_LSCALE=${EPBL_LSCALE:-500.E3,1000.E3,2000.E3,2000.E3,2000.E3}
+    ISEED_EPBL=$((CDATE*10000 + ${MEMBER#0}*100 + 13)),$((CDATE*10000 + ${MEMBER#0}*100 + 14)),$((CDATE*10000 + ${MEMBER#0}*100 + 15)),$((CDATE*10000 + ${MEMBER#0}*100 + 16)),$((CDATE*10000 + ${MEMBER#0}*100 + 17))
   fi
 
   # Link output files

diff --git a/ush/parsing_namelists_FV3.sh b/ush/parsing_namelists_FV3.sh
@@ -631,6 +631,24 @@ EOF
   use_zmtnblck = ${use_zmtnblck:-".true."}
 EOF
   fi
+
+  if [[ "${DO_OCN_SPPT:-NO}" == "YES" ]]; then
+    cat >> input.nml <<EOF
+  OCNSPPT=${OCNSPPT:-1.0}
+  OCNSPPT_LSCALE=${OCNSPPT_LSCALE:-500e3}
+  OCNSPPT_TAU=${OCNSPPT_TAU:-21600}
+  ISEED_OCNSPPT=${ISEED_OCNSPPT:-${ISEED}}
+EOF
+  fi
+
+  if [[ "${DO_OCN_PERT_EPBL:-NO}" == "YES" ]]; then
+    cat >> input.nml <<EOF
+  EPBL=${EPBL:-1.0}
+  EPBL_LSCALE=${EPBL_LSCALE:-500e3}
+  EPBL_TAU=${EPBL_TAU:-21600}
+  ISEED_EPBL=${ISEED_EPBL:-${ISEED}}
+EOF
+  fi
 
   cat >> input.nml << EOF
   ${nam_stochy_nml:-}

diff --git a/ush/parsing_namelists_MOM6.sh b/ush/parsing_namelists_MOM6.sh
@@ -25,24 +25,6 @@ EOF
 #  new_lscale=.true.
 #EOF
 
-if [[ "${DO_OCN_SPPT}" == "YES" ]]; then
-  cat >> input.nml <<EOF
-  OCNSPPT=${OCNSPPT:-1.0}
-  OCNSPPT_LSCALE=${OCNSPPT_LSCALE:-500e3}
-  OCNSPPT_TAU=${OCNSPPT_TAU:-21600}
-  ISEED_OCNSPPT=${ISEED_OCNSPPT:-${ISEED}}
-EOF
-fi
-
-if [[ "${DO_OCN_PERT_EPBL}" == "YES" ]]; then
-  cat >> input.nml <<EOF
-  EPBL=${EPBL:-1.0}
-  EPBL_LSCALE=${EPBL_LSCALE:-500e3}
-  EPBL_TAU=${EPBL_TAU:-21600}
-  ISEED_EPBL=${ISEED_EPBL:-${ISEED}}
-EOF
-fi
-
 #cat >> input.nml <<EOF
 #/
 #
@@ -85,12 +67,12 @@ else
   local MOM6_USE_WAVES="False"
 fi
 # == MOM_oda_incupd section ==
-local ODA_TEMPINC_VAR="Temp"
-local ODA_SALTINC_VAR="Salt"
-local ODA_THK_VAR="h"
+local ODA_TEMPINC_VAR=${ODA_TEMPINC_VAR:-"Temp"}
+local ODA_SALTINC_VAR=${ODA_SALTINC_VAR:-"Salt"}
+local ODA_THK_VAR=${ODA_THK_VAR:-"h"}
 local ODA_INCUPD_UV="True"
-local ODA_UINC_VAR="u"
-local ODA_VINC_VAR="v"
+local ODA_UINC_VAR=${ODA_UINC_VAR:-"u"}
+local ODA_VINC_VAR=${ODA_VINC_VAR:-"v"}
 # ODA_INCUPD
 # ODA_INCUPD_NHOURS
 # == MOM_surface_forcing section ==