Description:
Same as ctsm51_ctsm51d160_1deg_GSWP3V1_ABCDE_blk_A5BCD_1850AD but revert froot_leaf (C) and FUN_fracfixers (D) 50% of the way back to defaults.
Case directory:
Locally:
/glade/work/oleson/ctsm5.1.dev160/cime/scripts/ctsm51_ctsm51d160_1deg_GSWP3V1_ABpt5Cpt5DE_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
Output:
Output (if still available):
/glade/derecho/scratch/oleson/archive/ctsm51_ctsm51d160_1deg_GSWP3V1_ABpt5Cpt5DE_blk_A5BCD_1850AD/
Contacts:
@olyson
@wwieder
Extra details:
I've designated this combination simulation as ABpt5Cpt5DE_blk_A5BCD, where the first ABpt5Cpt5DE is:
A - SturmSnowtk (implemented via namelist snow_thermal_cond_method = 'Sturm1997')
B - decidphenotrig (implemented via src.clm/CNPhenologyMod.F90)
pt5C - back off frootleaf(11:12) to 50% of difference between default (1.5) and C (1.2) - so, 1.35 (implemented via ctsm51_params.c231117.norevert.nc.ABpt5Cpt5D_blk_A5BCD.nc)
pt5D - back off FUNfracfixers(11:12) to 50% of difference between default (0.25) and D (1.0) - so, 0.625 (implemented via ctsm51_params.c2311117.norevert.nc.ABpt5Cpt5D_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.ABpt5Cpt5D_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.ABpt5Cpt5D_blk_A5BCD.nc)
C - snw_rds_refrz = 1500 (implemented via ctsm51_params.c231117.norevert.nc.ABpt5Cpt5D_blk_A5BCD.nc)
D - fresh_snw_rds_max = 400._r8 (implemented via ctsm51_params.c231117.norevert.nc.ABpt5Cpt5D_blk_A5BCD.nc)
Description: Same as ctsm51_ctsm51d160_1deg_GSWP3V1_ABCDE_blk_A5BCD_1850AD but revert froot_leaf (C) and FUN_fracfixers (D) 50% of the way back to defaults.
Case directory: Locally: /glade/work/oleson/ctsm5.1.dev160/cime/scripts/ctsm51_ctsm51d160_1deg_GSWP3V1_ABpt5Cpt5DE_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.ABpt5Cpt5D_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:
https://webext.cgd.ucar.edu/I1850/ctsm51_ctsm51d160_1deg_GSWP3V1_ABpt5Cpt5DE_blk_A5BCD_1850AD/lnd/ctsm51_ctsm51d160_1deg_GSWP3V1_ABpt5Cpt5DE_blk_A5BCD_1850AD_311_320-ctsm51_ctsm51d160_1deg_GSWP3V1_ABCDE_blk_A5BCD_1850AD_311_320/setsIndex.html
Output: Output (if still available): /glade/derecho/scratch/oleson/archive/ctsm51_ctsm51d160_1deg_GSWP3V1_ABpt5Cpt5DE_blk_A5BCD_1850AD/
Contacts: @olyson @wwieder
Extra details: I've designated this combination simulation as ABpt5Cpt5DE_blk_A5BCD, where the first ABpt5Cpt5DE is: A - SturmSnowtk (implemented via namelist snow_thermal_cond_method = 'Sturm1997') B - decidphenotrig (implemented via src.clm/CNPhenologyMod.F90) pt5C - back off frootleaf(11:12) to 50% of difference between default (1.5) and C (1.2) - so, 1.35 (implemented via ctsm51_params.c231117.norevert.nc.ABpt5Cpt5D_blk_A5BCD.nc) pt5D - back off FUNfracfixers(11:12) to 50% of difference between default (0.25) and D (1.0) - so, 0.625 (implemented via ctsm51_params.c2311117.norevert.nc.ABpt5Cpt5D_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.ABpt5Cpt5D_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.ABpt5Cpt5D_blk_A5BCD.nc) C - snw_rds_refrz = 1500 (implemented via ctsm51_params.c231117.norevert.nc.ABpt5Cpt5D_blk_A5BCD.nc) D - fresh_snw_rds_max = 400._r8 (implemented via ctsm51_params.c231117.norevert.nc.ABpt5Cpt5D_blk_A5BCD.nc)
See discussion in https://github.com/NCAR/LMWG_dev/discussions/3 AD spunup to 2.79% disequilibrium, 320 years pAD spunup to X.XX% disequilibrium, X years