Wave field output
Introduction
The wavefield can be written in an hdf5 file, in order to visualize it in ParaView. To speed up the process, it is recommended to dedicate a few nodes to the writing tasks (see Asynchronous Output).
Refinement
iOutputMask
iOutputMask allows switching on and off the writing of SeisSol unknowns. The 6 first digits controls the components of the stress tensor (sigma_xx, sigma_yy, sigma_zz, sigma_xy, sigma_yz, and sigma_xz), and the 3 last digits the velocity components (u, v, w). When using poroelasticity, 4 more flags are added, for pore pressure (p) and fluid velocities (u_f, v_f, w_f).
iPlasticityMask
iPlasticityMask allows switching on and off the writing of plasticity variables. The 6 first digits controls the components of the off-fault plastic strain tensor (ep_xx, ep_yy, ep_zz, ep_xy, ep_yz, and ep_xz), and the last one the accumulated plastic strain (eta).
OutputRegionBounds
Using the OutputRegionBounds parameter, under the &Output heading, in the parameter.par file, the user can define the region for which the output is to be written. This region is provided in the following format:
OutputRegionBounds = xMin xMax yMin yMax zMin zMax
OutputGroups
Similar to the previous parameter, OutputGroups can be used to whitelist a set of mesh groups (as specified in the xdmf mesh file) that are included in the wavefield output. Cells whose group is not mentioned are not included in the output. This feature works with OutputRegionBounds, only cells that satisfy both criteria are included. It looks like this:
OutputGroups = 1 2 ! only include groups 1 and 2
Example
&Output
OutputFile = '/output/prefix'
iOutputMask = 0 0 0 0 0 0 1 1 1
iPlasticityMask = 0 0 0 0 0 0 1
OutputRegionBounds = -5e3 5e3 -10e3 10e3 -8e3 0e0
Format = 6 ! Format (6=hdf5, 10= no output)
TimeInterval = 5.0 ! Index of printed info at time
printIntervalCriterion = 2 ! Criterion for index of printed info: 1=timesteps,2=time,3=timesteps+time
refinement = 1
/