Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] acsh_v1_0.gz(81 Kbytes)
Manuscript Title: A finite element program package for magnetotelluric modelling.
Authors: E. Kisak, P. Silvester
Program title: E-PARALLEL FEMT-2D
Catalogue identifier: ACSH_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: 75K words
Word size: 32
Keywords: Geophysics, Magnetotelluric, Finite element, Laplace, Helmholtz, H-parallel, E-parallel, H-polarization, E-polarization, Cuthill-mckee, Sparse matrices.
Classification: 13.

Subprograms used:
Cat Id Title Reference
ACSJ_v1_0 ZFORMATS CPC 10(1975)421

Nature of problem:
Telluric and magnetotelluric response over two dimensional resistivity anomalies and topographic variations in the E-parallel case.

Solution method:
In addition to the features found in the program H-PARALLEL FEMT-2D a recursive approach has been developed to expand the region in space which will consist, a posterior, of a problem-dependent number of identical air strips. Field computations are restricted to nodes on or below the ground surface. Thus storage requirements are relatively not much larger than those of the H-parallel case.

Running time:
The same vertical fault treated in the H-parallel case has now 264 free nodes. It requires three air strips to yield a relative variation smaller than 1% in the values of the apparent resistivity. The total CPU time is 53 s.