Enable EnKF-only for atmosphere #4345
There are no files selected for viewing
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,30 @@ | ||
| experiment: | ||
| net: gfs | ||
| mode: cycled | ||
| pslot: {{ 'pslot' | getenv }} | ||
| app: ATM | ||
| resdetatmos: 96 | ||
| resensatmos: 48 | ||
| comroot: {{ 'RUNTESTS' | getenv }}/COMROOT | ||
| expdir: {{ 'RUNTESTS' | getenv }}/EXPDIR | ||
| icsdir: {{ 'ICSDIR_ROOT' | getenv }}/C96C48/20250808 | ||
| idate: 2024022318 | ||
| edate: 2024022406 | ||
| nens: 2 | ||
| interval: 0 | ||
| start: warm | ||
| yaml: {{ HOMEgfs }}/dev/ci/cases/yamls/gsienkf_atmDA_defaults.ci.yaml | ||
|
|
||
| skip_ci_on_hosts: | ||
| - gaeac5 | ||
| - orion | ||
| - hercules | ||
| - awsepicglobalworkflow | ||
|
|
||
| workflow: | ||
| engine: rocoto | ||
| rocoto: | ||
| maxtries: 2 | ||
| cyclethrottle: 3 | ||
| taskthrottle: 25 | ||
| verbosity: 2 |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,30 @@ | ||
| experiment: | ||
| net: gfs | ||
| mode: cycled | ||
| pslot: {{ 'pslot' | getenv }} | ||
| app: ATM | ||
| resdetatmos: 96 | ||
| resensatmos: 48 | ||
| comroot: {{ 'RUNTESTS' | getenv }}/COMROOT | ||
| expdir: {{ 'RUNTESTS' | getenv }}/EXPDIR | ||
| icsdir: {{ 'ICSDIR_ROOT' | getenv }}/C96C48/20250808 | ||
| idate: 2024022318 | ||
| edate: 2024022406 | ||
| nens: 2 | ||
| interval: 0 | ||
| start: warm | ||
| yaml: {{ HOMEgfs }}/dev/ci/cases/yamls/ufsenkf_atmDA_defaults.ci.yaml | ||
|
|
||
| skip_ci_on_hosts: | ||
| - gaeac5 | ||
| - orion | ||
| - hercules | ||
| - awsepicglobalworkflow | ||
|
|
||
| workflow: | ||
| engine: rocoto | ||
| rocoto: | ||
| maxtries: 2 | ||
| cyclethrottle: 3 | ||
| taskthrottle: 25 | ||
| verbosity: 2 |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,13 @@ | ||
| defaults: | ||
| !INC {{ HOMEgfs }}/dev/parm/config/gfs/yaml/defaults.yaml | ||
| base: | ||
| DOIAU: "NO" | ||
| DO_JEDIATMVAR: "NO" | ||
| DO_JEDIATMENS: "NO" | ||
| DOENKFONLY_ATM: "YES" | ||
| DO_TEST_MODE: "NO" | ||
| CASE_ANL: "C48" | ||
| CASE_HIST: "C96" | ||
| DONST: "NO" | ||
| DO_ARCHCOM: "YES" | ||
| USE_IODADIR: "NO" |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,13 @@ | ||
| defaults: | ||
| !INC {{ HOMEgfs }}/dev/parm/config/gfs/yaml/defaults.yaml | ||
| base: | ||
| DOIAU: "NO" | ||
| DO_JEDIATMVAR: "NO" | ||
| DO_JEDIATMENS: "YES" | ||
| DOENKFONLY_ATM: "YES" | ||
| DO_TEST_MODE: "NO" | ||
| CASE_ANL: "C48" | ||
| CASE_HIST: "C96" | ||
| DONST: "NO" | ||
| DO_ARCHCOM: "YES" | ||
| USE_IODADIR: "YES" |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,18 @@ | ||
| #! /usr/bin/env bash | ||
|
|
||
| ############################################################### | ||
| # Source UFSDA workflow modules | ||
| source "${HOMEgfs}/dev/ush/load_modules.sh" ufsda | ||
| status=$? | ||
| if [[ ${status} -ne 0 ]]; then | ||
| exit "${status}" | ||
| fi | ||
|
|
||
| export job="prepatmanlbias" | ||
| export jobid="${job}.$$" | ||
|
|
||
| ############################################################### | ||
| # Execute the JJOB | ||
| "${HOMEgfs}/dev/jobs/JGLOBAL_ATM_PREP_ANL_BIAS" | ||
| status=$? | ||
| exit "${status}" |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,137 @@ | ||
| #! /usr/bin/env bash | ||
| export WIPE_DATA="NO" | ||
| export DATA=${DATA:-${DATAROOT}/${RUN}fetch_${cyc}} | ||
| source "${HOMEgfs}/ush/jjob_header.sh" -e "prepatmanlbias" -c "base prepatmanlbias" | ||
|
|
||
| ############################################## | ||
| # Set variables used in the script | ||
| ############################################## | ||
|
|
||
| ############################################## | ||
| # Begin JOB SPECIFIC work | ||
| ############################################## | ||
| # Ignore possible spelling error (nothing is misspelled) | ||
| # shellcheck disable=SC2153 | ||
| GDATE=$(date --utc +%Y%m%d%H -d "${PDY} ${cyc} - ${assim_freq} hours") | ||
| export gPDY=${GDATE:0:8} | ||
| export gcyc=${GDATE:8:2} | ||
| export gyy=${GDATE:0:4} | ||
| export gmm=${GDATE:4:2} | ||
| export gdd=${GDATE:6:2} | ||
| export GDUMP="gdas" | ||
|
|
||
| # Generate COM variables from templates | ||
| MEMDIR="ensstat" RUN="enkf${GDUMP}" YMD=${gPDY} HH=${gcyc} declare_from_tmpl -rx \ | ||
| COMOUT_ATMOS_ANALYSIS_PREV:COM_ATMOS_ANALYSIS_TMPL | ||
|
|
||
| # Process bias corrections files grabbed from HPSS in task gdas_fetchatmanlbias | ||
| BIASDIR=${DATA}/${GDUMP}.${gPDY}/${gcyc}/atmos | ||
| BIASFILES="${GDUMP}.t${gcyc}z.abias ${GDUMP}.t${gcyc}z.abias_air ${GDUMP}.t${gcyc}z.abias_int ${GDUMP}.t${gcyc}z.abias_pc ${GDUMP}.t${gcyc}z.radstat" | ||
|
|
||
| ABIAS_SAT=${BIASDIR}/${GDUMP}.t${gcyc}z.abias | ||
| ABIASPC_SAT=${BIASDIR}/${GDUMP}.t${gcyc}z.abias_pc | ||
| ABIAS_SAT_JEDI_TAR=${COMOUT_ATMOS_ANALYSIS_PREV}/${GDUMP}.t${gcyc}z.rad_varbc_params.tar | ||
| ABIAS_AIR_JEDI_TAR=${COMOUT_ATMOS_ANALYSIS_PREV}/${GDUMP}.t${gcyc}z.air_varbc_params.tar | ||
|
|
||
| SATBIAS2IODAX=${SATBIAS2IODAX:-${EXECgfs}/satbias2ioda.x} | ||
| SATBIAS2IODAY=${HOMEgfs}/sorc/gdas.cd/ush/satbias_converter.yaml.tmpl | ||
| OUTPUT=${DATA}/output | ||
|
|
||
| mkdir -p "${COMOUT_ATMOS_ANALYSIS_PREV}" | ||
| mkdir -p "${OUTPUT}" | ||
|
|
||
| if [[ ${DO_JEDIATMENS} == "YES" ]]; then | ||
| # JEDI run | ||
| cd "${DATA}" || exit 1 | ||
| mkdir -p "testrun/varbc" | ||
| ln -sf "${ABIAS_SAT}" "./satbias_in" | ||
| ln -sf "${ABIASPC_SAT}" "./satbias_pc" | ||
| if grep -i "NaN" satbias_in > /dev/null 2>&1; then | ||
| ##export err=1 | ||
| ##err_chk "FATAL ERROR There are NaN in satbias_in." | ||
| echo "WARNING: There are NaN in satbias_in. Replace NaN with 0.000000" | ||
| cp satbias_in satbias_in_original | ||
| sed -i -e 's/ NaN/ 0.000000/g' \ | ||
| -e 's/ nan/ 0.000000/g' \ | ||
| -e 's/ NAN/ 0.000000/g' satbias_in | ||
| fi | ||
|
|
||
| # Get instruments from satbias_in | ||
| obsclass=$(grep "_" satbias_in | awk '{print $2}' | uniq) | ||
|
|
||
| # Loop over instruments. Covert GSI abias to JEDI format | ||
| for instrument in ${obsclass}; do | ||
| echo "${instrument}" | ||
| cp -f "${SATBIAS2IODAY}" satbias_converter.yaml | ||
| sed -i -e "s/INSTRUMENT/${instrument}/g" satbias_converter.yaml | ||
|
|
||
| ${SATBIAS2IODAX} satbias_converter.yaml | ||
|
|
||
| export err=$? | ||
| if [[ ${err} -ne 0 ]]; then | ||
| err_exit "FATAL ERROR satbias2iodas.x failed for ${instrument}, ABORT!" | ||
| fi | ||
|
|
||
| rm -f testrun/varbc/*nc | ||
| rm -f testrun/varbc/*txt | ||
|
|
||
| cd ./testrun/varbc/ | ||
| grep "${instrument}" ../../satbias_in | awk '{print $2" "$3" "$4}' > \ | ||
| "${OUTPUT}"/gdas.t"${gcyc}"z.radiance_"${instrument}".tlapse.txt | ||
| cpreq -p satbias_"${instrument}".nc4 "${OUTPUT}"/gdas.t"${gcyc}"z.radiance_"${instrument}".satbias.nc | ||
| mv satbias_"${instrument}".nc4 "${OUTPUT}"/gdas.t"${gcyc}"z.radiance_"${instrument}".satbias_cov.nc | ||
| cd "${DATA}" | ||
| rm -f satbias_converter.yaml | ||
| done | ||
|
|
||
| if [[ ${err} -eq 0 ]]; then | ||
| cd "${OUTPUT}" || exit | ||
| for file in ${BIASFILES}; do | ||
| cpreq "${BIASDIR}/${file}" "./enkf${file}.txt" | ||
| export err=$? | ||
| if [[ ${err} -ne 0 ]]; then | ||
| err_exit "Error copying operational anl bias correction files" | ||
| fi | ||
| done | ||
| # TODO: Temporary solution and need to process abias_air for JEDI when ready | ||
| tar -cvf "${ABIAS_AIR_JEDI_TAR}" "./enkf${GDUMP}.t${gcyc}z.abias_air.txt" | ||
| rm -f "./enkf${GDUMP}.t${gcyc}z.abias_air.txt" | ||
| tar -cvf "${ABIAS_SAT_JEDI_TAR}" ./*.nc ./*.txt | ||
| export err=$? | ||
| fi | ||
|
|
||
| else | ||
| # GSI run | ||
| for file in ${BIASFILES}; do | ||
| cpreq "${BIASDIR}/${file}" "${COMOUT_ATMOS_ANALYSIS_PREV}/enkf${file}.txt" | ||
| export err=$? | ||
| if [[ ${err} -ne 0 ]]; then | ||
| err_exit "FATAL ERROR copying operational anl bias correction files" | ||
| fi | ||
| done | ||
| fi | ||
|
|
||
| if [[ ${err} -ne 0 ]]; then | ||
| err_exit "FATAL ERROR executing ${EXSCRIPT}" | ||
| fi | ||
| set_trace | ||
|
|
||
| ############################################## | ||
| # End JOB SPECIFIC work | ||
| ############################################## | ||
|
|
||
| ############################################## | ||
| # Final processing | ||
| ############################################## | ||
| if [[ -e "${pgmout}" ]]; then | ||
| cat "${pgmout}" | ||
| fi | ||
|
|
||
| ################################### | ||
| # Remove temp directories | ||
| ################################### | ||
| if [[ "${KEEPDATA}" != "YES" ]]; then | ||
| rm -rf "${DATA}" | ||
| fi | ||
|
|
||
| exit 0 | ||
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
|
|
@@ -396,6 +396,10 @@ export BUILD_GSINFO_DIR="${PARMgfs}/gsinfo" | |
| export NMEM_ENS_GFS="{{ NMEM_ENS_GFS }}" | ||
| export NMEM_ENS_GFS_OFFSET="{{ NMEM_ENS_GFS_OFFSET }}" | ||
| export DO_CALC_INCREMENT_ENKF_GFS="NO" | ||
| export DOENKFONLY_ATM="{{ DOENKFONLY_ATM }}" | ||
| if [[ "${DOENKFONLY_ATM}" = "YES" && "${DO_JEDIATMENS}" = "YES" ]]; then | ||
| export DO_HISTORY_FILE_ON_NATIVE_GRID="YES" | ||
| fi | ||
|
|
||
| # EnKF output frequency | ||
| if [[ "${DOHYBVAR}" = "YES" ]]; then | ||
|
|
@@ -519,6 +523,17 @@ if [[ "${machine}" == "URSA" || "${machine}" == "GAEAC6" || "${machine}" == "ORI | |
| export DO_METP=NO | ||
| fi | ||
|
|
||
| if [[ "${DOENKFONLY_ATM:-NO}" == "YES" ]] ; then | ||
| export RECENTER_ENKF="NO" # Turn off recentering ensemble analysis | ||
| export DO_VERFOZN="NO" # Ozone data assimilation monitoring | ||
| export DO_VERFRAD="NO" # Radiance data assimilation monitoring | ||
| export DO_VMINMON="NO" # GSI minimization monitoring | ||
| export DO_METP="NO" # Run METPLUS jobs - set METPLUS settings in config.metp | ||
| export DO_FIT2OBS="NO" # Run fit to observations package | ||
| export DO_TRACKER="NO" # Hurricane track verification | ||
| export DO_GENESIS="NO" # Cyclone genesis verification | ||
| fi | ||
|
Comment on lines
+526
to
+535
Contributor
There was a problem hiding this comment. @aerorahul I wonder if this would be better to have in
Contributor
There was a problem hiding this comment. I think these can stay here for now. |
||
|
|
||
| # If starting ICs that are not at cycle hour | ||
| export OFFSET_START_HOUR=0 | ||
|
|
||
|
|
||
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Oops, something went wrong.
Oops, something went wrong.
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Uh oh!
There was an error while loading. Please reload this page.