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] aczt_v1_0.gz(6 Kbytes)
Manuscript Title: Algorithm for the first order hydrostatic ellipticity of a planet.
Authors: W. Moon
Program title: HYDELP
Catalogue identifier: ACZT_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 19(1980)63
Programming language: Fortran.
Computer: IBM 3033.
Operating system: OS.
RAM: 51K words
Word size: 8
Keywords: Geophysics, Hydrostatic ellipticity, Figure of planet.
Classification: 13.

Nature of problem:
In geophysics the first order ellipticity of a planet is often required for a given model. In this new method the hydrostatic ellipticity of a planet can be calculated as a function of the mean radius. The input data needed are a set of values of the density of the planet as a function of the mean radius. The density model may have several jump discontinuities.

Solution method:
The density model (with or without jump discontinuities) is interpolated by cubic Hermite basis functions and the Darwin-de Sitter-Radau transformed first order ellipticity equation is integrated exactly to obtain the solution. The numerical error is limited to that of the interpolation of the density model which is very small.

As far as the theory is valid there are no known restrictions. The number of model points and the jump discountinuities may be set by users.

Running time:
0.28 s for test run.