Description:
Same as ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCD_1850AD (issue #37 ) but adding on snw_rds_min=100
Motivation is to try to lower snow albedo further.
Just running an AD for now.
Note that this uses ctsm5.1.dev150.
NOTE: xdrdt=2 setting was not actually implemented due to a coding issue.
Case directory:
Locally:
/glade/work/oleson/ctsm5.1.dev150/cime/scripts/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCDE_1850AD
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/SoilTemperatureMod.F90 Sturm Snow tk
src.clm/CNPhenologyMod.F90 decreasing temperature thresholds to -3C, reducing the minimum critical daylength by 25%, doubling the snow threshold (I think to 20 cm)
src.clm/controlMod.F90 remove error associated with snicar_snobc_intmix
src.clm/AerosolMod.F90 change snw_rds_min from default of 54.526 to 100.
Output:
Output (if still available):
/glade/scratch/oleson/archive/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCDE_1850AD/
Contacts:
@olyson
@wwieder
Extra details:
My case name lengths are out of control, so I've designated this combination simulation as ABCDE_blk_ABCDE, where the first ABCDE is:
A - SturmSnowtk (implemented via src.clm/SoilTemperatureMod.F90)
B - decidphenotrig (implemented via src.clm/CNPhenologyMod.F90)
C - frootleaf (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc)
D - FUNfracfixers (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc)
E - new SNICAR
blk - snicar_snobc_intmix = .true. (implemented via namelist, per Issue #36 )
and the next ABCDE is
A - xdrdt = 2 (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc)
B - scvng_fct_mlt_sf = 0.5 (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc)
C - snw_rds_refrz = 1500 (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc)
D - fresh_snw_rds_max = 400._r8 (implemented via namelist)
E - snw_rds_min = 100._r8 (implemented via mod to AerosolMod.F90)
See discussion in https://github.com/NCAR/LMWG_dev/discussions/3 and Issue #36
AD spunup to 2.75% disequilibrium, 280 years
pAD spunup to X.XX% disequilibrium, X years
Complete for now.
Description: Same as ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCD_1850AD (issue #37 ) but adding on snw_rds_min=100 Motivation is to try to lower snow albedo further. Just running an AD for now. Note that this uses ctsm5.1.dev150. NOTE: xdrdt=2 setting was not actually implemented due to a coding issue.
Case directory: Locally: /glade/work/oleson/ctsm5.1.dev150/cime/scripts/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCDE_1850AD
Sandbox: Locally: /glade/work/oleson/ctsm5.1.dev150
usernl changes:
user_nl_clm: 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' use_c13 = .true. use_c14 = .true. use_c13_timeseries = .true. use_c14_bombspike = .true. fsurdat = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc' paramfile = '/glade/p/cgd/tss/people/oleson/modify_param/ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc'
snicar_snobc_intmix = .true. fresh_snw_rds_max = 400.d
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/SoilTemperatureMod.F90 Sturm Snow tk src.clm/CNPhenologyMod.F90 decreasing temperature thresholds to -3C, reducing the minimum critical daylength by 25%, doubling the snow threshold (I think to 20 cm) src.clm/controlMod.F90 remove error associated with snicar_snobc_intmix src.clm/AerosolMod.F90 change snw_rds_min from default of 54.526 to 100.
Diagnostics: https://webext.cgd.ucar.edu/I1850/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCDE_1850AD/lnd/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCDE_1850AD.191_200-ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCD_1850AD.191_200/setsIndex.html
Output: Output (if still available): /glade/scratch/oleson/archive/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_blk_ABCDE_1850AD/
Contacts: @olyson @wwieder
Extra details: My case name lengths are out of control, so I've designated this combination simulation as ABCDE_blk_ABCDE, where the first ABCDE is: A - SturmSnowtk (implemented via src.clm/SoilTemperatureMod.F90) B - decidphenotrig (implemented via src.clm/CNPhenologyMod.F90) C - frootleaf (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc) D - FUNfracfixers (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc) E - new SNICAR blk - snicar_snobc_intmix = .true. (implemented via namelist, per Issue #36 ) and the next ABCDE is A - xdrdt = 2 (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc) B - scvng_fct_mlt_sf = 0.5 (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc) C - snw_rds_refrz = 1500 (implemented via ctsm51_params.c211112.nc.ABCD_blk_ABCD.nc) D - fresh_snw_rds_max = 400._r8 (implemented via namelist) E - snw_rds_min = 100._r8 (implemented via mod to AerosolMod.F90)
See discussion in https://github.com/NCAR/LMWG_dev/discussions/3 and Issue #36 AD spunup to 2.75% disequilibrium, 280 years pAD spunup to X.XX% disequilibrium, X years Complete for now.