Changeset 1000 for LMDZ4/trunk/libf/dyn3dpar/fluxstokenc_p.F

Timestamp:

Oct 6, 2008, 10:43:22 AM (16 years ago)

Author:

Laurent Fairhead

Message:

• Modifs sur le parallelisme: masquage dans la physique
• Inclusion strato
• mise en coherence etat0
• le mode offline fonctionne maintenant en parallele,
• les fichiers de la dynamiques sont correctement sortis et peuvent etre reconstruit avec rebuild
• la version parallele de la dynamique peut s'executer sans MPI (sur 1 proc)
• L'OPENMP fonctionne maintenant sans la parallelisation MPI.

YM
LF

File:

• LMDZ4/trunk/libf/dyn3dpar/fluxstokenc_p.F (modified) (6 diffs)

LMDZ4/trunk/libf/dyn3dpar/fluxstokenc_p.F

@@ -27 +27 @@
       REAL phis(ip1jmp1)
 
-      REAL pbaruc(ip1jmp1,llm),pbarvc(ip1jm,llm)
+      REAL,SAVE :: pbaruc(ip1jmp1,llm),pbarvc(ip1jm,llm)
       REAL massem(ip1jmp1,llm),tetac(ip1jmp1,llm),phic(ip1jmp1,llm)
 
@@ -55 +55 @@
 cym      wg(:,:)       = 0.
 
+c$OMP MASTER
+
       if(first) then
 
@@ -98 +100 @@
         pbaruc(ijb:ije,1:llm)=0
         
-        if (pole_sud) ije=ij_end-iip1
+        IF (pole_sud) ije=ij_end-iip1
         pbarvc(ijb:ije,1:llm)=0
       ENDIF
@@ -134 +136 @@
       iadvtr   = iadvtr+1
 
-
+c$OMP END MASTER
+c$OMP BARRIER
 c   Test pour savoir si on advecte a ce pas de temps
       IF ( iadvtr.EQ.istdyn ) THEN
+c$OMP MASTER
 c    normalisation
       ijb=ij_begin
@@ -162 +166 @@
 c     1. calcul de w
 c     2. groupement des mailles pres du pole.
-
+c$OMP END MASTER
+c$OMP BARRIER 
         call Register_Hallo(pbaruc,ip1jmp1,llm,1,1,1,1,Req)
         call Register_Hallo(pbarvc,ip1jm,llm,1,1,1,1,Req)
         call SendRequest(Req)
+c$OMP BARRIER
         call WaitRequest(Req)
-
+c$OMP BARRIER
+c$OMP MASTER
         CALL groupe_p( massem, pbaruc,pbarvc, pbarug,pbarvg,wg )
         
@@ -226 +233 @@
         
 C
-
+c$OMP END MASTER
       ENDIF ! if iadvtr.EQ.istdyn