SUBROUTINE qinit(Info)

Defined in:

contrib/astro/i_qinit.f90

Inputs:

NodeInfo Info

Description:

qinit is called after the routine setprob to initialize the conditions on the grid before beginning integration. It is called only once, at the beginning of the run. First the ambient values of q are set to be 0, except as otherwise specified in domain.data. See i_setprob for more information.

Then the routines PlaceBackPlate PlaceBackPlate(Info) and PlaceSidePlates PlaceSidePlates(Info) are called, in the case of a laboratory simulation, with values set in domain.data

Next, the qinit or Place routine of the specific module desired is called, which should set all initial values of q for the grid.

Finally, in the case of an MHD simulation (lMHD == .TRUE.) with MaintainAuxArrays across AMR boundaries turned on, the Auxiliary array is initialized. Also, in the case that Entropy is conserved (set in physics.data), the iEntropy (iE+1) field of q is calculated from the pressure and density.

Note: Entropy conservation is not fully supported by AstroBEAR, and is deprecated. Currently i_setprob.f90 is set to stop the program if iProtect is set to 3 or 4 (for Entropy conservation). May be deleted soon.

Called In:

• BearEZ Bearez.f90 via the function Initialq and included from InfoFieldUtils, either by InfoFieldUtils.mpi_cpp.f90 or InfoFieldUtils.single_cpp.f90, depending upon whether mpibear or single xbear is run.

• MpiExec via the function EstimateError, and subroutines RestrictFields and InitFields; included from InfoFieldUtils

Modules Used:

None.

Files Included:

None.