code now has STRICT compliance with Fortran77