Gravity_accelOneRow(integer(2)(IN) :: pos, integer(IN) :: sweepDir, integer(IN) :: blockID, integer(IN) :: numCells, real(:)(OUT) :: grav, integer(IN),optional :: potentialIndex)
This routine computes the gravitational acceleration for a row of cells in a specified direction in a given block.
pos : Row indices transverse to the sweep direction sweepDir : The sweep direction: allowed values are SWEEP_X, SWEEP_Y, and SWEEP_Z. These values are defined in constants.h. blockID : The local identifier of the block to work on numCells : Number of cells to update in grav() grav() : Array to receive result potentialIndex : Variable # to take as potential if present ( this argument is applicable only in selfgravity formulated as Poisson's equation)