The computer program described in this report calculates the height anomaly, gravity anomaly, gravity disturbance, and the two components of the deflection of the vertical using fully normalized potential coefficients of a spherical harmonic expansion. The program is designed to calculate these quantities on a point to point basis although certain calculations are not repeated if the latitude of the point does not change. The point input consists of the geodetic latitude, longitude, and height above the reference ellipsoid. Expansions up to degree 180 have been tested with the program. The report first describes the theory to be implemented. The program is described with a set of results for five sample points computed with three different potential coefficient fields to degree 180. The computer time for a single point is 0.46 seconds per point on the Amdahl 470 V/8.