I ran it with your test files and it works, as you said.
However, it still fails when I try and use my fortran file - the file has worked with all my previous runs. Testing it in part, the run fails when I use the CALCOT routine to distribute layers and when I use SOURCE to apply head losses at certain layers.
What's odd is that even if I remove code to something simple .e.g
Which makes me think it's one of the numerical options in the *.cas file. I played around with the preconditionings but haven't had much luck. I additionally tested it with refined banks and had the same issue.