Buoyancy

\[\dmombuo\]

src/fluid/predict/ux.c

static int buoyancy (
    const domain_t * domain,
    const double * restrict t,
    double * restrict src
) {
  // impose it only when desired
  if (!param_add_buoyancy) {
    return 0;
  }
  const int isize = domain->mysizes[0];
  const int jsize = domain->mysizes[1];
#if NDIMS == 3
  const int ksize = domain->mysizes[2];
#endif
#if NDIMS == 2
  BEGIN
    src[cnt] +=
      + 0.5 * T(i-1, j  )
      + 0.5 * T(i  , j  );
  END
#else
  BEGIN
    src[cnt] +=
      + 0.5 * T(i-1, j  , k  )
      + 0.5 * T(i  , j  , k  );
  END
#endif
  return 0;
}