diff --git a/fix/gsi/gsiparm.anl.sh b/fix/gsi/gsiparm.anl.sh index 997bd7d1b..e5418bbcd 100755 --- a/fix/gsi/gsiparm.anl.sh +++ b/fix/gsi/gsiparm.anl.sh @@ -189,8 +189,8 @@ OBS_INPUT:: build_cloud_frac_p=0.50, clear_cloud_frac_p=0.10, iclean_hydro_withRef_allcol=1, - i_use_2mQ4B=${ii_use_2mq4b}, - i_use_2mT4B=${ii_use_2mt4b}, + i_use_2mQ4B=${i_use_2mQ4B}, + i_use_2mT4B=${i_use_2mT4B}, i_gsdcldanal_type=0, i_gsdsfc_uselist=1, i_lightpcp=1, @@ -199,8 +199,8 @@ OBS_INPUT:: i_gsdqc=2, / &CHEM - laeroana_fv3smoke=.true., - berror_fv3_cmaq_regional = .true., + laeroana_fv3smoke=${laeroana_fv3smoke}, + berror_fv3_cmaq_regional=${berror_fv3_cmaq_regional}, / &NST / diff --git a/scripts/exregional_run_analysis.sh b/scripts/exregional_run_analysis.sh index 5d9f06ac8..830fe95fc 100755 --- a/scripts/exregional_run_analysis.sh +++ b/scripts/exregional_run_analysis.sh @@ -348,8 +348,8 @@ niter2=50 lread_obs_save=.false. lread_obs_skip=.false. if_model_dbz=.false. -ii_use_2mq4b=2 -ii_use_2mt4b=1 +i_use_2mQ4B=2 +i_use_2mT4B=1 # Determine if hybrid option is available memname='atmf009' @@ -641,6 +641,10 @@ BERROR=${FIX_GSI}/${BERROR_FN} if [[ ${gsi_type} == "ANALYSIS" && ${ob_type} == "AERO" ]]; then + if [ ${BKTYPE} -eq 1 ]; then + echo "cold start, skip GSI SD DA" + exit 0 + fi ANAVINFO=${FIX_GSI}/${ANAVINFO_SD_FN} CONVINFO=${FIX_GSI}/${CONVINFO_SD_FN} BERROR=${FIX_GSI}/${BERROR_SD_FN} @@ -652,11 +656,13 @@ if [[ ${gsi_type} == "ANALYSIS" && ${ob_type} == "AERO" ]]; then l_hyb_ens=.false. nummem=0 beta1_inv=0.0 - ii_use_2mq4b=0 - ii_use_2mt4b=0 + i_use_2mQ4B=0 + i_use_2mT4B=0 netcdf_diag=.true. binary_diag=.false. usenewgfsberror=.false. + laeroana_fv3smoke=.true. + berror_fv3_cmaq_regional=.true. fi SATINFO=${FIX_GSI}/global_satinfo.txt @@ -858,14 +864,6 @@ cat << EOF > gsiparm.anl $gsi_namelist EOF -if [[ ${gsi_type} == "ANALYSIS" && ${ob_type} != "AERO" ]]; then -# base_e GSI can't handle the below two parameter for RRFS-SD DA. So delete them. - mv gsiparm.anl gsiparm.anl.sd - sed '/laeroana_fv3smoke/d' gsiparm.anl.sd > gsiparm.anl.1 - sed '/berror_fv3_cmaq_regional/d' gsiparm.anl.1 > gsiparm.anl - rm -fr gsiparm.anl.sd gsiparm.anl.1 -fi - # #----------------------------------------------------------------------- # @@ -908,39 +906,6 @@ fi #----------------------------------------------------------------------- # # comment out for testing -if [[ ${gsi_type} == "ANALYSIS" && ${ob_type} == "AERO" ]]; then - -if [ ${BKTYPE} -eq 1 ]; then - print_info_msg "$VERBOSE" " - Coldstart skip GSI SD DA" -else -# backup background - if [ "${IO_LAYOUT_Y}" == "1" ]; then - cp -fr ${bkpath}/fv_tracer.res.tile1.nc ${bkpath}/fv_tracer.res.tile1.nc.org - else - for ii in ${list_iolayout} - do - iii=`printf %4.4i $ii` - cp -fr ${bkpath}/fv_tracer.res.tile1.nc.${iii} ${bkpath}/fv_tracer.res.tile1.nc.${iii}.org - done - fi - - $APRUN ./gsi.x < gsiparm.anl > stdout 2>&1 || print_err_msg_exit "\ - Call to executable to run GSI returned with nonzero exit code." - -# copy updated fv3_tracer back to restart - if [ "${IO_LAYOUT_Y}" == "1" ]; then - cp -fr fv3_tracer ${bkpath}/fv_tracer.res.tile1.nc - else - for ii in ${list_iolayout} - do - iii=`printf %4.4i $ii` - cp -fr fv3_tracer.${iii} ${bkpath}/fv_tracer.res.tile1.nc.${iii} - done - fi -fi - -else $APRUN ./gsi.x < gsiparm.anl > stdout 2>&1 || print_err_msg_exit "\ Call to executable to run GSI returned with nonzero exit code." diff --git a/ush/config_defaults.sh b/ush/config_defaults.sh index 5b542ce76..0fdbdd508 100644 --- a/ush/config_defaults.sh +++ b/ush/config_defaults.sh @@ -695,6 +695,9 @@ i_use_2mT4B=0 i_T_Q_adjust=1 l_rtma3d=.false. i_precip_vertical_check=0 +# &CHEM +laeroana_fv3smoke=.false. +berror_fv3_cmaq_regional=.false. #----------------------------------------------------------------------- # HYBENSMEM_NMIN: # Minimum number of ensemble members required a hybrid GSI analysis