include 'lyu2.f' !the main program include 'GEOPACK.f' !package for common subroutines include 'EXTERN.f' !to call T89M.f include 'T89M.f' !the tsy89 magnetic field model