Original Fortran g77 code