gfortran and mpif77 version of Traub et al 2005 FORTRAN code.