NOTE: 07-20-16, this evolved into the diagnostics.py script. Maybe this page should contain some descriptions of that script and how one might use it?

Ideally these are done in the code

Sum over PFTs

check that all the C and N stocks and fluxes are summed correctly over all the pfts npp, gpp, nuptake, vegC,

e.g., GPPAll on the plots is leaf, stem, root. Need to take GPPAll for each PFT and sum that and make sure it is the same as the GPPAll across the whole ecosystem.

• more important* Litterfall N, LitterfallC, same as above, but more important --> linkage between vegetation and soil

Close C balance (start here)

check that

change in stocks = sum of fluxes

do at ecosystem level (integrate soil and veg) do separately for veg and soil

veg

delta vegC: (sum Veg C across (root, stem, leaves)) = NPP - litterfallC - veg fire emission - deadC

deadC is "dead burned biomass"

soil

delta soilC: (sum of mineral, shallow, deep, moss C) = litterfallC - Rh - soil fire emission

dead woody debris (dwd)

FIX (this is more complicated):

delta dwd = ?? dwd(timestep-1) / k - dwdRH

where k is some constant.

Close N balance

check that

change in stocks = sum of fluxes

do at ecosystem level (integrate soil and veg) do separately for veg and soil

veg

delta vegN: (sum Veg N across (root, stem, leaves)) = NUptake - litterfallN - veg fire emission - deadN

deadN is "dead burned biomass"

soil

delta soil org N: (sum of mineral, shallow, deep, moss N) = litterfallN + netNmin

delta available N: netNmin - Nuptake + Ninput - N lost

Method

Run model for 5 years, extract annual output for 5 years so that you can compute the deltas Then extract all these variables,