Warning: this ratfor does not correspond to the fortran code.
The latter has been changed by ?? and no longer corresponds.
This cost TH one day of misery in June 2015