Changeset 985 for LMDZ4/trunk/libf/dyn3dpar/qminimum_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/qminimum_p.F (modified) (5 diffs)

LMDZ4/trunk/libf/dyn3dpar/qminimum_p.F

@@ -34 +34 @@
       SAVE imprim
       DATA imprim /0/
+c$OMP THREADPRIVATE(imprim)
       INTEGER ijb,ije
+      INTEGER Index_pump(ip1jmp1)
+      INTEGER nb_pump
 c
 c Quand l'eau liquide est trop petite (ou negative), on prend
@@ -40 +43 @@
 c (sans changer la temperature !)
 c
+
       ijb=ij_begin
       ije=ij_end
-      
+
+c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)       
       DO 1000 k = 1, llm
       DO 1040 i = ijb, ije
@@ -50 +55 @@
  1040 CONTINUE
  1000 CONTINUE
+c$OMP END DO NOWAIT
+c$OMP BARRIER
+c --->  SYNCHRO OPENMP ICI
+
 c
 c Quand l'eau vapeur est trop faible (ou negative), on complete
@@ -58 +67 @@
       DO k = llm, 2, -1
 ccc      zx_abc = dpres(k) / dpres(k-1)
+c$OMP DO SCHEDULE(STATIC)
       DO i = ijb, ije
          zx_abc = deltap(i,k)/deltap(i,k-1)
@@ -64 +74 @@
          q(i,k,iq)   =  q(i,k,iq)   + zx_defau  
       ENDDO
+c$OMP END DO NOWAIT
       ENDDO
+c$OMP BARRIER
 c
 c Quand il s'agit de la premiere couche au-dessus du sol, on
 c doit imprimer un message d'avertissement (saturation possible).
 c
+      nb_pump=0
+c$OMP DO SCHEDULE(STATIC)
       DO i = ijb, ije
          zx_pump(i) = AMAX1( 0.0, seuil_vap - q(i,1,iq) )
          q(i,1,iq)  = AMAX1( q(i,1,iq), seuil_vap )
+         IF (zx_pump(i) > 0.0) THEN
+            nb_pump = nb_pump+1
+            Index_pump(nb_pump)=i
+         ENDIF
       ENDDO
-      pompe = SSUM(ije-ijb+1,zx_pump(ijb),1)
-      IF (imprim.LE.500 .AND. pompe.GT.0.0) THEN
-         WRITE(6,'(1x,"ATT!:on pompe de l eau au sol",e15.7)') pompe
-         DO i = ijb, ije
-            IF (zx_pump(i).GT.0.0) THEN
+c$OMP END DO  
+!      pompe = SSUM(ije-ijb+1,zx_pump(ijb),1)
+
+      IF (imprim.LE.100 .AND. nb_pump .GT. 0 ) THEN
+         PRINT *, 'ATT!:on pompe de l eau au sol'
+         DO i = 1, nb_pump
                imprim = imprim + 1
-               PRINT*,'  en ',i,zx_pump(i)
-            ENDIF
+               PRINT*,'  en ',index_pump(i),zx_pump(index_pump(i))
          ENDDO
       ENDIF