Skip to content

Commit c902645

Browse files
Update jgdas enkf update job with COMIN or COMOUT prefix instead of COM (NOAA-EMC#3333)
NCO has requested that each COM variable specify whether it is an input or an output. This completes that process for the global JGDAS_ENKF_UPDATE job and associated bash and python scripts Refs NOAA-EMC#2451
1 parent 6a287b7 commit c902645

File tree

2 files changed

+29
-26
lines changed

2 files changed

+29
-26
lines changed

jobs/JGDAS_ENKF_UPDATE

Lines changed: 6 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -25,10 +25,13 @@ export APREFIX="${RUN}.t${cyc}z."
2525
export GPREFIX="${GDUMP_ENS}.t${gcyc}z."
2626

2727
MEMDIR="ensstat" YMD=${PDY} HH=${cyc} declare_from_tmpl -rx \
28-
COM_ATMOS_ANALYSIS_STAT:COM_ATMOS_ANALYSIS_TMPL
28+
COMIN_ATMOS_ANALYSIS_STAT:COM_ATMOS_ANALYSIS_TMPL
29+
30+
MEMDIR="ensstat" YMD=${PDY} HH=${cyc} declare_from_tmpl -rx \
31+
COMOUT_ATMOS_ANALYSIS_STAT:COM_ATMOS_ANALYSIS_TMPL
2932

3033
MEMDIR="ensstat" RUN="enkfgdas" YMD=${gPDY} HH=${gcyc} declare_from_tmpl -rx \
31-
COM_ATMOS_HISTORY_STAT_PREV:COM_ATMOS_HISTORY_TMPL
34+
COMIN_ATMOS_HISTORY_STAT_PREV:COM_ATMOS_HISTORY_TMPL
3235

3336

3437
###############################################################
@@ -43,7 +46,7 @@ status=$?
4346
# Send Alerts
4447
##############################################
4548
if [ ${SENDDBN} = YES ] ; then
46-
"${DBNROOT}/bin/dbn_alert" "MODEL" "ENKF1_MSC_enkfstat" "${job}" "${COM_ATMOS_ANALYSIS_STAT}/${APREFIX}enkfstat"
49+
"${DBNROOT}/bin/dbn_alert" "MODEL" "ENKF1_MSC_enkfstat" "${job}" "${COMOUT_ATMOS_ANALYSIS_STAT}/${APREFIX}enkfstat"
4750
fi
4851

4952

scripts/exgdas_enkf_update.sh

Lines changed: 23 additions & 23 deletions
Original file line numberDiff line numberDiff line change
@@ -95,7 +95,7 @@ GSI_SOILANAL=${GSI_SOILANAL:-"NO"}
9595

9696
################################################################################
9797

98-
ATMGES_ENSMEAN="${COM_ATMOS_HISTORY_STAT_PREV}/${GPREFIX}atmf006.ensmean.nc"
98+
ATMGES_ENSMEAN="${COMIN_ATMOS_HISTORY_STAT_PREV}/${GPREFIX}atmf006.ensmean.nc"
9999
LONB_ENKF=${LONB_ENKF:-$($NCLEN $ATMGES_ENSMEAN grid_xt)} # get LONB_ENKF
100100
LATB_ENKF=${LATB_ENKF:-$($NCLEN $ATMGES_ENSMEAN grid_yt)} # get LATB_ENFK
101101
LEVS_ENKF=${LEVS_ENKF:-$($NCLEN $ATMGES_ENSMEAN pfull)} # get LEVS_ENFK
@@ -142,7 +142,7 @@ $NLN $ANAVINFO anavinfo
142142
$NLN $VLOCALEIG vlocal_eig.dat
143143

144144
# Bias correction coefficients based on the ensemble mean
145-
${NLN} "${COM_ATMOS_ANALYSIS_STAT}/${GBIASe}" "satbias_in"
145+
${NLN} "${COMIN_ATMOS_ANALYSIS_STAT}/${GBIASe}" "satbias_in"
146146

147147
################################################################################
148148

@@ -152,13 +152,13 @@ if [ $USE_CFP = "YES" ]; then
152152
cat > $DATA/untar.sh << EOFuntar
153153
#!/bin/sh
154154
memchar=\$1
155-
COM_ATMOS_ANALYSIS=\$2
155+
COMOUT_ATMOS_ANALYSIS=\$2
156156
flist="$CNVSTAT $OZNSTAT $RADSTAT"
157157
for ftype in \$flist; do
158158
if [ \$memchar = "ensmean" ]; then
159-
fname=\${COM_ATMOS_ANALYSIS}/\${ftype}.ensmean
159+
fname=\${COMOUT_ATMOS_ANALYSIS}/\${ftype}.ensmean
160160
else
161-
fname=\${COM_ATMOS_ANALYSIS}/\${ftype}
161+
fname=\${COMOUT_ATMOS_ANALYSIS}/\${ftype}
162162
fi
163163
tar -xvf \$fname
164164
done
@@ -171,13 +171,13 @@ fi
171171

172172
flist="$CNVSTAT $OZNSTAT $RADSTAT"
173173
if [ $USE_CFP = "YES" ]; then
174-
echo "${nm} ${DATA}/untar.sh ensmean ${COM_ATMOS_ANALYSIS_STAT}" | tee -a "${DATA}/mp_untar.sh"
174+
echo "${nm} ${DATA}/untar.sh ensmean ${COMIN_ATMOS_ANALYSIS_STAT}" | tee -a "${DATA}/mp_untar.sh"
175175
if [ ${CFP_MP:-"NO"} = "YES" ]; then
176176
nm=$((nm+1))
177177
fi
178178
else
179179
for ftype in $flist; do
180-
fname="${COM_ATMOS_ANALYSIS_STAT}/${ftype}.ensmean"
180+
fname="${COMIN_ATMOS_ANALYSIS_STAT}/${ftype}.ensmean"
181181
tar -xvf $fname
182182
done
183183
fi
@@ -191,55 +191,55 @@ for imem in $(seq 1 $NMEM_ENS); do
191191
memchar="mem"$(printf %03i $imem)
192192

193193
MEMDIR=${gmemchar} RUN=${GDUMP_ENS} YMD=${gPDY} HH=${gcyc} declare_from_tmpl -x \
194-
COM_ATMOS_HISTORY_MEM_PREV:COM_ATMOS_HISTORY_TMPL
194+
COMIN_ATMOS_HISTORY_MEM_PREV:COM_ATMOS_HISTORY_TMPL
195195

196196
MEMDIR=${memchar} YMD=${PDY} HH=${cyc} declare_from_tmpl -x \
197-
COM_ATMOS_ANALYSIS_MEM:COM_ATMOS_ANALYSIS_TMPL
197+
COMOUT_ATMOS_ANALYSIS_MEM:COM_ATMOS_ANALYSIS_TMPL
198198

199199
if [ $lobsdiag_forenkf = ".false." ]; then
200200
if [ $USE_CFP = "YES" ]; then
201-
echo "${nm} ${DATA}/untar.sh ${memchar} ${COM_ATMOS_ANALYSIS_MEM}" | tee -a "${DATA}/mp_untar.sh"
201+
echo "${nm} ${DATA}/untar.sh ${memchar} ${COMOUT_ATMOS_ANALYSIS_MEM}" | tee -a "${DATA}/mp_untar.sh"
202202
if [ ${CFP_MP:-"NO"} = "YES" ]; then
203203
nm=$((nm+1))
204204
fi
205205
else
206206
for ftype in $flist; do
207-
fname="${COM_ATMOS_ANALYSIS_MEM}/${ftype}"
207+
fname="${COMOUT_ATMOS_ANALYSIS_MEM}/${ftype}"
208208
tar -xvf $fname
209209
done
210210
fi
211211
fi
212-
mkdir -p "${COM_ATMOS_ANALYSIS_MEM}"
212+
mkdir -p "${COMOUT_ATMOS_ANALYSIS_MEM}"
213213
for FHR in $nfhrs; do
214-
${NLN} "${COM_ATMOS_HISTORY_MEM_PREV}/${GPREFIX}atmf00${FHR}${ENKF_SUFFIX}.nc" \
214+
${NLN} "${COMIN_ATMOS_HISTORY_MEM_PREV}/${GPREFIX}atmf00${FHR}${ENKF_SUFFIX}.nc" \
215215
"sfg_${PDY}${cyc}_fhr0${FHR}_${memchar}"
216216
if [ $GSI_SOILANAL = "YES" ]; then
217-
${NLN} "${COM_ATMOS_HISTORY_MEM_PREV}/${GPREFIX}sfcf00${FHR}${ENKF_SUFFIX}.nc" \
217+
${NLN} "${COMIN_ATMOS_HISTORY_MEM_PREV}/${GPREFIX}sfcf00${FHR}${ENKF_SUFFIX}.nc" \
218218
"bfg_${PDY}${cyc}_fhr0${FHR}_${memchar}"
219219
fi
220220
if [ $cnvw_option = ".true." ]; then
221-
${NLN} "${COM_ATMOS_HISTORY_MEM_PREV}/${GPREFIX}sfcf00${FHR}.nc" \
221+
${NLN} "${COMIN_ATMOS_HISTORY_MEM_PREV}/${GPREFIX}sfcf00${FHR}.nc" \
222222
"sfgsfc_${PDY}${cyc}_fhr0${FHR}_${memchar}"
223223
fi
224224
if [ $FHR -eq 6 ]; then
225225
if [ $DO_CALC_INCREMENT = "YES" ]; then
226-
${NLN} "${COM_ATMOS_ANALYSIS_MEM}/${APREFIX}atmanl.nc" \
226+
${NLN} "${COMOUT_ATMOS_ANALYSIS_MEM}/${APREFIX}atmanl.nc" \
227227
"sanl_${PDY}${cyc}_fhr0${FHR}_${memchar}"
228228
else
229-
${NLN} "${COM_ATMOS_ANALYSIS_MEM}/${APREFIX}atminc.nc" \
229+
${NLN} "${COMOUT_ATMOS_ANALYSIS_MEM}/${APREFIX}atminc.nc" \
230230
"incr_${PDY}${cyc}_fhr0${FHR}_${memchar}"
231231
fi
232232
else
233233
if [ $DO_CALC_INCREMENT = "YES" ]; then
234-
${NLN} "${COM_ATMOS_ANALYSIS_MEM}/${APREFIX}atma00${FHR}.nc" \
234+
${NLN} "${COMOUT_ATMOS_ANALYSIS_MEM}/${APREFIX}atma00${FHR}.nc" \
235235
"sanl_${PDY}${cyc}_fhr0${FHR}_${memchar}"
236236
else
237-
${NLN} "${COM_ATMOS_ANALYSIS_MEM}/${APREFIX}atmi00${FHR}.nc" \
237+
${NLN} "${COMOUT_ATMOS_ANALYSIS_MEM}/${APREFIX}atmi00${FHR}.nc" \
238238
"incr_${PDY}${cyc}_fhr0${FHR}_${memchar}"
239239
fi
240240
fi
241241
if [ $GSI_SOILANAL = "YES" ]; then
242-
${NLN} "${COM_ATMOS_ANALYSIS_MEM}/${APREFIX}sfci00${FHR}.nc" \
242+
${NLN} "${COMOUT_ATMOS_ANALYSIS_MEM}/${APREFIX}sfci00${FHR}.nc" \
243243
"sfcincr_${PDY}${cyc}_fhr0${FHR}_${memchar}"
244244
fi
245245
done
@@ -248,10 +248,10 @@ done
248248
# Ensemble mean guess
249249
for FHR in $nfhrs; do
250250

251-
${NLN} "${COM_ATMOS_HISTORY_STAT_PREV}/${GPREFIX}atmf00${FHR}.ensmean.nc" \
251+
${NLN} "${COMIN_ATMOS_HISTORY_STAT_PREV}/${GPREFIX}atmf00${FHR}.ensmean.nc" \
252252
"sfg_${PDY}${cyc}_fhr0${FHR}_ensmean"
253253
if [ $cnvw_option = ".true." ]; then
254-
${NLN} "${COM_ATMOS_HISTORY_STAT_PREV}/${GPREFIX}sfcf00${FHR}.ensmean.nc" \
254+
${NLN} "${COMIN_ATMOS_HISTORY_STAT_PREV}/${GPREFIX}sfcf00${FHR}.ensmean.nc" \
255255
"sfgsfc_${PDY}${cyc}_fhr0${FHR}_ensmean"
256256
fi
257257
done
@@ -412,7 +412,7 @@ $APRUN_ENKF ${DATA}/$(basename $ENKFEXEC) 1>stdout 2>stderr
412412
export err=$?; err_chk
413413

414414
# Cat runtime output files.
415-
cat stdout stderr > "${COM_ATMOS_ANALYSIS_STAT}/${ENKFSTAT}"
415+
cat stdout stderr > "${COMOUT_ATMOS_ANALYSIS_STAT}/${ENKFSTAT}"
416416

417417
################################################################################
418418
# Postprocessing

0 commit comments

Comments
 (0)