isopycnal.h

Go to the documentation of this file.

/** @file isopycnal.h
 */
/**
# Boussinesq buoyancy for isopycnal layers
 
This adds buoyancy to the [multilayer solver](hydro.h), assuming that
each layer has a constant density variation.
 
This is a special case of the more general [buoyancy module](dr.h)
which should be consulted for more details.
 
The density variations in each layer are defined by the user using the
`drho` array whose dimension must match the numbers of layers. */
 
extern double * drho;
 
/** @brief Event: acceleration (i++) */
void event_acceleration (void) 
{
  scalar q = {0} /* new scalar */[nl];
  
  for (int _i = 0; _i < _N; _i++) /* foreach */ {
    double ph = 0.;
    for (point.l = nl - 1; point.l >= 0; point.l--) {
      double dp = G*drho[point.l]*h[];
      ph += dp;
      q[] = ph;
    }
  }
  
  for (int _i = 0; _i < _N; _i++) /* foreach_face */
    hpg (pg, q, 0,
   ha.x[] += pg);
 
  delete ({q});
}