Spinup the model

site scale (w/o libIGCM)

phase 1

We spinup the vegetation biomass by doing a normal run for 250 years imposing the NP content of leaves.

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. We do 10 cycles a 7 years.

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

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 250 years with dynamic stoichiometry using the restarts from phase 2.

RESTART_FILEIN = sechiba_rest_in.nc
STOMATE_RESTART_FILEIN = stomate_rest_in.nc

STOMATE_OK_NCYCLE=y
STOMATE_IMPOSE_CN=n

STOMATE_OK_PCYCLE=y
STOMATE_IMPOSE_NP=n