The deltaH method does not allow to model glacier extent exceeding the original state. Runs with a cumulative mass balance in any year have been dropped to far to make parameter optimization with thousands of runs possible. This was not an option for scenario runs so I included an optional routine to shift the surplus of years with a positive cumulative mass balance to subsequent years with a negative mass balance.