Changeset 985 for LMDZ4/trunk/libf/dyn3dpar/exner_hyb_p.F

Timestamp:

Jul 30, 2008, 5:50:03 PM (16 years ago)

Author:

Laurent Fairhead

Message:

Mise a jour de dyn3dpar par rapport a dyn3d, inclusion OpenMP et filtre FFT YM
LF

File:

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

LMDZ4/trunk/libf/dyn3dpar/exner_hyb_p.F

@@ -51 +51 @@
       INTEGER ije,ijb,jje,jjb
 c
-c$OMP MASTER           
+c$OMP BARRIER           
       unpl2k    = 1.+ 2.* kappa
 c
@@ -57 +57 @@
       ije=ij_end
 
-
+c$OMP DO SCHEDULE(STATIC)
       DO   ij  = ijb, ije
         pks(ij) = cpp * ( ps(ij)/preff ) ** kappa
       ENDDO
+c$OMP ENDDO
+c Synchro OPENMP ici
 
+c$OMP MASTER
       if (pole_nord) then
         DO  ij   = 1, iim
@@ -83 +86 @@
         ENDDO
       endif
-
+c$OMP END MASTER
 c
 c
 c    .... Calcul des coeff. alpha et beta  pour la couche l = llm ..
 c
+c$OMP DO SCHEDULE(STATIC)
       DO     ij      = ijb,ije
        alpha(ij,llm) = 0.
        beta (ij,llm) = 1./ unpl2k
       ENDDO
+c$OMP ENDDO NOWAIT
 c
 c     ... Calcul des coeff. alpha et beta  pour l = llm-1  a l = 2 ...
@@ -97 +102 @@
       DO l = llm -1 , 2 , -1
 c
+c$OMP DO SCHEDULE(STATIC)
         DO ij = ijb, ije
         dellta = p(ij,l)* unpl2k + p(ij,l+1)* ( beta(ij,l+1)-unpl2k )
@@ -102 +108 @@
         beta (ij,l)  =   p(ij,l  ) / dellta   
         ENDDO
+c$OMP ENDDO NOWAIT
 c
       ENDDO
@@ -109 +116 @@
 c     .....  Calcul de pk pour la couche 1 , pres du sol  ....
 c
-
+c$OMP DO SCHEDULE(STATIC)
       DO   ij   = ijb, ije
        pk(ij,1) = ( p(ij,1)*pks(ij) - 0.5*alpha(ij,2)*p(ij,2) )  /
      *    (  p(ij,1)* (1.+kappa) + 0.5*( beta(ij,2)-unpl2k )* p(ij,2) )
       ENDDO
+c$OMP ENDDO NOWAIT
 c
 c    ..... Calcul de pk(ij,l) , pour l = 2 a l = llm  ........
 c
       DO l = 2, llm
+c$OMP DO SCHEDULE(STATIC)
         DO   ij   = ijb, ije
          pk(ij,l) = alpha(ij,l) + beta(ij,l) * pk(ij,l-1)
         ENDDO
+c$OMP ENDDO NOWAIT        
       ENDDO
 c
 c
 c      CALL SCOPY   ( ngrid * llm, pk, 1, pkf, 1 )
-      pkf(ijb:ije,1:llm)=pk(ijb:ije,1:llm)
-c$OMP END MASTER
+      DO l = 1, llm
+c$OMP DO SCHEDULE(STATIC)
+         DO   ij   = ijb, ije
+           pkf(ij,l)=pk(ij,l)
+         ENDDO
+c$OMP ENDDO NOWAIT             
+      ENDDO
+
 c$OMP BARRIER