RATPOE values (ratio of gravitational to total energy) are
uniform density values 3GM/5Rc**2 for Mercury, Venus, Mars,
Pluto, and Moon.  Earth is from PRL 36, March 15, 1976, p555.
The gas giants are based on exponential density profiles with
scale height set to give the right moment of inertia (from Allen).

for nshp=1,shape model is two dimensional
fourier series.(see comments in bodred)
internal controls necessary are:
tlat(1)=lower latitude on planet
tlat(2)=upper latitude on planet
nzone=20
ntess=10
nzone and ntess are needed for shppar to calc. four. partials

nshp=2, local shape model - altitude grid
grid points for appropriate range on planet represent
height in kilometers of the surface at the
lat and longitude corresponding to the grid point
(see bodred)
set defaults for lat & long range, grid spacing,
grid dimensions and total # of grid points