Description:
The purpose here is to establish a new baseline on Derecho using ctsm5.1.dev160 and including our ABCDE_blk_A5BCD changes. This will serve as a baseline with the most recent tag to compare against simulations where we plan to back off some of the ABCDE_blk_A5BCD changes. Also will serve as a baseline to compare against a similar set of simulations using CRUJRA forcing.
Note that we also turn on Meier here since it's inclusion on main is imminent.
We are also switching to 1deg resolution here for ease/speed of post-processing.
Case directory:
Locally:
/glade/work/oleson/ctsm5.1.dev160/cime/scripts/ctsm51_ctsm51d160_1deg_GSWP3V1_ABCDE_blk_A5BCD_1850AD
!----------------------------------------------------------------------------------
! Settings from cmip6_carbon_isotopes
! CAN ONLY USE THIS IF STARTING FROM AN INITIAL FILE THAT HAS ISOTOPES
!----------------------------------------------------------------------------------
src.clm/atm2lndType.F90: Add some additional history fields
src.clm/clm_varpar.F90: nlevurb=5 -> nlevurb=10
src.clm/surfrdMod.F90: Remove read of PFTDATA_MASK
src.clm/Wateratm2lndType.F90: Add some additional history fields
src.clm/CNPhenologyMod.F90 decreasing temperature thresholds to -3C, reducing the minimum critical daylength by 25%, doubling the snow threshold
src.clm/controlMod.F90 remove error associated with snicar_snobc_intmix
src.clm/SnowSnicarMod.F90 set flg_snoage_scl to .true. so that xdrdt change takes effect
Diagnostics:
Output:
Output (if still available):
/glade/derecho/scratch/oleson/archive/ctsm51_ctsm51d160_1deg_GSWP3V1_ABCDE_blk_A5BCD_1850AD/
Contacts:
@olyson
@wwieder
Extra details:
I've designated this combination simulation as ABCDE_blk_A5BCD, where the first ABCDE is:
A - SturmSnowtk (implemented via namelist snow_thermal_cond_method = 'Sturm1997')
B - decidphenotrig (implemented via src.clm/CNPhenologyMod.F90)
C - frootleaf (implemented via ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc)
D - FUNfracfixers (implemented via ctsm51_params.c2311117.norevert.nc.ABCD_blk_A5BCD.nc)
E - new SNICAR (on by default)
blk - snicar_snobc_intmix = .true. (implemented via namelist, per Issue #36 )
and the next ABCD is (per Issue #36 )
A5 - xdrdt = 5 (implemented via ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc and takes effect by code change to SnowSnicarMod.F90)
B - scvng_fct_mlt_sf = 0.5 (implemented via ctsm51_params.c2311117.norevert.nc.ABCD_blk_A5BCD.nc)
C - snw_rds_refrz = 1500 (implemented via ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc)
D - fresh_snw_rds_max = 400._r8 (implemented via ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc)
Description: The purpose here is to establish a new baseline on Derecho using ctsm5.1.dev160 and including our ABCDE_blk_A5BCD changes. This will serve as a baseline with the most recent tag to compare against simulations where we plan to back off some of the ABCDE_blk_A5BCD changes. Also will serve as a baseline to compare against a similar set of simulations using CRUJRA forcing. Note that we also turn on Meier here since it's inclusion on main is imminent. We are also switching to 1deg resolution here for ease/speed of post-processing.
Case directory: Locally: /glade/work/oleson/ctsm5.1.dev160/cime/scripts/ctsm51_ctsm51d160_1deg_GSWP3V1_ABCDE_blk_A5BCD_1850AD
Sandbox: Locally: /glade/work/oleson/ctsm5.1.dev160
usernl changes:
user_nl_clm:
! Full output hist_nhtfrq = 0 hist_mfilt = 1 hist_empty_htapes = .false. hist_fincl1 = 'FSDSND','FSDSVD','FSDSNI','FSDSVI','RAIN_FROM_ATM','SNOW_FROM_ATM','FLDS_NOT_DOWNSCALED','ZBOT','Tair_from_atm','Thair_from_atm','QBOT_NOT_DOWNSCALED','PBOT_NOT_DOWNSCALED','Rho_from_atm','PCO2','UWIND','VWIND','BCPHIDRY','BCPHODRY','BCPHIWET','OCPHIDRY','OCPHODRY','OCPHIWET','DSTWET1','DSTDRY1','DSTWET2','DSTDRY2','DSTWET3','DSTDRY3','DSTWET4','DSTDRY4','ATM_TOPO'
!---------------------------------------------------------------------------------- ! Settings from cmip6_carbon_isotopes ! CAN ONLY USE THIS IF STARTING FROM AN INITIAL FILE THAT HAS ISOTOPES !----------------------------------------------------------------------------------
use_c13 = .true. use_c14 = .true. use_c13_timeseries = .true. use_c14_bombspike = .true. fsurdat = '/glade/campaign/cesm/cesmdata/cseg/inputdata/lnd/clm2/surfdata_esmf/ctsm5.2.0/surfdata_0.9x1.25_hist_78pfts_CMIP6_1850_c230517.nc' paramfile = '/glade/campaign/cgd/tss/people/oleson/modify_param/ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc' snicar_snobc_intmix = .true. snow_thermal_cond_method = 'Sturm1997' z0param_method = 'Meier2022'
SourceMods:
src.clm/atm2lndType.F90: Add some additional history fields src.clm/clm_varpar.F90: nlevurb=5 -> nlevurb=10 src.clm/surfrdMod.F90: Remove read of PFTDATA_MASK src.clm/Wateratm2lndType.F90: Add some additional history fields src.clm/CNPhenologyMod.F90 decreasing temperature thresholds to -3C, reducing the minimum critical daylength by 25%, doubling the snow threshold src.clm/controlMod.F90 remove error associated with snicar_snobc_intmix src.clm/SnowSnicarMod.F90 set flg_snoage_scl to .true. so that xdrdt change takes effect
Diagnostics:
Output: Output (if still available): /glade/derecho/scratch/oleson/archive/ctsm51_ctsm51d160_1deg_GSWP3V1_ABCDE_blk_A5BCD_1850AD/
Contacts: @olyson @wwieder
Extra details: I've designated this combination simulation as ABCDE_blk_A5BCD, where the first ABCDE is: A - SturmSnowtk (implemented via namelist snow_thermal_cond_method = 'Sturm1997') B - decidphenotrig (implemented via src.clm/CNPhenologyMod.F90) C - frootleaf (implemented via ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc) D - FUNfracfixers (implemented via ctsm51_params.c2311117.norevert.nc.ABCD_blk_A5BCD.nc) E - new SNICAR (on by default) blk - snicar_snobc_intmix = .true. (implemented via namelist, per Issue #36 ) and the next ABCD is (per Issue #36 ) A5 - xdrdt = 5 (implemented via ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc and takes effect by code change to SnowSnicarMod.F90) B - scvng_fct_mlt_sf = 0.5 (implemented via ctsm51_params.c2311117.norevert.nc.ABCD_blk_A5BCD.nc) C - snw_rds_refrz = 1500 (implemented via ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc) D - fresh_snw_rds_max = 400._r8 (implemented via ctsm51_params.c231117.norevert.nc.ABCD_blk_A5BCD.nc)
See discussion in https://github.com/NCAR/LMWG_dev/discussions/3 AD spunup to 2.81% disequilibrium, 320 years pAD spunup to X.XX% disequilibrium, X years