Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

DERECHO + Sturm_vegMods + back off froot_leaf 50% + back_off FUN_fracfixers 50% + new SNICAR + SNICAR param changes w/xdrdt=5 (ctsm51_ctsm51d160_1deg_GSWP3V1_ABpt5Cpt5DE_blk_A5BCD_1850AD) #50

Closed
olyson opened this issue Jan 4, 2024 · 0 comments

Comments

@olyson
Copy link
Collaborator

olyson commented Jan 4, 2024

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


user_nl_ 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 #3
AD spunup to 2.79% disequilibrium, 320 years
pAD spunup to X.XX% disequilibrium, X years

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant