From 55b13fb8954c526469eaf26fbeecfadeea4d575d Mon Sep 17 00:00:00 2001
From: Grzegorz Kowal <grzegorz@amuncode.org>
Date: Tue, 29 Aug 2023 17:09:48 -0300
Subject: [PATCH] EVOLUTION: Add missing variable update after adding
 turbulence forcing.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
---
 sources/evolution.F90 | 17 ++++++++---------
 1 file changed, 8 insertions(+), 9 deletions(-)

diff --git a/sources/evolution.F90 b/sources/evolution.F90
index 215a3ea..a847cf5 100644
--- a/sources/evolution.F90
+++ b/sources/evolution.F90
@@ -1088,30 +1088,29 @@ module evolution
 !
   subroutine advance(status)
 
-! references
-!
     use blocks     , only : set_blocks_update
     use coordinates, only : toplev
     use forcing    , only : update_forcing, forcing_enabled
     use mesh       , only : update_mesh
 
-! local variables are not implicit by default
-!
     implicit none
 
-! input variables
-!
     integer, intent(out) :: status
-!
+
 !-------------------------------------------------------------------------------
 !
 ! advance the solution using the selected method
 !
     call evolve()
 
-! add forcing contribution
+! add forcing contribution, requires the boundary and primitive variable update
 !
-    if (forcing_enabled) call update_forcing(time, dt)
+    if (forcing_enabled) then
+      call update_forcing(time, dt)
+
+      call update_variables(time + dt, 0.0d+00, status)
+      if (status /= 0) go to 100
+    end if
 
 ! check if we need to perform the refinement step
 !