diff options
Diffstat (limited to 'src/AHFinder_output.F')
-rw-r--r-- | src/AHFinder_output.F | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/AHFinder_output.F b/src/AHFinder_output.F index c5a492e..3f5b339 100644 --- a/src/AHFinder_output.F +++ b/src/AHFinder_output.F @@ -106,7 +106,7 @@ c some compilers cannot trim an empty string, so we add at least one char end if if (status.and.report) then - out_mass = 0.141047396D0 * dsqrt(intarea_h) + out_mass = 0.141047396D0 * sqrt(intarea_h) else out_mass = 0.0 end if @@ -218,9 +218,9 @@ c some compilers cannot trim an empty string, so we add at least one char do l=1,lmax do m=1,l if (mod(m,2).ne.0) then - out_asymx = out_asymx + dabs(cc(l,m)) + out_asymx = out_asymx + abs(cc(l,m)) else - out_asymx = out_asymx + dabs(cs(l,m)) + out_asymx = out_asymx + abs(cs(l,m)) end if end do end do @@ -244,7 +244,7 @@ c some compilers cannot trim an empty string, so we add at least one char if (status.and.report.and..not.refy) then do l=1,lmax do m=1,l - out_asymy = out_asymy + dabs(cs(l,m)) + out_asymy = out_asymy + abs(cs(l,m)) end do end do out_asymy = out_asymy / c0(0) @@ -267,10 +267,10 @@ c some compilers cannot trim an empty string, so we add at least one char out_asymz = 0.0D0 if (status.and.report.and..not.refz) then do l=1,lmax - if (mod(l,2).ne.0) out_asymz = out_asymz + dabs(c0(l)) + if (mod(l,2).ne.0) out_asymz = out_asymz + abs(c0(l)) do m=1,l if (mod(l-m,2).ne.0) then - out_asymz = out_asymz + dabs(cc(l,m)) + dabs(cs(l,m)) + out_asymz = out_asymz + abs(cc(l,m)) + abs(cs(l,m)) end if end do end do |