Add an option to use the maximal lapse from libteukolskydata

2 files changed, 12 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index 7c3b340..cd81760 100644
--- a/param.ccl
+++ b/param.ccl
@@ -5,6 +5,11 @@ EXTENDS KEYWORD initial_data
   "teukolskydata" :: "Teukolsky wave initial data"
 }
 
+EXTENDS KEYWORD initial_lapse
+{
+  "teukolsky_max" :: "Teukolsky maximal slicing lapse"
+}
+
 RESTRICTED:
 CCTK_REAL amplitude "Wave amplitude A."
 {
diff --git a/src/teukolsky.c b/src/teukolsky.c
index 60cef99..234fedb 100644
--- a/src/teukolsky.c
+++ b/src/teukolsky.c
@@ -144,6 +144,13 @@ void teukolsky_data(CCTK_ARGUMENTS)
         }
     }
 
+    if (!strcmp(initial_lapse, "teukolsky_max")) {
+        ret = td_eval_lapse(td, cctk_lsh[2] * cctk_lsh[0], r_val, theta_val,
+                            (const unsigned int [2]){ 0, 0}, alp);
+        if (ret)
+            CCTK_WARN(0, "Error evaluating the lapse\n");
+    }
+
     free(r_val);
     free(theta_val);
     free(psi_val);