• update Alan's FORTRAN code