Changeset 295 for codes/icosagcm/trunk/src/pression.f90

Timestamp:

10/31/14 14:52:01 (10 years ago)

Author:

ymipsl

Message:

Merging OpenMP parallisme mode : by subdomain and on vertical level.
This feature is actually experimental but may be retro-compatible with the last method based only on subdomain

YM

File:

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

codes/icosagcm/trunk/src/pression.f90

@@ -13 +13 @@
     INTEGER :: ind
 
+!$OMP BARRIER
     DO ind=1,ndomain
       IF (.NOT. assigned_domain(ind)) CYCLE
@@ -21 +22 @@
       CALL compute_pression(ps, p,0)
     ENDDO
+!$OMP BARRIER
   
   END SUBROUTINE pression
@@ -27 +29 @@
   USE icosa
   USE disvert_mod
+  USE omp_para
   IMPLICIT NONE
     REAL(rstd),INTENT(IN) :: ps(iim*jjm)
@@ -34 +37 @@
 
     IF(ap_bp_present) THEN
-    DO    l    = 1, llm+1
-      DO j=jj_begin-offset,jj_end+offset
-        DO i=ii_begin-offset,ii_end+offset
-          ij=(j-1)*iim+i
-          p(ij,l) = ap(l) + bp(l) * ps(ij)
+      DO    l    = ll_begin, ll_endp1
+!      DO    l    = 1, llm + 1
+        DO j=jj_begin-offset,jj_end+offset
+          DO i=ii_begin-offset,ii_end+offset
+            ij=(j-1)*iim+i
+            p(ij,l) = ap(l) + bp(l) * ps(ij)
+          ENDDO
         ENDDO
       ENDDO
-    ENDDO
     END IF
+
   END SUBROUTINE compute_pression