Bxy has to be of type Field2D, so get from TokamakCoordinatesFactory …

…object
tomc271 · Nov 10, 2024 · 680498a · 680498a
1 parent 9417099
commit 680498a
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/examples/elm-pb-outerloop/elm_pb_outerloop.cxx b/examples/elm-pb-outerloop/elm_pb_outerloop.cxx
@@ -288,6 +288,8 @@ class ELMpb : public PhysicsModel {
   int damp_width;        // Width of inner damped region
   BoutReal damp_t_const; // Timescale of damping
 
+  TokamakCoordinatesFactory tokamak_coordinates_factory = TokamakCoordinatesFactory(*mesh);
+
   const BoutReal MU0 = 4.0e-7 * PI;
   const BoutReal Mi = 2.0 * 1.6726e-27; // Ion mass
   const BoutReal Me = 9.1094e-31;       // Electron mass
@@ -704,7 +706,6 @@ class ELMpb : public PhysicsModel {
       Dphi0 *= -1;
     }
 
-    auto tokamak_coordinates_factory = TokamakCoordinatesFactory(*mesh);
     const auto& metric = tokamak_coordinates_factory.make_tokamak_coordinates();
 
     V0 = -tokamak_coordinates_factory.get_Rxy() * tokamak_coordinates_factory.get_Bpxy() * Dphi0 / tokamak_coordinates_factory.get_Bxy();
@@ -1522,7 +1523,7 @@ class ELMpb : public PhysicsModel {
     auto P0_acc = Field2DAccessor<>(P0);
     auto J0_acc = Field2DAccessor<>(J0);
     auto phi0_acc = Field2DAccessor<>(phi0);
-    FieldMetric Bxy = metric->Bxy();
+    const Field2D& Bxy = tokamak_coordinates_factory.get_Bxy();
     auto B0_acc = Field2DAccessor<>(Bxy);
 
     // Evolving fields