Programs in Physics & Physical Chemistry
|[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.|
Nature of problem:
Telluric and magnetotelluric response over two dimensional resistivity anomalies and topographic variations in the E-parallel case.
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.
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.
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|