From b3fd38367400a7e4929d2015ad647d0ebd35a600 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Max=20L=C3=BCbke?= <mluebke@gfz-potsdam.de>
Date: Wed, 10 Apr 2024 09:57:56 +0000
Subject: [PATCH] Fix diffusion and chem field update order in RunMasterLoop

---
 src/poet.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/poet.cpp b/src/poet.cpp
index 68c39592f..7dac258af 100644
--- a/src/poet.cpp
+++ b/src/poet.cpp
@@ -259,8 +259,6 @@ static Rcpp::List RunMasterLoop(const RuntimeParameters &params,
 
     chem.simulate(dt);
 
-    diffusion.getField().update(chem.getField());
-
     // MDL master_iteration_end just writes on disk state_T and
     // state_C after every iteration if the cmdline option
     // --ignore-results is not given (and thus the R variable
@@ -273,6 +271,8 @@ static Rcpp::List RunMasterLoop(const RuntimeParameters &params,
           *global_rt_setup, trans_field, chem_field);
     }
 
+    diffusion.getField().update(chem.getField());
+
     MSG("End of *coupling* iteration " + std::to_string(iter) + "/" +
         std::to_string(maxiter));
     MSG();