diff options
| -rw-r--r-- | param.ccl | 2 | ||||
| -rw-r--r-- | schedule.ccl | 2 | ||||
| -rw-r--r-- | src/TwoPunctures.c | 1 |
3 files changed, 3 insertions, 2 deletions
@@ -89,7 +89,7 @@ REAL TP_Extend_Radius "Radius of an extended spacetime instead of the puncture" 0:* :: "anything positive, should be smaller than the horizon" } 0.0 -REAL par_b "x coordinate of the m+ puncture" +REAL par_b "x coordinate of the m+ puncture" STEERABLE=always { (0.0:*) :: "" } 1.0 diff --git a/schedule.ccl b/schedule.ccl index a603446..4353c66 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -22,7 +22,7 @@ if (CCTK_Equals(initial_data, "twopunctures")) } else { - SCHEDULE TwoPunctures AT Initial AFTER ADMBase_InitialData BEFORE ADMBase_PostInitial AFTER HydroBase_Initial + SCHEDULE TwoPunctures AT Initial AFTER ADMBase_InitialData BEFORE ADMBase_PostInitial AFTER HydroBase_Initial before SetTmunu before HydroBase_Prim2ConInitial { LANG: C STORAGE: puncture_u bare_mass diff --git a/src/TwoPunctures.c b/src/TwoPunctures.c index d961815..e9284d4 100644 --- a/src/TwoPunctures.c +++ b/src/TwoPunctures.c @@ -229,6 +229,7 @@ TwoPunctures (CCTK_ARGUMENTS) } else { CCTK_INFO ("Solving puncture equation for BH-BH system"); } + CCTK_VInfo (CCTK_THORNSTRING, "b = %g", par_b); /* initialise to 0 */ for (int j = 0; j < ntotal; j++) |