Fortran package RMATRIX
Subroutine package to calculate the frequency-dependent
transmission coefficients written by Colin J. Thomson.
Copyright (c) 1996 C. J. Thomson.
All rights reserved by the author
RMATRIX and associated programs may be used freely by University
researchers and teachers.
Commercial parties interested in using
the program should contact C. J. Thomson.
The author can also provide information on
new or different versions of the program and how to obtain them
via anonymous ftp.
E-mail:
thomson@geol.queensu.ca
Compiling and linking programs of package MODEL
Used Fortran 90 features not present in Fortran 77
- Names longer than 6 characters.
- Underscores in variable names.
- Variable declarations REAL*8 and COMPLEX*16.
- COMPLEX*16 constants like (1.d0,0.d0).
- DO loops without labels like DO I=1,6 ... END DO.
Intrinsic functions of the VAX/IBM extensions present in the source code
- CDABS
- CDEXP
- DCONJG
- DCMPLX
- CDSQRT
- DREAL
- DIMAG
Changes in the source code made by L.Klimes
- 'rmatrix1.f'
- The file has been split into main program 'rmatrix0.f' and subroutines
'rmatrix1.f'.
- Originally the second line of the input data (read in subroutine INPUT)
is now read as the first line of the input data in the main program.
- Input logical unit number 10 has been changed to LU in subroutine INPUT.
- Dummy arguments PX,PY,OMEGA of subroutine INPUT have been discarded and
replaced with new dummy argument LU.
- Prints to the console (device *) have been disabled.
- Calls to output subroutines WRITE4, J1SYMTST and J2SYMTST
have been disabled.
- 'rmatrix2.f'
- All four calls to output subroutine WRITE1 have been disabled.
- Limits for reporting of the rounding errors in subroutine LAYER
have been moderated.
- 'isotroc.f'
- The correction for the normal incidence has been supplemented.
The changes are mostly denoted by string 'c-LK'.
Program 'rmatrix.for'
Main file 'rmatrix.for' contains Fortran 90
INCLUDE statements to include the files with the main
program RMATRIX and the subroutines required.
It is thus sufficient to compile and link single file 'rmatrix.for'.
Note the above mentioned changes in the input data for program
RMATRIX!
Running demo program RMATRIX
Demo input data are read from file 'rmatrix.inp'.
Program RMATRIX generates output file 'rmatrix.out'.
Note the above mentioned limitation of the output of program
RMATRIX.
File list of package RMATRIX
- (A) Documentation:
- rmatrix.htm (this file)
- (B) Demo program:
- rmatrix.for File to be compiled.
- rmatrix0.f Main program RMATRIX.
- (C) Subroutines to calculate the transmission coefficients:
- rmatrix1.f
- rmatrix2.f
- isotroc.f
- aniso6cg.f
- aniso6c1.f
- aniso3c.f
- deigv.f
- cg.f
- (D) Sample input data for program RMATRIX:
- rmatrix.inp