Modify data['vars'] in place. Use Temp(t-1) for calculating productivity fraction.
Change magnitude of constant in additive damages.
For decade rather than year? Dunno why but the output seems to be correct this way.