Introduction

This is a comparison between mineos and earth-flattened wavenumber integration synthetics (Computer Programs in Seismology - 3.30 ) for the AK-135f continental model

I use the same model for both sets of synthetics, except that the wavenumber integration model does not have an inner core [The reason for this is that my formulation for wavenumber integration only permits fluid layers at the top or at the bottom of the elastic stack but not sandwiched between two solid layers.

The program used for the conversion is conv.f and the starting model for an AK135-F continental is ak135-f. Running the program as 

g77 conv.f
a.out < ak135-f

creates the files tak135-f.txt for mineos and tak135sph.mod for the CPS programs.

This reformatting program does the following:

The model used for the wavenumber integration code is tak135sph.mod and the model used for mineos is tak135-f.txt[Note on your browser, you may have to right-click and save-as for the files ending with .mod, since the browser may think that the file is a MOD music file].

The earth flattening mapping used in the CPS programs starts with the tak135sph.mod model and then internally performs the  following steps, using the symbolism  a for the radius of the Earth and r for the distance from the center of the Earth to a spherical shell.

David Harkrider assisted with the density mapping for the P-SV problem.  He demonstrated that the 2.275 was appropriate for the fundamental mode.  This exercise focused on the entire waveform, which effectively tests the appropriateness for higher modes.

Program Installation

Tests

The tests performed consisted of comparing travel times, free-oscillation synthetics to wavenumber-integration synthetics that used Earth flattening, and a comparison of phase and group velocity dispersion between the free oscillation and plane-layer codes.





Last changed February 2, 2008