Migrated from use of 'INCLUDE mpif.h' to 'use mpi' as that is preferred for modern Fortran; Trimmed commented old code