Changeset 295 for codes/icosagcm/trunk/src/vertical_interp.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/vertical_interp.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/vertical_interp.f90

@@ -21 +21 @@
   USE icosa
   USE pression_mod
+  USE omp_para
   IMPLICIT NONE
     TYPE(t_field),POINTER :: f_ps(:)
@@ -48 +49 @@
 
   SUBROUTINE compute_vertical_interp(p,in,out,pval)
+  USE omp_para
   IMPLICIT NONE
     REAL(rstd),INTENT(IN) :: p(iim*jjm,llm+1)
@@ -56 +58 @@
     INTEGER :: i,j,ij,l
         
-    DO j=jj_begin-1,jj_end+1
-      DO i=ii_begin-1,ii_end+1
-        ij=(j-1)*iim+i
-        l=llm-1
-        DO WHILE(0.5*(p(ij,l)+p(ij,l+1))<pval .AND. l>1)
-          l=l-1
+!$OMP BARRIER    
+    IF (is_omp_level_master) THEN
+    
+      DO j=jj_begin-1,jj_end+1
+        DO i=ii_begin-1,ii_end+1
+          ij=(j-1)*iim+i
+          l=llm-1
+          DO WHILE(0.5*(p(ij,l)+p(ij,l+1))<pval .AND. l>1)
+            l=l-1
+          ENDDO
+          pmid=0.5*(p(ij,l)+p(ij,l+1))
+          pmidp1=0.5*(p(ij,l+1)+p(ij,l+2))
+
+          coeff=(pval-pmid)/(pmid-pmidp1)
+        
+          out(ij)=in(ij,l)+coeff*(in(ij,l)-in(ij,l+1))
         ENDDO
-        pmid=0.5*(p(ij,l)+p(ij,l+1))
-        pmidp1=0.5*(p(ij,l+1)+p(ij,l+2))
-
-        coeff=(pval-pmid)/(pmid-pmidp1)
-        
-        out(ij)=in(ij,l)+coeff*(in(ij,l)-in(ij,l+1))
       ENDDO
-    ENDDO
+   
+    ENDIF
+!$OMP BARRIER    
         
   END SUBROUTINE compute_vertical_interp