NAME
Hydro_detectShock
SYNOPSIS
call Hydro_detectShock(real(IN),pointer :: solnData(:,:,:,:),
real(OUT) :: shock(:,:,:),
integer(IN) :: blkLimits(2,MDIM),
integer(IN) :: blkLimitsGC(2,MDIM),
integer(IN) :: guardCells(MDIM),
real(IN) :: primaryCoord(:),
real(IN) :: secondCoord(:),
real(IN) :: thirdCoord(:))
DESCRIPTION
Hydro_detectShock is a multidimensional shock detection algorithm. This is currently used by the burner to cut off burning in shocks if desired. This is helpful for detonations in some circumstances.
ARGUMENTS
solnData -- pointer to a block of solution data to inspect
shock -- an array indicating if there is a shock in a zone (=1) =0 otherwise
blkLimits -- index limits of block, interior cells only; see Grid_getBlkIndexLimits
blkLimitsGC -- index limits of block, including guard cells.
guardCells -- number of layers of guard cells which output is requested in addition
to interior cells
primaryCoord -- x coordinate of solnData, i.e., coordinate in IAXIS direction, including guardcells
secondCoord -- y coordinate of solnData, i.e., coordinate in JAXIS direction, including guardcells
thirdCoord -- z coordinate of solnData, i.e., coordinate in KAXIS direction, including guardcells
NOTE
The guard cells need to be filled before calling this routine.
The pressure needs to be updated before calling this routine
{For future reference : The stub of this routine used to include
a final argument called integer, intent(IN) :: igc)
if =1, include the first row of guardcells in the detection
However, the implementation does not! May need to be addressed in
future}
SEE ALSO
Grid_getBlkIndexLimits