Use ieee intrinsics to avoid trapping ieee_inexact signals

Author: peterdschwartz

components/elm/src/biogeophys/SoilStateType.F90

@@ -331,6 +331,7 @@ subroutine InitCold(this, bounds)
     use SharedParamsMod   , only : ParamsShareInst
     use FuncPedotransferMod , only : pedotransf, get_ipedof
     use RootBiophysMod      , only : init_vegrootfr
+    use, intrinsic :: ieee_exceptions
     !
     ! !ARGUMENTS:
     class(soilstate_type) :: this
@@ -468,10 +469,11 @@ subroutine InitCold(this, bounds)
     ! Try to read soil information from the file.
     call ncd_io(ncid=ncid, varname='ZSOI', flag='read', data=zsoifl, dim1name=grlnd, readvar=readvar)
     if (.not. readvar ) then
-       !    Variable ZSOI not found, use the ELM parameters.
-       do j = 1, nlevsoifl
-          zsoifl(j) = scalez*(exp(zecoeff*(j-0.5_r8))-1._r8)    !node depths
-       end do
+       call ieee_set_flag(ieee_all,.false.)
+       call ieee_set_halting_mode(ieee_inexact, .false.)
+       do j = 1, nlevgrnd
+          zsoi(j) = scalez*(exp(zecoeff*(dble(j)-0.5_r8))-1._r8)    !node depths
+       enddo
     end if
 
     dzsoifl(1) = 0.5_r8*(zsoifl(1)+zsoifl(2))             !thickness b/n two interfaces

components/elm/src/main/initVerticalMod.F90

@@ -44,6 +44,7 @@ module initVerticalMod
 
   !------------------------------------------------------------------------
   subroutine initVertical(bounds, snow_depth, thick_wall, thick_roof)
+    use, intrinsic :: ieee_exceptions
     !
     ! !ARGUMENTS:
     type(bounds_type)   , intent(in)    :: bounds
@@ -170,8 +171,11 @@ subroutine initVertical(bounds, snow_depth, thick_wall, thick_roof)
        !    Soil layers not available from the input, and no additional layers needed. Use the
        ! default soil thickness settings.
        ! -----------------------------------------------------------------
+       call ieee_set_flag(ieee_all,.false.)
+       call ieee_set_halting_mode(ieee_inexact, .false.)
+
        do j = 1, nlevgrnd
-          zsoi(j) = scalez*(exp(zecoeff*(j-0.5_r8))-1._r8)    !node depths
+          zsoi(j) = scalez*(exp(zecoeff*(dble(j)-0.5_r8))-1._r8)    !node depths
        enddo
     end if
     deallocate(zsoi_in)