cmanager.f

Go to the documentation of this file.

#if defined (KEKA) || defined (KEKB)
#define DOMPI
#endif
!     ****************************************************************
!     * cmanager: Managing Cosmos Simulation
!     ****************************************************************
      subroutine cmanager
      implicit none
#include "Zevhnp.h"


#if defined (DOMPI)
#include "mpif.h"
#include "Zmaxdef.h"
#include "Ztrack.h"
#include "Zmpi.h"
#include "Zmanagerp.h"

      integer err, intdata
      character*120 parampath  ! absolute param path given by job script 
                      ! via envrionment

      call mpi_init(err)
      call mpi_comm_size(mpi_comm_world, mpisize, err)
      call mpi_comm_rank(mpi_comm_world, mpirank, err)
!        this is to treat exising name list as it is.
!         no delim="apstrophe" is needed
#if defined (KEKB) || defined (KEKA)      
      call setrteopts("namelist=old")
#endif
!!!      call kgetenv2("PARAMPATH", parampath) ! this is OK ??
      err = kgetenv2("PARAMPATH", parampath)
!           to avoid simultaneous access to HD, we
!           read namelist from rank0, rank1,... 
      if(mpirank .eq. 0) then
         write(0,*) ' mpisize=',mpisize
         open(tempdev, file=parampath, action="read")
         call creadparam(tempdev)
         close(tempdev)
         if( mpisize .gt. 1) then
            call mpi_send(mpirank, 1, mpi_integer, 1, 1,
     *       mpi_comm_world, err)
         endif
      else
         call mpi_recv(intdata, 1, mpi_integer, mpirank-1, 1,
     *      mpi_comm_world, mpistat, err)
         open(tempdev, file=parampath, action="read")
         call creadparam(tempdev)
         close(tempdev)
!              mpisize=5 and mpirank=4  ==> no more
!                            mpirank=3  ==> send
         if( mpirank .lt. mpisize-1 ) then
            call mpi_send(mpirank, 1, mpi_integer, mpirank+1, 1,
     *       mpi_comm_world, err)
         endif
      endif
!      write(0,*)' rank=', mpirank, ' Job =', Job
!//////////////
#else
#if defined (MacIFC)
      write(0,*) ' '
#endif
      call creadparam(5)        ! read execution conditions from stdin
#endif
#ifdef NEXT486
      if( index(intmodel, 'fritiof1.6') .gt. 0 ) then
         call cerrormsg('fritiof1.6 cannot be used'//
     *  ' for NEXT-Absoft Fortran', 0)
      endif
#endif
      call cbeginrun            ! initialize the simulation
      call ceventloop           ! begin simulation and enter  event Loop
      call cendrun              ! close the simulation
      end