diff --git a/GravityField/realfield.f90 b/GravityField/realfield.f90 index d43d72e..66c99c8 100644 --- a/GravityField/realfield.f90 +++ b/GravityField/realfield.f90 @@ -102,12 +102,14 @@ function compute_gravity(fx, fy, body) rx = fx - body%posx ry = fy - body%posy - dist = sqrt( (rx*rx) + (ry*ry) ) + ! ??? dist = sqrt( (rx*rx) + (ry*ry) ) + dist = (rx*rx) + (ry*ry) if (dist .LT. 0.08) then ! write (0, *) "dist too small ", dist compute_gravity = 0e0 else - compute_gravity = body%mass / (dist ** 2) + ! ??? compute_gravity = body%mass / (dist ** 2) + compute_gravity = body%mass / dist endif end function