Hi @tom, **co2conc** in current implementation is a constant value. What I meant is to pass a new time-series instead of MaunaLoaCO2.txt (for example to run 2021-2100 crop production).

I found another issue with the code in **solution.py: line 296**. Instead of `for ii in range(compi):`

it should be `for ii in range(compi+1):`

, so you need to add **+1**. Example: in python if compi = 1, iteration would go only one time with ii = 0, but we need to have it twice with ii = 0 and ii = 1 (as in MATLAB). Because of this, the second compartment (ii = 1) will have thfcAdj = 0, which leads to miscalculation of water balance and hence everything else. If you need more explanations, reach me by email