Spinup the model

site scale (w/o libIGCM)

phase 1

We spinup the vegetation biomass by doing a normal run for 1k years imposing the NP content of leaves. Figure below shows total vegetation biomass during (1k yrs a split into three simulations 1 (red), 1b (green), 1c (blue)) .

config.card

#-- leap(gregorian), noleap, 360d
CalendarType=noleap
#-- Start and End of Job
#-- "YYYY-MM-DD"
# 1k year
DateBegin=1949-01-01
DateEnd=2948-12-31

# Forcing data between 1901 and 1920; depending on your forcing data
CyclicBegin=1949
CyclicEnd=1968

PARAM/run.def:

SPINUP_ANALYTIC= n

STOMATE_OK_NCYCLE=y
STOMATE_IMPOSE_CN=y

STOMATE_OK_PCYCLE=y
STOMATE_IMPOSE_NP=y


## minimize IO ##

# SECHIBA history output level (0..10)
# default = 5
SECHIBA_HISTLEVEL = 1

# SECHIBA history 2 output level (0..10)
# default = 1
SECHIBA_HISTLEVEL2 = 1

# STOMATE history output level (0..10)
# default = 10
STOMATE_HISTLEVEL = 1

# Writefrequency in seconds in sechiba_history.nc
# default = 86400.0
#WRITE_STEP = 1800.

# Writefrequency in seconds sechiba_out_2.nc
# default = 1800.0
#WRITE_STEP2 = 10800.0

# Writefrequency in days in stomate_history.nc
# default = 10.
#STOMATE_HIST_DT = 1

phase 2

We do an analytic spinup using the restarts from phase 1. With 10 cycles a 7 years. Figure 1 below shows total vegetation biomass during analytic spinup phase using restarts from after 250 yr (green) and 1k yr (red) of phase 1, respectively. Figure 2 below shows the evolution of the slowest pool (SOM passive) for an earlier setting.

RESTART_FILEIN = sechiba_rest_in.nc
STOMATE_RESTART_FILEIN = stomate_rest_in.nc

SPINUP_ANALYTIC= y
SPINUP_PERIOD= 3
TIME_LENGTH=20Y

STOMATE_OK_NCYCLE=y
STOMATE_IMPOSE_CN=y

STOMATE_OK_PCYCLE=y
STOMATE_IMPOSE_NP=y

phase 3-C: final simulations for simulations w/o nutrient effects

We do a normal run for 250 years imposing the NP content of leaves using the restarts from phase 2.

RESTART_FILEIN = sechiba_rest_in.nc
STOMATE_RESTART_FILEIN = stomate_rest_in.nc

SPINUP_ANALYTIC= n

TIME_LENGTH=250Y

STOMATE_OK_NCYCLE=y
STOMATE_IMPOSE_CN=y

STOMATE_OK_PCYCLE=y
STOMATE_IMPOSE_NP=y

phase 3: intermediate step to introduce nutrient effects

We do a normal run for 1k years with dynamic stoichiometry using the restarts from phase 2.

RESTART_FILEIN = sechiba_rest_in.nc
STOMATE_RESTART_FILEIN = stomate_rest_in.nc

SPINUP_ANALYTIC= n

TIME_LENGTH=1000Y
STOMATE_OK_NCYCLE=y
STOMATE_IMPOSE_CN=n

STOMATE_OK_PCYCLE=y
STOMATE_IMPOSE_NP=n