add RTMA-specific cloud analysis parameters, update RTMA config files…

…, bug fix (NOAA-EMC#413)
EricJames-NOAA · Aug 3, 2022 · e51a45f · e51a45f
1 parent 84ed595
commit e51a45f
Show file tree

Hide file tree

Showing 5 changed files with 45 additions and 10 deletions.
diff --git a/scripts/exregional_nonvarcldanl.sh b/scripts/exregional_nonvarcldanl.sh
@@ -316,6 +316,9 @@ cat << EOF > gsiparm.anl
  l_saturate_bkCloud=.true.,
  l_qnr_from_qr=${l_qnr_from_qr},
  n0_rain=100000000.0
+ i_T_Q_adjust=${i_T_Q_adjust},
+ l_rtma3d=${l_rtma3d},
+ i_precip_vertical_check=${i_precip_vertical_check},
  /
 EOF
 

diff --git a/ush/config.sh.RTMA_CONUS_3km b/ush/config.sh.RTMA_CONUS_3km
@@ -3,6 +3,7 @@ source fix_rrfs_locations.sh #determine current platform
 MACHINE=${PLATFORM}
 ACCOUNT="nrtrr"
 RESERVATION="rrfsdet"
+RESERVATION_POST="rrfsdet"
 EXPT_BASEDIR="/home/rtrr/RTMA"
 EXPT_SUBDIR="RTMA_CONUS"
 
@@ -14,13 +15,17 @@ if [[ -n $RESERVATION ]] ; then
  PARTITION_GRAPHICS=kjet
  PARTITION_ANALYSIS=kjet
  PARTITION_WGRIB2=kjet
+ PARTITION_POST=kjet
 
  QUEUE_DEFAULT="rt"
  QUEUE_FCST="rt"
  QUEUE_ANALYSIS="rt"
  QUEUE_WGRIB2="rt"
  QUEUE_GRAPHICS="rt"
  QUEUE_HPSS="rt"
+ if [[ -n ${RESERVATION_POST} ]] ; then
+ QUEUE_POST="rt"
+ fi
 fi
 
 OBSPATH_NSSLMOSIAC=/public/data/radar/nssl/mrms/conus
@@ -85,16 +90,21 @@ readin_localization=.true. #if true, it overwrites the "beta1_inv/ens_h/ens_
 beta1_inv=0.5 #beata_inv is 1-ensemble_wgt
 ens_h=20 #110
 ens_v=1 #3
-regional_ensemble_option=1 #1 for GDAS
+regional_ensemble_option=5 #1 for GDAS, 5 for RRFSDAS
 grid_ratio_fv3=2.0 #fv3 resolution 3km, so analysis=3*2=6km
 grid_ratio_ens=3 #analysis 6km, so ensemble=6*3=18km. GDAS ensemble is 20km
+##MEMO1: regional_ensemble_option=5 will be changed to 1 if there are not enough RRFS ensembles at the runtime.
+##MEMO2: grid_ratio_ens will be hardwired to 1 for RRFSDAS, so grid_ratio_ens here only affects GDAS.
 i_en_perts_io=1 #0 or 1: original file 3: pre-processed ensembles
 
 # &RAPIDREFRESH_CLDSURF
 l_PBL_pseudo_SurfobsT=.false.
 l_PBL_pseudo_SurfobsQ=.false
 i_use_2mQ4B=0
 i_use_2mT4B=0
+i_T_Q_adjust=1
+l_rtma3d=.true.
+i_precip_vertical_check=3
 #-------------------------------------------------------------------------------------
 #
 BERROR_FN="berror.rtma" #under $FIX_GSI
@@ -116,8 +126,8 @@ EXTRN_MDL_LBCS_SEARCH_OFFSET_HRS="0"
 BOUNDARY_LEN_HRS="1"
 BOUNDARY_LONG_LEN_HRS="1"
 
-DATE_FIRST_CYCL="20220406"
-DATE_LAST_CYCL="20220626"
+DATE_FIRST_CYCL="20220801"
+DATE_LAST_CYCL="20221231"
 CYCL_HRS=( "00" "12" )
 CYCL_HRS_SPINSTART=("03" "15")
 CYCL_HRS_PRODSTART=("99") # RTMA always warm-starts
@@ -144,6 +154,9 @@ for i in {0..23}; do FCST_LEN_HRS_CYCLES[$i]=0.0003; done
 DA_CYCLE_INTERV="1" #need to be 1 to have hourly DA
 RESTART_INTERVAL="1"
 
+USE_RRFSE_ENS="TRUE"
+CYCL_HRS_HYB_FV3LAM_ENS=("19" "20" "21" "22" "23" "00")
+
 NCORES_RUN_ANAL=240
 WTIME_RUN_ANAL="01:00:00"
 HYBENSMEM_NMIN=66
@@ -221,11 +234,16 @@ OBSPATH="/lfs4/BMC/wrfruc/gge/rtma/wcoss_obs"
 STMP="/lfs4/BMC/nrtrr/NCO_dirs/stmp" # Path to directory STMP that mostly contains input files.
 PTMP="/lfs4/BMC/nrtrr/NCO_dirs/ptmp" # Path to directory STMP that mostly contains input files.
 NWGES="/lfs4/BMC/nrtrr/NCO_dirs/nwges" # Path to directory NWGES that save boundary, cold initial, restart files
+if [[ ${regional_ensemble_option} == "5" ]]; then
+ RRFSE_NWGES="/lfs4/BMC/wrfruc/RRFSE/NCO_dirs/nwges" # Path to RRFSE directory NWGES that mostly contains ensemble restart files for GSI hybrid.
+ NUM_ENS_MEMBERS=30 # FV3LAM ensemble size for GSI hybrid analysis
+fi
 
 # change settings if it is wcoss_dell_p3
 if [[ "$MACHINE" == "wcoss_dell_p3" ]]; then
  ACCOUNT="RTMA-T2O"
  RESERVATION=""
+ RESERVATION_POST=""
  EXPT_BASEDIR="/gpfs/dell2/emc/modeling/save/${USER}/RTMA-CONUS"
  SERVICE_ACCOUNT=""
  PARTITION_DEFAULT=""
@@ -238,6 +256,10 @@ if [[ "$MACHINE" == "wcoss_dell_p3" ]]; then
  OBSPATH_NSSLMOSIAC="/gpfs/dell1/nco/ops/dcom/prod/ldmdata/obs/upperair/mrms/conus/MergedReflectivityQC"
  AIRCRAFT_REJECT="/gpfs/dell2/emc/modeling/noscrub/Edward.Colon/FixData/obsuselist/amdar_reject_lists"
  SFCOBS_USELIST="/gpfs/dell2/emc/modeling/noscrub/Edward.Colon/FixData/obsuselist/mesonet_uselists"
+
+ regional_ensemble_option=1
+ USE_RRFSE_ENS="FALSE"
+ CYCL_HRS_HYB_FV3LAM_ENS=("99")
 
  WTIME_RUN_ANAL="00:30:00"
  NNODES_RUN_ANAL="28"
@@ -262,7 +284,4 @@ if [[ "$MACHINE" == "wcoss_dell_p3" ]]; then
  PTMP="/gpfs/dell2/ptmp/${USER}/ptmp"
  NWGES="/gpfs/dell2/ptmp/${USER}/nwges"
 fi
-# Wait longer for rrfsdas
-#if [[ "$MACHINE" == "jet" ]] && [[ $regional_ensemble_option -eq 5 ]] ; then
-# START_TIME_ANAL="5:00:00"
-#fi
+
diff --git a/ush/config.sh.RTMA_NA_3km b/ush/config.sh.RTMA_NA_3km
@@ -14,6 +14,7 @@ if [[ -n $RESERVATION ]] ; then
  PARTITION_GRAPHICS=kjet
  PARTITION_ANALYSIS=kjet
  PARTITION_WGRIB2=kjet
+ PARTITION_POST=kjet
 
  QUEUE_DEFAULT="rt"
  QUEUE_FCST="rt"
@@ -88,16 +89,21 @@ readin_localization=.true. #if true, it overwrites the "beta1_inv/ens_h/ens_
 beta1_inv=0.5 #beata_inv is 1-ensemble_wgt
 ens_h=20 #110
 ens_v=1 #3
-regional_ensemble_option=1 #1 for GDAS
+regional_ensemble_option=1 #1 for GDAS, 5 for RRFSDAS
 grid_ratio_fv3=2.0 #fv3 resolution 3km, so analysis=3*2=6km
 grid_ratio_ens=3 #analysis 6km, so ensemble=6*3=18km. GDAS ensemble is 20km
+##MEMO1: regional_ensemble_option=5 will be changed to 1 if there are not enough RRFS ensembles at the runtime.
+##MEMO2: grid_ratio_ens will be hardwired to 1 for RRFSDAS, so grid_ratio_ens here only affects GDAS.
 i_en_perts_io=1 #0 or 1: original file 3: pre-processed ensembles
 
 # &RAPIDREFRESH_CLDSURF
 l_PBL_pseudo_SurfobsT=.false.
 l_PBL_pseudo_SurfobsQ=.false
 i_use_2mQ4B=0
 i_use_2mT4B=0
+i_T_Q_adjust=1
+l_rtma3d=.true.
+i_precip_vertical_check=3
 #-------------------------------------------------------------------------------------
 #
 BERROR_FN="berror.rtma" #under $FIX_GSI

diff --git a/ush/config_defaults.sh b/ush/config_defaults.sh
@@ -633,6 +633,9 @@ l_PBL_pseudo_SurfobsT=.false.
 l_PBL_pseudo_SurfobsQ=.false.
 i_use_2mQ4B=0
 i_use_2mT4B=0
+i_T_Q_adjust=1
+l_rtma3d=.false.
+i_precip_vertical_check=0
 #-----------------------------------------------------------------------
 # HYBENSMEM_NMIN:
 # Minimum number of ensemble members required a hybrid GSI analysis 

diff --git a/ush/templates/FV3LAM_wflow.xml b/ush/templates/FV3LAM_wflow.xml
@@ -1130,7 +1130,7 @@ MODULES_RUN_TASK_FP script.
  <taskdep task="&RUN_ENKFUPDT_TN;"/>
  {%- endif %}
  </and>
- {%- else %}
+ {%- elif do_envar_radar_ref %}
  <or>
  <and>
  <taskdep task="&ANAL_GSI_TN;_#type#"/>
@@ -1147,7 +1147,10 @@ MODULES_RUN_TASK_FP script.
  </or> 
  </and>
  </or>
+ {%- else %}
+ <taskdep task="&ANAL_GSI_TN;_#type#"/>
  {%- endif %}
+ </and>
  </dependency>
 
  </task>
@@ -1234,6 +1237,7 @@ MODULES_RUN_TASK_FP script.
 ************************************************************************
 ************************************************************************
 -->
+{% if not is_rtma %}
  <metatask name="&SAVE_RESTART_TN;_prod">
 
  <var name="fhr"> {{ restart_hrs_prod }} </var>
@@ -1272,7 +1276,7 @@ MODULES_RUN_TASK_FP script.
  </task>
 
  </metatask>
-
+{%- endif %}
 <!--
 ************************************************************************
 ************************************************************************