diff options
Diffstat (limited to 'src/AHFinder.F')
-rw-r--r-- | src/AHFinder.F | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/src/AHFinder.F b/src/AHFinder.F index 364e301..8bd32ff 100644 --- a/src/AHFinder.F +++ b/src/AHFinder.F @@ -969,8 +969,8 @@ aux = intexp2_h/intarea_h - (intexp_h/intarea_h)**2 if (aux.ge.zero) then - aux = dsqrt(aux) - if (dabs(intexp_h/intarea_h)-aux.gt.zero) horizon = .false. + aux = sqrt(aux) + if (abs(intexp_h/intarea_h)-aux.gt.zero) horizon = .false. else horizon = .false. end if @@ -1070,18 +1070,18 @@ ! ! M = sqrt( A / (16 pi) ) = 0.141047396 sqrt(A) - ahfmass(mfind) = 0.141047396D0*dsqrt(intarea_h) + ahfmass(mfind) = 0.141047396D0*sqrt(intarea_h) write(*,*) write(*,*) 'Surface found, details below.' write(*,*) write(*,"(A21,ES14.6)") ' Surface area =',intarea_h write(*,"(A21,ES14.6)") ' Surface mass =', - . 0.141047396D0*dsqrt(intarea_h) + . 0.141047396D0*sqrt(intarea_h) write(*,"(A21,ES14.6)") ' Mean value of H =',aux*intexp_h write(*,"(A21,ES14.6)") ' Mean value of H^2 =',aux*intexp2_h write(*,"(A21,ES14.6)") ' Standard deviation =', - . dsqrt(dabs(aux*intexp2_h - (aux*intexp_h)**2)) + . sqrt(abs(aux*intexp2_h - (aux*intexp_h)**2)) write(*,*) if (ahf_gaussout.ne.0) then write(*,"(A11,ES14.6)") ' circ_eq =',circ_eq @@ -1096,11 +1096,11 @@ write(1,*) write(1,"(A21,ES14.6)") ' Surface area =',intarea_h write(1,"(A21,ES14.6)") ' Surface mass =', - . 0.141047396D0*dsqrt(intarea_h) + . 0.141047396D0*sqrt(intarea_h) write(1,"(A21,ES14.6)") ' Mean value of H =',aux*intexp_h write(1,"(A21,ES14.6)") ' Mean value of H^2 =',aux*intexp2_h write(1,"(A21,ES14.6)") ' Standard deviation =', - . dsqrt(dabs(aux*intexp2_h - (aux*intexp_h)**2)) + . sqrt(abs(aux*intexp2_h - (aux*intexp_h)**2)) write(1,*) if (ahf_gaussout.ne.0) then write(1,"(A11,ES14.6)") ' circ_eq =',circ_eq |