Hello All,
My name is Kyle Hussey, a current undergraduate student at The College of New Jersey, and I am doing research under Dr. Shannon Graham. We would like to use the Green’s Functions that are created when running step 7 in the subduction 3D, and use them in another inversion program, namely TDefnode. To get a sense of our project and goal, here is a link to a paper that goes through a similar process
https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2014GL062505
We will be using it to study the subduction zone in Mexico. The output files for the green’s functions are in .h5 format. However, TDefnode just takes ASCII files. So we would like to convert or read the output of the .h5 file so that we can just get a list of the station’s positions after the impulse. Any help would be greatly appreaciated and if you need more information from me, please let me know.
Thank you,
Kyle Hussey
Hi Kyle,
Sorry for the slow response. I have done what you’re talking about, primarily with an older version of PyLith (v. 2.x). I have written a package that should allow this with either PyLith 2.x or 3.x+, however. The git repository is part of the following paper: https://academic.oup.com/gji/article/240/2/886/7906614. If I’m remembering correctly, the package uses existing TDEFNODE Green’s functions as templates and then integrates the PyLith Green’s functions to produce TDEFNODE-type Green’s functions. Let me know if you’re able to access the repository, and if the code works OK for you. As I mentioned, I did this with an older version of PyLith, so it hasn’t really been tested with newer versions.
Cheers,
Charles
Thank you so much for this resource! I have access to the repository and will be trying it out. I’ll let you know if we run into any issues and if we are able to solve them or not. We will try incorporating your code with Pylith 4.2.1. If we need to fall back on an older version, that shouldn’t be an issue. But I will talk with my advisor about the best course of action if we run into too many problems