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

1 files changed, 27 insertions, 0 deletions

diff --git a/ML_BSSN_Helper/schedule.ccl b/ML_BSSN_Helper/schedule.ccl
index a2134a1..48bea28 100644
--- a/ML_BSSN_Helper/schedule.ccl
+++ b/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"