On Fri, Oct 05, 2012 at 05:30:18PM -0400, Petr Tsatsin wrote:
Hi, I did a run with TmunuBase::support_old_CalcTmunu_mechanism = yes and it works.
Hi Petr,
Can you please try the attached patch to ADMMacros with TmunuBase::support_old_CalcTmunu_mechanism = no and see if that still works and gives the same results? You will have to add inheritance of some thorns to TmunuBase though, or maybe add that to ADMMacros (which is not done by the patch).
However at the beginning several iterations i got 'poison values' after the values ok.
This probably means that past timelevels haven't been correctly initialized by the initial data routines.
Another thing is that Surface Integrals Values osculate from 0 to correct value for the ADM mass. Volume Integrated ADM Mass looks ok.
My best guess would again be that something isn't initializing past timelevels correctly. This would show in correct values whenever refinement levels 'exist' at the time of integration. '0' is however a pretty strange value. Something else must be wrong as well.
You could try to see whether you see something strange using the following parameter:
ADMMass_Debug = "yes"
Frank