Changeset 764 for LMDZ4/trunk/libf/dyn3dpar/integrd_p.F

Timestamp:

Jun 4, 2007, 4:13:10 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Merge entre la version V3_conv et le HEAD
YM, JG, LF

File:

• LMDZ4/trunk/libf/dyn3dpar/integrd_p.F (modified) (15 diffs)

LMDZ4/trunk/libf/dyn3dpar/integrd_p.F

@@ -4 +4 @@
       SUBROUTINE integrd_p
      $  (  nq,vcovm1,ucovm1,tetam1,psm1,massem1,
-     $     dv,du,dteta,dq,dp,vcov,ucov,teta,q,ps,masse,phis,finvmaold )
+     $     dv,du,dteta,dq,dp,vcov,ucov,teta,q,ps0,masse,phis,finvmaold)
       USE parallel
       IMPLICIT NONE
@@ -41 +41 @@
       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
       REAL q(ip1jmp1,llm,nq)
-      REAL ps(ip1jmp1),masse(ip1jmp1,llm),phis(ip1jmp1)
+      REAL ps0(ip1jmp1),masse(ip1jmp1,llm),phis(ip1jmp1)
 
       REAL vcovm1(ip1jm,llm),ucovm1(ip1jmp1,llm)
@@ -55 +55 @@
       REAL vscr( ip1jm ),uscr( ip1jmp1 ),hscr( ip1jmp1 ),pscr(ip1jmp1)
       REAL massescr( ip1jmp1,llm ), finvmasse(ip1jmp1,llm)
-      REAL p(ip1jmp1,llmp1)
+      REAL,SAVE :: p(ip1jmp1,llmp1)
       REAL tpn,tps,tppn(iim),tpps(iim)
       REAL qpn,qps,qppn(iim),qpps(iim)
@@ -62 +62 @@
       INTEGER  l,ij,iq
 
-      EXTERNAL  filtreg,massdair,pression
-      EXTERNAL  SCOPY 
       REAL SSUM
       EXTERNAL SSUM
       INTEGER ijb,ije,jjb,jje
+      REAL,SAVE :: ps(ip1jmp1)
 c-----------------------------------------------------------------------
+     
       if (pole_nord) THEN
-      
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
         DO  l = 1,llm
           DO  ij = 1,iip1
@@ -76 +76 @@
            ENDDO
         ENDDO
-      
+c$OMP END DO NOWAIT        
       ENDIF
 
       if (pole_sud) THEN
-      
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
         DO  l = 1,llm
           DO  ij = 1,iip1
@@ -87 +87 @@
           ENDDO
         ENDDO
-      
+c$OMP END DO NOWAIT      
       ENDIF
 
@@ -96 +96 @@
       ijb=ij_begin
       ije=ij_end
-      massescr(ijb:ije,:)=masse(ijb:ije,:)
-      
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 
+      DO  l = 1,llm
+        massescr(ijb:ije,l)=masse(ijb:ije,l)
+      ENDDO
+c$OMP END DO NOWAIT 
+
+c$OMP MASTER      
       DO 2 ij = ijb,ije
-       pscr (ij)    = ps(ij)
+       pscr (ij)    = ps0(ij)
        ps (ij)      = psm1(ij) + dt * dp(ij)
    2  CONTINUE
@@ -133 +138 @@
       
       ENDIF
+c$OMP END MASTER
+c$OMP BARRIER
 c
 c  ... Calcul  de la nouvelle masse d'air au dernier temps integre t+1 ...
 c
+
       CALL pression_p ( ip1jmp1, ap, bp, ps, p )
+c$OMP BARRIER
       CALL massdair_p (     p  , masse         )
 
@@ -142 +151 @@
       ijb=ij_begin
       ije=ij_end
-      finvmasse(ijb:ije,:)=masse(ijb:ije,:)
+      
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 
+      DO  l = 1,llm
+        finvmasse(ijb:ije,l)=masse(ijb:ije,l)
+      ENDDO
+c$OMP END DO NOWAIT
 
       jjb=jj_begin
@@ -151 +165 @@
 c    ............   integration  de  ucov, vcov,  h     ..............
 
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 
       DO 10 l = 1,llm
       
@@ -225 +240 @@
 
   10  CONTINUE
-
+c$OMP END DO NOWAIT
 
 c
@@ -233 +248 @@
       ije=ij_end
       
-
+c$OMP MASTER
          DO l = 1, llm
           DO ij = ijb, ije
@@ -282 +297 @@
       
       ENDIF
+      
+c$OMP END MASTER
+c$OMP BARRIER
 
 c         CALL  SCOPY( ijp1llm , finvmasse, 1, finvmaold, 1 )
-      
-      finvmaold(ijb:ije,:)=finvmasse(ijb:ije,:)        
+
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
+      DO l = 1, llm      
+        finvmaold(ijb:ije,l)=finvmasse(ijb:ije,l)        
+      ENDDO
+c$OMP END DO NOWAIT
 c
 c
@@ -292 +314 @@
 15    continue
 
+c$OMP MASTER
+        ps0(ijb:ije)=ps(ijb:ije)
+c$OMP END MASTER
 c    .................................................................
 
@@ -298 +323 @@
 c       CALL SCOPY (    ip1jmp1 ,  pscr   , 1,   psm1  , 1 )
 c       CALL SCOPY ( ip1jmp1*llm, massescr, 1,  massem1, 1 )
+c$OMP MASTER
         psm1(ijb:ije)=pscr(ijb:ije)
-        massem1(ijb:ije,:)=massescr(ijb:ije,:)
+c$OMP END MASTER
+
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
+          DO l = 1, llm
+            massem1(ijb:ije,l)=massescr(ijb:ije,l)
+          ENDDO
+c$OMP END DO NOWAIT       
       END IF