Update model and settings for Prototype P8b (#681)

Updates model and setting defaults for prototype 8b.

The UFS version is updated to a recent version of UFS (tag Prototype-P8b). As part of this, the print_esmf option is removed from the model configure file and replaced with a new setting in nems.configure. We have not made this option user-configurable at this point, but that may come in the future.

The default CCPP suites have been updated for all modes to either FV3_GFS_v17_p8 or FV3_GFS_v17_coupled_p8. This involves changing the microphysics to Thompson and changes the gravity wave drag version (knob_ugwp_version) to 0. This change breaks cycled mode for the time being. A future PR will revert the defaults to ones that allow cycled to run along with a more robust system for changing settings.

Tiled fix files are now used for all modes instead of just coupled.

Cellular automata now defaults on ON.

There are also other miscellaneous settings that have had their defaults changed.

Fixes #641, #687

modulefiles/module_base.hera

@@ -27,18 +27,17 @@ module load png/1.6.35
 module load hdf5/1.10.6
 module load netcdf/4.7.4
 module load pio/2.5.2
-module load esmf/8_2_0_beta_snapshot_14
+module load esmf/8.2.1b04 
 module load fms/2021.03
 
 module load bacio/2.4.1
-module load g2/3.4.1
-module load g2tmpl/1.9.1
+module load g2/3.4.2
+module load g2tmpl/1.10.0
 module load ip/3.3.3
 module load nemsio/2.5.2
 module load sp/2.3.3
 module load w3emc/2.7.3
 module load w3nco/2.4.1
-module load upp/10.0.8
 
 module load wgrib2/2.0.8
 setenv WGRIB2 wgrib2
@@ -49,4 +48,3 @@ module load anaconda/2.3.0
 
 # waveprep
 module load cdo/1.9.5 
-

modulefiles/module_base.orion

@@ -26,18 +26,17 @@ module load png/1.6.35
 module load hdf5/1.10.6
 module load netcdf/4.7.4
 module load pio/2.5.2
-module load esmf/8_2_0_beta_snapshot_14
+module load esmf/8.2.1b04
 module load fms/2021.03
 
 module load bacio/2.4.1
-module load g2/3.4.1
-module load g2tmpl/1.9.1
+module load g2/3.4.2
+module load g2tmpl/1.10.0
 module load ip/3.3.3
 module load nemsio/2.5.2
 module load sp/2.3.3
 module load w3emc/2.7.3
 module load w3nco/2.4.1
-module load upp/10.0.8
 
 module load wgrib/2.0.8
 setenv WGRIB2 wgrib2

modulefiles/module_base.wcoss_dell_p3

@@ -34,18 +34,17 @@ module load png/1.6.35
 module load hdf5/1.10.6
 module load netcdf/4.7.4
 module load pio/2.5.2
-module load esmf/8_2_0_beta_snapshot_14
+module load esmf/8.2.1b04
 module load fms/2021.03
 
 module load bacio/2.4.1
-module load g2/3.4.1
-module load g2tmpl/1.9.1
+module load g2/3.4.2
+module load g2tmpl/1.10.0
 module load ip/3.3.3
 module load nemsio/2.5.2
 module load sp/2.3.3
 module load w3emc/2.7.3
 module load w3nco/2.4.1
-module load upp/10.0.8
 
 module load wgrib2/2.0.8
 setenv WGRIB2 wgrib2

parm/chem/CA2G_instance_CA.bc.rc

@@ -7,6 +7,9 @@ nbins:  2
 aerosol_radBands_optics_file: ExtData/optics/opticsBands_BC.v1_3.RRTMG.nc
 aerosol_monochromatic_optics_file: ExtData/monochromatic/optics_BC.v1_3.nc
 
+# Aircraft emission factor: convert input unit to kg C
+aircraft_fuel_emission_factor: 1.0000
+
 # Heights [m] of LTO, CDS and CRS aviation emissions layers
 aviation_vertical_layers: 0.0 100.0 9.0e3 10.0e3
 
@@ -35,3 +38,4 @@ sigma: 2.0  2.0
 
 pressure_lid_in_hPa: 0.01
 
+point_emissions_srcfilen: /dev/null

parm/chem/CA2G_instance_CA.br.rc

@@ -5,6 +5,9 @@
 aerosol_radBands_optics_file: ExtData/optics/opticsBands_BRC.v1_5.RRTMG.nc
 aerosol_monochromatic_optics_file: ExtData/monochromatic/optics_BRC.v1_5.nc
 
+# Aircraft emission factor: convert input unit to kg C
+aircraft_fuel_emission_factor: 1.0000
+
 # Heights [m] of LTO, CDS and CRS aviation emissions layers
 aviation_vertical_layers: 0.0 100.0 9.0e3 10.0e3
 
@@ -41,3 +44,5 @@ sigma: 2.20  2.20
 nbins: 2
 
 pressure_lid_in_hPa: 0.01
+
+point_emissions_srcfilen: /dev/null

parm/chem/CA2G_instance_CA.oc.rc

@@ -5,6 +5,9 @@
 aerosol_radBands_optics_file: ExtData/optics/opticsBands_OC.v1_3.RRTMG.nc
 aerosol_monochromatic_optics_file: ExtData/monochromatic/optics_OC.v1_3.nc
 
+# Aircraft emission factor: convert input unit to kg C
+aircraft_fuel_emission_factor: 1.0000
+
 # Heights [m] of LTO, CDS and CRS aviation emissions layers
 aviation_vertical_layers: 0.0 100.0 9.0e3 10.0e3
 
@@ -41,3 +44,5 @@ sigma: 2.20  2.20
 pressure_lid_in_hPa: 0.01
 
 nbins: 2 
+
+point_emissions_srcfilen: /dev/null

parm/chem/ExtData.other

@@ -72,6 +72,9 @@ OC_ANTEOC2 NA  N Y %y4-%m2-%d2t12:00:00 none none OC_elev  ExtData/nexus/CEDS/v2
 # 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
@@ -94,6 +97,9 @@ BC_ANTEBC2 NA  N Y %y4-%m2-%d2t12:00:00 none none BC_elev  ExtData/nexus/CEDS/v2
 # 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
@@ -118,6 +124,9 @@ 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

parm/config/config.base.emc.dyn

@@ -150,7 +150,7 @@ export DO_WAVE="NO"
 export DO_OCN="NO"
 export DO_ICE="NO"
 export DO_AERO=@DO_AERO@
-export CCPP_SUITE="FV3_GFS_v16"
+export CCPP_SUITE="FV3_GFS_v17_p8"
 export WAVE_CDUMP="" # When to include wave suite: gdas, gfs, or both
 export cplwav2atm=".false."
 
@@ -176,15 +176,15 @@ case "${APP}" in
     export DO_COUPLED="YES"
     export DO_OCN="YES"
     export DO_ICE="YES"
-    export CCPP_SUITE="FV3_GFS_v16_coupled_nsstNoahmpUGWPv1"
+    export CCPP_SUITE="FV3_GFS_v17_coupled_p8"
     export confignamevarfornems="cpld${aero_nems_string}"
     ;;
   S2SW)
     export DO_COUPLED="YES"
     export DO_WAVE="YES"
     export DO_OCN="YES"
     export DO_ICE="YES"
-    export CCPP_SUITE="FV3_GFS_v16_coupled_nsstNoahmpUGWPv1"
+    export CCPP_SUITE="FV3_GFS_v17_coupled_p8"
     export WAVE_CDUMP="both"
     export cplwav2atm=".true."
     export confignamevarfornems="cpld${aero_nems_string}_wave"
@@ -299,7 +299,7 @@ export gldas_cyc=00
 # fi
 
 # Microphysics Options: 99-ZhaoCarr, 8-Thompson; 6-WSM6, 10-MG, 11-GFDL
-export imp_physics=11
+export imp_physics=8
 
 # Shared parameters
 # Hybrid related

parm/config/config.base.nco.static

@@ -212,7 +212,7 @@ export DO_WAVE=YES
 export WAVE_CDUMP="both"
 
 # Microphysics Options: 99-ZhaoCarr, 8-Thompson; 6-WSM6, 10-MG, 11-GFDL
-export imp_physics=11
+export imp_physics=8
 
 # Shared parameters
 # Hybrid related

parm/config/config.defaults.s2sw

@@ -12,20 +12,17 @@ FHOUT_GFS=6
 FHOUT_HF_GFS=-1
 
 # config.fcst
-IALB=2
-IEMS=2
-DO_CA="YES"
-min_seaice="1e-6"
+min_seaice="1.0e-6"
 use_cice_alb=".true."
-export iopt_sfc="3"
 
 # config.fv3
 DELTIM=300
-layout_x_gfs=8
-layout_y_gfs=12
+layout_x_gfs=24
+layout_y_gfs=16
 WRITE_GROUP_GFS=1
 WRTTASK_PER_GROUP_GFS=80
 WRTIOBUF="32M"
+MEDPETS=300
 
 # config.wave
 waveGRD='gwes_30m'

parm/config/config.fcst

@@ -92,14 +92,59 @@ export h2o_phys=".true."
 # Options of stratosphere O3 physics reaction coefficients
 export new_o3forc="YES"
 
-# export launch_level=$(echo "$LEVS/2.35" |bc)
+export gwd_opt=2  
+
+# --GFS.v16 uGWD.v0, used for suite FV3_GFS_v16 and UFS p6 etc
+# do_ugwp=T: use unified CGWD and OGWD, and turbulent orographic form drag (TOFD)
+# do_ugwp=F: use unified CGWD but old OGWD, TOFD is not uded.
+if [ $gwd_opt -eq 1 ]; then
+    export knob_ugwp_version=0
+    export do_ugwp=".false."
+    export do_tofd=".false."
+    export launch_level=$(echo "$LEVS/2.35" |bc)
+fi
+
+
+# -- uGWD.v1, for suite FV3_GFS_v17 and FV3_GFS_v17p8b etc
+if [ $gwd_opt -eq 2 ]; then
+
+    #--used for UFS p7 and p8a
+    #export  knob_ugwp_version=1
+    #export  do_ugwp=".false."
+    #export  do_tofd=".false."
+    #export  do_ugwp_v0=".false."
+    #export  do_ugwp_v1=".true."
+    #export  do_ugwp_v0_orog_only=".false."
+    #export  do_ugwp_v0_nst_only=".false."
+    #export  do_gsl_drag_ls_bl=".true."
+    #export  do_gsl_drag_ss=".true."
+    #export  do_gsl_drag_tofd=".true."
+    #export  do_ugwp_v1_orog_only=".false."
+
+    #--used for UFS p8b 
+    export knob_ugwp_version=0
+    export do_ugwp=".false."
+    export do_tofd=".false."
+    export do_ugwp_v0=".true."
+    export do_ugwp_v1=".false."
+    export do_ugwp_v0_orog_only=".false."
+    export do_ugwp_v0_nst_only=".false."
+    export do_gsl_drag_ls_bl=".false."
+    export do_gsl_drag_ss=".true."
+    export do_gsl_drag_tofd=".true."
+    export do_ugwp_v1_orog_only=".false."
+    export launch_level=$(echo "$LEVS/2.35" |bc)
+fi
+
+
+
 
 # Sponge layer settings
 export tau=10.0
 export rf_cutoff=7.5e2
 export d2_bg_k1=0.20
 export d2_bg_k2=0.04
-export dz_min=2
+export dz_min=6
 export n_sponge=42
 if [ $LEVS = "128" -a "$CDUMP" = "gdas" ]; then
    export tau=5.0
@@ -124,10 +169,7 @@ export icliq_sw=2   ; #cloud optical coeffs from AER's newer version v3.9-v4.0 f
 export isubc_sw=2
 export isubc_lw=2
 
-export IALB=${IALB:-1}
-export IEMS=${IEMS:-1}
-export ISOL=2
-export ICO2=2
+export iopt_sfc="3"
 
 # Microphysics configuration
 export dnats=0
@@ -147,9 +189,23 @@ elif [ $imp_physics -eq 6 ]; then # WSM6
 
 elif [ $imp_physics -eq 8 ]; then # Thompson
     export ncld=2
-    export FIELD_TABLE="$HOMEgfs/parm/parm_fv3diag/field_table_thompson${tbf}"
+    export FIELD_TABLE="$HOMEgfs/parm/parm_fv3diag/field_table_thompson_noaero_tke"
     export nwat=6
 
+    export cal_pre=".false."
+    export random_clds=".false."
+    export effr_in=".true."
+    export ltaerosol=".false."
+    export lradar=".false."
+    export ttendlim="-999"
+    export dt_inner=$((DELTIM/2))
+
+    export hord_mt_nh_nonmono=5
+    export hord_xx_nh_nonmono=5
+    export vtdm4_nh_nonmono=0.02
+    export nord=2
+    export dddmp=0.1
+    export d4_bg=0.12
 elif [ $imp_physics -eq 11 ]; then # GFDL
     export ncld=5
     export FIELD_TABLE="$HOMEgfs/parm/parm_fv3diag/field_table_gfdl${tbf}"
@@ -179,7 +235,7 @@ export DO_SPPT=${DO_SPPT:-"NO"}
 export DO_SKEB=${DO_SKEB:-"NO"}
 export DO_SHUM=${DO_SHUM:-"NO"}
 export DO_LAND_PERT=${DO_LAND_PERT:-"NO"}
-export DO_CA=${DO_CA:-"NO"}
+export DO_CA=${DO_CA:-"YES"}
 export DO_OCN_SPPT=${DO_OCN_SPPT:-"NO"}
 export DO_OCN_PERT_EPBL=${DO_OCN_PERT_EPBL:-"NO"}
 
@@ -192,7 +248,8 @@ export min_lakeice="0.15"
 export min_seaice=${min_seaice:-"0.15"}
 export use_cice_alb=${use_cice_alb:-".false."}
 
-export FSICL="99999"
+export FSICL="0"
+export FSICS="0"
 
 #---------------------------------------------------------------------
 

parm/config/config.fv3

@@ -87,14 +87,14 @@ case $case_in in
         export WRTIOBUF="8M"
         ;;
     "C384")
-        export DELTIM=${DELTIM:-240}
+        export DELTIM=${DELTIM:-300}
         export layout_x=6
         export layout_y=8
-        export layout_x_gfs=${layout_x_gfs:-6}
-        export layout_y_gfs=${layout_y_gfs:-8}
+        export layout_x_gfs=${layout_x_gfs:-8}
+        export layout_y_gfs=${layout_y_gfs:-12}
         export nth_fv3=2
         export nth_fv3_gfs=${nth_fv3_gfs:-2}
-        export cdmbgwd="1.0,2.2,1.0,1.0"  # mountain blocking, ogwd, cgwd, cgwd src scaling
+        export cdmbgwd="1.1,0.72,1.0,1.0"  # mountain blocking, ogwd, cgwd, cgwd src scaling
         export WRITE_GROUP=1
         export WRTTASK_PER_GROUP=$npe_node_max
         export WRITE_GROUP_GFS=${WRITE_GROUP_GFS:-2}

parm/config/config.resources

@@ -198,8 +198,8 @@ elif [ $step = "fcst" ]; then
 
         (( ATMPETS = layout_x * layout_y * 6 ))
 
-        # Mediator only uses the atm model PETS
-        export MEDPETS=$ATMPETS
+        # Mediator only uses the atm model PETS or less
+        export MEDPETS=${MEDPETS:-ATMPETS}
 
         if [[ $DO_AERO == "YES" ]]; then
             # Aerosol model only uses the atm model PETS

parm/mom6/MOM_input_template_025

@@ -98,7 +98,7 @@ WRITE_GEOM = 2                  ! default = 1
                                 ! If =0, never write the geometry and vertical grid files. If =1, write the
                                 ! geometry and vertical grid files only for a new simulation. If =2, always
                                 ! write the geometry and vertical grid files. Other values are invalid.
-SAVE_INITIAL_CONDS = True      !   [Boolean] default = False
+SAVE_INITIAL_CONDS = False      !   [Boolean] default = False
                                 ! If true, write the initial conditions to a file given by IC_OUTPUT_FILE.
 
 ! === module MOM_hor_index ===
@@ -406,6 +406,9 @@ GILL_EQUATORIAL_LD = True       !   [Boolean] default = False
                                 ! radius, otherwise, if false, use Pedlosky's definition. These definitions
                                 ! differ by a factor of 2 in front of the beta term in the denominator. Gill's
                                 ! is the more appropriate definition.
+INTERNAL_WAVE_SPEED_BETTER_EST = False !   [Boolean] default = True 
+                                ! If true, use a more robust estimate of the first mode wave speed as the
+                                ! starting point for iterations.
 
 ! === module MOM_set_visc ===
 CHANNEL_DRAG = True             !   [Boolean] default = False
@@ -728,6 +731,11 @@ NSTAR = 0.06                    !   [nondim] default = 0.2
                                 ! The portion of the buoyant potential energy imparted by surface fluxes that is
                                 ! available to drive entrainment at the base of mixed layer when that energy is
                                 ! positive.
+EPBL_MLD_BISECTION = True       !   [Boolean] default = False 
+                                ! If true, use bisection with the iterative determination of the self-consistent
+                                ! mixed layer depth.  Otherwise use the false position after a maximum and
+                                ! minimum bound have been evaluated and the returned value or bisection before
+                                ! this.
 MSTAR_CONV_ADJ = 0.667          !   [nondim] default = 0.0
                                 ! Coefficient used for reducing mstar during convection due to reduction of
                                 ! stable density gradient.