#103 rotation

Closed
asteinbe wants to merge 2 commits from rotation_quick into master

Allows to set stored_quantity: rotation in the gf_store config when using qssp2017 and retrieves the rotation components instead of displacement components.

Allows to set stored_quantity: rotation in the gf_store config when using qssp2017 and retrieves the rotation components instead of displacement components.
emolch commented 2 months ago
Owner

Cool!

I’ll rebase it in a second.

Cool! I'll rebase it in a second.
hvasbath commented 2 months ago
Collaborator

Awesome! Andreas bekommt den Orden der Arbeit.

Awesome! Andreas bekommt den Orden der Arbeit.
emolch added the
Add a Test
label 2 months ago
emolch commented 2 months ago
Owner

In what units are the rotations given?

In what units are the rotations given?
asteinbe commented 2 months ago
Poster

same as displacement, there is an option to also output the roation rate, but for further analysis the roation seems to be used in the community.

same as displacement, there is an option to also output the roation rate, but for further analysis the roation seems to be used in the community.
emolch commented 2 months ago
Owner

Dimensionless as curl(u), right?

Dimensionless as curl(u), right?
asteinbe commented 2 months ago
Poster

yes, the seismogram than has the moments units.
https://library.seg.org/doi/epdf/10.1190/geo2017-0492.1

regarding the test: we should have in general a qssp test like we have for qseis (test_gf_qseis) to compare output with the direct qssp outputs. There we could also test the rotations output in comparison. I tried to implement that similar in test_gf_qssp but failed miserably at having the qssp runner outputs retrieved (At line 132 of file qpgrnspec.f (unit = 21)
Fortran runtime error: Cannot open file ‘/tmp/gfstorekfiv27fw/qssp_green/U_TEMP_10km_0km’: No such file or directory).

yes, the seismogram than has the moments units. https://library.seg.org/doi/epdf/10.1190/geo2017-0492.1 regarding the test: we should have in general a qssp test like we have for qseis (test_gf_qseis) to compare output with the direct qssp outputs. There we could also test the rotations output in comparison. I tried to implement that similar in test_gf_qssp but failed miserably at having the qssp runner outputs retrieved (At line 132 of file qpgrnspec.f (unit = 21) Fortran runtime error: Cannot open file '/tmp/gfstorekfiv27fw/qssp_green/U_TEMP_10km_0km': No such file or directory).
emolch commented 2 months ago
Owner

Test is working (with and without rotations): a86cb736f9

Not a very deep test but it checks that all is technically operational. Would be nice to compare with rotations computed numerically from the displacement field.

Test is working (with and without rotations): https://git.pyrocko.org/pyrocko/pyrocko/commit/a86cb736f9bcad084555c2bf0020d8bef6784adc Not a very deep test but it checks that all is technically operational. Would be nice to compare with rotations computed numerically from the displacement field.
emolch commented 2 months ago
Owner
Minimal doc entry is in now in the QSSP section of Fomosto: https://data.pyrocko.org/builds/a81f41ad99c17eb2ff090a2a47588ad8fade5105/docs/apps/fomosto/backends.html#the-qssp-backend Ready for master?
asteinbe commented 2 months ago
Poster

sure, cool, fridays fun projects.

sure, cool, fridays fun projects.
emolch removed the
Add a Test
label 2 months ago
emolch commented 2 months ago
Owner

Still on the wishlist:

  • add thorough test
  • add docs example
Still on the wishlist: * add thorough test * add docs example
emolch commented 2 months ago
Owner

I’ll merge and close now and open an issue with the wishlist... If you are already working on rotational seismograms, you may want to take care of these?

I'll merge and close now and open an issue with the wishlist... If you are already working on rotational seismograms, you may want to take care of these?
emolch commented 2 months ago
Owner

Manually rebased and pushed to master.

Manually rebased and pushed to master.
emolch closed this pull request 2 months ago
All checks were successful
continuous-integration/drone/push Build is passing
Please reopen this pull request to perform a merge.
Sign in to join this conversation.
No reviewers
No Milestone
No Assignees
3 Participants
Notifications
Due Date

No due date set.

Dependencies

This pull request currently doesn't have any dependencies.

Loading…
There is no content yet.