Changeset 59

Timestamp:

07/31/12 18:39:34 (12 years ago)

Author:

ymipsl

Message:

Correct non compliant operation in MPI_Allreduce

YM

File:

• codes/icosagcm/trunk/src/dissip_gcm.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/dissip_gcm.f90

@@ -61 +61 @@
   REAL(rstd),POINTER    :: theta(:)
   REAL(rstd),POINTER    :: dtheta(:)
-  REAL(rstd)            :: dumax
-  REAL(rstd)            :: dthetamax
+  REAL(rstd)            :: dumax,dumax1
+  REAL(rstd)            :: dthetamax,dthetamax1
   REAL(rstd)            :: r
   REAL(rstd)            :: tau
@@ -197 +197 @@
       ENDDO
 
-      IF (using_mpi) CALL MPI_ALLREDUCE(dumax,dumax,1,MPI_REAL8,MPI_MAX,comm_icosa,ierr)
-
+      IF (using_mpi) THEN 
+                          CALL MPI_ALLREDUCE(dumax,dumax1,1,MPI_REAL8,MPI_MAX,comm_icosa,ierr)
+        dumax=dumax1
+      ENDIF
+                        
       DO ind=1,ndomain
         CALL swap_dimensions(ind)
@@ -267 +270 @@
       ENDDO
 
-      IF (using_mpi) CALL MPI_ALLREDUCE(dumax,dumax,1,MPI_REAL8,MPI_MAX,comm_icosa,ierr)
+      IF (using_mpi) THEN
+        CALL MPI_ALLREDUCE(dumax,dumax1,1,MPI_REAL8,MPI_MAX,comm_icosa,ierr)
+        dumax=dumax1
+      ENDIF
       
       DO ind=1,ndomain
@@ -333 +339 @@
         ENDDO
       ENDDO
-      IF (using_mpi) CALL MPI_ALLREDUCE(dthetamax,dthetamax,1,MPI_REAL8,MPI_MAX,comm_icosa,ierr)
+      IF (using_mpi) THEN
+        CALL MPI_ALLREDUCE(dthetamax,dthetamax1,1,MPI_REAL8,MPI_MAX,comm_icosa,ierr)
+        dthetamax=dthetamax1
+      ENDIF
       PRINT *,"divgrad : it :",it ,": dthetamax",dthetamax