NCAR / LMWG_dev

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

Split between snow and surface water (ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_[1850AD,1850pAD,1850) #43

Open olyson opened 9 months ago

olyson commented 9 months ago

Description: Same as GSWP3V1 control (ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17GSWP3V1[1850AD,1850pAD,1850], Issue #1) but "split overlap between snow and surface water rather than reducing only surface water" per @swensosc . Purpose is to assess whether this has any effect on vegetation viability at high northern latitudes.

Case directory: Locally: /glade/work/oleson/cesm2_3_alpha16b/cime/scripts/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1h2osfc(1850AD, 1850pAD, 1850)

Sandbox: Locally: /glade/work/oleson/cesm2_3_alpha16b

user_nl_clm changes:

! 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' 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'

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/SurfaceWaterMod.F90: Implement split between fsno and surface water

Diagnostics:

(https://webext.cgd.ucar.edu/I1850/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_h2osfc_1850AD/lnd/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_h2osfc_1850AD.71_80-ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_fullout_1850AD.71_80/setsIndex.html)

Output: /glade/scratch/oleson/archive/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1h2osfc(1850AD, 1850pAD, 1850)

Contacts: @olyson @swensosc

Extra details: See discussion in https://github.com/NCAR/LMWG_dev/discussions/3 AD spunup to 20.63 % disequilibrium, 120 years pAD spunup to X.XX % disequilibrium, X years Complete for now.

First attempt at this run stalled at year 37 for no reason that I can see. Resubmitting with more frequent restart files. Second attempt was successful.