diff options
Diffstat (limited to 'src/AHFinder_mask.F')
-rw-r--r-- | src/AHFinder_mask.F | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/src/AHFinder_mask.F b/src/AHFinder_mask.F index 194f13c..cef8bbb 100644 --- a/src/AHFinder_mask.F +++ b/src/AHFinder_mask.F @@ -84,7 +84,7 @@ dd = max(dx,dy,dz) aux = ahf_maskbuffer*dd - if (dabs(buffer).lt.aux) buffer = - aux + if (abs(buffer).lt.aux) buffer = - aux ! ***************************************** @@ -101,7 +101,7 @@ do i=1,nx if (ahfgrid(i,j,k).ge.buffer) then - aux = dsqrt((x(i,j,k)-xc)**2 + (y(i,j,k)-yc)**2 + aux = sqrt((x(i,j,k)-xc)**2 + (y(i,j,k)-yc)**2 . + (z(i,j,k)-zc)**2) if (rhor.gt.aux) rhor = aux end if @@ -167,9 +167,9 @@ aux = 0.57D0*rhor - xa = dabs(x(i,j,k)-xc) - ya = dabs(y(i,j,k)-yc) - za = dabs(z(i,j,k)-zc) + xa = abs(x(i,j,k)-xc) + ya = abs(y(i,j,k)-yc) + za = abs(z(i,j,k)-zc) if ((xa.lt.aux).and.(ya.lt.aux).and.(za.lt.aux)) then ahmask(i,j,k) = zero @@ -202,9 +202,9 @@ aux = 0.7D0*rhor - xa = dabs(x(i,j,k)-xc) - ya = dabs(y(i,j,k)-yc) - za = dabs(z(i,j,k)-zc) + xa = abs(x(i,j,k)-xc) + ya = abs(y(i,j,k)-yc) + za = abs(z(i,j,k)-zc) if ((xa.lt.aux).and.(ya.lt.aux).and.(za.lt.aux).and. . ((xa+ya+za).lt.2.0D0*aux)) then |