General Hartree-Fock program. Nature of problem: Numerical non-relativistic Hartree-Fock wavefunctions are determined for atoms in a bound state. These are used to predict a variety of atomic parameters. Solution method: The self-consistent field method of solution is employed. Two methods of solution of the differential equations are used, each improving only a single function at a time (Comp. Phys. Commun. 15(1977)221), and rotations are introduced explicitly to find an energy stationary with respect to a rotation of the orbital basis. Otherwise the procedures are essentially the same as those in MCHF77 (Comp. Phys. Commun. 14(1978)145).