modified: src/makefile

	modified:   src/parallel.f90

src/makefile

@@ -15,9 +15,9 @@ INC =
 # gcc:
 CMP = mpifort
 #O3: optimization
-#FLAGS   = -O3 -mcmodel=medium
+FLAGS   = -O3 -mcmodel=medium
 #debugging:
-FLAGS   = -g3 -O0 -fbounds-check -mcmodel=medium -fimplicit-none -fcheck=all -fbacktrace -floop-nest-optimize -ffpe-trap=invalid,zero,overflow -Wconversion -fno-tree-vectorize
+#FLAGS   = -g3 -O0 -fbounds-check -mcmodel=medium -fimplicit-none -fcheck=all -fbacktrace -floop-nest-optimize -ffpe-trap=invalid,zero,overflow -Wconversion -fno-tree-vectorize
 FFLAGS  = -c $(FLAGS)
 
 

src/parallel.f90

@@ -298,6 +298,7 @@
 !      
 !      myvar=myvar/nrprocs
 
+      syncSum=myvar
       if (myid.ne.0) then
         call mpi_isend(myvar,1,mpi_pr,0,0,mpi_comm_world,rs,mpierr)
         call mpi_wait(rs,mpistatus,mpierr)
@@ -307,7 +308,7 @@
       do 1 proc=1,nrprocs-1
         call mpi_irecv(myvar2,1,mpi_pr,proc,0,mpi_comm_world,rr,mpierr)
         call mpi_wait(rr,mpistatus,mpierr)
-        syncSum=myvar+myvar2
+        syncSum=syncSum+myvar2
 1     enddo
       endif
 
@@ -358,7 +359,7 @@
       do 1 proc=1,nrprocs-1
         call mpi_irecv(myvar2,1,mpi_integer,proc,0,mpi_comm_world,rr,mpierr)
         call mpi_wait(rr,mpistatus,mpierr)
-        syncSumI=myvar+myvar2
+        syncSumI=syncSumI+myvar2
 1     enddo
       endif
 
@@ -398,6 +399,7 @@
 !      
 !      myvar=myvar/nrprocs
 
+      syncSumI=myvar
       if (myid.ne.0) then
         call mpi_isend(myvar,1,mpi_pr,0,0,mpi_comm_world,rs,mpierr)
         call mpi_wait(rs,mpistatus,mpierr)
@@ -407,7 +409,7 @@
       do 1 proc=1,nrprocs-1
         call mpi_irecv(myvar2,1,mpi_pr,proc,0,mpi_comm_world,rr,mpierr)
         call mpi_wait(rr,mpistatus,mpierr)
-        syncAv=myvar+myvar2
+        syncAv=syncAv+myvar2
 1     enddo
       syncAv=syncAv/nrprocs
       endif