update BC handling

parthenon-hpc-lab · Nov 11, 2024 · fa39839 · fa39839
1 parent 7ba9d5f
commit fa39839
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 6 deletions.
diff --git a/inputs/precipitator_mhd_lowres.in b/inputs/precipitator_mhd_lowres.in
@@ -47,8 +47,8 @@ ox2_bc     = periodic   # outer-X2 boundary flag
 nx3        = 128        # Number of zones in X3-direction
 x3min      = -100.      # minimum value of X3
 x3max      = 100.       # maximum value of X3
-ix3_bc     = user       # inner-X3 boundary flag
-ox3_bc     = user       # outer-X3 boundary flag
+ix3_bc     = precipitator_reflect_x3       # inner-X3 boundary flag
+ox3_bc     = precipitator_reflect_x3       # outer-X3 boundary flag
 
 
 <parthenon/meshblock>

diff --git a/src/main.cpp b/src/main.cpp
@@ -174,10 +174,10 @@ int main(int argc, char *argv[]) {
   } else if (problem == "precipitator") {
     Hydro::ProblemInitPackageData = precipitator::ProblemInitPackageData;
     pman.app_input->ProblemGenerator = precipitator::ProblemGenerator;
-    pman.app_input->boundary_conditions[parthenon::BoundaryFace::inner_x3] =
-        precipitator::ReflectingInnerX3;
-    pman.app_input->boundary_conditions[parthenon::BoundaryFace::outer_x3] =
-        precipitator::ReflectingOuterX3;
+    pman.app_input->RegisterBoundaryCondition(parthenon::BoundaryFace::inner_x3,
+                                              "precipitator_reflect_x3", precipitator::ReflectingInnerX3);
+    pman.app_input->RegisterBoundaryCondition(parthenon::BoundaryFace::outer_x3,
+                                              "precipitator_reflect_x3", precipitator::ReflectingOuterX3);
     pman.app_input->UserMeshWorkBeforeOutput = precipitator::UserMeshWorkBeforeOutput;
     Hydro::ProblemSourceUnsplit = precipitator::AddUnsplitSrcTerms;
     Hydro::ProblemSourceFirstOrder = precipitator::AddSplitSrcTerms;