diff options
Diffstat (limited to 'src/findnormals.F90')
| -rw-r--r-- | src/findnormals.F90 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/findnormals.F90 b/src/findnormals.F90 index 4c963cf..518679c 100644 --- a/src/findnormals.F90 +++ b/src/findnormals.F90 @@ -82,7 +82,7 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk) ! Find magnitude of (sx,sy,sz). - smag = dsqrt(sx**2 + sy**2 + sz**2) + smag = sqrt(sx**2 + sy**2 + sz**2) ! Normalize (sx,sy,sz). @@ -125,7 +125,7 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk) vz = dble(kk) dot = (sx*vx + sy*vy + sz*vz) & - / dsqrt(vx**2 + vy**2 + vz**2) + / sqrt(vx**2 + vy**2 + vz**2) ! If the new product is larger than the ! largest so far, redifine the normal. |