Don’t these lines imply that you’re modeling the log of sea level as quadratic?
rate = q*T^2 + a*T + c # a quadratic model on the rate of sea level
slr = pt.cumsum(rate)
Otherwise I would expect slr = (1 + rate).cumprod()
Don’t these lines imply that you’re modeling the log of sea level as quadratic?
rate = q*T^2 + a*T + c # a quadratic model on the rate of sea level
slr = pt.cumsum(rate)
Otherwise I would expect slr = (1 + rate).cumprod()