A couple of issues came up when trying to run with Gfortran:
- the subroutine trace1 has to be commented out. See r5107.
- could not run with M7. The code crashes with a segmentation fault in m7_cumor.