Dear all

Is there a way to add different gauge conditions for the lapse? I'm trying to add a gauge shock avoid lapse for the Bona-Masso family of slicing conditions

dot alpha = -alpha^2 f(alpha) K

with f(alpha)=1+constant/alpha^2

I don't know if I should add a thorn that complements the McLachlan thorn, or what is the correct procedure to implement the gauge condition? Furthermore, this kind of gauge condition allows to the lapse become negative which I know is a serious drawback when the lapse collapses to zero, however I want to do some numerical experiments with this kind of lapse.

Any help is very welcome. Thanks

Erik