call Particles_sinkAccelGasOnSinks(OPTIONAL,integer(IN) :: accelProps(MDIM))
Computes gas -> sinks gravitational accelerations by direct summation over all sink particles and grid cells. For cosmology, will also want to get contribution from PDE (mapped DM delegate particle density).
accelProps : optionally give the indices of the sink particle properties into which the gas-on-sink accelerations should be stored. Default is ACCX_PART_PROP, ACCY_PART_PROP, ACCZ_PART_PROP.
written by Robi Banerjee, 2007-2008 modified by Christoph Federrath, 2008-2012 ported to FLASH3.3/4 by Chalence Safranek-Shrader, 2010-2012 modified by Nathan Goldbaum, 2012 refactored for FLASH4 by John Bachan, 2012 debugged and renamed to reflect symmetry with Particles_sinkAccelSinksOnGas (Christoph Federrath, 2013) If accelProps is given but contains an invalid index (e.g., 0), the routine returns without updating any accelerations, but pt_sinkGatherGlobal will still have been called.