search

NCAR / LMWG_dev

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

CRUJRA(w/FLDS) + d166deadveg branch + init tsoisno=272K + 1940-1949 atm forcing (ctsm51_ctsm51d166deadveg_1deg_CRUJRA1940-1949_FLDS_ABsnoCDE_blk_A5BCD_1850AD, 1850pAD) #55

Open olyson opened 7 months ago

olyson commented 7 months ago

Description: Same as #52 but using 1940-1949 forcing. Purpose is to see how different this simulation is from #52 (which uses our default years of 1901-1920 forcing). This will help assess if ERA5 forcing, which is only available starting in 1940, is adequate for our use in spinup and historical simulations. The AD and pAD simulations here will be followed by a historical that starts by looping over 1940-1949 until we get to the year 1940. Include init tsoisno=272K (Issue #2338)

Case directory: Locally: /glade/work/oleson/ctsm_deadveg/cime/scripts/ctsm51_ctsm51d166deadveg_1deg_CRUJRA1940-1949_FLDS_ABsnoCDE_blk_A5BCD_1850AD

Sandbox: Locally: /glade/work/oleson/ctsm_deadveg

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'

user_nl_datm_streams (note "datafiles" strings have been shortened below)

CLMCRUNCEPv7.Solar:tintalgo = coszen
CLMCRUNCEPv7.Solar:mapalgo = bilinear
CLMCRUNCEPv7.Solar:taxmode = cycle
CLMCRUNCEPv7.Solar:datavars = FSDS_DIRECT Faxa_swdndr,FSDS_DIFFUSE Faxa_swdndf
CLMCRUNCEPv7.Solar:meshfile = /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/lnd_mesh_cdf5.nc
CLMCRUNCEPv7.Solar:datafiles = /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/clmforc.TRENDY.c2023_0.5x0.5.Solr.1940_cdf5.nc,...,/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/clmforc.TRENDY.c2023_0.5x0.5.Solr.1949_cdf5.nc
CLMCRUNCEPv7.Precip:tintalgo = nearest
CLMCRUNCEPv7.Precip:mapalgo = bilinear
CLMCRUNCEPv7.Precip:taxmode = cycle
CLMCRUNCEPv7.Precip:datavars = PRECTmms Faxa_precn
CLMCRUNCEPv7.Precip:meshfile = /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/lnd_mesh_cdf5.nc
CLMCRUNCEPv7.Precip:datafiles = /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/clmforc.TRENDY.c2023_0.5x0.5.Prec.1940_cdf5.nc,...,/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/clmforc.TRENDY.c2023_0.5x0.5.Prec.1949_cdf5.nc
CLMCRUNCEPv7.TPQW:tintalgo = linear
CLMCRUNCEPv7.TPQW:mapalgo = bilinear
CLMCRUNCEPv7.TPQW:taxmode = cycle
CLMCRUNCEPv7.TPQW:datavars = QBOT Sa_shum,PSRF Sa_pbot,TBOT Sa_tbot,WIND Sa_wind, FLDS Faxa_lwdn
CLMCRUNCEPv7.TPQW:meshfile = /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/lnd_mesh_cdf5.nc
CLMCRUNCEPv7.TPQW:datafiles = /glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/clmforc.TRENDY.c2023_0.5x0.5.TPQWL.1940_cdf5.nc,...,/glade/campaign/cesm/cesmdata/cseg/inputdata/atm/datm7/atm_forcing.datm7.CRUJRA.0.5d.c2023/TRENDY_cdf5/clmforc.TRENDY.c2023_0.5x0.5.TPQWL.1949_cdf5.nc

SourceMods:

src.clm/clm_varpar.F90: nlevurb=5 -> nlevurb=10 src.clm/surfrdMod.F90: Remove read of PFTDATA_MASK src.clm/TemperatureType.F90: Set initial tsoisno to 272K

Diagnostics:

A comparison of the AD simulations (1940-1949 spin versus 1901-1920 spin) (without comment for now):

https://webext.cgd.ucar.edu/I1850/ctsm51_ctsm51d166deadveg_1deg_CRUJRA1940-1949_FLDS_ABsnoCDE_blk_A5BCD_1850AD/lnd/ctsm51_ctsm51d166deadveg_1deg_CRUJRA1940-1949_FLDS_ABsnoCDE_blk_A5BCD_1850AD_301_320-ctsm51_ctsm51d166deadveg_1deg_CRUJRA_FLDS_ABsnoCDE_blk_A5BCD_1850AD_301_320/setsIndex.html

Output: Output (if still available): /glade/derecho/scratch/oleson/archive/ctsm51_ctsm51d166deadveg_1deg_CRUJRA1940-1949_FLDS_ABsnoCDE_blk_A5BCD_1850AD/

Contacts: @olyson @wwieder

Extra details: I've designated this combination simulation as ABsnoCDE_blk_A5BCD, where the first ABsnoCDE is: A - SturmSnowtk (implemented in the default lnd_in) Bsno - snow5d_thresh_for_onset = 0.2 (implemented by default in ctsm51_params.c240208.nc, other CNPhenologyMod.F90 mods back to defaults) C - frootleaf (implemented by default in ctsm51_params.c240208.nc) D - FUN_fracfixers (implemented by default in ctsm51_params.c240208.nc) E - new SNICAR (implemented in the default lnd_in) blk - snicar_snobc_intmix = .true. (implemented in the default lnd_in) and the next A5BCD is A5 - xdrdt = 5 (implemented by default in ctsm51_params.c240208.nc) B - scvng_fct_mlt_sf = 0.5 (implemented by default in ctsm51_params.c240208.nc) C - snw_rds_refrz = 1500 (implemented by default in ctsm51_params.c240208.nc) D - fresh_snw_rds_max = 400._r8 (implemented by default in ctsm51_params.c240208.nc)

See discussion in https://github.com/NCAR/LMWG_dev/discussions/3 AD spunup to 2.71% disequilibrium, 320 years pAD spunup to 5.33% disequilibrium, 1160 years Complete for now.