P wave type (version 2)

[krisvanneste]

This feature branch implements inversion of P-wave spectra:

• Added 'rpp' configuration parameter for average P-wave radiation coefficient in sourcespec.conf.
• Start and end time of the P-window are computed, making sure it does not overlap with the S phase, and stored in the trace headers (in ssp_process_traces.py).
• Selection of P/S-window based on config.wave_type for computing spectra (in ssp_process_traces.py and ssp_build_spectra.py) and plotting the traces (ssp_plot_traces.py).
• Selection of Z-component only for P-wave spectra (in ssp_build_spectra.py).
• Selection of average P/S radiation coefficient or of P/S takeoff angle (when using focal mechanism) based on config.wave_type (in ssp_radiation_pattern.py).
• Selection of P/S velocities based on config.wave_type to convert displacement to moment (in ssp_build_spectra.py).
• Selection of P/S velocities based on config.wave_type to compute quality factor (in ssp_inversion.py) and to determine bounds on t_star from Qo bounds (in ssp_data_types.py).

[claudiodsf]

Hi @krisvanneste, I added some more commits and it's now good to go for me.

Please note that some configuration options have been renamed.

Could you please double-check?

[krisvanneste]

Hi @krisvanneste, I added some more commits and it's now good to go for me.

Please note that some configuration options have been renamed.

Could you please double-check?

Claudio,

I don't get the same results as before. I need some time to figure out where the differences come from.

[krisvanneste]

OK, my problem was due to incorrect setting of the 'rp_from_focal_mechanism' option. I now get exactly the same results for Mw, fc and t_star. Only the Qo values are different. I checked and this is due to a problem in the ssp_wave_arrival._wave_arrival function, where travel time calculation for constant velocity fails if vp_tt or vs_tt is None. However, this has not strictly something to do with this branch, I can open another issue for that.

[claudiodsf]

Ok, so I'll merge this one and wait for your issue on vp_tt and vs_tt.