Temporarily disable FPE halting modes in diffusion solver to prevent …

…underflow
MintonGroup · Feb 21, 2022 · d67a569 · d67a569
1 parent 7f768c1
commit d67a569
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/src/util/util_diffusion_solver.f90 b/src/util/util_diffusion_solver.f90
@@ -22,6 +22,7 @@
 subroutine util_diffusion_solver(user,surf,N,indarray,kdiff,cumulative_elchange,maxhits)
    use module_globals
    use module_util, EXCEPT_THIS_ONE => util_diffusion_solver
+   use, intrinsic :: ieee_exceptions
    implicit none
 
    ! Arguments
@@ -40,7 +41,12 @@ subroutine util_diffusion_solver(user,surf,N,indarray,kdiff,cumulative_elchange,
    integer(I4B) :: mi,mip1,mim1,mj,mjp1,mjm1,nloops
    real(DP),dimension(N,N) :: elchange
    integer(I4B),dimension(6,N,N) :: indarray_extended
+   logical, dimension(size(IEEE_ALL))   :: fpe_halting_modes, fpe_quiet_modes
+   logical, dimension(size(IEEE_USUAL)) :: fpe_flag 
 
+   call ieee_get_halting_mode(IEEE_ALL,fpe_halting_modes)  ! Save the current halting modes so we can turn them off temporarily
+   fpe_quiet_modes(:) = .false.
+   call ieee_set_halting_mode(IEEE_ALL,fpe_quiet_modes)
 
    cumulative_elchange = 0.0_DP
    kdiffmax = maxval(kdiff)
@@ -122,6 +128,8 @@ subroutine util_diffusion_solver(user,surf,N,indarray,kdiff,cumulative_elchange,
       cumulative_elchange = cumulative_elchange + elchange 
    end do
 
+   call ieee_set_halting_mode(IEEE_ALL,fpe_halting_modes)  ! Restore the original FPE halting modes
+
 return
 end subroutine util_diffusion_solver