Comment (by Roland Haas):

As far as I can tell the parfile in the examples repo, which matches the parfile in the paper repo, never worked as advertised. Certainly not with the Curie release (the release version at the time the paper was prepared) and given that Maxwell (the release just after the paper was published) fails with NaN I cannot say if it possibly ever worked with Maxwell.

My own guess would be that when the data for the paper was produced, a modified version of McLachlan was used that had gauge evolution disabled. No version of McLachlan of any release has lapse_evolution_method appear anywhere other than in param.ccl so no version could have respected it being set to static. Since ML_BSSN e.g, has its own ML_lapse group, even having ADMBase reset alp to 1 after each step would not avoid the metric blowing up since alp is no used for input by ML_BSSN.

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2373/kasner-example-from-et-paper-cannot-be