Add new RHS boundary condition NewRad, implemented in the thorn LSUThorns/NewRad.

4 files changed, 67 insertions, 2 deletions

diff --git a/m/prototype/ML_BSSN_Helper/interface.ccl b/m/prototype/ML_BSSN_Helper/interface.ccl
index 8ad9ba2..4c57f29 100644
--- a/m/prototype/ML_BSSN_Helper/interface.ccl
+++ b/m/prototype/ML_BSSN_Helper/interface.ccl
@@ -6,5 +6,32 @@ USES INCLUDE: Symmetry.h
 
 
 
-CCTK_INT FUNCTION Boundary_SelectGroupForBC(CCTK_POINTER_TO_CONST IN GH, CCTK_INT IN faces, CCTK_INT IN boundary_width, CCTK_INT IN table_handle, CCTK_STRING IN group_name, CCTK_STRING IN bc_name)
+CCTK_INT FUNCTION MoLRegisterConstrainedGroup \
+    (CCTK_INT IN ConstrainedIndex)
+USES FUNCTION MoLRegisterConstrainedGroup
+
+CCTK_INT FUNCTION                         \
+    Boundary_SelectGroupForBC             \
+        (CCTK_POINTER_TO_CONST IN cctkGH, \
+          CCTK_INT IN faces,              \
+          CCTK_INT IN boundary_width,     \
+          CCTK_INT IN table_handle,       \
+          CCTK_STRING IN group_name,      \
+          CCTK_STRING IN bc_name)
 USES FUNCTION Boundary_SelectGroupForBC
+
+CCTK_INT FUNCTION                         \
+    ExtrapolateGammas                     \
+        (CCTK_POINTER_TO_CONST IN cctkGH, \
+         CCTK_REAL ARRAY INOUT var)
+USES FUNCTION ExtrapolateGammas
+
+CCTK_INT FUNCTION                         \
+    NewRad_Apply                          \
+        (CCTK_POINTER_TO_CONST IN cctkGH, \
+         CCTK_REAL ARRAY IN var,          \
+         CCTK_REAL ARRAY INOUT rhs,       \
+         CCTK_REAL IN var0,               \
+         CCTK_REAL IN v0,                 \
+         CCTK_INT IN radpower)
+USES FUNCTION NewRad_Apply
diff --git a/m/prototype/ML_BSSN_Helper/param.ccl b/m/prototype/ML_BSSN_Helper/param.ccl
index 605ea7a..338d0a2 100644
--- a/m/prototype/ML_BSSN_Helper/param.ccl
+++ b/m/prototype/ML_BSSN_Helper/param.ccl
@@ -1,3 +1,7 @@
+SHARES: Boundary
+
+USES CCTK_INT radpower
+
 SHARES: MethodofLines
 
 USES CCTK_INT MoL_Num_Constrained_Vars
@@ -15,6 +19,13 @@ SHARES: ML_BSSN
 USES CCTK_INT timelevels
 USES KEYWORD calculate_ADMBase_variables_at
 
+USES CCTK_REAL harmonicF
+
+EXTENDS KEYWORD my_rhs_boundary_condition
+{
+  "NewRad" :: "NewRad"
+}
+
 PRIVATE:
 
 CCTK_INT ML_BSSN_MaxNumConstrainedVars "Number of constrained variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Constrained_Vars
diff --git a/m/prototype/ML_BSSN_Helper/schedule.ccl b/m/prototype/ML_BSSN_Helper/schedule.ccl
index a2134a1..48bea28 100644
--- a/m/prototype/ML_BSSN_Helper/schedule.ccl
+++ b/m/prototype/ML_BSSN_Helper/schedule.ccl
@@ -62,6 +62,33 @@ if (CCTK_EQUALS (evolution_method, "ML_BSSN")) {
   
   
   
+  if (CCTK_EQUALS (my_rhs_boundary_condition, "NewRad"))
+  {
+    SCHEDULE ML_BSSN_ExtrapolateGammas AT initial AFTER ML_BSSN_convertFromADMBaseGamma
+    {
+      LANG: C
+      SYNC: ML_Gamma
+      SYNC: ML_dtlapse
+      SYNC: ML_dtshift
+    } "Extrapolate Gammas and time derivatives of lapse and shift"
+    
+    SCHEDULE ML_BSSN_NewRad IN ML_BSSN_evolCalcGroup AFTER ML_BSSN_RHS
+    {
+      LANG: C
+      SYNC: ML_curvrhs
+      SYNC: ML_dtlapserhs
+      SYNC: ML_dtshiftrhs
+      SYNC: ML_Gammarhs
+      SYNC: ML_lapserhs
+      SYNC: ML_log_confacrhs
+      SYNC: ML_metricrhs
+      SYNC: ML_shiftrhs
+      SYNC: ML_trace_curvrhs
+    } "Apply NewRad boundary conditions to RHS"
+  }
+  
+  
+  
   SCHEDULE GROUP ML_BSSN_convertToADMBaseGroup IN ML_BSSN_convertToADMBaseGroupWrapper
   {
   } "Calculate ADM variables"
diff --git a/m/prototype/ML_BSSN_Helper/src/make.code.defn b/m/prototype/ML_BSSN_Helper/src/make.code.defn
index dcfe66e..7cc4e72 100644
--- a/m/prototype/ML_BSSN_Helper/src/make.code.defn
+++ b/m/prototype/ML_BSSN_Helper/src/make.code.defn
@@ -1,3 +1,3 @@
 # -*-Makefile-*-
-SRCS = RegisterConstrained.c RegisterSlicing.c SelectBCsADMBase.c UnsetCheckpointTags.c
+SRCS = ExtrapolateGammas.c NewRad.c RegisterConstrained.c RegisterSlicing.c SelectBCsADMBase.c UnsetCheckpointTags.c
 #  CopyADMBase.c