Piece of classic FORTRAN code