compilable using gfortran; use Make nompi to compile only no mpi fisource!!!