Enable EP4a experiments (DO NOT MERGE)

jobs/JGDAS_ENKF_FCST

@@ -8,6 +8,7 @@ source "${HOMEgfs}/ush/jjob_header.sh" -e "efcs" -c "base fcst efcs"
 # Set variables used in the script
 ##############################################
 export CDUMP=${RUN/enkf}
+[[ ${DO_WAVE} == "YES" ]] && declare -rx RUNwave="${RUN}wave"
 
 ##############################################
 # Begin JOB SPECIFIC work
@@ -30,7 +31,6 @@ fi
 export ENSEND=$((NMEM_EFCSGRP * 10#${ENSGRP}))
 export ENSBEG=$((ENSEND - NMEM_EFCSGRP + 1))
 
-
 ###############################################################
 # Run relevant script
 

parm/config/gefs/config.base.emc.dyn

@@ -128,6 +128,9 @@ export SENDDBN_NTC=${SENDDBN_NTC:-"NO"}
 export SENDDBN=${SENDDBN:-"NO"}
 export DBNROOT=${DBNROOT:-${UTILROOT}/fakedbn}
 
+# PROTO settings
+export PROTO=@PROTO@
+
 # APP settings
 export APP=@APP@
 
@@ -151,6 +154,19 @@ export OPS_RES="C768" # Do not change
 export LEVS=128
 export CASE="@CASECTL@"
 export CASE_ENS="@CASEENS@"
+#
+# Setup prototype cases
+case "${PROTO}" in
+  EP4A)
+    APP="S2SWA"
+    CASE="C384"
+    CASE_ENS="C384"
+    ;;
+  *)
+    ;;
+esac
+
+
 # TODO: This should not depend on $CASE or $CASE_ENS
 # These are the currently available grid-combinations
 case "${CASE}" in
@@ -322,4 +338,17 @@ export ARCH_FCSTICFREQ=1     # Archive frequency in days for gdas and gfs foreca
 
 export DELETE_COM_IN_ARCHIVE_JOB="YES"   # NO=retain ROTDIR.  YES default in arch.sh and earc.sh.
 
+# Setup prototype cases
+case "${PROTO}" in
+  EP4A)
+    APP="S2SWA"
+    CASE="C384"
+    CASE_ENS="C384"
+    OCNRES=025
+    waveGRD='glo_025'
+    ;;
+  *)
+    ;;
+esac
+
 echo "END: config.base"

parm/config/gefs/config.efcs

@@ -31,8 +31,14 @@ source "${EXPDIR}/config.ufs" ${string}
 source "${EXPDIR}/config.resources" efcs
 
 # Use serial I/O for ensemble (lustre?)
-export OUTPUT_FILETYPE_ATM="netcdf"
-export OUTPUT_FILETYPE_SFC="netcdf"
+case "${PROTO:-}" in
+  EP4A)
+  ;;
+    *)
+    export OUTPUT_FILETYPE_ATM="netcdf"
+    export OUTPUT_FILETYPE_SFC="netcdf"
+  ;;
+esac
 
 # Number of enkf members per fcst job
 export NMEM_EFCSGRP=1
@@ -42,23 +48,46 @@ export RERUN_EFCSGRP="NO"
 export WRITE_DOPOST=".true."
 
 # Stochastic physics parameters (only for ensemble forecasts)
-export DO_SKEB="YES"
-export SKEB=0.3
-export SKEB_TAU=21600.
-export SKEB_LSCALE=250000.
-export SKEBNORM=0
-export SKEB_NPASS=30
-export SKEB_VDOF=5
-export DO_SHUM="YES"
-export SHUM=0.005
-export SHUM_TAU=21600.
-export SHUM_LSCALE=500000.
-export DO_SPPT="YES"
-export SPPT=0.5
-export SPPT_TAU=21600.
-export SPPT_LSCALE=500000.
-export SPPT_LOGIT=".true."
-export SPPT_SFCLIMIT=".true."
+case "${PROTO:-}" in
+  EP4A)
+    export MOM6_RESTART_SETTING=n
+    export IEMS_ENKF=2
+    export FHCYC_ENKF=24
+    export DO_SKEB="YES"
+    export SKEB=0.8
+    export SKEB_TAU=21600.
+    export SKEB_LSCALE=500000.
+    export SKEBNORM=1
+    export SKEB_NPASS=30
+    export SKEB_VDOF=5
+    export DO_SHUM="NO"
+    export DO_SPPT="YES"
+    export SPPT=0.56
+    export SPPT_TAU=21600.
+    export SPPT_LSCALE=500000.
+    export SPPT_LOGIT=".true."
+    export SPPT_SFCLIMIT=".true."
+  ;;
+    *)
+    export DO_SKEB="YES"
+    export SKEB=0.3
+    export SKEB_TAU=21600.
+    export SKEB_LSCALE=250000.
+    export SKEBNORM=0
+    export SKEB_NPASS=30
+    export SKEB_VDOF=5
+    export DO_SHUM="YES"
+    export SHUM=0.005
+    export SHUM_TAU=21600.
+    export SHUM_LSCALE=500000.
+    export DO_SPPT="YES"
+    export SPPT=0.5
+    export SPPT_TAU=21600.
+    export SPPT_LSCALE=500000.
+    export SPPT_LOGIT=".true."
+    export SPPT_SFCLIMIT=".true."
+  ;;
+esac
 
 export restart_interval=${restart_interval:-6}
 

parm/config/gefs/config.fcst

@@ -298,6 +298,30 @@ else
   export io_layout="1,1"
 fi
 
+# Use prototype settings
+case "${PROTO}" in
+   EP4A)
+      export tau=0.0
+      export rf_cutoff=10.
+      export d2_bg_k1=0.15
+      export d2_bg_k2=0
+      export IAER=2011
+      export io_layout="16,16"
+      export MOM6_RESTART_SETTING=n
+      export ODA_INCUPD_NHOURS=0
+      AERO_CONFIG_DIR=${HOMEgfs}/parm/ufs/ep4a
+      AERO_DIAG_TABLE=${AERO_CONFIG_DIR}/diag_table.aero
+      AERO_FIELD_TABLE=${AERO_CONFIG_DIR}/field_table.aero
+      aero_conv_scav_factors="'*:0.3','so2:0.0','msa:0.0','dms:0.0','bc1:0.6','bc2:0.6','oc1:0.4','oc2:0.4','dust1:0.6','dust2:0.6', 'dust3:0.6','dust4:0.6','dust5:0.6','seas1:0.5','seas2:0.5','seas3:0.5','seas4:0.5','seas5:0.5'"
+      rst_wav=$(( restart_interval_gfs > 0 ? restart_interval_gfs : FHMAX_WAV ))
+      export DT_1_RST_WAV=0
+      export DT_2_RST_WAV=$(( 3600 * rst_wav ))
+      export RSTIOFF_WAV=0
+      ;;
+   *)
+      ;;
+esac
+
 if [[ "${DO_AERO}" = "YES" ]]; then # temporary settings for aerosol coupling
     export AERO_DIAG_TABLE="${AERO_DIAG_TABLE:-${HOMEgfs}/parm/ufs/fv3/diag_table.aero}"
     export AERO_FIELD_TABLE="${AERO_FIELD_TABLE:-${HOMEgfs}/parm/ufs/fv3/field_table.aero}"
@@ -309,6 +333,7 @@ if [[ "${DO_AERO}" = "YES" ]]; then # temporary settings for aerosol coupling
     export dnats_aero="${aero_diag_tracers:-0}"
 fi
 
+
 # Remember config.efcs will over-ride these values for ensemble forecasts
 # if these variables are re-defined there.
 # Otherwise, the ensemble forecast will inherit from config.fcst

parm/config/gefs/config.ufs

@@ -381,4 +381,33 @@ if [[ "${skip_ww3}" == "false" ]]; then
   export ntasks_ww3 nthreads_ww3
 fi
 
+case "${PROTO}" in
+  EP4A)
+    export layout_x=16
+    export layout_y=16
+    export nthreads_fv3=1
+    export nthreads_fv3_gfs=1
+    (( ntasks_fv3 = layout_x * layout_y * 6 ))
+    (( ntasks_fv3_gfs = layout_x_gfs * layout_y_gfs * 6 ))
+    export ntasks_fv3
+    export ntasks_fv3_gfs
+    export WRITE_GROUP=1
+    export WRITE_GROUP_GFS=1
+    (( ntasks_quilt = WRITE_GROUP * WRTTASK_PER_GROUP_PER_THREAD ))
+    (( ntasks_quilt_gfs = WRITE_GROUP_GFS * WRTTASK_PER_GROUP_PER_THREAD_GFS ))
+    export ntasks_quilt
+    export ntasks_quilt_gfs
+    export ntasks_cice6=72
+    export nthreads_cice6=1
+    export ntasks_mom6=130
+    export nthreads_mom6=1
+    export ntasks_ww3=262
+    export nthreads_ww3=1
+    export OUTPUT_FILETYPE_ATM="netcdf_parallel"
+    export OUTPUT_FILETYPE_SFC="netcdf_parallel"
+    ;;
+     *)
+    ;;
+esac
+
 echo "END: config.ufs"

parm/ufs/ep4a/AERO.rc

@@ -0,0 +1,10 @@
+NX: 4
+NY: 24
+
+# Atmospheric Model Configuration Parameters
+# ------------------------------------------
+IOSERVER_NODES: 0
+
+DYCORE: NONE
+
+NUM_BANDS: 30

parm/ufs/ep4a/AERO_ExtData.rc

@@ -0,0 +1,160 @@
+PrimaryExports%%
+# -------------|-------|-------|--------|----------------------|--------|--------|-------------|----------|
+#  Import      |       |       | Regrid |        Refresh       | OffSet | Scale  | Variable On |   File   |
+#  Name        | Units | Clim  | Method |     Time Template    | Factor | Factor |     File    | Template |
+# -------------|-------|-------|--------|----------------------|--------|--------|-------------|----------|
+
+#====== Atmospheric Parameters =======================================
+TROPP            'Pa'    Y       N                  -            0.0      1.0      TROPP         /dev/null:10000.
+
+#====== Dust Imports =================================================
+# Ginoux input files
+DU_SRC            NA  N Y - none none du_src     ExtData/Dust/gocart.dust_source.v5a.x1152_y721.nc
+
+# FENGSHA input files. Note: regridding should be N or E - Use files with _FillValue != NaN 
+DU_CLAY           '1'  Y E           -           none none clayfrac    ExtData/Dust/FENGSHA_p81_10km_inputs.nc
+DU_SAND           '1'  Y E           -           none none sandfrac    ExtData/Dust/FENGSHA_p81_10km_inputs.nc
+DU_SILT           '1'  Y E           -           none none siltfrac    /dev/null
+DU_SSM            '1'  Y E           -           none none ssm         /dev/null:1.0
+DU_RDRAG          '1'  Y E %y4-%m2-%d2t12:00:00  none none albedo_drag ExtData/Dust/FENGSHA_p81_10km_inputs.nc
+DU_UTHRES         '1'  Y E           -           none none uthres      ExtData/Dust/FENGSHA_p81_10km_inputs.nc
+
+#====== Sulfate Sources =================================================
+# Anthropogenic (BF & FF) emissions -- allowed to input as two layers
+SU_ANTHROL1 NA  N Y %y4-%m2-%d2t12:00:00 none none SO2      ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+SU_ANTHROL2 NA  N Y %y4-%m2-%d2t12:00:00 none none SO2_elev ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+
+# Ship emissions
+SU_SHIPSO2 NA  N Y %y4-%m2-%d2t12:00:00 none none SO2_ship  ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+SU_SHIPSO4 NA  N Y %y4-%m2-%d2t12:00:00 none none SO4_ship  ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+
+# Aircraft fuel consumption
+SU_AIRCRAFT NA  Y Y %y4-%m2-%d2t12:00:00 none none none /dev/null
+
+# DMS concentration
+SU_DMSO NA  Y Y %y4-%m2-%d2t12:00:00 none none conc         ExtData/MERRA2/sfc/DMSclim_sfcconcentration.x360_y181_t12.Lana2011.nc4
+
+# Aviation emissions during the three phases of flight
+SU_AVIATION_LTO NA  Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_so2.aviation_lto.x3600_y1800_t12.2010.nc4
+SU_AVIATION_CDS NA  Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_so2.aviation_cds.x3600_y1800_t12.2010.nc4
+SU_AVIATION_CRS NA  Y Y %y4-%m2-%d2t12:00:00 none none so2_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_so2.aviation_crs.x3600_y1800_t12.2010.nc4
+
+# H2O2, OH and NO3 mixing ratios
+# --------------------------------------------------------------
+# If using 64 levels please replace this section with the correct values (ie replace 127 with 64)
+
+SU_H2O2 NA  N Y %y4-%m2-%d2t12:00:00 none none h2o2 ExtData/PIESA/L127/A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.%y4.nc
+SU_OH   NA  N Y %y4-%m2-%d2t12:00:00 none none oh   ExtData/PIESA/L127/A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.%y4.nc
+SU_NO3  NA  N Y %y4-%m2-%d2t12:00:00 none none no3  ExtData/PIESA/L127/A2_ACCMIP_gmic_MERRA_oh_h2o2_no3.x144_y91_z127_t14.%y4.nc
+#---------------------------------------------------------------
+
+# Production of SO2 from OCS oxidation
+pSO2_OCS NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null
+
+#SU_regionMask   NA  N v - none none REGION_MASK /scratch1/NCEPDEV/nems/Raffaele.Montuoro/data/NASA/ExtData/PIESA/sfc/ARCTAS.region_mask.x540_y361.2008.nc
+
+#=========== Carbonaceous aerosol sources ===========================================
+# ORGANIC CARBON
+# ---------------
+
+# # VOCs - OFFLINE MEGAN BIOG
+OC_ISOPRENE NA  N Y %y4-%m2-%d2t12:00:00 none none isoprene ExtData/nexus/MEGAN_OFFLINE_BVOC/v2019-10/%y4/MEGAN.OFFLINE.BIOVOC.%y4.emis.%y4%m2%d2.nc
+OC_LIMO     NA  N Y %y4-%m2-%d2t12:00:00 none none limo     ExtData/nexus/MEGAN_OFFLINE_BVOC/v2019-10/%y4/MEGAN.OFFLINE.BIOVOC.%y4.emis.%y4%m2%d2.nc
+OC_MTPA     NA  N Y %y4-%m2-%d2t12:00:00 none none mtpa     ExtData/nexus/MEGAN_OFFLINE_BVOC/v2019-10/%y4/MEGAN.OFFLINE.BIOVOC.%y4.emis.%y4%m2%d2.nc
+OC_MTPO     NA  N Y %y4-%m2-%d2t12:00:00 none none mtpo     ExtData/nexus/MEGAN_OFFLINE_BVOC/v2019-10/%y4/MEGAN.OFFLINE.BIOVOC.%y4.emis.%y4%m2%d2.nc
+
+# Biofuel Source -- Included in AeroCom anthropogenic emissions
+OC_BIOFUEL NA  Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null
+
+# Anthropogenic (BF & FF) emissions -- allowed to input as two layers
+OC_ANTEOC1 NA  N Y %y4-%m2-%d2t12:00:00 none none OC       ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+OC_ANTEOC2 NA  N Y %y4-%m2-%d2t12:00:00 none none OC_elev  ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+
+# EDGAR based ship emissions
+OC_SHIP    NA  N Y %y4-%m2-%d2t12:00:00 none none OC_ship  ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+
+# Aircraft fuel consumption
+OC_AIRCRAFT NA  N Y %y4-%m2-%d2t12:00:00 none none oc_aviation /dev/null
+
+# Aviation emissions during the three phases of flight
+OC_AVIATION_LTO NA  Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_oc.aviation_lto.x3600_y1800_t12.2010.nc4
+OC_AVIATION_CDS NA  Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_oc.aviation_cds.x3600_y1800_t12.2010.nc4
+OC_AVIATION_CRS NA  Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_oc.aviation_crs.x3600_y1800_t12.2010.nc4
+
+# SOA production
+pSOA_ANTHRO_VOC NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null
+
+#============================================================================================================
+# BLACK CARBON
+# ------------
+
+# Biofuel Source -- Included in AeroCom anthropogenic emissions
+BC_BIOFUEL NA  Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null
+
+# Anthropogenic (BF & FF) emissions -- allowed to input as two layers
+BC_ANTEBC1 NA  N Y %y4-%m2-%d2t12:00:00 none none BC       ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+BC_ANTEBC2 NA  N Y %y4-%m2-%d2t12:00:00 none none BC_elev  ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+
+# EDGAR based ship emissions
+BC_SHIP NA  N Y %y4-%m2-%d2t12:00:00 none none BC_ship     ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+
+# Aircraft fuel consumption
+BC_AIRCRAFT NA  N Y %y4-%m2-%d2t12:00:00 none none bc_aviation /dev/null
+
+# Aviation emissions during the LTO, SDC and CRS phases of flight
+BC_AVIATION_LTO NA  Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_bc.aviation_lto.x3600_y1800_t12.2010.nc4
+BC_AVIATION_CDS NA  Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_bc.aviation_cds.x3600_y1800_t12.2010.nc4
+BC_AVIATION_CRS NA  Y Y %y4-%m2-%d2t12:00:00 none none bc_aviation ExtData/PIESA/sfc/HTAP/v2.2/htap-v2.2.emis_bc.aviation_crs.x3600_y1800_t12.2010.nc4
+
+#============================================================================================================
+# BROWN CARBON
+# ------------
+# Biomass burning -- QFED-v2.x
+BRC_BIOMASS NA  N Y %y4-%m2-%d2t12:00:00 none none biomass /dev/null
+
+# Terpene emission
+BRC_TERPENE NA  Y Y %y4-%m2-%d2t12:00:00 none none terpene /dev/null
+
+# Biofuel Source -- Included in AeroCom anthropogenic emissions
+BRC_BIOFUEL NA  Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null
+
+# Anthropogenic (BF & FF) emissions -- allowed to input as two layers
+BRC_ANTEBRC1 NA  Y Y %y4-%m2-%d2t12:00:00 none none anteoc1 /dev/null
+BRC_ANTEBRC2 NA  Y Y %y4-%m2-%d2t12:00:00 none none anteoc2 /dev/null
+
+# EDGAR based ship emissions
+BRC_SHIP    NA  Y Y %y4-%m2-%d2t12:00:00 none none oc_ship /dev/null
+
+# Aircraft fuel consumption
+BRC_AIRCRAFT NA  N Y %y4-%m2-%d2t12:00:00 none none none /dev/null
+
+# Aviation emissions during the three phases of flight
+BRC_AVIATION_LTO NA  Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation /dev/null
+BRC_AVIATION_CDS NA  Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation /dev/null
+BRC_AVIATION_CRS NA  Y Y %y4-%m2-%d2t12:00:00 none none oc_aviation /dev/null
+
+# SOA production
+pSOA_BIOB_VOC NA Y Y %y4-%m2-%d2t12:00:00 none none biofuel /dev/null
+
+# ======= Nitrate Sources ========
+EMI_NH3_AG   'kg m-2 s-1'  N Y      %y4-%m2-%d2T12:00:00    none     none     NH3_ag      ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+EMI_NH3_EN   'kg m-2 s-1'  N Y      %y4-%m2-%d2T12:00:00    none     none     NH3_en      /dev/null
+EMI_NH3_IN   'kg m-2 s-1'  N Y      %y4-%m2-%d2T12:00:00    none     none     NH3_in      ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+EMI_NH3_RE   'kg m-2 s-1'  N Y      %y4-%m2-%d2T12:00:00    none     none     NH3_re      ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+EMI_NH3_TR   'kg m-2 s-1'  N Y      %y4-%m2-%d2T12:00:00    none     none     NH3_tr      ExtData/nexus/CEDS/v2019/%y4/CEDS.2019.emis.%y4%m2%d2.nc
+EMI_NH3_OC   'kg m-2 s-1'  Y Y      %y4-%m2-%d2T12:00:00    none     none     emiss_ocn   ExtData/PIESA/sfc/GEIA.emis_NH3.ocean.x576_y361.t12.20080715_12z.nc4
+
+# --------------------------------------------------------------
+# If using 64 levels please replace this section with the correct values (ie replace 127 with 64)
+NITRATE_HNO3 'mol mol-1'   Y N      %y4-%m2-%d2T12:00:00    none     0.20     hno3        ExtData/PIESA/L127/GMI.vmr_HNO3.x144_y91.t12.2006.nc4
+# --------------------------------------------------------------
+NI_regionMask NA           Y V      -                       none     none     REGION_MASK ExtData/PIESA/sfc/ARCTAS.region_mask.x540_y361.2008.nc
+#====== BIOMASS BURNING EMISSIONS =======================================
+
+# QFED
+#--------------------------------------------------------------------------------------------------------------------------------
+SU_BIOMASS NA  N Y %y4-%m2-%d2t12:00:00 none none biomass ExtData/nexus/QFED/%y4/%m2/qfed2.emis_so2.006.%y4%m2%d2.nc4
+OC_BIOMASS NA  N Y %y4-%m2-%d2t12:00:00 none none biomass ExtData/nexus/QFED/%y4/%m2/qfed2.emis_oc.006.%y4%m2%d2.nc4
+BC_BIOMASS NA  N Y %y4-%m2-%d2t12:00:00 none none biomass ExtData/nexus/QFED/%y4/%m2/qfed2.emis_bc.006.%y4%m2%d2.nc4
+EMI_NH3_BB NA  N Y %y4-%m2-%d2t12:00:00 none none biomass ExtData/nexus/QFED/%y4/%m2/qfed2.emis_nh3.006.%y4%m2%d2.nc4
+%%