COSMOS v7.655  COSMOSv7655
(AirShowerMC)
inc_gatherNrfai.f
Go to the documentation of this file.
1 ! we note the last index may be actual max used
2 !
3  do k = 1, ansites
4  do j = 1, 4
5  do l = 1,nfai
6  do i = 1, nrbin
7  nrfaireca(i, l, j, k)=0.
8  enddo
9  enddo
10  enddo
11  enddo
12  do k = 1, maxnoofsites
13  do j = 1, 4
14  do l = 1,nfai
15  do i = 1, nrbin
16  nrfaialla(i, l, j, k)=0.
17  enddo
18  enddo
19  enddo
20  enddo
21  do k = 1, maxnoofsites
22  do l = 1,nfai
23  do i = 1, nrbin
24  derfaia(i, l, k) = 0.
25  enddo
26  enddo
27  enddo
28 ! if(mpirank .eq. 0 ) then
29  call mpi_reduce(nrfairec, nrfaireca, ansites*4*nfai*nrbin,
30  * mpi_real, mpi_sum, 0, mpi_comm_world, icon)
31 ! call mpi_reduce(MPI_IN_PLACE, nrfairec, ansites*4*nfai*nrbin,
32 ! * MPI_REAL, MPI_SUM, 0, MPI_COMM_WORLD, icon)
33 
34  call mpi_reduce(nrfaiall, nrfaialla, maxnoofsites*4*nfai*nrbin,
35  * mpi_real, mpi_sum, 0, mpi_comm_world, icon)
36 
37 ! call mpi_reduce(MPI_IN_PLACE, nrfaiAll,
38 ! * MaxNoOfSites*4*nfai*nrbin,
39 ! * MPI_REAL, MPI_SUM, 0, MPI_COMM_WORLD, icon)
40 
41  call mpi_reduce(derfai, derfaia, maxnoofsites*nfai*nrbin,
42  * mpi_real, mpi_sum, 0, mpi_comm_world, icon)
43 ! call mpi_reduce(MPI_IN_PLACE, dErfai, MaxNoOfSites*nfai*nrbin,
44 ! * MPI_REAL, MPI_SUM, 0, MPI_COMM_WORLD, icon)
45 ! else
46 
47 ! call mpi_reduce(nrfaiAll, nrfaiAll, MaxNoOfSites*4*nfai*nrbin,
48 ! * MPI_REAL, MPI_SUM, 0, MPI_COMM_WORLD, icon)
49 
50 ! call mpi_reduce(dErfai, dErfai, MaxNoOfSites*nfai*nrbin,
51 ! * MPI_REAL, MPI_SUM, 0, MPI_COMM_WORLD, icon)
52 
53 ! endif
maxnoofsites
const int maxnoofsites
Definition: Zobs.h:7
i
nodes i
Definition: cstdatmos0.fNew2.h:93
nrbin
! timing nrbin
Definition: Zprivate2.h:12
j
********************block data cblkHeavy ********************integer j data *HeavyG2symbol *data *HeavyG2code kiron data j
Definition: cblkHeavy.h:36
nfai
! timing nfai
Definition: Zprivate2.h:12