gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

EVOLUTION: Use ieee_is_nan() instead of isnan() in IBM BG/P.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2014-08-13 07:31:05 -03:00
parent 02a67e86b2
commit 39b560065d
1 changed files with 12 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
src/evolution.F90
View File

@ -843,6 +843,9 @@ module evolution
! !
use coordinates , only : im, jm, km use coordinates , only : im, jm, km
use equations , only : nv, pvars, cvars use equations , only : nv, pvars, cvars
#ifdef IBM
use, intrinsic :: ieee_arithmetic
#endif /* IBM */
! include external variables ! include external variables
! !
@ -882,12 +885,21 @@ module evolution
do j = 1, jm do j = 1, jm
do i = 1, im do i = 1, im
do p = 1, nv do p = 1, nv
#ifdef IBM
if (ieee_is_nan(pdata%u(p,i,j,k))) then
print *, 'U NaN:', cvars(p), pdata%meta%id, i, j, k
end if
if (ieee_is_isnan(pdata%q(p,i,j,k))) then
print *, 'Q NaN:', pvars(p), pdata%meta%id, i, j, k
end if
#else /* IBM */
if (isnan(pdata%u(p,i,j,k))) then if (isnan(pdata%u(p,i,j,k))) then
print *, 'U NaN:', cvars(p), pdata%meta%id, i, j, k print *, 'U NaN:', cvars(p), pdata%meta%id, i, j, k
end if end if
if (isnan(pdata%q(p,i,j,k))) then if (isnan(pdata%q(p,i,j,k))) then
print *, 'Q NaN:', pvars(p), pdata%meta%id, i, j, k print *, 'Q NaN:', pvars(p), pdata%meta%id, i, j, k
end if end if
#endif /* IBM */
end do ! p = 1, nv end do ! p = 1, nv
end do ! i = 1, im end do ! i = 1, im
end do ! j = 1, jm end do ! j = 1, jm