anchor.diff:             vel(2) = ymom(i,j,k) * rhoInv
              vel(3) = zmom(i,j,k) * rhoInv
   
  -           ! Account for symmetric boundaries.
  -
  -           vel = merge(ZERO, vel, physbc_lo(:) .eq. Symmetry)
  -           vel = merge(ZERO, vel, physbc_hi(:) .eq. Symmetry)
  -
              vel = inertial_velocity(pos, vel, time)
   
              g(1) = gravx(i,j,k)
  
positive.code:                                         vol, vo_lo, vo_hi, &
                                          lo, hi, dx, time, Qtt) bind(C,name='quadrupole_tensor_double_dot')
  
    use bl_constants_module, only: ZERO, THIRD, HALF, ONE, TWO
    use prob_params_module, only: center, problo, probhi, physbc_lo, physbc_hi, Symmetry
    use wdmerger_util_module, only: inertial_rotation, inertial_velocity
    use castro_util_module, only: position
  
  
negative.code:   use fundamental_constants_module, only: Gconst, c_light, parsec
    use probdata_module, only: gw_dist, axis_1, axis_2, axis_3
    use meth_params_module, only: rot_axis
    use prob_params_module, only: dim
  
    implicit none
  
  