Skip to content

Commit

Permalink
fix no exit face termination with ifort?
Browse files Browse the repository at this point in the history
  • Loading branch information
wpbonelli committed Dec 5, 2024
1 parent 83e6671 commit a20e8ff
Showing 1 changed file with 12 additions and 7 deletions.
19 changes: 12 additions & 7 deletions src/Solution/ParticleTracker/MethodDis.f90
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand All @@ -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

Expand Down

0 comments on commit a20e8ff

Please sign in to comment.