DAVINT crash in ifort12
Occurs with O2 flag. Resolved with O1 flag.
Reason still to be investigated. Compiler bug?
Manifests itself as spurious modification of kthr array in mod_saturation. Can be bizarrely resolved by seemingly random interventions, e.g., adding write statements (for debugging purposes) within the subroutine