From a20e8ff0ac13cd797b4011a8f66125e5ab58c19f Mon Sep 17 00:00:00 2001
From: wpbonelli <wbonelli@ucar.edu>
Date: Thu, 5 Dec 2024 10:46:22 -0500
Subject: [PATCH] fix no exit face termination with ifort?

---
 src/Solution/ParticleTracker/MethodDis.f90 | 19 ++++++++++++-------
 1 file changed, 12 insertions(+), 7 deletions(-)

diff --git a/src/Solution/ParticleTracker/MethodDis.f90 b/src/Solution/ParticleTracker/MethodDis.f90
index 61c67430f14..3700e3e80ce 100644
--- a/src/Solution/ParticleTracker/MethodDis.f90
+++ b/src/Solution/ParticleTracker/MethodDis.f90
@@ -100,6 +100,15 @@ subroutine load_cell(this, ic, cell)
       cell%zOrigin = cell%defn%bot
       cell%ipvOrigin = 1
 
+      factor = DONE / cell%defn%retfactor
+      factor = factor / cell%defn%porosity
+
+      areaz = dx * dy
+      term = factor / areaz
+
+      cell%vz1 = cell%defn%faceflow(6) * term
+      cell%vz2 = -cell%defn%faceflow(7) * term
+
       if (this%fmi%ibdgwfsat0(ic) == 0) then
         cell%vx1 = DZERO
         cell%vx2 = DZERO
@@ -111,19 +120,15 @@ subroutine load_cell(this, ic, cell)
       end if
 
       areax = dy * dz
-      areay = dx * dz
-      areaz = dx * dy
-      factor = DONE / cell%defn%retfactor
-      factor = factor / cell%defn%porosity
       term = factor / areax
       cell%vx1 = cell%defn%faceflow(1) * term
       cell%vx2 = -cell%defn%faceflow(3) * term
+
+      areay = dx * dz
       term = factor / areay
       cell%vy1 = cell%defn%faceflow(4) * term
       cell%vy2 = -cell%defn%faceflow(2) * term
-      term = factor / areaz
-      cell%vz1 = cell%defn%faceflow(6) * term
-      cell%vz2 = -cell%defn%faceflow(7) * term
+
     end select
   end subroutine load_cell