Restore mesh->getCoordinates() calls

since TokamakCoordinates function now doesn't return Coordinates.

Author: tomchapman

examples/6field-simple/elm_6f.cxx

@@ -378,6 +378,9 @@ class Elm_6f : public PhysicsModel {
 
     bool noshear;
 
+    // Get the metric tensor
+    Coordinates* coord = mesh->getCoordinates();
+
     output.write("Solving high-beta flute reduced equations\n");
     output.write("\tFile    : {:s}\n", __FILE__);
     output.write("\tCompiled: {:s} at {:s}\n", __DATE__, __TIME__);

examples/conducting-wall-mode/cwm.cxx

@@ -62,6 +62,9 @@ class CWM : public PhysicsModel {
     // Load 2D profiles (set to zero if not found)
     GRID_LOAD(Ni0, Te0);
 
+    coord = mesh->getCoordinates();
+
+
     // Load normalisation values
     GRID_LOAD(Te_x, Ni_x, bmag);
 

examples/constraints/alfven-wave/alfven.cxx

@@ -160,6 +160,9 @@ class Alfven : public PhysicsModel {
 
   void LoadMetric(BoutReal Lnorm, BoutReal Bnorm) {
 
+    // Get the coordinates object
+    Coordinates* coord = mesh->getCoordinates();
+
     bool noshear;
 
     // Check type of parallel transform

examples/dalf3/dalf3.cxx

@@ -154,6 +154,8 @@ class DALF3 : public PhysicsModel {
       return 1;
     }
 
+    Coordinates* coord = mesh->getCoordinates();
+
     // SHIFTED RADIAL COORDINATES
 
     // Check type of parallel transform

examples/elm-pb-outerloop/elm_pb_outerloop.cxx

@@ -357,6 +357,7 @@ class ELMpb : public PhysicsModel {
   // Note: The rhs() function needs to be public so that RAJA can use CUDA
 
   int init(bool restarting) override {
+    Coordinates* metric = mesh->getCoordinates();
 
     output.write("Solving high-beta flute reduced equations\n");
     output.write("\tFile    : {:s}\n", __FILE__);

examples/elm-pb/elm_pb.cxx

@@ -300,7 +300,7 @@ class ELMpb : public PhysicsModel {
 
         bool noshear;
 
-        Coordinates *metric = mesh->getCoordinates();
+        Coordinates* metric = mesh->getCoordinates();
 
         output.write("Solving high-beta flute reduced equations\n");
         output.write("\tFile    : {:s}\n", __FILE__);

examples/gyro-gem/gem.cxx

@@ -345,6 +345,7 @@ class GEM : public PhysicsModel {
     output << "\tNormalised rho_e = " << rho_e << endl;
     output << "\tNormalised rho_i = " << rho_i << endl;
 
+    coord = mesh->getCoordinates();
     // Set B field vector
 
     B0vec.covariant = false;

examples/laplacexy/laplace_perp/test.cxx

@@ -18,6 +18,7 @@ int main(int argc, char** argv) {
     set_tokamak_coordinates_on_mesh(tokamak_coordinates, *mesh, true, 1.0, 1.0);
   }
 
+    Coordinates* coord = mesh->getCoordinates();
   ///////////////////////////////////////
 
   // Read an analytic input

tests/integrated/test-drift-instability/2fluid.cxx

@@ -86,6 +86,7 @@ class TwoFluid : public PhysicsModel {
     b0xcv.covariant = false;  // Read contravariant components
     mesh->get(b0xcv, "bxcv"); // b0xkappa terms
 
+    coord = mesh->getCoordinates();
     // Load normalisation values
     GRID_LOAD(Te_x);
     GRID_LOAD(Ti_x);

tests/integrated/test-interchange-instability/2fluid.cxx

@@ -54,6 +54,7 @@ class Interchange : public PhysicsModel {
 
     b0xcv *= -1.0; // NOTE: THIS IS FOR 'OLD' GRID FILES ONLY
 
+    coord = mesh->getCoordinates();
     // Load normalisation values
     GRID_LOAD(Te_x);
     GRID_LOAD(Ti_x);