New in version 2.1.
For relevant parameters, please also see Shock Finding Parameters.
Shock finding is accomplished using one of four methods. The primary method uses a coordinate-unsplit temperature jump (method 1), as described in Skillman et. al. 2008 with the exception that instead of searching across multiple grids for the pre- and post-shock cells, we terminate the search at the edge of the ghost zones within each grid.
Shock finding is controlled by the ShockMethod parameter, which can take the following values:
0 - Off
1 - Unsplit Temperature Jumps
2 - Dimensionally Split Temperature Jumps
3 - Unsplit Velocity Jumps
4 - Dimensionally Split Velocity Jumps
When ShockMethod nonzero, this will create a “Mach” field in the output files.
Note: Method 1 has been used the most by the developer, and therefore is the primary method. Method 2 has been tested quite a bit, but the downsides of using a dimensionally split method are outlined in the above paper. Methods 3 and 4 are more experimental and will run, but results may vary.
Additional Shock Finding Parameters:
ShockTemperatureFloor - When calculating the mach number using temperature jumps, set the temperature floor in the calculation to this value.
StorePreShockFields - Optionally store the Pre-shock Density and Temperature during data output.