From 4ff2021a1f2f78e133d7c961ed740f65bbbe8b4e Mon Sep 17 00:00:00 2001 From: Andy Ridgwell Date: Thu, 15 Feb 2024 12:06:31 -0800 Subject: [PATCH] Add files via upload --- .../published/vervoortetal.2024/README.txt | 88 +++++++++++++ .../vervoortetal.2024/fkl_np30_EXP0_c.4Ma | 123 +++++++++++++++++ .../vervoortetal.2024/fkl_np30_EXP1_c.4Ma | 124 ++++++++++++++++++ .../vervoortetal.2024/fkl_np30_EXP2_c.4Ma | 121 +++++++++++++++++ .../vervoortetal.2024/fkl_np30_EXP3_c.4Ma | 121 +++++++++++++++++ .../vervoortetal.2024/fkl_np30_EXP4_c.4Ma | 120 +++++++++++++++++ .../vervoortetal.2024/fkl_np30_c.SPIN1 | 122 +++++++++++++++++ .../vervoortetal.2024/fkl_np30_c.SPIN2gl | 116 ++++++++++++++++ .../vervoortetal.2024/fkl_pp30_EXP0_c.4Ma | 123 +++++++++++++++++ .../vervoortetal.2024/fkl_pp30_EXP1_c.4Ma | 124 ++++++++++++++++++ .../vervoortetal.2024/fkl_pp30_EXP2_c.4Ma | 121 +++++++++++++++++ .../vervoortetal.2024/fkl_pp30_EXP3_c.4Ma | 121 +++++++++++++++++ .../vervoortetal.2024/fkl_pp30_EXP4_c.4Ma | 120 +++++++++++++++++ .../vervoortetal.2024/fkl_pp30_c.SPIN1 | 122 +++++++++++++++++ .../vervoortetal.2024/fkl_pp30_c.SPIN2gl | 116 ++++++++++++++++ 15 files changed, 1782 insertions(+) create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/README.txt create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP0_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP1_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP2_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP3_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP4_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_c.SPIN1 create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_c.SPIN2gl create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP0_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP1_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP2_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP3_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP4_c.4Ma create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_c.SPIN1 create mode 100644 genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_c.SPIN2gl diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/README.txt b/genie-userconfigs/PUBS/published/vervoortetal.2024/README.txt new file mode 100644 index 000000000..c7b1b3fb8 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/README.txt @@ -0,0 +1,88 @@ +################################################################ +### README.txt ################################################# +################################################################ + +For: +'Earth System Model Analysis of how Astronomical Forcing is Imprinted onto the Marine Geological Record: +The Role of the Inorganic (Carbonate) Carbon Cycle and Feedbacks' +P. Vervoort, S. Kirtland Turner, F. Rochholz, A. Ridgwell + +################################################################ +06/15/2020 -- README.txt file creation (PV) +07/01/2021 -- update (AR) +07/23/2021 -- updated with new experiment names (PV) +09/11/2023 -- add note (L.28); update experiment names +14/02/2024 -- updated with reference directory (PV) +################################################################ + +Provided is the code used to create the model experiments presented in the paper. +Also given are the configuration files necessary to run these experiments. + +The intention is to provide an opportunity to question the assumptions and +interpretation through re-analysis and the creation of new and different experiments. +(Plus, to provide a means to replicate results.) + +All experiments are run from: +$HOME/cgenie.muffin/genie-main +(unless a different installation directory has been used) + +### NOTE ######################################################### +% 09/11/2023 - The manuscript has had some updates after publication +% in order to correct for an error in the insolation code. Commands +% below are used to run the updated simulations with corrected code +################################################################# + +### model experiments -- spinups ################################ + +The commands to run the spinups are listed as follows: + +fkl_pp30 => hemispherically symmetric continental configuration with one continent from pole-to-pole +fkl_np30 => hemispherically asymmetric continental configuration with one continent in the NH + +(1a) INITIAL SPINUP + +The initial, 1st-stage closed system spin-up needs 20 kyr to reach equilibrium: + +./runmuffin.sh muffin.CBSG.fkl_pp30.BASES_PV PUBS/published/vervoortetal.2024 fkl_pp30_c.SPIN1 20000 +./runmuffin.sh muffin.CBSG.fkl_np30.BASES_PV PUBS/published/vervoortetal.2024 fkl_np30_c.SPIN1 20000 + +(1b) Second Stage SPINUP + +The follow-on, 2nd-stage (accelerated) open system spin-up requires a long spinup time of 1.5 Myr to reach equilibrium + +./runmuffin.sh muffin.CBSG.fkl_pp30.BASES_PV PUBS/published/vervoortetal.2024 fkl_pp30_c.SPIN2gl 1500000 fkl_pp30_c.SPIN1 +./runmuffin.sh muffin.CBSG.fkl_np30.BASES_PV PUBS/published/vervoortetal.2024 fkl_np30_c.SPIN2gl 1500000 fkl_np30_c.SPIN1 + +### model experiments -- main ensemble ########################## + +Below the commands to run the experiments with astronomical forcing for 4 Myr +Experiment 0-5 correspond to the numbering as described in the manuscript + +(2a) For EXPERIMENT 0, no carbon cycle feedbacks: + +./runmuffin.sh muffin.CBSG.fkl_pp30.BASES_PV PUBS/published/vervoortetal.2024 fkl_pp30_EXP0_c.4Ma 4000000 fkl_pp30_c.SPIN2gl +./runmuffin.sh muffin.CBSG.fkl_np30.BASES_PV PUBS/published/vervoortetal.2024 fkl_np30_EXP0_c.4Ma 4000000 fkl_np30_c.SPIN2gl + +(2b) For EXPERIMENT 1, climate + ocean circulation + CO2-ocean chemistry feedbacks only: + +./runmuffin.sh muffin.CBSG.fkl_pp30.BASES_PV PUBS/published/vervoortetal.2024 fkl_pp30_EXP1_c.4Ma 4000000 fkl_pp30_c.SPIN2gl +./runmuffin.sh muffin.CBSG.fkl_np30.BASES_PV PUBS/published/vervoortetal.2024 fkl_np30_EXP1_c.4Ma 4000000 fkl_np30_c.SPIN2gl + +(2c) For EXPERIMENT 2, EXPERIMENT 1 + primary productivity feedbacks: + +./runmuffin.sh muffin.CBSG.fkl_pp30.BASES_PV PUBS/published/vervoortetal.2024 fkl_pp30_EXP2_c.4Ma 4000000 fkl_pp30_c.SPIN2gl +./runmuffin.sh muffin.CBSG.fkl_np30.BASES_PV PUBS/published/vervoortetal.2024 fkl_np30_EXP2_c.4Ma 4000000 fkl_np30_c.SPIN2gl + +(2d) For EXPERIMENT 3, EXPERIMENT 2 + CaCO3 sedimentation feedbacks: + +./runmuffin.sh muffin.CBSG.fkl_pp30.BASES_PV PUBS/published/vervoortetal.2024 fkl_pp30_EXP3_c.4Ma 4000000 fkl_pp30_c.SPIN2gl +./runmuffin.sh muffin.CBSG.fkl_np30.BASES_PV PUBS/published/vervoortetal.2024 fkl_np30_EXP3_c.4Ma 4000000 fkl_np30_c.SPIN2gl + +(2e) For EXPERIMENT 4, EXPERIMENT 3 + weathering feedbacks: + +./runmuffin.sh muffin.CBSG.fkl_pp30.BASES_PV PUBS/published/vervoortetal.2024 fkl_pp30_EXP4_c.4Ma 4000000 fkl_pp30_c.SPIN2gl +./runmuffin.sh muffin.CBSG.fkl_np30.BASES_PV PUBS/published/vervoortetal.2024 fkl_np30_EXP4_c.4Ma 4000000 fkl_np30_c.SPIN2gl + +################################################################ +################################################################ +################################################################ \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP0_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP0_c.4Ma new file mode 100644 index 000000000..b8ed3b28a --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP0_c.4Ma @@ -0,0 +1,123 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set NO CO2-climate feedback +ea_36=n +# set radiative forcing to 3x278ppm = 834 ppm (2.99x278 = 831 ppm) +ea_radfor_scl_co2=2.99 +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# FIX biological production in surface ocean w/ forcing file +bg_par_bio_prodopt="bio_POCflux" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> defined in base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set a CLOSED system +bg_ctrl_force_sed_closedsystem=.true. +# set NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.false. +rg_opt_weather_T_Si=.false. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=19.47 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=6.1537e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +# derived from diag_weather_Ca.res of spin2 +rg_par_weather_CaSiO3=6.1537e+12 +rg_par_weather_CaCO3=6.1537e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.77 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify POC forcing field +bg_par_forcing_name="pyyyyz.FPOC_np30.vervoortetal2021" +# re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP1_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP1_c.4Ma new file mode 100644 index 000000000..2d45e5788 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP1_c.4Ma @@ -0,0 +1,124 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# FIX biological production in surface ocean w/ forcing file +bg_par_bio_prodopt="bio_POCflux" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> defined in base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# ******* Corg burial ******* +# N/A +# +# --- WEATHERING ------------------------------------------------------ +# +# set a CLOSED system +bg_ctrl_force_sed_closedsystem=.true. +# set NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.false. +rg_opt_weather_T_Si=.false. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=19.47 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=6.1537e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +# derived from diag_weather_Ca.res of spin2 +rg_par_weather_CaSiO3=6.1537e+12 +rg_par_weather_CaCO3=6.1537e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.77 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify POC forcing field +bg_par_forcing_name="pyyyyz.FPOC_np30.vervoortetal2021" +# re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP2_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP2_c.4Ma new file mode 100644 index 000000000..0eeb73ef0 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP2_c.4Ma @@ -0,0 +1,121 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> defined in base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set an OPEN system +bg_ctrl_force_sed_closedsystem=.true. +# set NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.false. +rg_opt_weather_T_Si=.false. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=19.47 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=6.1537e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +# derived from diag_weather_Ca.res of spin2 +rg_par_weather_CaSiO3=6.1537e+12 +rg_par_weather_CaCO3=6.1537e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.77 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +# re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP3_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP3_c.4Ma new file mode 100644 index 000000000..db71dd56a --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP3_c.4Ma @@ -0,0 +1,121 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> defined in base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set an OPEN system +bg_ctrl_force_sed_closedsystem=.false. +# set NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.false. +rg_opt_weather_T_Si=.false. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=19.47 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=6.1537e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +# derived from diag_weather_Ca.res of spin2 +rg_par_weather_CaSiO3=6.1537e+12 +rg_par_weather_CaCO3=6.1537e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.77 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +# re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP4_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP4_c.4Ma new file mode 100644 index 000000000..b2b93b4c1 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_EXP4_c.4Ma @@ -0,0 +1,120 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> defined in base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set an OPEN system +bg_ctrl_force_sed_closedsystem=.false. +# CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.true. +rg_opt_weather_T_Si=.true. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=19.47 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=7.0e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +rg_par_weather_CaSiO3=5.6e+12 +rg_par_weather_CaCO3=7.0e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.77 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +# re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_c.SPIN1 b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_c.SPIN1 new file mode 100644 index 000000000..c005ff8ad --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_c.SPIN1 @@ -0,0 +1,122 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> defined in base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.false. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set an OPEN system +bg_ctrl_force_sed_closedsystem=true. +# NO CaCO3 and CaSiO3 weathering-temperature feedback +#rg_opt_weather_T_Ca=.true. +#rg_opt_weather_T_Si=.true. +# weathering reference mean global land surface temperature (C) +#rg_par_ref_T0=19.45 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=7.0e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +rg_par_weather_CaSiO3=7.0e+12 +rg_par_weather_CaCO3=7.0e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=4.161 +# +# --- GEOCHEM ACCELERATION ------------------------------------- +# +#gl_ctrl_update_pCO2=.true. +#ma_gem_notyr_min=10 +#ma_gem_notyr_max=10 +#ma_gem_yr_min=90 +#ma_gem_yr_max=90 +#ma_gem_dyr=0 +#ma_gem_adapt_auto=.false. +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=834.0E-06 +bg_par_atm_force_scale_val_4=-4.9 +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# choose average orbit configuration +# mean ecc, mean obl (in rad), equinox at perihelion/aphelion so hemispheres receive the same amount of sun +ea_par_orbit_osce=0.026735 +ea_par_orbit_oscobl=0.40420 +ea_par_orbit_osclonperi=3.1416 +# +# --- MISC ----------------------------------------------------- +# +# Set initial ocean alkalinity +bg_ocn_init_12=1.975E-03 +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_c.SPIN2gl b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_c.SPIN2gl new file mode 100644 index 000000000..c52a05f96 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_np30_c.SPIN2gl @@ -0,0 +1,116 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> defined in base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set an OPEN system +bg_ctrl_force_sed_closedsystem=.false. +# NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.true. +rg_opt_weather_T_Si=.true. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=19.47 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=7.0e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +rg_par_weather_CaSiO3=5.6e+12 +rg_par_weather_CaCO3=7.0e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.77 +# +# --- GEOCHEM ACCELERATION ------------------------------------- +# +gl_ctrl_update_pCO2=.true. +ma_gem_notyr_min=10 +ma_gem_notyr_max=10 +ma_gem_yr_min=90 +ma_gem_yr_max=90 +ma_gem_dyr=0 +ma_gem_adapt_auto=.false. +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# choose average orbit configuration +# mean ecc, mean obl (in rad), equinox at perihelion/aphelion so hemispheres receive the same amount of sun +ea_par_orbit_osce=0.026735 +ea_par_orbit_oscobl=0.40420 +ea_par_orbit_osclonperi=3.1416 +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP0_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP0_c.4Ma new file mode 100644 index 000000000..98a8c59c4 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP0_c.4Ma @@ -0,0 +1,123 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set NO CO2-climate feedback +ea_36=n +# set radiative forcing to 3x278ppm = 834 ppm (3.00x278 = 834 ppm) +ea_radfor_scl_co2=3.00 +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# FIX biological production in surface ocean w/ forcing file +bg_par_bio_prodopt="bio_POCflux" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> also see base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set a CLOSED system +bg_ctrl_force_sed_closedsystem=.true. +# set NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.false. +rg_opt_weather_T_Si=.false. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=23.6645 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=6.9988e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +# derived from diag_weather_Ca.res of spin2 +rg_par_weather_CaSiO3=6.9988e+12 +rg_par_weather_CaCO3=6.9988e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.88 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify POC forcing field +bg_par_forcing_name="pyyyyz.FPOC_pp30.vervoortetal2021" +# re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP1_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP1_c.4Ma new file mode 100644 index 000000000..7d2e879ad --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP1_c.4Ma @@ -0,0 +1,124 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# FIX biological production in surface ocean w/ forcing file +bg_par_bio_prodopt="bio_POCflux" +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> also see base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# ******* Corg burial ******* +# N/A +# +# --- WEATHERING ------------------------------------------------------ +# +# set a CLOSED system +bg_ctrl_force_sed_closedsystem=.true. +# set NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.false. +rg_opt_weather_T_Si=.false. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=23.6645 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=6.9988e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +# derived from diag_weather_Ca.res of spin2 +rg_par_weather_CaSiO3=6.9988e+12 +rg_par_weather_CaCO3=6.9988e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.88 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify POC forcing field +bg_par_forcing_name="pyyyyz.FPOC_pp30.vervoortetal2021" +# re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP2_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP2_c.4Ma new file mode 100644 index 000000000..c550941a6 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP2_c.4Ma @@ -0,0 +1,121 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> also see base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set an OPEN system +bg_ctrl_force_sed_closedsystem=.true. +# set NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.false. +rg_opt_weather_T_Si=.false. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=23.6645 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=6.9988e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +# derived from diag_weather_Ca.res of spin2 +rg_par_weather_CaSiO3=6.9988e+12 +rg_par_weather_CaCO3=6.9988e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.88 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP3_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP3_c.4Ma new file mode 100644 index 000000000..50244b3fb --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP3_c.4Ma @@ -0,0 +1,121 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> also see base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set an OPEN system +bg_ctrl_force_sed_closedsystem=.false. +# set NO CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.false. +rg_opt_weather_T_Si=.false. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=23.6645 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=6.9988e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +# derived from diag_weather_Ca.res of spin2 +rg_par_weather_CaSiO3=6.9988e+12 +rg_par_weather_CaCO3=6.9988e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.88 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP4_c.4Ma b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP4_c.4Ma new file mode 100644 index 000000000..c822a509d --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_EXP4_c.4Ma @@ -0,0 +1,120 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> also see base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set an OPEN system +bg_ctrl_force_sed_closedsystem=.false. +# CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.true. +rg_opt_weather_T_Si=.true. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=23.6645 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=7.0e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +rg_par_weather_CaSiO3=6.9418e+12 +rg_par_weather_CaCO3=7.0e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.88 +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +#(saves annual averages, NOT discrete values in time) +bg_par_infile_sig_name='save_timeseries_EVERY001000_orbits.dat' +bg_par_infile_slice_name='save_timeslice_EVERY002500_orbits.dat' +# representative NH latitude and season to sample orbital variation of insolation +bg_par_t_sig_count_N=12 +bg_par_sig_j_N=17 +# SH latitude/season +bg_par_t_sig_count_S=12 +bg_par_sig_j_S=2 +# +# --- FORCINGS -------------------------------------------------------- +# +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# Call orbit_radfor +ea_38="y" +# Specify the type of orbital forcing default (0), time-varying (1), alternative config (2) +ea_39=1 +# Number of data points in orbits file +ea_40=5000 +# Interval between data points in goldstein time steps +ea_41=48000 +# filename for orbital parameters (must be in genie-embm/data/input) +ea_42="orbits_La2004_57.93_52.91Ma_new.dat" +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_c.SPIN1 b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_c.SPIN1 new file mode 100644 index 000000000..1650ab2e6 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_c.SPIN1 @@ -0,0 +1,122 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> also see base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.false. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set a CLOSED system during 1st spin-up +bg_ctrl_force_sed_closedsystem=.true. +# CaCO3 and CaSiO3 weathering-temperature feedback +#rg_opt_weather_T_Ca=.true. +#rg_opt_weather_T_Si=.true. +# weathering reference mean global land surface temperature (C) +#rg_par_ref_T0=23.785 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=7.0e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +rg_par_weather_CaSiO3=7.0e+12 +rg_par_weather_CaCO3=7.0e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=4.161 +# +# --- GEOCHEM ACCELERATION ------------------------------------- +# +#gl_ctrl_update_pCO2=.true. +#ma_gem_notyr_min=10 +#ma_gem_notyr_max=10 +#ma_gem_yr_min=90 +#ma_gem_yr_max=90 +#ma_gem_dyr=0 +#ma_gem_adapt_auto=.false. +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +# specify forcings +bg_par_forcing_name="pyyyyz.RpCO2_Rp13CO2" +bg_par_atm_force_scale_val_3=834.0E-06 +bg_par_atm_force_scale_val_4=-4.9 +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# choose average orbit configuration +# mean ecc, mean obl (in rad), equinox at perihelion/aphelion so hemispheres receive the same amount of sun +ea_par_orbit_osce=0.026735 +ea_par_orbit_oscobl=0.40420 +ea_par_orbit_osclonperi=3.1416 +# +# --- MISC ----------------------------------------------------- +# +# Set initial ocean alkalinity +bg_ocn_init_12=1.975E-03 +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file diff --git a/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_c.SPIN2gl b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_c.SPIN2gl new file mode 100644 index 000000000..5153c0000 --- /dev/null +++ b/genie-userconfigs/PUBS/published/vervoortetal.2024/fkl_pp30_c.SPIN2gl @@ -0,0 +1,116 @@ +# +# +# --- CLIMATE --------------------------------------------------------- +# +# set climate feedback +ea_36=y +# set seasonal cycle +ea_dosc=.true. +go_dosc=.true. +gs_dosc=.true. +# +# --- BIOLOGICAL NEW PRODUCTION --------------------------------------- +# +# maximum rate of conversion of dissolved PO4 into organic matter by phytoplankton (mol kg-1 yr-1) +bg_par_bio_k0_PO4=8.9876e-006 +# [PO4] M-M half-sat value (mol kg-1) +bg_par_bio_c0_PO4=8.9369e-007 +# +# --- ORGANIC MATTER EXPORT RATIOS ------------------------------------ +# +# production fraction of dissolved organic matter +bg_par_bio_red_DOMfrac=0.66 +# +# --- INORGANIC MATTER EXPORT RATIOS ---------------------------------- +# +# fixed CaCO3:POC +bg_opt_bio_CaCO3toPOCrainratio='FIXED' +# underlying export CaCO3 as a proportion of particulate organic matter (i.e., CaCO3/POC) +bg_par_bio_red_POC_CaCO3=0.200 +# +# --- REMINERALIZATION ------------------------------------------------ +# +# DOC lifetime (yrs) +bg_par_bio_remin_DOMlifetime=0.5 +# initial fractional abundance of POC component #2 +bg_par_bio_remin_POC_frac2=0.0557 +# depth of remineralization or particulate organic matter +bg_par_bio_remin_POC_eL1=589.9451 +# remineralization length #2 for POC +bg_par_bio_remin_POC_eL2=1000000.0 +# initial fractional abundance of CaCO3 component #2 +bg_par_bio_remin_CaCO3_frac2=0.45 +# depth of remineralization or CaCO3 +bg_par_bio_remin_CaCO3_eL1=1.8905e+003 +# remineralization length #2 for CaCO3 +bg_par_bio_remin_CaCO3_eL2=1000000.0 +# +# --- SEDIMENTS ---------------------------------------------------------- +# +# sediment grid options ==> also see base-config +sg_par_sed_Dmax_neritic=500.0 +# sediment diagenesis options +sg_par_sed_diagen_CaCO3opt="ridgwell2001lookup" +sg_ctrl_sed_bioturb=.true. +sg_ctrl_sed_bioturb_Archer=.false. +sg_par_n_sed_mix=20 +# set invarient surface bioturbation mixing rate +sg_par_sed_mix_k_sur_max=0.15 +sg_par_sed_mix_k_sur_min=0.15 +# additional detrital flux (g cm-2 kyr-1) +sg_par_sed_fdet=0.180 +# +# --- WEATHERING ------------------------------------------------------ +# +# set a CLOSED system during 1st spin-up +bg_ctrl_force_sed_closedsystem=.false. +# CaCO3 and CaSiO3 weathering-temperature feedback +rg_opt_weather_T_Ca=.true. +rg_opt_weather_T_Si=.true. +# weathering reference mean global land surface temperature (C) +rg_par_ref_T0=23.6645 +#CO2 outgassing rate (mol C yr-1) +rg_par_outgas_CO2=7.0e+12 +# global carbonate and silicate weathering rate (mol Ca2+ yr-1) +rg_par_weather_CaSiO3=6.9418e+12 +rg_par_weather_CaCO3=7.0e+12 +# d13C values +rg_par_outgas_CO2_d13C=-6.0 +rg_par_weather_CaCO3_d13C=13.88 +# +# --- GEOCHEM ACCELERATION ------------------------------------- +# +gl_ctrl_update_pCO2=.true. +ma_gem_notyr_min=10 +ma_gem_notyr_max=10 +ma_gem_yr_min=90 +ma_gem_yr_max=90 +ma_gem_dyr=0 +ma_gem_adapt_auto=.false. +# +# --- DATA SAVING ----------------------------------------------------- +# +# date saving & reporting options +bg_par_data_save_level=8 +bg_ctrl_debug_lvl0=.true. +ma_debug_loop=1 +# +# --- FORCINGS -------------------------------------------------------- +# +#re-scale gas transfer coeff to give ~0.058 mol m-2 yr-1 uatm-1 global mean air-sea coefficient (original: 0.310) +bg_par_gastransfer_a=0.715813093980993 +# select correct insolation code +ea_opt_orbit_old=.false. +# choose average orbit configuration +# mean ecc, mean obl (in rad), equinox at perihelion/aphelion so hemispheres receive the same amount of sun +ea_par_orbit_osce=0.026735 +ea_par_orbit_oscobl=0.40420 +ea_par_orbit_osclonperi=3.1416 +# +# --- MISC ----------------------------------------------------- +# +# add ventillation tracers +bg_ctrl_force_ocn_age1=.true. +# +# --- END ------------------------------------------------------ +# \ No newline at end of file