allow multiple frequency/period bands inputs for adjoint source creation
Closed this issue · 1 comments
bch0w commented
This issue was sparked by discussions with groups at Colorado School of Mines and UToronto.
During a workflow it can be useful isolate certain phases by systematically filtering and windowing for a collection of non-continuous period bands (e.g. 17-40s + 40 - 110s). This generates a set of adjoint sources, which can then be summed up to give to final adjoint source that is fed back into the numerical solver for an adjoint simulation.
This capability is not currently available in Pyatoa - at the moment Users simply set a minimum and maximum period range, effectively limiting the period band to one range. To implement this into Pyatoa, one option could be:
- Input a range of period bands into Config
- Manager class breaks a Stream into multiple traces (one for each period band)
- Each trace is filtered, windowed and measured (adjoint source) individually (using existing framework)
- At the write adjoint source stage, all adjoint sources for a given component are summed and written