bcast_inputs.f:      call MPI_BCAST(ipr, 1, integer_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f:      call MPI_BCAST(inorm, 1, integer_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f:      call MPI_BCAST(itmax, 1, integer_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f:      call MPI_BCAST(dt, 1, dp_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f:      call MPI_BCAST(omega, 1, dp_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f:      call MPI_BCAST(tolrsd, 5, dp_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f:      call MPI_BCAST(nx0, 1, dp_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f:      call MPI_BCAST(ny0, 1, dp_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f:      call MPI_BCAST(nz0, 1, dp_type, root, MPI_COMM_WORLD, ierr)
bcast_inputs.f-
bcast_inputs.f-      return
bcast_inputs.f-      end
bcast_inputs.f-
--
error.f:      call MPI_ALLREDUCE( dummy,
error.f-     >                    errnm,
error.f-     >                    5,
error.f-     >                    dp_type,
error.f:     >                    MPI_SUM,
error.f:     >                    MPI_COMM_WORLD,
error.f-     >                    IERROR )
error.f-
error.f-      do m = 1, 5
error.f-         errnm(m) = sqrt ( errnm(m) / ( (nx0-2)*(ny0-2)*(nz0-2) ) )
--
exchange_1.f:      integer STATUS(MPI_STATUS_SIZE)
exchange_1.f-      integer IERROR
exchange_1.f-
exchange_1.f-
exchange_1.f-
--
exchange_1.f:              call MPI_RECV( dum1(1,jst),
exchange_1.f-     >                       5*(jend-jst+1),
exchange_1.f-     >                       dp_type,
exchange_1.f-     >                       north,
exchange_1.f-     >                       from_n,
exchange_1.f:     >                       MPI_COMM_WORLD,
exchange_1.f-     >                       status,
exchange_1.f-     >                       IERROR )
exchange_1.f-              do j=jst,jend
exchange_1.f-                  g(1,0,j,k) = dum1(1,j)
--
exchange_1.f:              call MPI_RECV( dum1(1,ist),
exchange_1.f-     >                       5*(iend-ist+1),
exchange_1.f-     >                       dp_type,
exchange_1.f-     >                       west,
exchange_1.f-     >                       from_w,
exchange_1.f:     >                       MPI_COMM_WORLD,
exchange_1.f-     >                       status,
exchange_1.f-     >                       IERROR )
exchange_1.f-              do i=ist,iend
exchange_1.f-                  g(1,i,0,k) = dum1(1,i)
--
exchange_1.f:              call MPI_RECV( dum1(1,jst),
exchange_1.f-     >                       5*(jend-jst+1),
exchange_1.f-     >                       dp_type,
exchange_1.f-     >                       south,
exchange_1.f-     >                       from_s,
exchange_1.f:     >                       MPI_COMM_WORLD,
exchange_1.f-     >                       status,
exchange_1.f-     >                       IERROR )
exchange_1.f-              do j=jst,jend
exchange_1.f-                  g(1,nx+1,j,k) = dum1(1,j)
--
exchange_1.f:              call MPI_RECV( dum1(1,ist),
exchange_1.f-     >                       5*(iend-ist+1),
exchange_1.f-     >                       dp_type,
exchange_1.f-     >                       east,
exchange_1.f-     >                       from_e,
exchange_1.f:     >                       MPI_COMM_WORLD,
exchange_1.f-     >                       status,
exchange_1.f-     >                       IERROR )
exchange_1.f-              do i=ist,iend
exchange_1.f-                  g(1,i,ny+1,k) = dum1(1,i)
--
exchange_1.f:              call MPI_SEND( dum(1,jst), 
exchange_1.f-     >                       5*(jend-jst+1), 
exchange_1.f-     >                       dp_type, 
exchange_1.f-     >                       south, 
exchange_1.f-     >                       from_n, 
exchange_1.f:     >                       MPI_COMM_WORLD, 
exchange_1.f-     >                       IERROR )
exchange_1.f-          endif
exchange_1.f-
exchange_1.f-          if( east .ne. -1 ) then
--
exchange_1.f:              call MPI_SEND( dum(1,ist), 
exchange_1.f-     >                       5*(iend-ist+1), 
exchange_1.f-     >                       dp_type, 
exchange_1.f-     >                       east, 
exchange_1.f-     >                       from_w, 
exchange_1.f:     >                       MPI_COMM_WORLD, 
exchange_1.f-     >                       IERROR )
exchange_1.f-          endif
exchange_1.f-
exchange_1.f-      else
--
exchange_1.f:              call MPI_SEND( dum(1,jst), 
exchange_1.f-     >                       5*(jend-jst+1), 
exchange_1.f-     >                       dp_type, 
exchange_1.f-     >                       north, 
exchange_1.f-     >                       from_s, 
exchange_1.f:     >                       MPI_COMM_WORLD, 
exchange_1.f-     >                       IERROR )
exchange_1.f-          endif
exchange_1.f-
exchange_1.f-          if( west .ne. -1 ) then
--
exchange_1.f:              call MPI_SEND( dum(1,ist), 
exchange_1.f-     >                       5*(iend-ist+1), 
exchange_1.f-     >                       dp_type, 
exchange_1.f-     >                       west, 
exchange_1.f-     >                       from_e, 
exchange_1.f:     >                       MPI_COMM_WORLD, 
exchange_1.f-     >                       IERROR )
exchange_1.f-          endif
exchange_1.f-
exchange_1.f-      endif
--
exchange_3.f:      integer STATUS(MPI_STATUS_SIZE)
exchange_3.f-      integer IERROR
exchange_3.f-
exchange_3.f-
exchange_3.f-
--
exchange_3.f:          call MPI_IRECV( buf1,
exchange_3.f-     >                    10*ny*nz,
exchange_3.f-     >                    dp_type,
exchange_3.f:     >                    MPI_ANY_SOURCE,
exchange_3.f-     >                    from_n,
exchange_3.f:     >                    MPI_COMM_WORLD,
exchange_3.f-     >                    mid,
exchange_3.f-     >                    IERROR )
exchange_3.f-      end if
exchange_3.f-
--
exchange_3.f:          call MPI_SEND( buf,
exchange_3.f-     >                   10*ny*nz,
exchange_3.f-     >                   dp_type,
exchange_3.f-     >                   south,
exchange_3.f-     >                   from_n,
exchange_3.f:     >                   MPI_COMM_WORLD,
exchange_3.f-     >                   IERROR )
exchange_3.f-        end if
exchange_3.f-
exchange_3.f-c---------------------------------------------------------------------
--
exchange_3.f:          call MPI_WAIT( mid, STATUS, IERROR )
exchange_3.f-
exchange_3.f-          do k = 1,nz
exchange_3.f-            do j = 1,ny
exchange_3.f-              ipos1 = (k-1)*ny + j
--
exchange_3.f:          call MPI_IRECV( buf1,
exchange_3.f-     >                    10*ny*nz,
exchange_3.f-     >                    dp_type,
exchange_3.f:     >                    MPI_ANY_SOURCE,
exchange_3.f-     >                    from_s,
exchange_3.f:     >                    MPI_COMM_WORLD,
exchange_3.f-     >                    mid,
exchange_3.f-     >                    IERROR )
exchange_3.f-      end if
exchange_3.f-
--
exchange_3.f:          call MPI_SEND( buf,
exchange_3.f-     >                   10*ny*nz,
exchange_3.f-     >                   dp_type,
exchange_3.f-     >                   north,
exchange_3.f-     >                   from_s,
exchange_3.f:     >                   MPI_COMM_WORLD,
exchange_3.f-     >                   IERROR )
exchange_3.f-        end if
exchange_3.f-
exchange_3.f-c---------------------------------------------------------------------
--
exchange_3.f:          call MPI_WAIT( mid, STATUS, IERROR )
exchange_3.f-
exchange_3.f-          do k = 1,nz
exchange_3.f-            do j = 1,ny
exchange_3.f-              ipos1 = (k-1)*ny + j
--
exchange_3.f:          call MPI_IRECV( buf1,
exchange_3.f-     >                    10*nx*nz,
exchange_3.f-     >                    dp_type,
exchange_3.f:     >                    MPI_ANY_SOURCE,
exchange_3.f-     >                    from_w,
exchange_3.f:     >                    MPI_COMM_WORLD,
exchange_3.f-     >                    mid,
exchange_3.f-     >                    IERROR )
exchange_3.f-      end if
exchange_3.f-
--
exchange_3.f:          call MPI_SEND( buf,
exchange_3.f-     >                   10*nx*nz,
exchange_3.f-     >                   dp_type,
exchange_3.f-     >                   east,
exchange_3.f-     >                   from_w,
exchange_3.f:     >                   MPI_COMM_WORLD,
exchange_3.f-     >                   IERROR )
exchange_3.f-        end if
exchange_3.f-
exchange_3.f-c---------------------------------------------------------------------
--
exchange_3.f:          call MPI_WAIT( mid, STATUS, IERROR )
exchange_3.f-
exchange_3.f-          do k = 1,nz
exchange_3.f-            do i = 1,nx
exchange_3.f-              ipos1 = (k-1)*nx + i
--
exchange_3.f:          call MPI_IRECV( buf1,
exchange_3.f-     >                    10*nx*nz,
exchange_3.f-     >                    dp_type,
exchange_3.f:     >                    MPI_ANY_SOURCE,
exchange_3.f-     >                    from_e,
exchange_3.f:     >                    MPI_COMM_WORLD,
exchange_3.f-     >                    mid,
exchange_3.f-     >                    IERROR )
exchange_3.f-      end if
exchange_3.f-
--
exchange_3.f:          call MPI_SEND( buf,
exchange_3.f-     >                   10*nx*nz,
exchange_3.f-     >                   dp_type,
exchange_3.f-     >                   west,
exchange_3.f-     >                   from_e,
exchange_3.f:     >                   MPI_COMM_WORLD,
exchange_3.f-     >                   IERROR )
exchange_3.f-        end if
exchange_3.f-
exchange_3.f-c---------------------------------------------------------------------
--
exchange_3.f:          call MPI_WAIT( mid, STATUS, IERROR )
exchange_3.f-
exchange_3.f-          do k = 1,nz
exchange_3.f-            do i = 1,nx
exchange_3.f-              ipos1 = (k-1)*nx + i
--
exchange_4.f:      integer STATUS(MPI_STATUS_SIZE)
exchange_4.f-      integer IERROR
exchange_4.f-
exchange_4.f-
exchange_4.f-
--
exchange_4.f:        call MPI_IRECV( dum,
exchange_4.f-     >                  2*nx,
exchange_4.f-     >                  dp_type,
exchange_4.f:     >                  MPI_ANY_SOURCE,
exchange_4.f-     >                  from_e,
exchange_4.f:     >                  MPI_COMM_WORLD,
exchange_4.f-     >                  msgid3,
exchange_4.f-     >                  IERROR )
exchange_4.f-
exchange_4.f:        call MPI_WAIT( msgid3, STATUS, IERROR )
exchange_4.f-
exchange_4.f-        do i = 1,nx
exchange_4.f-          g(i,ny+1) = dum(i)
exchange_4.f-          h(i,ny+1) = dum(i+nx)
--
exchange_4.f:        call MPI_SEND( dum,
exchange_4.f-     >                 2*nx,
exchange_4.f-     >                 dp_type,
exchange_4.f-     >                 west,
exchange_4.f-     >                 from_e,
exchange_4.f:     >                 MPI_COMM_WORLD,
exchange_4.f-     >                 IERROR )
exchange_4.f-
exchange_4.f-      end if
exchange_4.f-
--
exchange_4.f:        call MPI_IRECV( dum,
exchange_4.f-     >                  2*ny2,
exchange_4.f-     >                  dp_type,
exchange_4.f:     >                  MPI_ANY_SOURCE,
exchange_4.f-     >                  from_s,
exchange_4.f:     >                  MPI_COMM_WORLD,
exchange_4.f-     >                  msgid1,
exchange_4.f-     >                  IERROR )
exchange_4.f-
exchange_4.f:        call MPI_WAIT( msgid1, STATUS, IERROR )
exchange_4.f-
exchange_4.f-        do j = 0,ny+1
exchange_4.f-          g(nx+1,j) = dum(j+1)
exchange_4.f-          h(nx+1,j) = dum(j+ny2+1)
--
exchange_4.f:        call MPI_SEND( dum,
exchange_4.f-     >                 2*ny2,
exchange_4.f-     >                 dp_type,
exchange_4.f-     >                 north,
exchange_4.f-     >                 from_s,
exchange_4.f:     >                 MPI_COMM_WORLD,
exchange_4.f-     >                 IERROR )
exchange_4.f-
exchange_4.f-      end if
exchange_4.f-
--
exchange_5.f:      integer STATUS(MPI_STATUS_SIZE)
exchange_5.f-      integer IERROR
exchange_5.f-
exchange_5.f-
exchange_5.f-
--
exchange_5.f:        call MPI_IRECV( dum,
exchange_5.f-     >                  nz,
exchange_5.f-     >                  dp_type,
exchange_5.f:     >                  MPI_ANY_SOURCE,
exchange_5.f-     >                  from_s,
exchange_5.f:     >                  MPI_COMM_WORLD,
exchange_5.f-     >                  msgid1,
exchange_5.f-     >                  IERROR )
exchange_5.f-
exchange_5.f:        call MPI_WAIT( msgid1, STATUS, IERROR )
exchange_5.f-
exchange_5.f-        do k = 1,nz
exchange_5.f-          g(nx+1,k) = dum(k)
exchange_5.f-        end do
--
exchange_5.f:        call MPI_SEND( dum,
exchange_5.f-     >                 nz,
exchange_5.f-     >                 dp_type,
exchange_5.f-     >                 north,
exchange_5.f-     >                 from_s,
exchange_5.f:     >                 MPI_COMM_WORLD,
exchange_5.f-     >                 IERROR )
exchange_5.f-
exchange_5.f-      end if
exchange_5.f-
--
exchange_6.f:      integer STATUS(MPI_STATUS_SIZE)
exchange_6.f-      integer IERROR
exchange_6.f-
exchange_6.f-
exchange_6.f-
--
exchange_6.f:        call MPI_IRECV( dum,
exchange_6.f-     >                  nz,
exchange_6.f-     >                  dp_type,
exchange_6.f:     >                  MPI_ANY_SOURCE,
exchange_6.f-     >                  from_e,
exchange_6.f:     >                  MPI_COMM_WORLD,
exchange_6.f-     >                  msgid3,
exchange_6.f-     >                  IERROR )
exchange_6.f-
exchange_6.f:        call MPI_WAIT( msgid3, STATUS, IERROR )
exchange_6.f-
exchange_6.f-        do k = 1,nz
exchange_6.f-          g(ny+1,k) = dum(k)
exchange_6.f-        end do
--
exchange_6.f:        call MPI_SEND( dum,
exchange_6.f-     >                 nz,
exchange_6.f-     >                 dp_type,
exchange_6.f-     >                 west,
exchange_6.f-     >                 from_e,
exchange_6.f:     >                 MPI_COMM_WORLD,
exchange_6.f-     >                 IERROR )
exchange_6.f-
exchange_6.f-      end if
exchange_6.f-
--
init_comm.f:      call MPI_INIT( IERROR )
init_comm.f-
init_comm.f-c---------------------------------------------------------------------
init_comm.f-c   establish the global rank of this process
init_comm.f-c---------------------------------------------------------------------
init_comm.f:      call MPI_COMM_RANK( MPI_COMM_WORLD,
init_comm.f-     >                     id,
init_comm.f-     >                     IERROR )
init_comm.f-
init_comm.f-c---------------------------------------------------------------------
--
init_comm.f:      call MPI_COMM_SIZE( MPI_COMM_WORLD,
init_comm.f-     >                     num,
init_comm.f-     >                     IERROR )
init_comm.f-
init_comm.f-      ndim   = nodedim(num)
--
l2norm.f:      call MPI_ALLREDUCE( dummy,
l2norm.f-     >                    sum,
l2norm.f-     >                    5,
l2norm.f-     >                    dp_type,
l2norm.f:     >                    MPI_SUM,
l2norm.f:     >                    MPI_COMM_WORLD,
l2norm.f-     >                    IERROR )
l2norm.f-
l2norm.f-      do m = 1, 5
l2norm.f-         sum(m) = sqrt ( sum(m) / ( (nx0-2)*(ny0-2)*(nz0-2) ) )
--
lu.f:      call mpi_finalize(ierr)
lu.f-      end
lu.f-
lu.f-
--
mpifrag.f:        dp_type = MPI_DOUBLE_PRECISION
mpifrag.f:        integer_type = MPI_INTEGER
mpifrag.f:        dc_type = MPI_DOUBLE_COMPLEX
--
pintgr.f:      call MPI_ALLREDUCE( dummy,
pintgr.f-     >                    frc1,
pintgr.f-     >                    1,
pintgr.f-     >                    dp_type,
pintgr.f:     >                    MPI_SUM,
pintgr.f:     >                    MPI_COMM_WORLD,
pintgr.f-     >                    IERROR )
pintgr.f-
pintgr.f-      frc1 = dxi * deta * frc1
pintgr.f-
--
pintgr.f:      call MPI_ALLREDUCE( dummy,
pintgr.f-     >                    frc2,
pintgr.f-     >                    1,
pintgr.f-     >                    dp_type,
pintgr.f:     >                    MPI_SUM,
pintgr.f:     >                    MPI_COMM_WORLD,
pintgr.f-     >                    IERROR )
pintgr.f-
pintgr.f-      frc2 = dxi * dzeta * frc2
pintgr.f-
--
pintgr.f:      call MPI_ALLREDUCE( dummy,
pintgr.f-     >                    frc3,
pintgr.f-     >                    1,
pintgr.f-     >                    dp_type,
pintgr.f:     >                    MPI_SUM,
pintgr.f:     >                    MPI_COMM_WORLD,
pintgr.f-     >                    IERROR )
pintgr.f-
pintgr.f-      frc3 = deta * dzeta * frc3
pintgr.f-      frc = 0.25d+00 * ( frc1 + frc2 + frc3 )
--
read_input.f:         call MPI_COMM_SIZE(MPI_COMM_WORLD, nnodes, ierror)
read_input.f-         if (nnodes .ne. nnodes_compiled) then
read_input.f-            write (*, 2000) nnodes, nnodes_compiled
read_input.f- 2000       format (5x,'Warning: program is running on',i3,' processors'
read_input.f-     >             /5x,'but was compiled for ', i3)
--
read_input.f:            CALL MPI_ABORT( MPI_COMM_WORLD, MPI_ERR_OTHER, IERROR )
read_input.f-
read_input.f-         end if
read_input.f-
read_input.f-         if ( ( nx0 .gt. isiz01 ) .or.
--
read_input.f:            CALL MPI_ABORT( MPI_COMM_WORLD, MPI_ERR_OTHER, IERROR )
read_input.f-
read_input.f-         end if
read_input.f-
read_input.f-
--
ssor.f:      call MPI_BARRIER( MPI_COMM_WORLD, IERROR )
ssor.f- 
ssor.f:      wtime = MPI_WTIME()
ssor.f- 
ssor.f-c---------------------------------------------------------------------
ssor.f-c   the timestep loop
ssor.f-c---------------------------------------------------------------------
--
ssor.f:      wtime = MPI_WTIME() - wtime
ssor.f- 
ssor.f:      call MPI_ALLREDUCE( wtime, 
ssor.f-     >                    maxtime, 
ssor.f-     >                    1, 
ssor.f:     >                    MPI_DOUBLE_PRECISION, 
ssor.f:     >                    MPI_MAX, 
ssor.f:     >                    MPI_COMM_WORLD,
ssor.f-     >                    IERROR )
ssor.f- 
ssor.f-
ssor.f-
--
subdomain.f:          CALL MPI_ABORT( MPI_COMM_WORLD,
subdomain.f-     >                    ERRORCODE,
subdomain.f-     >                    IERROR )
subdomain.f-      end if
subdomain.f-
--
subdomain.f:          CALL MPI_ABORT( MPI_COMM_WORLD,
subdomain.f-     >                    ERRORCODE,
subdomain.f-     >                    IERROR )
subdomain.f-      end if
subdomain.f-
