NCAR / LMWG_dev

Repository to track LMWG development simulations
3 stars 0 forks source link

Sturm_vegMods + new SNICAR (ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_1850AD, 1850pAD, 1850) #33

Open olyson opened 11 months ago

olyson commented 11 months ago

Description: Same as ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1ABCD* but adding on new SNICAR optical properties. Motivation is that we expect higher snow albedos with the new SNICAR and we want to see how much that might offset the improvement in high latitude vegetation we gained in Sturm_vegMods. Note that this uses ctsm5.1.dev150.


Case directory: Locally: /glade/work/oleson/ctsm5.1.dev150/cime/scripts/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1ABCDE(1850AD, 1850pAD, 1850)


Sandbox: Locally: /glade/work/oleson/ctsm5.1.dev150


usernl changes: hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS', 'H2OSNO','NDEP_TO_SMINN' 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' finidat = '/glade/p/cgd/tss/people/oleson/CLM5_restarts/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_SturmSnowtk_1850pAD.clm2.r.0401-01-01-00000.nc' paramfile = '/glade/p/cgd/tss/people/oleson/modify_param/ctsm51_params.c211112.nc.ABCD.nc'


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)


Diagnostics: Compared to the GSWP3 control: https://webext.cgd.ucar.edu/I1850/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_1850/lnd/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_1850.11_20-ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_1850.11_20/setsIndex.html Compared to GSWP3 ABCD: https://webext.cgd.ucar.edu/I1850/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_1850/lnd/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_1850.11_20-ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_ABCD_1850.11_20/setsIndex.html


Output: Output (if still available): /glade/scratch/oleson/archive/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_1850AD/ /glade/scratch/oleson/archive/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_1850pAD/ /glade/scratch/oleson/archive/ctsm51_ctsm51d150_ne30pg3ne30pg3mg17_GSWP3V1_ABCDE_1850/


Contacts: @olyson @wwieder


Extra details: My case name lengths are out of control, so I've designated this combination simulation as ABCDE, where 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.nc) D - FUNfracfixers (implemented via ctsm51_params.c211112.nc.ABCD.nc) E - new SNICAR (but note that this simulation is using a new CTSM tag, so other changes are involved here)

See discussion in https://github.com/NCAR/LMWG_dev/discussions/3 AD spunup to 2.42% disequilibrium, 280 years pAD spunup to 5.77% disequilibrium, 640 years Complete for now