! A driver coupler code for MPH version 3
! Written by Yun (Helen) He and Chris Ding, NERSC/LBNL.

      program main
      use MPH_module
      implicit none
      include 'mpif.h'

      integer main_world, comm
      integer myProc, myproc_exe, myProc_global

      integer alpha
      real beta
      character*80 dynamics

      call MPH_debug (1)
      main_World = MPH_components (name1="ocean", name2="atmosphere",
     &                             name3="coupler")
      if (MPH_global_proc_id() == 0) call MPH_help ('on')

      call MPI_COMM_RANK (main_World, myProc_exe, ierr)
      call MPI_COMM_RANK (MPH_Global_World, myProc_global, ierr)
      write(*,*)'local_id=', myProc_exe, ', global_id=', myProc_global

      if (Proc_in_component("ocean", comm)) then
         call MPI_COMM_RANK (comm, myProc, ierr)
         write(*,*)'component=ocean, local_id=', myProc, ', exe_pid=',
     & myProc_exe, ', global_id=', myProc_global
         call MPH_get_argument("alpha", alpha)
         write(*,*)'alpha=', alpha
         call MPH_get_argument("beta", beta)
         write(*,*)'beta=', beta
      endif

      if (Proc_in_component("atmosphere", comm)) then
         call MPI_COMM_RANK (comm, myProc, ierr)
         write(*,*)'component=atmosphere, local_id=', myProc,
     & ', exe_pid=', myProc_exe, ', global_id=', myProc_global
         call MPH_get_argument("alpha", alpha)
         write(*,*)'alpha=', alpha
         call MPH_get_argument("beta", beta)
         write(*,*)'beta=', beta
      endif

      if (Proc_in_component("coupler", comm)) then
         call MPI_COMM_RANK (comm, myProc, ierr)
         write(*,*)'component=coupler, local_id=', myProc,
     & ', exe_pid=', myProc_exe, ', global_id=', myProc_global
         call MPH_get_argument("dynamics", dynamics)
         write(*,*)'dynamics=', dynamics
         call MPH_get_argument("beta", beta)
         write(*,*)'beta=', beta
      endif

      call MPI_FINALIZE (ierr)
      end program