QINIT(APP)

(in application module)

Each application module has a subroutine (typically named something like qinitJet or qinitEnvelope) which sets up the initial conditions for the domain. By iterating over every point in the domain, all of the components of the state vector q must be specified.

The elements of q can vary, depending on geometry and dimensionality as specified in bear.data, and the use of MHD as specified in domain.data.

Here are the possible structures of the state vector q:

In 2D, Cartesian geometry: qCartesian2D

In 2D, cylindrical geometry with iCyl = 1: qCylindrical2D

In 2D, cylindrical geometry with iCyl = 2: qCylindrical2D2 (Note, this was a special tag assigned to handle angular momentum for hydro only, but is still in 2D)

Finally, in 3D, Cartesian geometry: qCartesian3D