Gravity_computeDt(integer(IN) :: blockID, integer(IN) :: MyPE, real (OUT) :: dt_grav, integer(:)(INOUT) :: dt_minloc(5))
Compute the timestep limiter due to the gravitational solver.
MyPE: local processor number dt_grav: Will Return the limiting timestep. Should be set to a large value (1.D99) on input. dt_minloc(5): An array to receive information about which processor, block, and zone was responsible for setting the limiting timestep. The order is i, j, k, b, p, where (i,j,k) = zone indices, b = local block ID, and p = PE #. This routine should only modify these values if it changes dt_grav. blockID: The local ID of the block to compute the limiter on.