program main
  use mpi
  implicit none
  ! include 'mpif.h'
  integer :: ptot, rank, ierr

  call MPI_Init(ierr)
  call MPI_Comm_Size(mpi_comm_world,ptot,ierr)
  call MPI_Comm_Rank(mpi_comm_world,rank,ierr)

  print"(a,2i8)",'Hello World FORTRAN',rank,ptot

  call MPI_Finalize(ierr)
end program main