Added preliminary version of the file using parts of old FORTRAN code