Changeset 709 for LMDZ4/branches/V3_test/libf/dyn3dpar/groupe_p.F

Timestamp:

Sep 20, 2006, 12:12:39 PM (18 years ago)

Author:

Laurent Fairhead

Message:

Nouvelles versions de la dynamique YM
LF

File:

• LMDZ4/branches/V3_test/libf/dyn3dpar/groupe_p.F (modified) (7 diffs)

LMDZ4/branches/V3_test/libf/dyn3dpar/groupe_p.F

@@ -30 +30 @@
       real wm(iip1,jjp1,llm)
 
-      real zconvm(iip1,jjp1,llm),zconvmm(iip1,jjp1,llm)
+      real,save :: zconvm(iip1,jjp1,llm)
+      real,save :: zconvmm(iip1,jjp1,llm)
 
       real uu
@@ -38 +39 @@
       logical firstcall
       save firstcall
+c$OMP THREADPRIVATE(firstcall)
 
       data firstcall/.true./
@@ -57 +59 @@
       jjb=jj_begin
       jje=jj_end
-      zconvmm(:,jjb:jje,:)=zconvm(:,jjb:jje,:)
+
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
+      do l=1,llm
+        zconvmm(:,jjb:jje,l)=zconvm(:,jjb:jje,l)
+      enddo
+c$OMP END DO NOWAIT
+
       call groupeun_p(jjp1,llm,jjb,jje,zconvmm)
       
@@ -64 +72 @@
       if (pole_nord) jjb=jj_begin
       if (pole_sud)  jje=jj_end-1
-      pbarvm(:,jjb:jje,:)=pbarv(:,jjb:jje,:)
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
+      do l=1,llm
+        pbarvm(:,jjb:jje,l)=pbarv(:,jjb:jje,l)
+      enddo
+c$OMP END DO NOWAIT
+
       call groupeun_p(jjm,llm,jjb,jje,pbarvm)
 
@@ -74 +87 @@
       if (pole_sud)  jje=jj_end-1
       
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
       do l=1,llm
          do j=jjb,jje
@@ -86 +100 @@
          enddo
       enddo
-
+c$OMP END DO NOWAIT
 c    integration de la convergence de masse de haut  en bas ......
    
       jjb=jj_begin
       jje=jj_end
-      
-      do l=1,llm
-         do j=jjb,jje
-            do i=1,iip1
-               zconvmm(i,j,l)=zconvmm(i,j,l)
-            enddo
-         enddo
-      enddo
-      
+
+c$OMP BARRIER
+c$OMP MASTER      
       do  l = llm-1,1,-1
           do j=jjb,jje
@@ -107 +115 @@
           enddo
       enddo
-      
+
       if (.not. pole_sud) then
         zconvmm(:,jj_end+1,:)=0
-        wm(:,jj_end+1,:)=0
+cym     wm(:,jj_end+1,:)=0
       endif
+      
+c$OMP END MASTER
+c$OMP BARRIER      
+
       CALL vitvert_p(zconvmm(1,1,1),wm(1,1,1))