Programs in Physics & Physical Chemistry
|[Licence| Download | New Version Template] acsg_v1_0.gz(58 Kbytes)|
|Manuscript Title: A finite element program package for magnetotelluric modelling.|
|Authors: E. Kisak, P. Silvester|
|Program title: H-PARALLEL FEMT-2D|
|Catalogue identifier: ACSG_v1_0|
Distribution format: gz
|Journal reference: Comput. Phys. Commun. 10(1975)421|
|Programming language: Fortran.|
|Computer: IBM 360/75.|
|Operating system: OS/360 WITH HASP II.|
|Program overlaid: yes|
|RAM: 50K words|
|Word size: 32|
|Keywords: Geophysics, Magnetotelluric, Finite element, Laplace, Helmholtz, H-parallel, H-polarization, E-parallel, E-polarization, Cuthill-mckee, Sparse matrices.|
Nature of problem:
Telluric and magnetotelluric response over two dimensional resistivity anomalies and topographic variation in the H-parallel case.
The high-order triangular finite element method is used with polynomials of degree up to six. The resultant linear system of equations is not solved by an iterative process; rather, the complex symmetric positive- definite coefficient matrix is factorized into LDL T. The generation of nodes is mostly automatic and is followed by reordering according to the reverse Cuthill-Mckee algorithm to improve matrix profile and thus greatly reduce storage requirements. Surface electric field values are computed with the matrix representation of the normal derivative operator.
The usual assumptions of magnetotelluric modelling are assumed. Only geological structures having infinite length in one direction are treated. It is assumed that all boundaries within the problem region can be well approximated by straightline segments.
Data in the BLOCK DATA BLOCK1 of the program are punched in Z format, a feature of the IBM FORTRAN IV. This BLOCK DATA can be generated in other formats using the prorgam ZFORMATS described in this paper and which must be suitably modified.
The execution times are greatly problem-dependent. The sample problem of a vertical fault described in the long write-up and involving the solution of a complex linear system of 235 unknowns required 25 s in total CPU time on the IBM 360/75 computer.
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|