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

{{{
SPINUP_ANALYTIC= n

TIME_LENGTH=1000Y

STOMATE_OK_NCYCLE=y
STOMATE_IMPOSE_CN=y

STOMATE_OK_PCYCLE=y
STOMATE_IMPOSE_NP=y
}}}

[[Image(totalMass.jpg)]]

=== 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

}}}
[[Image(totC.png)]]

=== 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
}}}