Ryan Modrak
Ryan Modrak
This was a very useful report that identifies a rounding or numerical issue that arises when strike and dip are close to 90 degrees. Hoping it can be investigated...
To plot observed and synthetic waveforms, we could do something like ``` from mtuq import MomentTensor from mtuq.graphics import plot_waveforms1 mt = MomentTensor([mrr,mtt,mpp,mrt,mrp,mtp]) plot_waveforms1(filename,data, greens.get_synthetics(mt), stations, origin, header=None) ``` Agreed,...
Here is how one could modify any of the _examples/_ scripts: ``` mt_dict = { "Mrr": -5800435174432632.0, "Mtt": 789342714264084.5, "Mpp": 5011092460168549.0, "Mrt": 3034156721905327.0, "Mrp": 2110077787951300.2, "Mtp": 2602039428905879.0, } mt_array =...
Alternatively, one could generate the synthetics first and pass them directly ``` from mtuq.graphics import plot_waveforms2 synthetics_bw = greens_bw.get_synthetics(mt) synthetics_sw = greens_sw.get_synthetics(mt) plot_waveforms2(event_id+'DC_waveforms.png', data_bw, data_sw, synthetics_bw, synthetics_sw, stations, origin) ````
Attached is a complete script illustrating the above [direct_evaluation.py](https://github.com/user-attachments/files/22963526/direct_evaluation.py)
a "misfit function evaluation" page could be added here? https://mtuqorg.github.io/mtuq/user_guide/06.html
Hi @Jinyin-Hu , Many thanks for the report. Stepping back, below are some very general comments that I hope might be relevant. Feedback Robert or others might also be helpful....
Drawing from private functions such as `_get_greens` may be error prone? Perhaps try the following to generate unprocessed MTUQ synthetics ``` db = open_db(path_cps, format='CPS') greens = db.get_greens_tensors(stations, origin) synthetics...
Hi @Jinyin-Hu , Perhaps I could defer to others for those questions, and focus for now on a related issue-- From speaking with Herrmann's and Ammon's students, my understanding is...
Moving forward, I think we could have an MTUQ example based on CPS Green's functions that exactly reproduces the current examples based on FK. I believe Robert and I independently...