Hi all

I am trying to run a BH-NS simulation using at least one of some (and sometimes more) unofficial thors.

That would be Kadath for the initial data and possibly THC for hydro evolution.

I'm trying to run the simulation either with GRHydro or with THC.

In both cases I believe that the initial data gets imported correctly - I have looked at some simple visuals which I'm attaching to this mail.

However, as the simulation proceeds, huge hamiltonian constraint violations occur and e.x. the lapse gets driven to 0 globally. This can also be seen.

How can I help this? Is it because the puncture is not tracked correctly? Trying to set up a PunctureTracker at the location of the black hole yields a failure very soon, as the shift at puncture is said to be (nan, nan, nan).

I'm sending the .par files, outputs (.out), errors (.err) and visualization (.png).

Best regards