https://9pdf.info/document/z1dl5w3v-fortran-code-subroutine-published-thesis-a.html