[two-scale-heat-conduction] Fix parallel macro: Write data only for interior elements (#665)

* [two-scale-heat-conduction] Write data only for interior elements

* Formatting

---------

Co-authored-by: Ishaan Desai <ishaandesai@gmail.com>

Author: mathiskelm

two-scale-heat-conduction/macro-dumux/appl/main.cc

@@ -289,8 +289,15 @@ int main(int argc, char **argv)
     // linearize & solve
     nonLinearSolver.solve(x, *timeLoop);
 
-    for (int solIdx = 0; solIdx < numberOfElements; ++solIdx)
-      temperatures[solIdx] = x[solIdx][problem->returnTemperatureIdx()];
+    int solIdx = 0;
+    for (const auto &element : elements(leafGridView, Dune::Partitions::interior)) {
+      auto fvGeometry = localView(*gridGeometry);
+      fvGeometry.bindElement(element);
+      for (const auto &scv : scvs(fvGeometry)) {
+        temperatures[solIdx++] =
+            x[scv.elementIndex()][problem->returnTemperatureIdx()];
+      }
+    }
 
     if (getParam<bool>("Precice.RunWithCoupling") == true) {
       couplingParticipant.writeQuantityVector(meshName,