Merge pull request #415 from RussTreadon-NOAA/develop

GitHub Issue #266. remove global_cycle from DA jobs and scripts
NOAA-EMC · Jun 23, 2022 · 2570e90 · 2570e90
2 parents 483ff97 + 009b74d
commit 2570e90
Show file tree

Hide file tree

Showing 4 changed files with 1 addition and 116 deletions.
diff --git a/jobs/JGDAS_ATMOS_ANALYSIS_DIAG b/jobs/JGDAS_ATMOS_ANALYSIS_DIAG
@@ -110,15 +110,6 @@ if [ $DOHYBVAR = "YES" ]; then
 fi
 
 
-
-# Update surface fields with global_cycle
-export DOGCYCLE=${DOGCYCLE:-"YES"}
-
-
-# Generate Gaussian surface analysis
-export DOGAUSFCANL=${DOGAUSFCANL:-"YES"}
-
-
 ###############################################################
 # Run relevant script
 env

diff --git a/jobs/JGLOBAL_ATMOS_ANALYSIS b/jobs/JGLOBAL_ATMOS_ANALYSIS
@@ -131,14 +131,6 @@ export TCVITL="${COMOUT}/${OPREFIX}syndata.tcvitals.tm00"
 [[ $DONST = "YES" ]] && export NSSTBF="${COMIN_OBS}/${OPREFIX}nsstbufr"
 
 
-# Update surface fields with global_cycle
-export DOGCYCLE=${DOGCYCLE:-"YES"}
-
-
-# Generate Gaussian surface analysis
-export DOGAUSFCANL=${DOGAUSFCANL:-"YES"}
-
-
 # Copy fix file for obsproc
 if [ $RUN = "gfs" ]; then
     mkdir -p $ROTDIR/fix

diff --git a/jobs/JGLOBAL_ATMOS_ANALYSIS_CALC b/jobs/JGLOBAL_ATMOS_ANALYSIS_CALC
@@ -121,11 +121,6 @@ if [ $DOHYBVAR = "YES" ]; then
 fi
 
 
-
-# Update surface fields with global_cycle
-export DOGCYCLE=${DOGCYCLE:-"YES"}
-
-
 # Generate Gaussian surface analysis
 export DOGAUSFCANL=${DOGAUSFCANL:-"YES"}
 

diff --git a/scripts/exglobal_atmos_analysis.sh b/scripts/exglobal_atmos_analysis.sh
@@ -53,29 +53,6 @@ APRUNCFP=${APRUNCFP:-""}
 APRUN_GSI=${APRUN_GSI:-${APRUN:-""}}
 NTHREADS_GSI=${NTHREADS_GSI:-${NTHREADS:-1}}
 
-# Surface cycle related parameters
-DOGCYCLE=${DOGCYCLE:-"NO"}
-CYCLESH=${CYCLESH:-$HOMEgfs/ush/global_cycle.sh}
-export CYCLEXEC=${CYCLEXEC:-$HOMEgfs/exec/global_cycle}
-NTHREADS_CYCLE=${NTHREADS_CYCLE:-24}
-APRUN_CYCLE=${APRUN_CYCLE:-${APRUN:-""}}
-export SNOW_NUDGE_COEFF=${SNOW_NUDGE_COEFF:-'-2.'}
-export CYCLVARS=${CYCLVARS:-""}
-export FHOUR=${FHOUR:-0}
-export DELTSFC=${DELTSFC:-6}
-export FIXgsm=${FIXgsm:-$HOMEgfs/fix/fix_am}
-export FIXfv3=${FIXfv3:-$HOMEgfs/fix/fix_fv3_gmted2010}
-
-DOGAUSFCANL=${DOGAUSFCANL-"NO"}
-GAUSFCANLSH=${GAUSFCANLSH:-$HOMEgfs/ush/gaussian_sfcanl.sh}
-export GAUSFCANLEXE=${GAUSFCANLEXE:-$HOMEgfs/exec/gaussian_sfcanl.exe}
-NTHREADS_GAUSFCANL=${NTHREADS_GAUSFCANL:-1}
-APRUN_GAUSFCANL=${APRUN_GAUSFCANL:-${APRUN:-""}}
-
-# FV3 specific info (required for global_cycle)
-export CASE=${CASE:-"C384"}
-ntiles=${ntiles:-6}
-
 # Microphysics in the model; 99:ZC, 11:GFDLMP
 export imp_physics=${imp_physics:-99}
 lupp=${lupp:-".true."}
@@ -639,7 +616,6 @@ fi
 
 ##############################################################
 # Output files
-# $SFCANL is no longer created here since global_cycle is not called
 $NLN $ATMANL siganl
 $NLN $ATMINC siginc.nc
 if [ $DOHYBVAR = "YES" -a $l4densvar = ".true." -a $lwrite4danl = ".true." ]; then
@@ -990,75 +966,6 @@ if [ $DO_CALC_INCREMENT = "YES" ]; then
   export err=$?; err_chk
 fi
 
-##############################################################
-# Update surface fields in the FV3 restart's using global_cycle
-if [ $DOGCYCLE = "YES" ]; then
-
-    mkdir -p $COMOUT/RESTART
-
-    # Global cycle requires these files
-    export FNTSFA=${FNTSFA:-$COMIN_OBS/${OPREFIX}rtgssthr.grb}
-    export FNACNA=${FNACNA:-$COMIN/${OPREFIX}seaice.5min.blend.grb}
-    export FNSNOA=${FNSNOA:-$COMIN/${OPREFIX}snogrb_t${JCAP_CASE}.${LONB_CASE}.${LATB_CASE}}
-    [[ ! -f $FNSNOA ]] && export FNSNOA="$COMIN/${OPREFIX}snogrb_t1534.3072.1536"
-    FNSNOG=${FNSNOG:-$COMIN_GES/${GPREFIX}snogrb_t${JCAP_CASE}.${LONB_CASE}.${LATB_CASE}}
-    [[ ! -f $FNSNOG ]] && FNSNOG="$COMIN_GES/${GPREFIX}snogrb_t1534.3072.1536"
-
-    # Set CYCLVARS by checking grib date of current snogrb vs that of prev cycle
-    if [ $RUN_GETGES = "YES" ]; then
-        snoprv=$($GETGESSH -q -t snogrb_$JCAP_CASE -e $gesenvir -n $GDUMP -v $GDATE)
-    else
-        snoprv=${snoprv:-$FNSNOG}
-    fi
-
-    if [ $($WGRIB -4yr $FNSNOA 2>/dev/null | grep -i snowc | awk -F: '{print $3}' | awk -F= '{print $2}') -le \
-         $($WGRIB -4yr $snoprv 2>/dev/null | grep -i snowc | awk -F: '{print $3}' | awk -F= '{print $2}') ] ; then
-        export FNSNOA=" "
-        export CYCLVARS="FSNOL=99999.,FSNOS=99999.,"
-    else
-        export SNOW_NUDGE_COEFF=${SNOW_NUDGE_COEFF:-0.}
-        export CYCLVARS="FSNOL=${SNOW_NUDGE_COEFF},$CYCLVARS"
-    fi
-
-    if [ $DONST = "YES" ]; then
-        export NST_FILE=${NST_FILE:-$COMOUT/${APREFIX}dtfanl.nc}
-    else
-        export NST_FILE="NULL"
-    fi
-
-    if [ $DOIAU = "YES" ]; then
-        # update surface restarts at the beginning of the window, if IAU
-        # For now assume/hold dtfanl.nc valid at beginning of window
-        for n in $(seq 1 $ntiles); do
-            $NLN $COMIN_GES/RESTART/$bPDY.${bcyc}0000.sfc_data.tile${n}.nc $DATA/fnbgsi.00$n
-            $NLN $COMOUT/RESTART/$bPDY.${bcyc}0000.sfcanl_data.tile${n}.nc $DATA/fnbgso.00$n
-            $NLN $FIXfv3/$CASE/${CASE}_grid.tile${n}.nc                    $DATA/fngrid.00$n
-            $NLN $FIXfv3/$CASE/${CASE}_oro_data.tile${n}.nc                $DATA/fnorog.00$n
-        done
-
-        export APRUNCY=$APRUN_CYCLE
-        export OMP_NUM_THREADS_CY=$NTHREADS_CYCLE
-        export MAX_TASKS_CY=$ntiles
-
-        $CYCLESH
-        export err=$?; err_chk
-    fi
-    # update surface restarts at middle of window
-    for n in $(seq 1 $ntiles); do
-        $NLN $COMIN_GES/RESTART/$PDY.${cyc}0000.sfc_data.tile${n}.nc $DATA/fnbgsi.00$n
-        $NLN $COMOUT/RESTART/$PDY.${cyc}0000.sfcanl_data.tile${n}.nc $DATA/fnbgso.00$n
-        $NLN $FIXfv3/$CASE/${CASE}_grid.tile${n}.nc                  $DATA/fngrid.00$n
-        $NLN $FIXfv3/$CASE/${CASE}_oro_data.tile${n}.nc              $DATA/fnorog.00$n
-    done
-
-    export APRUNCY=$APRUN_CYCLE
-    export OMP_NUM_THREADS_CY=$NTHREADS_CYCLE
-    export MAX_TASKS_CY=$ntiles
-
-    $CYCLESH
-    export err=$?; err_chk
-fi
-
 
 ##############################################################
 # For eupd
@@ -1105,7 +1012,7 @@ cd $pwd
 if [ $SENDECF = "YES" -a "$RUN" != "enkf" ]; then
    ecflow_client --event release_fcst
 fi
-echo "$CDUMP $CDATE atminc and tiled sfcanl done at `date`" > $COMOUT/${APREFIX}loginc.txt
+echo "$CDUMP $CDATE atminc done at `date`" > $COMOUT/${APREFIX}loginc.txt
 
 ################################################################################
 set +x