Package ANRAY, version 4.01

Ivan Psencik

The new version of the package ANRAY does not differ dramatically from the version 4.00. Nevertheless, limited modifications and corrections should make it, hopefully, more reliable and friendly to a user.

The most important changes are the following:

a) The formulae for the calculation of radiation patterns of waves generated by point sources situated in inhomogeneous anisotropic media were modified according to the corrections of the paper by Psencik & Teles (1996), see the Appendix. Modifications based on the formulae proposed by Klimes (1997) for the evaluation of the KMAH index in inhomogeneous anisotropic media have not yet been finished in this version.

b) A new strategy was used in dealing with problems of shear wave singularities. In the earlier versions, various modifications were used, which were often unreliable. The problem will require substantial modifications of the program in near future. In this version, whenever a ray gets to a region of a shear wave singularity, either inside a layer or at an interface, its calculation stops. In this way, calculation of a ray of a multiply reflected unconverted qP wave can be stopped if at some interface, generated qS waves would propagate into a singular direction. Due to this measure, the wavefield at some receivers, although illuminated by the considered wave, is not computed. The code, seems to work now more reliably in such situations.

c) Substantial changes have been made in the specification of input and output files. They should lead to a simpler manipulation with input and output files. In connection with these changes, minor changes in input data (practically only in the first lines) were also made.

d) A possibility to specify a moment tensor source by the elements of the tensor was added to the existing specifications of point sources.

e) Extensive debugging in programs as well as in their documentation was performed. Many confusing statements in the documentation were removed or substituted by, hopefully, more clear explanations.

Hardcopy of the paper

The paper is available in PostScript (68 kB) and GZIPped PostScript (20 kB).


In: Seismic Waves in Complex 3-D Structures, Report 6, pp. 179-181, Dep. Geophys., Charles Univ., Prague, 1997.
SW3D - main page of consortium Seismic Waves in Complex 3-D Structures .