Things I learned yesterday and today.
Python doesn't understand quadruple-precision real numbers. f2py will happily compile Fortran code containing real(16)s into a form that Python can talk to, but Python will segfault as soon as it receives one of those real(16)s.
Also, Fortran's built-in fraction function doesn't do what I thought it did.