call Hydro_shockStrength(real(IN),pointer :: solnData(:,:,:,:), real(INOUT) :: shock(:,:,:), integer(IN) :: blkLimits(2,MDIM), integer(IN) :: blkLimitsGC(2,MDIM), integer(IN) :: guardCells(MDIM), real(IN) :: primaryCoord(blkLimitsGC(LOW,IAXIS):blkLimitsGC(HIGH,IAXIS)), real(IN) :: secondCoord(blkLimitsGC(LOW,IAXIS):blkLimitsGC(HIGH,IAXIS)), real(IN) :: thirdCoord(blkLimitsGC(LOW,IAXIS):blkLimitsGC(HIGH,IAXIS)), real(IN) :: threshold, integer(IN) :: mode)
Hydro_shockStrength computes a measure of shock strength for each cell in a block. Hydro_shockStrength includes the same multidimensional shock detection algorithm as Hydro_detectShhock, from which it is derived.
solnData -- pointer to a block of solution data to inspect shock -- an array returning a measure of shock strenght for cells where a shock (of sufficient strength, see threshold) is detected. A combination of current shock strenght with the previous value may be returned if mode > 1, see mode below. 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 threshold - threshold value for shock strength mode - 1 to forget previous values in shock array, 2 to add to prev value, 3 for max of previous and current value NOTE The guard cells need to be filled before calling this routine. The pressure needs to be updated before calling this routine