Fix Field3D::setBoundaryTo for FCI methods

src/field/field3d.cxx

@@ -491,8 +491,28 @@ void Field3D::setBoundaryTo(const Field3D& f3d) {
 
   allocate(); // Make sure data allocated
 
-  /// Loop over boundary regions
+  if (isFci()) {
+    // Set yup/ydown using midpoint values from f3d
+    ASSERT1(f3d.hasParallelSlices());
+    ASSERT1(hasParallelSlices());
+
+    for (auto& region : fieldmesh->getBoundariesPar()) {
+      for (const auto& pnt : *region) {
+        // Interpolate midpoint value in f3d
+        const BoutReal val = pnt.interpolate_sheath_o1(f3d);
+        // Set the same boundary value in this field
+        pnt.dirichlet_o1(*this, val);
+      }
+    }
+  }
+
+  // Non-FCI.
+  // Transform to field-aligned coordinates?
+  // Loop over boundary regions
   for (const auto& reg : fieldmesh->getBoundaries()) {
+    if (isFci() && reg->by != 0) {
+      continue;
+    }
     /// Loop within each region
     for (reg->first(); !reg->isDone(); reg->next()) {
       for (int z = 0; z < nz; z++) {

tests/integrated/test-fci-mpi/CMakeLists.txt

@@ -3,7 +3,7 @@ bout_add_mms_test(test-fci-mpi
   USE_RUNTEST
   USE_DATA_BOUT_INP
   PROCESSORS 6
-  DOWNLOAD https://zenodo.org/record/7614499/files/W7X-conf4-36x8x128.fci.nc?download=1
+  DOWNLOAD https://zenodo.org/records/14221309/files/W7X-conf0-36x8x128.fci.nc?download=1
   DOWNLOAD_NAME grid.fci.nc
   REQUIRES BOUT_HAS_PETSC
 )

tests/integrated/test-laplace-petsc3d/CMakeLists.txt

@@ -6,5 +6,5 @@ bout_add_integrated_test(test-laplace-petsc3d
     data_slab_core/BOUT.inp
     data_slab_sol/BOUT.inp
   USE_RUNTEST
-  REQUIRES BOUT_HAS_PETSC
+  REQUIRES BOUT_HAS_PETSC BOUT_ENABLE_ALL_TESTS
   )

tests/unit/include/bout/test_single_index_ops.cxx

@@ -276,6 +276,9 @@ TEST_F(SingleIndexOpsTest, Div_par) {
 
   // Need parallel derivatives of input
   input.calcParallelSlices();
+  // and of coordinates
+  input.getMesh()->getCoordinates()->J.calcParallelSlices();
+  input.getMesh()->getCoordinates()->g_22.calcParallelSlices();
 
   // Differentiate whole field
   Field3D difops = Div_par(input);