Fortran routines in the directory taken from cmlib (public domain)
