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