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)) .

SPINUP_ANALYTIC= n

TIME_LENGTH=1000Y

STOMATE_OK_NCYCLE=y
STOMATE_IMPOSE_CN=y

STOMATE_OK_PCYCLE=y
STOMATE_IMPOSE_NP=y

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).

RESTART_FILEIN = sechiba_rest_in.nc
STOMATE_RESTART_FILEIN = stomate_rest_in.nc

SPINUP_ANALYTIC= y
SPINUP_PERIOD= 7
TIME_LENGTH=71Y

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 500 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=500Y
STOMATE_OK_NCYCLE=y
STOMATE_IMPOSE_CN=n

STOMATE_OK_PCYCLE=y
STOMATE_IMPOSE_NP=n