Changeset 869 for trunk/NEMO/LIM_SRC_3/limitd_th.F90

Timestamp:

2008-03-26T10:21:54+01:00 (16 years ago)

Author:

rblod

Message:

Parallelisation of LIM3. This commit seems to ensure the reproducibility mono/mpp. See ticket #77.

File:

• trunk/NEMO/LIM_SRC_3/limitd_th.F90 (modified) (7 diffs)

trunk/NEMO/LIM_SRC_3/limitd_th.F90

@@ -27 +27 @@
    USE limcons
    USE prtctl           ! Print control
+   USE lib_mpp 
  
    IMPLICIT NONE
@@ -1062 +1063 @@
          jl                 ! category index
 
-      LOGICAL ::   &  !:
+      INTEGER ::   &  !:
          zshiftflag          ! = .true. if ice must be shifted
 
@@ -1149 +1150 @@
       ! identify thicknesses that are too big
       !---------------------------------------
-         zshiftflag = .false.
+         zshiftflag = 0
 
          DO jj = 1, jpj 
             DO ji = 1, jpi 
                IF (a_i(ji,jj,jl) .GT. zeps .AND. ht_i(ji,jj,jl) .GT. hi_max(jl) ) THEN 
-                  zshiftflag        = .true.
+                  zshiftflag        = 1
                   zdonor(ji,jj,jl)  = jl 
                   zdaice(ji,jj,jl)  = a_i(ji,jj,jl)
@@ -1161 +1162 @@
             END DO                 ! ji
          END DO                 ! jj
-
-         IF (zshiftflag) THEN
+         IF( lk_mpp ) CALL mpp_max(zshiftflag)
+
+         IF ( zshiftflag == 1 ) THEN
 
       !------------------------------
@@ -1193 +1195 @@
       ! Identify thicknesses that are too small
       !-----------------------------------------
-         zshiftflag = .false.
+         zshiftflag = 0
 
          DO jj = 1, jpj
@@ -1200 +1202 @@
                   ht_i(ji,jj,jl+1) .LE. hi_max(jl)) THEN
 
-                  zshiftflag = .true.
+                  zshiftflag = 1
                   zdonor(ji,jj,jl) = jl + 1
                   zdaice(ji,jj,jl) = a_i(ji,jj,jl+1) 
@@ -1208 +1210 @@
          END DO                 ! jj
 
-         IF (zshiftflag) THEN
+         IF(lk_mpp) CALL mpp_max(zshiftflag)
+         IF (zshiftflag==1) THEN
 
       !------------------------------