Bugs of TidalStrain

We released the Fortran program "TidalStrain" that calculates the theoretical tidal response of the solid earth on an arbitrary fault plane at an arbitrary depth from the earth’s surface to the core-mantle boundary (~2891 km) on 27 August 2019.

After that, some bugs of "TidalStrain" were founded.

This page reports bugs of "TidalStrain" and the revised version.

Bugs

• "TidalStrain" included bugs in GOTIC2 (ver. 2001.05.16). Please see "13. Version history" in Matsumoto [2004] for the details of bugs.
• We found errors in calculating for an event of shallower than 10 km in deep.
• Depending on computing environment, an error case appears.

Corrections / Improvements

• We reconstructed the program "TidalStrain.2" based on GOTIC2 (ver. 2004.10.25) for computation of oceanic tidal loading effect. We also corrected above other bugs.
  Please see following figure for the comparison between old and new.
• We also enclosed calculation codes of tidal phase angle, tidal strain/stress level, and principal stress.

*Acknowledgments: Our program "TidalStrain.2" uses a most part of subroutines in "GOTIC2 (version 2004.10.15) [Matsumoto et al., 2001; Matsumoto, 2004]" to calculate ocean tide loading effects.

Note that "source_mod2020/*_mod2020.f" indicate subroutines that are different from original subroutines in "GOTIC2".

References
Matsumoto, K., T. Sato, T. Takanezawa, and M. Ooe, 2001, GOTIC2: A program for computation of oceanic tidal loading effect, J. Geod. Soc. Japan, 47, 243-248.
Matsumoto, K., 2004, https://www.miz.nao.ac.jp/staffs/nao99/, we browsed it on June 3, 2020.