The alpha equation in this model has been changed to a d2dt2 scheme, rather than a hard-coded correction on a ddt scheme that relies on the old-time value remaining unchanged from the previous timestep. This modification makes the model compatible with some forms of mesh change.