Skip to content

Commit b60d23e

Browse files
ADCollardADCollard
authored
DA changes from GFSv16.3.25 (NOAA-EMC#3498)
Script changes to turn on the CADS cloud detection code and the `optconv` modification to improve GSI minimization. These changes are included in the ops package for GFSv16.3.25. Refs NOAA-EMC#2987 Refs NOAA-EMC#2934
1 parent 0a00b98 commit b60d23e

4 files changed

Lines changed: 26 additions & 9 deletions

File tree

‎scripts/exgdas_enkf_update.sh‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -364,6 +364,7 @@ cat > enkf.nml << EOFnml
364364
sattypes_rad(77)= 'viirs-m_j2', dsis(77)= 'viirs-m_j2',
365365
sattypes_rad(78)= 'atms_n21', dsis(78)= 'atms_n21',
366366
sattypes_rad(79)= 'cris-fsr_n21', dsis(79)= 'cris-fsr_n21',
367+
sattypes_rad(80)= 'abi_g19', dsis(80)= 'abi_g19',
367368
${SATOBS_ENKF}
368369
/
369370
&ozobs_enkf

‎scripts/exglobal_atmos_analysis.sh‎

Lines changed: 23 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -50,6 +50,10 @@ cnvw_option=${cnvw_option:-".false."}
5050
# Observation usage options
5151
cao_check=${cao_check:-".true."}
5252
ta2tb=${ta2tb:-".true."}
53+
optconv=${optconv:-0.06}
54+
AIRS_CADS=${AIRS_CADS:-".false."}
55+
IASI_CADS=${IASI_CADS:-".false."}
56+
CRIS_CADS=${CRIS_CADS:-".false."}
5357

5458
# Diagnostic files options
5559
netcdf_diag=${netcdf_diag:-".true."}
@@ -115,7 +119,7 @@ ESIASI=${ESIASI:-${COMIN_OBS}/${OPREFIX}esiasi.tm00.bufr_d${OSUFFIX}}
115119
IASIDB=${IASIDB:-${COMIN_OBS}/${OPREFIX}iasidb.tm00.bufr_d${OSUFFIX}}
116120
AMSREBF=${AMSREBF:-${COMIN_OBS}/${OPREFIX}amsre.tm00.bufr_d${OSUFFIX}}
117121
AMSR2BF=${AMSR2BF:-${COMIN_OBS}/${OPREFIX}amsr2.tm00.bufr_d${OSUFFIX}}
118-
GMI1CRBF=${GMI1CRBF:-${COMIN_OBS}/${OPREFIX}gmi1cr.tm00.bufr_d${OSUFFIX}} # GMI temporarily disabled due to array overflow.
122+
GMI1CRBF=${GMI1CRBF:-${COMIN_OBS}/${OPREFIX}gmi1cr.tm00.bufr_d${OSUFFIX}}
119123
SAPHIRBF=${SAPHIRBF:-${COMIN_OBS}/${OPREFIX}saphir.tm00.bufr_d${OSUFFIX}}
120124
SEVIRIBF=${SEVIRIBF:-${COMIN_OBS}/${OPREFIX}sevcsr.tm00.bufr_d${OSUFFIX}}
121125
AHIBF=${AHIBF:-${COMIN_OBS}/${OPREFIX}ahicsr.tm00.bufr_d${OSUFFIX}}
@@ -146,6 +150,8 @@ TCVITL=${TCVITL:-${COMIN_OBS}/${OPREFIX}syndata.tcvitals.tm00}
146150
B1AVHAM=${B1AVHAM:-${COMIN_OBS}/${OPREFIX}avcsam.tm00.bufr_d${OSUFFIX}}
147151
B1AVHPM=${B1AVHPM:-${COMIN_OBS}/${OPREFIX}avcspm.tm00.bufr_d${OSUFFIX}}
148152
HDOB=${HDOB:-${COMIN_OBS}/${OPREFIX}hdob.tm00.bufr_d${OSUFFIX}}
153+
SAILDRONE=${SAILDRONE:-${COMIN_OBS}/${OPREFIX}saldrn.tm00.bufr_d${OSUFFIX}}
154+
GSBBF=${GSBBF:-${COMIN_OBS}/${OPREFIX}gsbprf.tm00.bufr_d${OSUFFIX}}
149155

150156
# Guess files
151157
GPREFIX=${GPREFIX:-""}
@@ -453,10 +459,10 @@ ${NLN} ${B1AMUB} amsubbufr
453459
${NLN} ${B1MHS} mhsbufr
454460
${NLN} ${ESAMUA} amsuabufrears
455461
${NLN} ${ESAMUB} amsubbufrears
456-
#$NLN $ESMHS mhsbufrears
462+
#$NLN $ESMHS mhsbufrears
457463
${NLN} ${AMUADB} amsuabufr_db
458464
${NLN} ${AMUBDB} amsubbufr_db
459-
#$NLN $MHSDB mhsbufr_db
465+
#$NLN $MHSDB mhsbufr_db
460466
${NLN} ${SBUVBF} sbuvbufr
461467
${NLN} ${OMPSNPBF} ompsnpbufr
462468
${NLN} ${OMPSLPBF} ompslpbufr
@@ -472,7 +478,7 @@ ${NLN} ${ESIASI} iasibufrears
472478
${NLN} ${IASIDB} iasibufr_db
473479
${NLN} ${AMSREBF} amsrebufr
474480
${NLN} ${AMSR2BF} amsr2bufr
475-
#${NLN} ${GMI1CRBF} gmibufr # GMI temporarily disabled due to array overflow.
481+
${NLN} ${GMI1CRBF} gmibufr
476482
${NLN} ${SAPHIRBF} saphirbufr
477483
${NLN} ${SEVIRIBF} seviribufr
478484
${NLN} ${CRISBF} crisbufr
@@ -494,6 +500,8 @@ ${NLN} ${AHIBF} ahibufr
494500
${NLN} ${ABIBF} abibufr
495501
${NLN} ${HDOB} hdobbufr
496502
${NLN} ${SSTVIIRS} sstviirs
503+
${NLN} ${SAILDRONE} sdbufr
504+
${NLN} ${GSBBF} wbbufr
497505

498506
if [[ "${DONST}" == "YES" ]]; then
499507
${NLN} "${NSSTBF}" nsstbufr
@@ -752,7 +760,7 @@ cat > gsiparm.anl << EOF
752760
lobsdiag_forenkf=${lobsdiag_forenkf},
753761
write_fv3_incr=${write_fv3_increment},
754762
nhr_anal=${IAUFHRS},
755-
ta2tb=${ta2tb},
763+
ta2tb=${ta2tb},optconv=${optconv},
756764
${WRITE_INCR_ZERO}
757765
${WRITE_ZERO_STRAT}
758766
${WRITE_STRAT_EFOLD}
@@ -787,7 +795,7 @@ cat > gsiparm.anl << EOF
787795
/
788796
&OBSQC
789797
dfact=0.75,dfact1=3.0,noiqc=.true.,oberrflg=.false.,c_varqc=0.02,
790-
use_poq7=.true.,qc_noirjaco3_pole=.true.,vqc=.false.,nvqc=.true.,
798+
use_poq7=.true.,qc_noirjaco3_pole=.false.,vqc=.false.,nvqc=.true.,
791799
aircraft_t_bc=.true.,biaspredt=1.0e5,upd_aircraft=.true.,cleanup_tail=.true.,
792800
tcp_width=70.0,tcp_ermax=7.35,airs_cads=${AIRS_CADS},cris_cads=${CRIS_CADS},
793801
iasi_cads=${IASI_CADS},blacklst=.true.,
@@ -810,6 +818,13 @@ OBS_INPUT::
810818
prepbufr pw null pw 0.0 0 0
811819
prepbufr uv null uv 0.0 0 0
812820
prepbufr_profl uv null uv 0.0 0 0
821+
wbbufr t null t 0.0 0 0
822+
wbbufr q null q 0.0 0 0
823+
wbbufr uv null uv 0.0 0 0
824+
sdbufr ps null ps 0.0 0 0
825+
sdbufr t null t 0.0 0 0
826+
sdbufr q null q 0.0 0 0
827+
sdbufr uv null uv 0.0 0 0
813828
satwndbufr uv null uv 0.0 0 0
814829
hdobbufr uv null uv 0.0 0 0
815830
prepbufr spd null spd 0.0 0 0
@@ -896,6 +911,7 @@ OBS_INPUT::
896911
abibufr abi g16 abi_g16 0.0 1 0
897912
abibufr abi g17 abi_g17 0.0 1 0
898913
abibufr abi g18 abi_g18 0.0 1 0
914+
abibufr abi g19 abi_g19 0.0 1 0
899915
rapidscatbufr uv null uv 0.0 0 0
900916
ompsnpbufr ompsnp npp ompsnp_npp 0.0 0 0
901917
ompslpbufr ompslp npp ompslp_npp 0.0 0 0
@@ -1008,7 +1024,7 @@ fi
10081024
# Send alerts
10091025
if [[ ${SENDDBN} == "YES" ]]; then
10101026
if [[ ${RUN} == "gfs" ]]; then
1011-
${DBNROOT}/bin/dbn_alert MODEL GFS_abias ${job} ${ABIAS}
1027+
"${DBNROOT}/bin/dbn_alert" MODEL GFS_abias "${job}" "${ABIAS}"
10121028
fi
10131029
fi
10141030

‎scripts/exglobal_diag.sh‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -84,7 +84,7 @@ if [[ "${GENDIAG}" == "YES" ]] ; then
8484
diagtype[0]="conv conv_gps conv_ps conv_pw conv_q conv_sst conv_t conv_tcp conv_uv conv_spd"
8585
diagtype[1]="pcp_ssmi_dmsp pcp_tmi_trmm"
8686
diagtype[2]="sbuv2_n16 sbuv2_n17 sbuv2_n18 sbuv2_n19 gome_metop-a gome_metop-b omi_aura mls30_aura ompsnp_npp ompstc8_npp ompstc8_n20 ompsnp_n20 ompstc8_n21 ompsnp_n21 ompslp_npp gome_metop-c"
87-
diagtype[3]="msu_n14 sndr_g08 sndr_g11 sndr_g12 sndr_g13 sndr_g08_prep sndr_g11_prep sndr_g12_prep sndr_g13_prep sndrd1_g11 sndrd2_g11 sndrd3_g11 sndrd4_g11 sndrd1_g12 sndrd2_g12 sndrd3_g12 sndrd4_g12 sndrd1_g13 sndrd2_g13 sndrd3_g13 sndrd4_g13 sndrd1_g14 sndrd2_g14 sndrd3_g14 sndrd4_g14 sndrd1_g15 sndrd2_g15 sndrd3_g15 sndrd4_g15 amsua_n15 amsua_n16 amsua_n17 amsub_n15 amsub_n16 amsub_n17 hsb_aqua airs_aqua amsua_aqua imgr_g08 imgr_g11 imgr_g12 imgr_g14 imgr_g15 ssmi_f13 ssmi_f15 amsua_n18 amsua_metop-a mhs_n18 mhs_metop-a amsre_low_aqua amsre_mid_aqua amsre_hig_aqua ssmis_f16 ssmis_f17 ssmis_f18 ssmis_f19 ssmis_f20 iasi_metop-a amsua_n19 mhs_n19 seviri_m08 seviri_m09 seviri_m10 seviri_m11 cris_npp cris-fsr_npp cris-fsr_n20 atms_npp atms_n20 amsua_metop-b mhs_metop-b iasi_metop-b avhrr_metop-b avhrr_n18 avhrr_n19 avhrr_metop-a amsr2_gcom-w1 gmi_gpm saphir_meghat ahi_himawari8 abi_g16 abi_g17 amsua_metop-c mhs_metop-c iasi_metop-c avhrr_metop-c viirs-m_npp viirs-m_j1 abi_g18 ahi_himawari9 viirs-m_j2 cris-fsr_n21 atms_n21"
87+
diagtype[3]="msu_n14 sndr_g08 sndr_g11 sndr_g12 sndr_g13 sndr_g08_prep sndr_g11_prep sndr_g12_prep sndr_g13_prep sndrd1_g11 sndrd2_g11 sndrd3_g11 sndrd4_g11 sndrd1_g12 sndrd2_g12 sndrd3_g12 sndrd4_g12 sndrd1_g13 sndrd2_g13 sndrd3_g13 sndrd4_g13 sndrd1_g14 sndrd2_g14 sndrd3_g14 sndrd4_g14 sndrd1_g15 sndrd2_g15 sndrd3_g15 sndrd4_g15 amsua_n15 amsua_n16 amsua_n17 amsub_n15 amsub_n16 amsub_n17 hsb_aqua airs_aqua amsua_aqua imgr_g08 imgr_g11 imgr_g12 imgr_g14 imgr_g15 ssmi_f13 ssmi_f15 amsua_n18 amsua_metop-a mhs_n18 mhs_metop-a amsre_low_aqua amsre_mid_aqua amsre_hig_aqua ssmis_f16 ssmis_f17 ssmis_f18 ssmis_f19 ssmis_f20 iasi_metop-a amsua_n19 mhs_n19 seviri_m08 seviri_m09 seviri_m10 seviri_m11 cris_npp cris-fsr_npp cris-fsr_n20 atms_npp atms_n20 amsua_metop-b mhs_metop-b iasi_metop-b avhrr_metop-b avhrr_n18 avhrr_n19 avhrr_metop-a amsr2_gcom-w1 gmi_gpm saphir_meghat ahi_himawari8 abi_g16 abi_g17 amsua_metop-c mhs_metop-c iasi_metop-c avhrr_metop-c viirs-m_npp viirs-m_j1 abi_g18 ahi_himawari9 viirs-m_j2 cris-fsr_n21 atms_n21 abi_g19"
8888

8989
diaglist[0]=listcnv
9090
diaglist[1]=listpcp

‎sorc/gsi_enkf.fd‎

Submodule gsi_enkf.fd updated 67 files

0 commit comments

Comments
 (0)