diff --git a/scripts/exgdas_efsoi.sh b/scripts/exgdas_efsoi.sh
new file mode 100755
index 0000000000..f2763ff311
--- /dev/null
+++ b/scripts/exgdas_efsoi.sh
@@ -0,0 +1,388 @@
+#!/bin/ksh
+################################################################################
+####  UNIX Script Documentation Block
+#                      .                                             .
+# Script name:         exgdas_efsoi.sh
+# Script description:  Runs efsoi to make observation sensitivities
+#
+# Author:        Liaofan Lin/Andrew Eichmann          Date: 2020-12-03
+#
+# Abstract: This script runs the efsoi executable
+#
+# $Id$
+#
+# Attributes:
+#   Language: POSIX shell
+#   Machine: Hera
+#
+################################################################################
+
+# Set environment.
+VERBOSE=${VERBOSE:-"YES"}
+if [ $VERBOSE = "YES" ] ; then
+   echo $(date) EXECUTING $0 $* >&2
+   set -x
+fi
+
+# Directories.
+pwd=$(pwd)
+
+# Utilities
+NCP=${NCP:-"/bin/cp -p"}
+NLN=${NLN:-"/bin/ln -sf"}
+ERRSCRIPT=${ERRSCRIPT:-'eval [[ $err = 0 ]]'}
+NEMSIOGET=${NEMSIOGET:-$NWPROD/utils/exec/nemsio_get}
+NCLEN=${NCLEN:-$HOMEgfs/ush/getncdimlen}
+USE_CFP=${USE_CFP:-"NO"}
+CFP_MP=${CFP_MP:-"NO"}
+nm=""
+if [ $CFP_MP = "YES" ]; then
+    nm=0
+fi
+APRUNCFP=${APRUNCFP:-""}
+APRUN_ENKF=${APRUN_ENKF:-${APRUN:-""}}
+NTHREADS_ENKF=${NTHREADS_ENKF:-${NTHREADS:-1}}
+
+# Executables
+EFSOIEXEC=${EFSOIEXEC:-$HOMEgfs/exec/efsoi.x}
+
+# Cycling and forecast hour specific parameters
+CDATE=${CDATE:-"2001010100"}
+
+# Filenames.
+GPREFIX=${GPREFIX:-""}
+GSUFFIX=${GSUFFIX:-$SUFFIX}
+APREFIX=${APREFIX:-""}
+ASUFFIX=${ASUFFIX:-$SUFFIX}
+
+SMOOTH_ENKF=${SMOOTH_ENKF:-"YES"}
+
+GBIASe=${GBIASe:-${APREFIX}abias_int.ensmean}
+CNVSTAT=${CNVSTAT:-${APREFIX}cnvstat} # AFE not needed?
+OZNSTAT=${OZNSTAT:-${APREFIX}oznstat} # AFE not needed?
+RADSTAT=${RADSTAT:-${APREFIX}radstat} # AFE not needed?
+#ENKFSTAT=${ENKFSTAT:-${APREFIX}enkfstat}
+EFSOISTAT=${EFSOISTAT:-${APREFIX}efsoistat}
+
+#AFE for EFSOI
+#VERFANL=${VERFANL:-${APREFIX}atmanl.nemsio}
+VERFANL=${VERFANL:-${APREFIX}atmanl.ensres.nc}
+FCSTLONG=${GPREFIX}atmf030.ensmean.nc
+FCSTSHORT=${APREFIX}atmf024.ensmean.nc
+FCST6HRS=${PPREFIX}atmf006.ensmean.nc
+OSENSEIN=osense_${CDATE}_init.dat
+OSENSEOUT=osense_${CDATE}_final.dat
+
+# this needs to be set manually because params in enkf will default to fhr03
+fgfileprefixes=sfg_${CDATE}_fhr06_
+
+
+
+
+#analysise Namelst parameters
+USE_CORRELATED_OBERRS=${USE_CORRELATED_OBERRS:-"NO"}
+NMEM_ENKF=${NMEM_ENKF:-80}
+NAM_ENKF=${NAM_ENKF:-""}
+SATOBS_ENKF=${SATOBS_ENKF:-""}
+OZOBS_ENKF=${OZOBS_ENKF:-""}
+use_correlated_oberrs=${use_correlated_oberrs:-".false."}
+if [ $USE_CORRELATED_OBERRS == "YES" ]; then
+   use_correlated_oberrs=".true."
+fi
+imp_physics=${imp_physics:-"99"}
+#lupp=${lupp:-".true."}
+lupp=${lupp:-".false."} # AFE to match old EFSOI
+corrlength=${corrlength:-1250}
+lnsigcutoff=${lnsigcutoff:-2.5}
+analpertwt=${analpertwt:-0.85}
+#readin_localization_enkf=${readin_localization_enkf:-".true."}
+readin_localization_enkf=${readin_localization_enkf:-".false."} # AFE
+reducedgrid=${reducedgrid:-".true."}
+letkf_flag=${letkf_flag:-".false."}
+getkf=${getkf:-".false."}
+denkf=${denkf:-".false."}
+nobsl_max=${nobsl_max:-10000}
+lobsdiag_forenkf=${lobsdiag_forenkf:-".false."}
+write_spread_diag=${write_spread_diag:-".false."}
+cnvw_option=${cnvw_option:-".false."}
+netcdf_diag=${netcdf_diag:-".true."}
+modelspace_vloc=${modelspace_vloc:-".false."} # if true, 'vlocal_eig.dat' is needed
+IAUFHRS_ENKF=${IAUFHRS_ENKF:-6}
+DO_CALC_INCREMENT=${DO_CALC_INCREMENT:-"NO"}
+INCREMENTS_TO_ZERO=${INCREMENTS_TO_ZERO:-"'NONE'"}
+
+################################################################################
+#ATMGES_ENSMEAN=$COMIN_GES_ENS/${GPREFIX}atmf006.ensmean${GSUFFIX}
+ATMGES_ENSMEAN=$COMIN_ANL/$VERFANL
+if [ $SUFFIX = ".nc" ]; then
+   LONB_ENKF=${LONB_ENKF:-$($NCLEN $ATMGES_ENSMEAN grid_xt)} # get LONB_ENKF
+   LATB_ENKF=${LATB_ENKF:-$($NCLEN $ATMGES_ENSMEAN grid_yt)} # get LATB_ENFK
+   LEVS_ENKF=${LEVS_ENKF:-$($NCLEN $ATMGES_ENSMEAN pfull)} # get LEVS_ENFK
+   use_gfs_ncio=".true."
+   use_gfs_nemsio=".false."
+   paranc=${paranc:-".true."}
+   if [ $DO_CALC_INCREMENT = "YES" ]; then
+      write_fv3_incr=".false."
+   else
+      write_fv3_incr=".true."
+      WRITE_INCR_ZERO="incvars_to_zero= $INCREMENTS_TO_ZERO,"
+   fi
+else
+   LEVS_ENKF=${LEVS_ENKF:-$($NEMSIOGET $ATMGES_ENSMEAN dimz | awk '{print $2}')}
+   LATB_ENKF=${LATB_ENKF:-$($NEMSIOGET $ATMGES_ENSMEAN dimy | awk '{print $2}')}
+   LONB_ENKF=${LONB_ENKF:-$($NEMSIOGET $ATMGES_ENSMEAN dimx | awk '{print $2}')}
+   use_gfs_ncio=".false."
+   use_gfs_nemsio=".true."
+   paranc=${paranc:-".false."}
+fi
+LATA_ENKF=${LATA_ENKF:-$LATB_ENKF}
+LONA_ENKF=${LONA_ENKF:-$LONB_ENKF}
+
+SATANGL=${SATANGL:-${FIXgsi}/global_satangbias.txt}
+SATINFO=${SATINFO:-${FIXgsi}/global_satinfo.txt}
+CONVINFO=${CONVINFO:-${FIXgsi}/global_convinfo.txt}
+OZINFO=${OZINFO:-${FIXgsi}/global_ozinfo.txt}
+SCANINFO=${SCANINFO:-${FIXgsi}/global_scaninfo.txt}
+HYBENSINFO=${HYBENSINFO:-${FIXgsi}/global_hybens_info.l${LEVS_ENKF}.txt}
+ANAVINFO=${ANAVINFO:-${FIXgsi}/global_anavinfo.l${LEVS_ENKF}.txt}
+VLOCALEIG=${VLOCALEIG:-${FIXgsi}/vlocal_eig_l${LEVS_ENKF}.dat}
+
+ENKF_SUFFIX="s"
+[[ $SMOOTH_ENKF = "NO" ]] && ENKF_SUFFIX=""
+
+################################################################################
+# Preprocessing
+mkdata=NO
+if [ ! -d $DATA ]; then
+   mkdata=YES
+   mkdir -p $DATA
+fi
+cd $DATA || exit 99
+
+################################################################################
+# Fixed files
+$NLN $SATANGL    satbias_angle
+$NLN $SATINFO    satinfo
+$NLN $SCANINFO   scaninfo
+$NLN $CONVINFO   convinfo
+$NLN $OZINFO     ozinfo
+$NLN $HYBENSINFO hybens_info
+$NLN $ANAVINFO   anavinfo
+$NLN $VLOCALEIG  vlocal_eig.dat
+
+# Bias correction coefficients based on the ensemble mean
+$NLN $COMOUT_ANL_ENSFSOI/$GBIASe satbias_in
+
+################################################################################
+# Ensemble guess, observational data and analyses/increments
+
+nfhrs=`echo $IAUFHRS_ENKF | sed 's/,/ /g'`
+for imem in $(seq 1 $NMEM_ENKF); do
+   memchar="mem"$(printf %03i $imem)
+   mkdir ${memchar}
+   $NLN $COMOUT_ANL_ENSFSOI/$memchar/${APREFIX}atmf024.nc ${memchar}
+done
+
+$NLN $COMIN_GES_ENS/${GPREFIX}atmf006.ensmean${GSUFFIX} sfg_${CDATE}_fhr06_ensmean
+$NLN $COMIN_GES_ENS/${GPREFIX}atmf006.ensmean${GSUFFIX} sfg_${CDATE}_fhr03_ensmean
+
+# verifying analysis
+$NLN ${COMIN_ANL}/${VERFANL} .
+# the above way is the proper one, but this bit of subterfuge
+# is necessary because ensmean is hard coded into GSI - needs
+# to changed
+$NLN ${COMIN_ANL}/${VERFANL} ${APREFIX}atmanl.ensmean.nc
+
+# forecasts
+$NLN $COMIN_GES_ENS/$FCSTLONG .
+$NLN $COMOUT_ANL_ENSFSOI/$FCSTSHORT .
+
+# inital osense file
+# efsoi.x will read then clobber this
+$NCP $COMOUT_ANL_ENSFSOI/$OSENSEIN osense_${CDATE}.dat
+
+
+
+
+if [ $USE_CFP = "YES" ]; then
+   chmod 755 $DATA/mp_untar.sh
+   ncmd=$(cat $DATA/mp_untar.sh | wc -l)
+   if [ $ncmd -gt 0 ]; then
+      ncmd_max=$((ncmd < npe_node_max ? ncmd : npe_node_max))
+      APRUNCFP=$(eval echo $APRUNCFP)
+      $APRUNCFP $DATA/mp_untar.sh
+      export ERR=$?
+      export err=$ERR
+      $ERRSCRIPT || exit 3
+   fi
+fi
+
+################################################################################
+# Create global_enkf namelist
+# AFE changed from original:
+# gdatehr, datehr, andataname added
+# analpertwt and lnsigcutoff changed upstream
+# numiter from 0 to 1
+#   fso_cycling=.true.,
+#   efsoi_flag=.true.,
+#   wmoist=1.0,adrate=0.75
+
+
+
+cat > enkf.nml << EOFnml
+&nam_enkf
+   datestring="$CDATE",datapath="$DATA/",
+   gdatehr=$gcyc,
+   datehr=$cyc,
+   fgfileprefixes=$fgfileprefixes
+   andataname="$VERFANL",
+   analpertwtnh=${analpertwt},analpertwtsh=${analpertwt},analpertwttr=${analpertwt},
+   covinflatemax=1.e2,covinflatemin=1,pseudo_rh=.true.,iassim_order=0,
+   corrlengthnh=${corrlength},corrlengthsh=${corrlength},corrlengthtr=${corrlength},
+   lnsigcutoffnh=${lnsigcutoff},lnsigcutoffsh=${lnsigcutoff},lnsigcutofftr=${lnsigcutoff},
+   lnsigcutoffpsnh=${lnsigcutoff},lnsigcutoffpssh=${lnsigcutoff},lnsigcutoffpstr=${lnsigcutoff},
+   lnsigcutoffsatnh=${lnsigcutoff},lnsigcutoffsatsh=${lnsigcutoff},lnsigcutoffsattr=${lnsigcutoff},
+   obtimelnh=1.e30,obtimelsh=1.e30,obtimeltr=1.e30,
+   saterrfact=1.0,numiter=1,
+   sprd_tol=1.e30,paoverpb_thresh=0.98,
+   nlons=$LONA_ENKF,nlats=$LATA_ENKF,nlevs=$LEVS_ENKF,nanals=$NMEM_ENKF,
+   deterministic=.true.,sortinc=.true.,lupd_satbiasc=.false.,
+   reducedgrid=${reducedgrid},readin_localization=${readin_localization_enkf}.,
+   use_gfs_nemsio=${use_gfs_nemsio},use_gfs_ncio=${use_gfs_ncio},imp_physics=$imp_physics,lupp=$lupp,
+   univaroz=.false.,adp_anglebc=.true.,angord=4,use_edges=.false.,emiss_bc=.true.,
+   letkf_flag=${letkf_flag},nobsl_max=${nobsl_max},denkf=${denkf},getkf=${getkf}.,
+   nhr_anal=${IAUFHRS_ENKF},nhr_state=${IAUFHRS_ENKF},use_qsatensmean=.true.,
+   lobsdiag_forenkf=$lobsdiag_forenkf,
+   write_spread_diag=$write_spread_diag,
+   modelspace_vloc=$modelspace_vloc,
+   use_correlated_oberrs=${use_correlated_oberrs},
+   netcdf_diag=$netcdf_diag,cnvw_option=$cnvw_option,
+   paranc=$paranc,write_fv3_incr=$write_fv3_incr,
+   efsoi_cycling=.true.,
+   efsoi_flag=.true.,
+   wmoist=1.0,adrate=0.75
+   $WRITE_INCR_ZERO
+   $NAM_ENKF  
+/
+&satobs_enkf
+   sattypes_rad(1) = 'amsua_n15',     dsis(1) = 'amsua_n15',
+   sattypes_rad(2) = 'amsua_n18',     dsis(2) = 'amsua_n18',
+   sattypes_rad(3) = 'amsua_n19',     dsis(3) = 'amsua_n19',
+   sattypes_rad(4) = 'amsub_n16',     dsis(4) = 'amsub_n16',
+   sattypes_rad(5) = 'amsub_n17',     dsis(5) = 'amsub_n17',
+   sattypes_rad(6) = 'amsua_aqua',    dsis(6) = 'amsua_aqua',
+   sattypes_rad(7) = 'amsua_metop-a', dsis(7) = 'amsua_metop-a',
+   sattypes_rad(8) = 'airs_aqua',     dsis(8) = 'airs_aqua',
+   sattypes_rad(9) = 'hirs3_n17',     dsis(9) = 'hirs3_n17',
+   sattypes_rad(10)= 'hirs4_n19',     dsis(10)= 'hirs4_n19',
+   sattypes_rad(11)= 'hirs4_metop-a', dsis(11)= 'hirs4_metop-a',
+   sattypes_rad(12)= 'mhs_n18',       dsis(12)= 'mhs_n18',
+   sattypes_rad(13)= 'mhs_n19',       dsis(13)= 'mhs_n19',
+   sattypes_rad(14)= 'mhs_metop-a',   dsis(14)= 'mhs_metop-a',
+   sattypes_rad(15)= 'goes_img_g11',  dsis(15)= 'imgr_g11',
+   sattypes_rad(16)= 'goes_img_g12',  dsis(16)= 'imgr_g12',
+   sattypes_rad(17)= 'goes_img_g13',  dsis(17)= 'imgr_g13',
+   sattypes_rad(18)= 'goes_img_g14',  dsis(18)= 'imgr_g14',
+   sattypes_rad(19)= 'goes_img_g15',  dsis(19)= 'imgr_g15',
+   sattypes_rad(20)= 'avhrr_n18',     dsis(20)= 'avhrr3_n18',
+   sattypes_rad(21)= 'avhrr_metop-a', dsis(21)= 'avhrr3_metop-a',
+   sattypes_rad(22)= 'avhrr_n19',     dsis(22)= 'avhrr3_n19',
+   sattypes_rad(23)= 'amsre_aqua',    dsis(23)= 'amsre_aqua',
+   sattypes_rad(24)= 'ssmis_f16',     dsis(24)= 'ssmis_f16',
+   sattypes_rad(25)= 'ssmis_f17',     dsis(25)= 'ssmis_f17',
+   sattypes_rad(26)= 'ssmis_f18',     dsis(26)= 'ssmis_f18',
+   sattypes_rad(27)= 'ssmis_f19',     dsis(27)= 'ssmis_f19',
+   sattypes_rad(28)= 'ssmis_f20',     dsis(28)= 'ssmis_f20',
+   sattypes_rad(29)= 'sndrd1_g11',    dsis(29)= 'sndrD1_g11',
+   sattypes_rad(30)= 'sndrd2_g11',    dsis(30)= 'sndrD2_g11',
+   sattypes_rad(31)= 'sndrd3_g11',    dsis(31)= 'sndrD3_g11',
+   sattypes_rad(32)= 'sndrd4_g11',    dsis(32)= 'sndrD4_g11',
+   sattypes_rad(33)= 'sndrd1_g12',    dsis(33)= 'sndrD1_g12',
+   sattypes_rad(34)= 'sndrd2_g12',    dsis(34)= 'sndrD2_g12',
+   sattypes_rad(35)= 'sndrd3_g12',    dsis(35)= 'sndrD3_g12',
+   sattypes_rad(36)= 'sndrd4_g12',    dsis(36)= 'sndrD4_g12',
+   sattypes_rad(37)= 'sndrd1_g13',    dsis(37)= 'sndrD1_g13',
+   sattypes_rad(38)= 'sndrd2_g13',    dsis(38)= 'sndrD2_g13',
+   sattypes_rad(39)= 'sndrd3_g13',    dsis(39)= 'sndrD3_g13',
+   sattypes_rad(40)= 'sndrd4_g13',    dsis(40)= 'sndrD4_g13',
+   sattypes_rad(41)= 'sndrd1_g14',    dsis(41)= 'sndrD1_g14',
+   sattypes_rad(42)= 'sndrd2_g14',    dsis(42)= 'sndrD2_g14',
+   sattypes_rad(43)= 'sndrd3_g14',    dsis(43)= 'sndrD3_g14',
+   sattypes_rad(44)= 'sndrd4_g14',    dsis(44)= 'sndrD4_g14',
+   sattypes_rad(45)= 'sndrd1_g15',    dsis(45)= 'sndrD1_g15',
+   sattypes_rad(46)= 'sndrd2_g15',    dsis(46)= 'sndrD2_g15',
+   sattypes_rad(47)= 'sndrd3_g15',    dsis(47)= 'sndrD3_g15',
+   sattypes_rad(48)= 'sndrd4_g15',    dsis(48)= 'sndrD4_g15',
+   sattypes_rad(49)= 'iasi_metop-a',  dsis(49)= 'iasi_metop-a',
+   sattypes_rad(50)= 'seviri_m08',    dsis(50)= 'seviri_m08',
+   sattypes_rad(51)= 'seviri_m09',    dsis(51)= 'seviri_m09',
+   sattypes_rad(52)= 'seviri_m10',    dsis(52)= 'seviri_m10',
+   sattypes_rad(53)= 'seviri_m11',    dsis(53)= 'seviri_m11',
+   sattypes_rad(54)= 'amsua_metop-b', dsis(54)= 'amsua_metop-b',
+   sattypes_rad(55)= 'hirs4_metop-b', dsis(55)= 'hirs4_metop-b',
+   sattypes_rad(56)= 'mhs_metop-b',   dsis(56)= 'mhs_metop-b',
+   sattypes_rad(57)= 'iasi_metop-b',  dsis(57)= 'iasi_metop-b',
+   sattypes_rad(58)= 'avhrr_metop-b', dsis(58)= 'avhrr3_metop-b',
+   sattypes_rad(59)= 'atms_npp',      dsis(59)= 'atms_npp',
+   sattypes_rad(60)= 'atms_n20',      dsis(60)= 'atms_n20',
+   sattypes_rad(61)= 'cris_npp',      dsis(61)= 'cris_npp',
+   sattypes_rad(62)= 'cris-fsr_npp',  dsis(62)= 'cris-fsr_npp',
+   sattypes_rad(63)= 'cris-fsr_n20',  dsis(63)= 'cris-fsr_n20',
+   sattypes_rad(64)= 'gmi_gpm',       dsis(64)= 'gmi_gpm',
+   sattypes_rad(65)= 'saphir_meghat', dsis(65)= 'saphir_meghat',
+   sattypes_rad(66)= 'amsua_metop-c', dsis(66)= 'amsua_metop-c',
+   sattypes_rad(67)= 'mhs_metop-c',   dsis(67)= 'mhs_metop-c',
+   sattypes_rad(68)= 'ahi_himawari8', dsis(68)= 'ahi_himawari8',
+   sattypes_rad(69)= 'abi_g16',       dsis(69)= 'abi_g16',
+   sattypes_rad(70)= 'abi_g17',       dsis(70)= 'abi_g17',
+   $SATOBS_ENKF
+/
+&ozobs_enkf
+   sattypes_oz(1) = 'sbuv2_n16',
+   sattypes_oz(2) = 'sbuv2_n17',
+   sattypes_oz(3) = 'sbuv2_n18',
+   sattypes_oz(4) = 'sbuv2_n19',
+   sattypes_oz(5) = 'omi_aura',
+   sattypes_oz(6) = 'gome_metop-a',
+   sattypes_oz(7) = 'gome_metop-b',
+   sattypes_oz(8) = 'mls30_aura',
+   sattypes_oz(9) = 'ompsnp_npp',
+   sattypes_oz(10) = 'ompstc8_npp',
+   $OZOBS_ENKF
+/
+EOFnml
+
+################################################################################
+# Run enkf update
+
+export OMP_NUM_THREADS=$NTHREADS_ENKF
+#export pgm=$ENKFEXEC AFE
+export pgm=$EFSOIEXEC
+. prep_step
+
+#$NCP $ENKFEXEC $DATA AFE
+#$APRUN_ENKF ${DATA}/$(basename $ENKFEXEC) 1>stdout 2>stderr AFE
+$NCP $EFSOIEXEC $DATA
+$APRUN_ENKF ${DATA}/$(basename $EFSOIEXEC) 1>stdout 2>stderr
+rc=$?
+
+export ERR=$rc
+export err=$ERR
+$ERRSCRIPT || exit 2
+
+# Cat runtime output files.
+cat stdout stderr > $COMOUT_ANL_ENSFSOI/$EFSOISTAT
+
+$NCP osense_${CDATE}.dat $COMOUT_ANL_ENSFSOI/$OSENSEOUT
+$NCP osense_${CDATE}.dat $OSENSE_SAVE_DIR/$OSENSEOUT
+
+################################################################################
+#  Postprocessing
+
+cd $pwd
+[[ $mkdata = "YES" ]] && rm -rf $DATA
+set +x
+if [ $VERBOSE = "YES" ]; then
+   echo $(date) EXITING $0 with return code $err >&2
+fi
+exit $err
diff --git a/scripts/exgdas_efsoi_fcst.sh b/scripts/exgdas_efsoi_fcst.sh
new file mode 100755
index 0000000000..13f75cd625
--- /dev/null
+++ b/scripts/exgdas_efsoi_fcst.sh
@@ -0,0 +1,243 @@
+#!/bin/ksh
+################################################################################
+####  UNIX Script Documentation Block
+#                      .                                             .
+# Script name:         exgdas_efsoi_fcst.sh
+# Script description:  Run ensemble forecasts
+#
+# Author:        Rahul Mahajan      Org: NCEP/EMC     Date: 2017-03-02
+# Author:        Andrew Eichmann    Org: NCEP/EMC     Date: 2021-02-25
+#
+# Abstract: This script runs ensemble forecasts serially one-after-another
+#           for EFSOI.  Based on exgdas_enkf_fcst.sh
+#
+# $Id$
+#
+# Attributes:
+#   Language: POSIX shell
+#   Machine: WCOSS-Cray/Theia
+#
+####
+################################################################################
+
+# Set environment.
+export VERBOSE=${VERBOSE:-"YES"}
+if [ $VERBOSE = "YES" ] ; then
+   echo $(date) EXECUTING $0 $* >&2
+   set -x
+fi
+
+# Directories.
+pwd=$(pwd)
+export FIX_DIR=${FIX_DIR:-$HOMEgfs/fix}
+export FIX_AM=${FIX_AM:-$FIX_DIR/fix_am}
+
+# Utilities
+export NCP=${NCP:-"/bin/cp -p"}
+export NMV=${NMV:-"/bin/mv"}
+export NLN=${NLN:-"/bin/ln -sf"}
+export ERRSCRIPT=${ERRSCRIPT:-'eval [[ $err = 0 ]]'}
+export NDATE=${NDATE:-/$NWPROD/util/exec/ndate}
+
+# Scripts.
+FORECASTSH=${FORECASTSH:-$HOMEgfs/scripts/exglobal_forecast.sh}
+
+# Enemble group, begin and end
+ENSGRP=${ENSGRP:-1}
+ENSBEG=${ENSBEG:-1}
+ENSEND=${ENSEND:-1}
+
+# Model builds
+export FCSTEXECDIR=${FCSTEXECDIR:-$HOMEgfs/sorc/fv3gfs.fd/BUILD/bin}
+export FCSTEXEC=${FCSTEXEC:-fv3gfs.x}
+
+# Get DA specific diag table.
+export PARM_FV3DIAG=${PARM_FV3DIAG:-$HOMEgfs/parm/parm_fv3diag}
+export DIAG_TABLE=${DIAG_TABLE_ENKF:-${DIAG_TABLE:-$PARM_FV3DIAG/diag_table_da}}
+
+# Cycling and forecast hour specific parameters
+export CDATE=${CDATE:-"2001010100"}
+export CDUMP=${CDUMP:-"gdas"}
+
+# Re-run failed members, or entire group
+RERUN_EFCSGRP=${RERUN_EFCSGRP:-"YES"}
+
+# Recenter flag and increment file prefix
+RECENTER_ENKF=${RECENTER_ENKF:-"YES"}
+export PREFIX_ATMINC=${PREFIX_ATMINC:-""}
+
+# Ops related stuff
+SENDECF=${SENDECF:-"NO"}
+SENDDBN=${SENDDBN:-"NO"}
+GSUFFIX=${GSUFFIX:-$SUFFIX}
+
+################################################################################
+# Preprocessing
+mkdata=NO
+if [ ! -d $DATA ]; then
+   mkdata=YES
+   mkdir -p $DATA
+fi
+cd $DATA || exit 99
+DATATOP=$DATA
+
+################################################################################
+# Set output data
+cymd=$(echo $CDATE | cut -c1-8)
+chh=$(echo  $CDATE | cut -c9-10)
+EFCSGRP=$COMOUT/efcs.grp${ENSGRP}
+if [ -f $EFCSGRP ]; then
+   if [ $RERUN_EFCSGRP = "YES" ]; then
+      rm -f $EFCSGRP
+   else
+      echo "RERUN_EFCSGRP = $RERUN_EFCSGRP, will re-run FAILED members only!"
+      $NMV $EFCSGRP ${EFCSGRP}.fail
+   fi
+fi
+
+################################################################################
+# Set namelist/model config options common to all members once
+
+# There are many many model namelist options
+# Some are resolution (CASE) dependent, some depend on the model configuration
+# and will need to be added here before $FORECASTSH is called
+# For now assume that
+# 1. the ensemble and the deterministic are same resolution
+# 2. the ensemble runs with the same configuration as the deterministic
+
+# Model config option for Ensemble
+export TYPE=${TYPE_ENKF:-${TYPE:-nh}}                  # choices:  nh, hydro
+export MONO=${MONO_ENKF:-${MONO:-non-mono}}            # choices:  mono, non-mono
+
+# fv_core_nml
+export CASE=${CASE_ENKF:-${CASE:-C768}}
+export layout_x=${layout_x_ENKF:-${layout_x:-8}}
+export layout_y=${layout_y_ENKF:-${layout_y:-16}}
+export LEVS=${LEVS_ENKF:-${LEVS:-64}}
+
+# nggps_diag_nml
+export FHOUT=${FHOUT_ENKF:-3}
+
+# model_configure
+export DELTIM=${DELTIM_ENKF:-${DELTIM:-225}}
+export FHMAX=${FHMAX_ENKF:-9}
+export restart_interval=${restart_interval_ENKF:-${restart_interval:-6}}
+
+# gfs_physics_nml
+export FHSWR=${FHSWR_ENKF:-${FHSWR:-3600.}}
+export FHLWR=${FHLWR_ENKF:-${FHLWR:-3600.}}
+export IEMS=${IEMS_ENKF:-${IEMS:-1}}
+export ISOL=${ISOL_ENKF:-${ISOL:-2}}
+export IAER=${IAER_ENKF:-${IAER:-111}}
+export ICO2=${ICO2_ENKF:-${ICO2:-2}}
+export cdmbgwd=${cdmbgwd_ENKF:-${cdmbgwd:-"3.5,0.25"}}
+export dspheat=${dspheat_ENKF:-${dspheat:-".true."}}
+export shal_cnv=${shal_cnv_ENKF:-${shal_cnv:-".true."}}
+export FHZER=${FHZER_ENKF:-${FHZER:-6}}
+export FHCYC=${FHCYC_ENKF:-${FHCYC:-6}}
+
+# Set PREFIX_ATMINC to r when recentering on
+if [ $RECENTER_ENKF = "YES" ]; then
+   export PREFIX_ATMINC="r"
+fi
+
+# APRUN for different executables
+export APRUN_FV3=${APRUN_FV3:-${APRUN:-""}}
+export NTHREADS_FV3=${NTHREADS_FV3:-${NTHREADS:-1}}
+
+################################################################################
+# Run forecast for ensemble member
+rc=0
+for imem in $(seq $ENSBEG $ENSEND); do
+
+   cd $DATATOP
+
+   cmem=$(printf %03i $imem)
+   memchar="mem$cmem"
+
+   echo "Processing MEMBER: $cmem"
+
+   ra=0
+
+   skip_mem="NO"
+   if [ -f ${EFCSGRP}.fail ]; then
+      memstat=$(cat ${EFCSGRP}.fail | grep "MEMBER $cmem" | grep "PASS" | wc -l)
+      [[ $memstat -eq 1 ]] && skip_mem="YES"
+   fi
+
+   if [ $skip_mem = "NO" ]; then
+
+      ra=0
+
+      export MEMBER=$imem
+      export DATA=$DATATOP/$memchar
+      if [ -d $DATA ]; then rm -rf $DATA; fi
+      $FORECASTSH
+      ra=$?
+
+      # Notify a member forecast failed and abort
+      if [ $ra -ne 0 ]; then
+         msg="FATAL ERROR:  forecast of member $cmem FAILED.  Aborting job"
+         print $msg
+         export err=$ra
+         $ERRSCRIPT || exit 2
+      fi
+
+      ((rc+=ra))
+
+   fi
+
+   if [ $SENDDBN = YES ]; then
+     fhr=$FHOUT
+     while [ $fhr -le $FHMAX ]; do
+       FH3=$(printf %03i $fhr)
+       if [ $(expr $fhr % 3) -eq 0 ]; then
+         $DBNROOT/bin/dbn_alert MODEL GFS_ENKF $job $COMOUT/$memchar/${CDUMP}.t${cyc}z.sfcf${FH3}${GSUFFIX}
+       fi
+       fhr=$((fhr+FHOUT))
+     done
+   fi
+
+   cd $DATATOP
+
+   if [ -s $EFCSGRP ]; then
+       $NCP $EFCSGRP log_old
+   fi
+   [[ -f log ]] && rm log
+   [[ -f log_new ]] && rm log_new
+   if [ $ra -ne 0 ]; then
+      echo "MEMBER $cmem : FAIL" > log
+   else
+      echo "MEMBER $cmem : PASS" > log
+   fi
+   if [ -s log_old ] ; then
+       cat log_old log > log_new
+   else
+       cat log > log_new
+   fi
+   $NCP log_new $EFCSGRP
+
+done
+
+################################################################################
+# Echo status of ensemble group
+cd $DATATOP
+echo "Status of ensemble members in group $ENSGRP:"
+cat $EFCSGRP
+[[ -f ${EFCSGRP}.fail ]] && rm ${EFCSGRP}.fail
+
+################################################################################
+# If any members failed, error out
+export ERR=$rc
+export err=$ERR
+$ERRSCRIPT || exit 2
+
+################################################################################
+#  Postprocessing
+cd $pwd
+[[ $mkdata = "YES" ]] && rm -rf $DATATOP
+set +x
+if [ $VERBOSE = "YES" ] ; then
+   echo $(date) EXITING $0 with return code $err >&2
+fi
+exit $err
diff --git a/scripts/exgdas_efsoi_update.sh b/scripts/exgdas_efsoi_update.sh
new file mode 100755
index 0000000000..73ece19b26
--- /dev/null
+++ b/scripts/exgdas_efsoi_update.sh
@@ -0,0 +1,411 @@
+#!/bin/ksh
+################################################################################
+####  UNIX Script Documentation Block
+#                      .                                             .
+# Script name:         exgdas_efsoi_update.sh
+# Script description:  Make global_enkf update for efsoi
+#
+# Author:        Rahul Mahajan      Org: NCEP/EMC     Date: 2017-03-02
+# Author:        Liaofan Lin/Andrew Eichmann          Date: 2020-12-03
+#
+# Abstract: This script runs the global_enkf update for efsoi
+#
+# $Id$
+#
+# Attributes:
+#   Language: POSIX shell
+#   Machine: WCOSS-Cray/Theia
+#
+################################################################################
+
+# Set environment.
+VERBOSE=${VERBOSE:-"YES"}
+if [ $VERBOSE = "YES" ] ; then
+   echo $(date) EXECUTING $0 $* >&2
+   set -x
+fi
+
+# Directories.
+pwd=$(pwd)
+
+# Utilities
+NCP=${NCP:-"/bin/cp -p"}
+NLN=${NLN:-"/bin/ln -sf"}
+ERRSCRIPT=${ERRSCRIPT:-'eval [[ $err = 0 ]]'}
+NEMSIOGET=${NEMSIOGET:-$NWPROD/utils/exec/nemsio_get}
+NCLEN=${NCLEN:-$HOMEgfs/ush/getncdimlen}
+USE_CFP=${USE_CFP:-"NO"}
+CFP_MP=${CFP_MP:-"NO"}
+nm=""
+if [ $CFP_MP = "YES" ]; then
+    nm=0
+fi
+APRUNCFP=${APRUNCFP:-""}
+APRUN_ENKF=${APRUN_ENKF:-${APRUN:-""}}
+NTHREADS_ENKF=${NTHREADS_ENKF:-${NTHREADS:-1}}
+
+# Executables
+ENKFEXEC=${ENKFEXEC:-$HOMEgfs/exec/global_enkf.x}
+
+# Cycling and forecast hour specific parameters
+CDATE=${CDATE:-"2001010100"}
+
+# Filenames.
+GPREFIX=${GPREFIX:-""}
+GSUFFIX=${GSUFFIX:-$SUFFIX}
+APREFIX=${APREFIX:-""}
+ASUFFIX=${ASUFFIX:-$SUFFIX}
+
+SMOOTH_ENKF=${SMOOTH_ENKF:-"YES"}
+
+GBIASe=${GBIASe:-${APREFIX}abias_int.ensmean}
+CNVSTAT=${CNVSTAT:-${APREFIX}cnvstat}
+OZNSTAT=${OZNSTAT:-${APREFIX}oznstat}
+RADSTAT=${RADSTAT:-${APREFIX}radstat}
+ENKFSTAT=${ENKFSTAT:-${APREFIX}enkfstat}
+
+# Namelist parameters
+USE_CORRELATED_OBERRS=${USE_CORRELATED_OBERRS:-"NO"}
+NMEM_ENKF=${NMEM_ENKF:-80}
+NAM_ENKF=${NAM_ENKF:-""}
+SATOBS_ENKF=${SATOBS_ENKF:-""}
+OZOBS_ENKF=${OZOBS_ENKF:-""}
+use_correlated_oberrs=${use_correlated_oberrs:-".false."}
+if [ $USE_CORRELATED_OBERRS == "YES" ]; then
+   use_correlated_oberrs=".true."
+fi
+imp_physics=${imp_physics:-"99"}
+lupp=${lupp:-".true."}
+corrlength=${corrlength:-1250}
+lnsigcutoff=${lnsigcutoff:-2.5}
+analpertwt=${analpertwt:-0.85}
+readin_localization_enkf=${readin_localization_enkf:-".true."}
+reducedgrid=${reducedgrid:-".true."}
+letkf_flag=${letkf_flag:-".false."}
+getkf=${getkf:-".false."}
+denkf=${denkf:-".false."}
+nobsl_max=${nobsl_max:-10000}
+lobsdiag_forenkf=${lobsdiag_forenkf:-".false."}
+write_spread_diag=${write_spread_diag:-".false."}
+cnvw_option=${cnvw_option:-".false."}
+netcdf_diag=${netcdf_diag:-".true."}
+modelspace_vloc=${modelspace_vloc:-".false."} # if true, 'vlocal_eig.dat' is needed
+IAUFHRS_ENKF=${IAUFHRS_ENKF:-6}
+DO_CALC_INCREMENT=${DO_CALC_INCREMENT:-"NO"}
+INCREMENTS_TO_ZERO=${INCREMENTS_TO_ZERO:-"'NONE'"}
+
+################################################################################
+ATMGES_ENSMEAN=$COMIN_GES_ENS/${GPREFIX}atmf006.ensmean${GSUFFIX}
+if [ $SUFFIX = ".nc" ]; then
+   LONB_ENKF=${LONB_ENKF:-$($NCLEN $ATMGES_ENSMEAN grid_xt)} # get LONB_ENKF
+   LATB_ENKF=${LATB_ENKF:-$($NCLEN $ATMGES_ENSMEAN grid_yt)} # get LATB_ENFK
+   LEVS_ENKF=${LEVS_ENKF:-$($NCLEN $ATMGES_ENSMEAN pfull)} # get LEVS_ENFK
+   use_gfs_ncio=".true."
+   use_gfs_nemsio=".false."
+   paranc=${paranc:-".true."}
+   if [ $DO_CALC_INCREMENT = "YES" ]; then
+      write_fv3_incr=".false."
+   else
+      write_fv3_incr=".true."
+      WRITE_INCR_ZERO="incvars_to_zero= $INCREMENTS_TO_ZERO,"
+   fi
+else
+   LEVS_ENKF=${LEVS_ENKF:-$($NEMSIOGET $ATMGES_ENSMEAN dimz | awk '{print $2}')}
+   LATB_ENKF=${LATB_ENKF:-$($NEMSIOGET $ATMGES_ENSMEAN dimy | awk '{print $2}')}
+   LONB_ENKF=${LONB_ENKF:-$($NEMSIOGET $ATMGES_ENSMEAN dimx | awk '{print $2}')}
+   use_gfs_ncio=".false."
+   use_gfs_nemsio=".true."
+   paranc=${paranc:-".false."}
+fi
+LATA_ENKF=${LATA_ENKF:-$LATB_ENKF}
+LONA_ENKF=${LONA_ENKF:-$LONB_ENKF}
+
+SATANGL=${SATANGL:-${FIXgsi}/global_satangbias.txt}
+SATINFO=${SATINFO:-${FIXgsi}/global_satinfo.txt}
+CONVINFO=${CONVINFO:-${FIXgsi}/global_convinfo.txt}
+OZINFO=${OZINFO:-${FIXgsi}/global_ozinfo.txt}
+SCANINFO=${SCANINFO:-${FIXgsi}/global_scaninfo.txt}
+HYBENSINFO=${HYBENSINFO:-${FIXgsi}/global_hybens_info.l${LEVS_ENKF}.txt}
+ANAVINFO=${ANAVINFO:-${FIXgsi}/global_anavinfo.l${LEVS_ENKF}.txt}
+VLOCALEIG=${VLOCALEIG:-${FIXgsi}/vlocal_eig_l${LEVS_ENKF}.dat}
+
+ENKF_SUFFIX="s"
+[[ $SMOOTH_ENKF = "NO" ]] && ENKF_SUFFIX=""
+
+################################################################################
+# Preprocessing
+mkdata=NO
+if [ ! -d $DATA ]; then
+   mkdata=YES
+   mkdir -p $DATA
+fi
+cd $DATA || exit 99
+
+################################################################################
+# Fixed files
+$NLN $SATANGL    satbias_angle
+$NLN $SATINFO    satinfo
+$NLN $SCANINFO   scaninfo
+$NLN $CONVINFO   convinfo
+$NLN $OZINFO     ozinfo
+$NLN $HYBENSINFO hybens_info
+$NLN $ANAVINFO   anavinfo
+$NLN $VLOCALEIG  vlocal_eig.dat
+
+# Bias correction coefficients based on the ensemble mean
+$NLN $COMOUT_ANL_ENS/$GBIASe satbias_in
+
+################################################################################
+
+if [ $USE_CFP = "YES" ]; then
+   [[ -f $DATA/untar.sh ]] && rm $DATA/untar.sh
+   [[ -f $DATA/mp_untar.sh ]] && rm $DATA/mp_untar.sh
+   set +x
+   cat > $DATA/untar.sh << EOFuntar
+#!/bin/sh
+memchar=\$1
+flist="$CNVSTAT $OZNSTAT $RADSTAT"
+for ftype in \$flist; do
+   if [ \$memchar = "ensmean" ]; then
+      fname=$COMOUT_ANL_ENS/\${ftype}.ensmean
+   else
+      fname=$COMOUT_ANL_ENS/\$memchar/\$ftype
+   fi
+   tar -xvf \$fname
+done
+EOFuntar
+   set -x
+   chmod 755 $DATA/untar.sh
+fi
+
+################################################################################
+# Ensemble guess, observational data and analyses/increments
+
+flist="$CNVSTAT $OZNSTAT $RADSTAT"
+if [ $USE_CFP = "YES" ]; then
+   echo "$nm $DATA/untar.sh ensmean" | tee -a $DATA/mp_untar.sh
+   if [ ${CFP_MP:-"NO"} = "YES" ]; then
+       nm=$((nm+1))
+   fi
+else
+   for ftype in $flist; do
+      fname=$COMOUT_ANL_ENS/${ftype}.ensmean
+      tar -xvf $fname
+   done
+fi
+nfhrs=`echo $IAUFHRS_ENKF | sed 's/,/ /g'`
+for imem in $(seq 1 $NMEM_ENKF); do
+   memchar="mem"$(printf %03i $imem)
+   if [ $lobsdiag_forenkf = ".false." ]; then
+      if [ $USE_CFP = "YES" ]; then
+         echo "$nm $DATA/untar.sh $memchar" | tee -a $DATA/mp_untar.sh
+         if [ ${CFP_MP:-"NO"} = "YES" ]; then
+             nm=$((nm+1))
+         fi
+      else
+         for ftype in $flist; do
+            fname=$COMOUT_ANL_ENS/$memchar/$ftype
+            tar -xvf $fname
+         done
+      fi
+   fi
+   mkdir -p $COMOUT_ANL_ENSFSOI/$memchar
+   for FHR in $nfhrs; do
+      $NLN $COMIN_GES_ENS/$memchar/${GPREFIX}atmf00${FHR}${ENKF_SUFFIX}${GSUFFIX}  sfg_${CDATE}_fhr0${FHR}_${memchar}
+      if [ $cnvw_option = ".true." ]; then
+         $NLN $COMIN_GES_ENS/$memchar/${GPREFIX}sfcf00${FHR}${GSUFFIX} sfgsfc_${CDATE}_fhr0${FHR}_${memchar}
+      fi
+      if [ $FHR -eq 6 ]; then
+         if [ $DO_CALC_INCREMENT = "YES" ]; then
+            #$NLN $COMOUT_ANL_ENS/$memchar/${APREFIX}atmanl${ASUFFIX}             sanl_${CDATE}_fhr0${FHR}_${memchar}
+            $NLN $COMOUT_ANL_ENSFSOI/$memchar/${APREFIX}atmanl${ASUFFIX}             sanl_${CDATE}_fhr0${FHR}_ni${memchar}
+         else
+            #$NLN $COMOUT_ANL_ENS/$memchar/${APREFIX}atminc${ASUFFIX}             incr_${CDATE}_fhr0${FHR}_${memchar}
+            $NLN $COMOUT_ANL_ENSFSOI/$memchar/${APREFIX}atminc${ASUFFIX}             incr_${CDATE}_fhr0${FHR}_ni${memchar}
+         fi
+      else
+         if [ $DO_CALC_INCREMENT = "YES" ]; then
+            #$NLN $COMOUT_ANL_ENS/$memchar/${APREFIX}atma00${FHR}${ASUFFIX}             sanl_${CDATE}_fhr0${FHR}_${memchar}
+            $NLN $COMOUT_ANL_ENSFSOI/$memchar/${APREFIX}atma00${FHR}${ASUFFIX}             sanl_${CDATE}_fhr0${FHR}_ni${memchar}
+         else
+            #$NLN $COMOUT_ANL_ENS/$memchar/${APREFIX}atmi00${FHR}${ASUFFIX}             incr_${CDATE}_fhr0${FHR}_${memchar}
+            $NLN $COMOUT_ANL_ENSFSOI/$memchar/${APREFIX}atmi00${FHR}${ASUFFIX}             incr_${CDATE}_fhr0${FHR}_ni${memchar}
+         fi
+      fi
+   done
+done
+
+# Ensemble mean guess
+for FHR in $nfhrs; do
+   $NLN $COMIN_GES_ENS/${GPREFIX}atmf00${FHR}.ensmean${GSUFFIX} sfg_${CDATE}_fhr0${FHR}_ensmean
+   if [ $cnvw_option = ".true." ]; then
+      $NLN $COMIN_GES_ENS/${GPREFIX}sfcf00${FHR}.ensmean${GSUFFIX} sfgsfc_${CDATE}_fhr0${FHR}_ensmean
+   fi
+done
+
+if [ $USE_CFP = "YES" ]; then
+   chmod 755 $DATA/mp_untar.sh
+   ncmd=$(cat $DATA/mp_untar.sh | wc -l)
+   if [ $ncmd -gt 0 ]; then
+      ncmd_max=$((ncmd < npe_node_max ? ncmd : npe_node_max))
+      APRUNCFP=$(eval echo $APRUNCFP)
+      $APRUNCFP $DATA/mp_untar.sh
+      export ERR=$?
+      export err=$ERR
+      $ERRSCRIPT || exit 3
+   fi
+fi
+
+# Keep osense data in comrot
+$NLN $COMOUT_ANL_ENSFSOI/osense_${CDATE}_init.dat  osense_${CDATE}.dat
+
+################################################################################
+# Create global_enkf namelist
+cat > enkf.nml << EOFnml
+&nam_enkf
+   datestring="$CDATE",datapath="$DATA/",
+   analpertwtnh=${analpertwt},analpertwtsh=${analpertwt},analpertwttr=${analpertwt},
+   covinflatemax=1.e2,covinflatemin=1,pseudo_rh=.true.,iassim_order=0,
+   corrlengthnh=${corrlength},corrlengthsh=${corrlength},corrlengthtr=${corrlength},
+   lnsigcutoffnh=${lnsigcutoff},lnsigcutoffsh=${lnsigcutoff},lnsigcutofftr=${lnsigcutoff},
+   lnsigcutoffpsnh=${lnsigcutoff},lnsigcutoffpssh=${lnsigcutoff},lnsigcutoffpstr=${lnsigcutoff},
+   lnsigcutoffsatnh=${lnsigcutoff},lnsigcutoffsatsh=${lnsigcutoff},lnsigcutoffsattr=${lnsigcutoff},
+   obtimelnh=1.e30,obtimelsh=1.e30,obtimeltr=1.e30,
+   saterrfact=1.0,numiter=0,
+   sprd_tol=1.e30,paoverpb_thresh=0.98,
+   nlons=$LONA_ENKF,nlats=$LATA_ENKF,nlevs=$LEVS_ENKF,nanals=$NMEM_ENKF,
+   deterministic=.true.,sortinc=.true.,lupd_satbiasc=.false.,
+   reducedgrid=${reducedgrid},readin_localization=${readin_localization_enkf}.,
+   use_gfs_nemsio=${use_gfs_nemsio},use_gfs_ncio=${use_gfs_ncio},imp_physics=$imp_physics,lupp=$lupp,
+   univaroz=.false.,adp_anglebc=.true.,angord=4,use_edges=.false.,emiss_bc=.true.,
+   letkf_flag=${letkf_flag},nobsl_max=${nobsl_max},denkf=${denkf},getkf=${getkf}.,
+   nhr_anal=${IAUFHRS_ENKF},nhr_state=${IAUFHRS_ENKF},use_qsatensmean=.true.,
+   lobsdiag_forenkf=$lobsdiag_forenkf,
+   write_spread_diag=$write_spread_diag,
+   modelspace_vloc=$modelspace_vloc,
+   use_correlated_oberrs=${use_correlated_oberrs},
+   netcdf_diag=$netcdf_diag,cnvw_option=$cnvw_option,
+   paranc=$paranc,write_fv3_incr=$write_fv3_incr,
+   efsoi_cycling=.true.,
+   $WRITE_INCR_ZERO
+   $NAM_ENKF  
+/
+&satobs_enkf
+   sattypes_rad(1) = 'amsua_n15',     dsis(1) = 'amsua_n15',
+   sattypes_rad(2) = 'amsua_n18',     dsis(2) = 'amsua_n18',
+   sattypes_rad(3) = 'amsua_n19',     dsis(3) = 'amsua_n19',
+   sattypes_rad(4) = 'amsub_n16',     dsis(4) = 'amsub_n16',
+   sattypes_rad(5) = 'amsub_n17',     dsis(5) = 'amsub_n17',
+   sattypes_rad(6) = 'amsua_aqua',    dsis(6) = 'amsua_aqua',
+   sattypes_rad(7) = 'amsua_metop-a', dsis(7) = 'amsua_metop-a',
+   sattypes_rad(8) = 'airs_aqua',     dsis(8) = 'airs_aqua',
+   sattypes_rad(9) = 'hirs3_n17',     dsis(9) = 'hirs3_n17',
+   sattypes_rad(10)= 'hirs4_n19',     dsis(10)= 'hirs4_n19',
+   sattypes_rad(11)= 'hirs4_metop-a', dsis(11)= 'hirs4_metop-a',
+   sattypes_rad(12)= 'mhs_n18',       dsis(12)= 'mhs_n18',
+   sattypes_rad(13)= 'mhs_n19',       dsis(13)= 'mhs_n19',
+   sattypes_rad(14)= 'mhs_metop-a',   dsis(14)= 'mhs_metop-a',
+   sattypes_rad(15)= 'goes_img_g11',  dsis(15)= 'imgr_g11',
+   sattypes_rad(16)= 'goes_img_g12',  dsis(16)= 'imgr_g12',
+   sattypes_rad(17)= 'goes_img_g13',  dsis(17)= 'imgr_g13',
+   sattypes_rad(18)= 'goes_img_g14',  dsis(18)= 'imgr_g14',
+   sattypes_rad(19)= 'goes_img_g15',  dsis(19)= 'imgr_g15',
+   sattypes_rad(20)= 'avhrr_n18',     dsis(20)= 'avhrr3_n18',
+   sattypes_rad(21)= 'avhrr_metop-a', dsis(21)= 'avhrr3_metop-a',
+   sattypes_rad(22)= 'avhrr_n19',     dsis(22)= 'avhrr3_n19',
+   sattypes_rad(23)= 'amsre_aqua',    dsis(23)= 'amsre_aqua',
+   sattypes_rad(24)= 'ssmis_f16',     dsis(24)= 'ssmis_f16',
+   sattypes_rad(25)= 'ssmis_f17',     dsis(25)= 'ssmis_f17',
+   sattypes_rad(26)= 'ssmis_f18',     dsis(26)= 'ssmis_f18',
+   sattypes_rad(27)= 'ssmis_f19',     dsis(27)= 'ssmis_f19',
+   sattypes_rad(28)= 'ssmis_f20',     dsis(28)= 'ssmis_f20',
+   sattypes_rad(29)= 'sndrd1_g11',    dsis(29)= 'sndrD1_g11',
+   sattypes_rad(30)= 'sndrd2_g11',    dsis(30)= 'sndrD2_g11',
+   sattypes_rad(31)= 'sndrd3_g11',    dsis(31)= 'sndrD3_g11',
+   sattypes_rad(32)= 'sndrd4_g11',    dsis(32)= 'sndrD4_g11',
+   sattypes_rad(33)= 'sndrd1_g12',    dsis(33)= 'sndrD1_g12',
+   sattypes_rad(34)= 'sndrd2_g12',    dsis(34)= 'sndrD2_g12',
+   sattypes_rad(35)= 'sndrd3_g12',    dsis(35)= 'sndrD3_g12',
+   sattypes_rad(36)= 'sndrd4_g12',    dsis(36)= 'sndrD4_g12',
+   sattypes_rad(37)= 'sndrd1_g13',    dsis(37)= 'sndrD1_g13',
+   sattypes_rad(38)= 'sndrd2_g13',    dsis(38)= 'sndrD2_g13',
+   sattypes_rad(39)= 'sndrd3_g13',    dsis(39)= 'sndrD3_g13',
+   sattypes_rad(40)= 'sndrd4_g13',    dsis(40)= 'sndrD4_g13',
+   sattypes_rad(41)= 'sndrd1_g14',    dsis(41)= 'sndrD1_g14',
+   sattypes_rad(42)= 'sndrd2_g14',    dsis(42)= 'sndrD2_g14',
+   sattypes_rad(43)= 'sndrd3_g14',    dsis(43)= 'sndrD3_g14',
+   sattypes_rad(44)= 'sndrd4_g14',    dsis(44)= 'sndrD4_g14',
+   sattypes_rad(45)= 'sndrd1_g15',    dsis(45)= 'sndrD1_g15',
+   sattypes_rad(46)= 'sndrd2_g15',    dsis(46)= 'sndrD2_g15',
+   sattypes_rad(47)= 'sndrd3_g15',    dsis(47)= 'sndrD3_g15',
+   sattypes_rad(48)= 'sndrd4_g15',    dsis(48)= 'sndrD4_g15',
+   sattypes_rad(49)= 'iasi_metop-a',  dsis(49)= 'iasi_metop-a',
+   sattypes_rad(50)= 'seviri_m08',    dsis(50)= 'seviri_m08',
+   sattypes_rad(51)= 'seviri_m09',    dsis(51)= 'seviri_m09',
+   sattypes_rad(52)= 'seviri_m10',    dsis(52)= 'seviri_m10',
+   sattypes_rad(53)= 'seviri_m11',    dsis(53)= 'seviri_m11',
+   sattypes_rad(54)= 'amsua_metop-b', dsis(54)= 'amsua_metop-b',
+   sattypes_rad(55)= 'hirs4_metop-b', dsis(55)= 'hirs4_metop-b',
+   sattypes_rad(56)= 'mhs_metop-b',   dsis(56)= 'mhs_metop-b',
+   sattypes_rad(57)= 'iasi_metop-b',  dsis(57)= 'iasi_metop-b',
+   sattypes_rad(58)= 'avhrr_metop-b', dsis(58)= 'avhrr3_metop-b',
+   sattypes_rad(59)= 'atms_npp',      dsis(59)= 'atms_npp',
+   sattypes_rad(60)= 'atms_n20',      dsis(60)= 'atms_n20',
+   sattypes_rad(61)= 'cris_npp',      dsis(61)= 'cris_npp',
+   sattypes_rad(62)= 'cris-fsr_npp',  dsis(62)= 'cris-fsr_npp',
+   sattypes_rad(63)= 'cris-fsr_n20',  dsis(63)= 'cris-fsr_n20',
+   sattypes_rad(64)= 'gmi_gpm',       dsis(64)= 'gmi_gpm',
+   sattypes_rad(65)= 'saphir_meghat', dsis(65)= 'saphir_meghat',
+   sattypes_rad(66)= 'amsua_metop-c', dsis(66)= 'amsua_metop-c',
+   sattypes_rad(67)= 'mhs_metop-c',   dsis(67)= 'mhs_metop-c',
+   sattypes_rad(68)= 'ahi_himawari8', dsis(68)= 'ahi_himawari8',
+   sattypes_rad(69)= 'abi_g16',       dsis(69)= 'abi_g16',
+   sattypes_rad(70)= 'abi_g17',       dsis(70)= 'abi_g17',
+   $SATOBS_ENKF
+/
+&ozobs_enkf
+   sattypes_oz(1) = 'sbuv2_n16',
+   sattypes_oz(2) = 'sbuv2_n17',
+   sattypes_oz(3) = 'sbuv2_n18',
+   sattypes_oz(4) = 'sbuv2_n19',
+   sattypes_oz(5) = 'omi_aura',
+   sattypes_oz(6) = 'gome_metop-a',
+   sattypes_oz(7) = 'gome_metop-b',
+   sattypes_oz(8) = 'mls30_aura',
+   sattypes_oz(9) = 'ompsnp_npp',
+   sattypes_oz(10) = 'ompstc8_npp',
+   $OZOBS_ENKF
+/
+EOFnml
+
+################################################################################
+# Run enkf update
+
+export OMP_NUM_THREADS=$NTHREADS_ENKF
+export pgm=$ENKFEXEC
+. prep_step
+
+$NCP $ENKFEXEC $DATA
+$APRUN_ENKF ${DATA}/$(basename $ENKFEXEC) 1>stdout 2>stderr
+rc=$?
+
+export ERR=$rc
+export err=$ERR
+$ERRSCRIPT || exit 2
+
+# save for EFSOI task (still needed?)
+$NCP $COMOUT_ANL_ENS/$GBIASe $COMOUT_ANL_ENSFSOI
+
+# Cat runtime output files.
+cat stdout stderr > $COMOUT_ANL_ENSFSOI/$ENKFSTAT
+
+################################################################################
+#  Postprocessing
+
+cd $pwd
+[[ $mkdata = "YES" ]] && rm -rf $DATA
+set +x
+if [ $VERBOSE = "YES" ]; then
+   echo $(date) EXITING $0 with return code $err >&2
+fi
+exit $err