diff --git a/src/dynamics/senh/dyn_comp.F90 b/src/dynamics/senh/dyn_comp.F90
index da9ed024ee..95931ff1f4 100644
--- a/src/dynamics/senh/dyn_comp.F90
+++ b/src/dynamics/senh/dyn_comp.F90
@@ -37,7 +37,7 @@ module dyn_comp
 use thread_mod,              only: nthreads, hthreads, vthreads, omp_get_max_threads, omp_get_thread_num
 use time_mod,                only: nsplit,tstep
 use time_manager,            only: is_first_step
-use ieee_arithmetic,         only: ieee_is_nan
+use shr_infnan_mod,          only: isnan => shr_infnan_isnan
 
 implicit none
 private
@@ -1943,7 +1943,7 @@ subroutine read_dyn_field_2d(fieldname, fh, dimname, buffer)
    ! to NaN.  In that case infld can return NaNs where the element GLL points
    ! are not "unique columns"
    ! Set NaNs or fillvalue points to zero
-   where (ieee_is_nan(buffer) .or. (buffer==fillvalue)) buffer = 0.0_r8
+   where (isnan(buffer) .or. (buffer==fillvalue)) buffer = 0.0_r8
 
 end subroutine read_dyn_field_2d
 
@@ -1973,7 +1973,7 @@ subroutine read_dyn_field_3d(fieldname, fh, dimname, buffer)
    ! to NaN.  In that case infld can return NaNs where the element GLL points
    ! are not "unique columns"
    ! Set NaNs or fillvalue points to zero
-   where (ieee_is_nan(buffer) .or. (buffer == fillvalue)) buffer = 0.0_r8
+   where (isnan(buffer) .or. (buffer == fillvalue)) buffer = 0.0_r8
 
 end subroutine read_dyn_field_3d