Pressure Seismogram

[int-ptr-ptr]

Is your feature request related to a problem? Please describe. A pressure seismogram is a feature in the original SPECFEM that would be beneficial to add in SPECFEM++.

Describe the solution you'd like One should be able to add -pressure to seismogram-type in the configuration file, and produce seismogram outputs for pressure at stations.

Additional context Add any other context or screenshots about the feature request here.

[int-ptr-ptr]

For this to work, as well as making it easier for other seismogram types to be added, the current compute_seismograms routine may need to be reworked, since additional shared views may be needed.

Moreover, there is currently no good way to isolate code between seismogram types that need more than just a "use this field, instead"-type difference in implementation.

For our pressure case, we need to compute the divergence of the displacement field (or for the acoustic case, just take the acceleration field). I can see two ways of doing it:

• As in commit 7e9d40bf949dc5708e1c8e418ddb8c61ba3f947d, we can add an auxiliary field within the receiver-seismogram parfor (requiring more scratch size), store the displacement field in there, then compute the divergence. This aux-field is used under the assumption that the receiver kernel's get_field() should be ignorant of the seismogram type.
• Alternatively, if we let get_field() access seismogram_type_l, then the divergence can be computed directly from the displacement field. Additionally, the acoustic kernel can just return the negative acceleration of $\chi$ and skip the $\chi$->displacement conversion.