Changeset 71 for trunk/SOURCES/mix-SIA-L1_mod.f90

Timestamp:

06/15/16 17:13:33 (8 years ago)

Author:

dumas

Message:

First parallelization instructions with OpenMP tested in debug : exactly the same results | isostasy called only every 50 years

File:

• trunk/SOURCES/mix-SIA-L1_mod.f90 (modified) (4 diffs)

trunk/SOURCES/mix-SIA-L1_mod.f90

@@ -53 +53 @@
 subroutine mix_SIA_L1
 
+!$ USE OMP_LIB
+! variables parallelisation openMP
+!!$  integer :: rang ,nb_taches
 
 ! subroutine qui associe SIA et L1
@@ -59 +62 @@
 
 ! use module_choix
-    if (itracebug.eq.1)  call tracebug(' routine mix_SIA_L1')
+if (itracebug.eq.1)  call tracebug(' routine mix_SIA_L1')
 
 debug_3D(:,:,9)=0.
@@ -73 +76 @@
    ! que celles du L1. On garde le L1 comme vitesse basale.
    ! version utilisee pour les runs MIS11 de Cairns
-
+        !$OMP PARALLEL
+   !$OMP DO
    do j=2,ny
       do i=2,nx
@@ -95 +99 @@
       end do
    end do
+   !$OMP END DO
+        !$OMP END PARALLEL
 
 else if (i_resolmeca.eq.2) then
 
 ! addition systematique
-
+!       !$OMP PARALLEL PRIVATE(rang,nb_taches)
+!       !$ rang=OMP_GET_THREAD_NUM()
+!       !$ nb_taches=OMP_GET_NUM_THREADS()
+!       !$OMP DO SCHEDULE(STATIC,NY/nb_taches)
+        !$OMP PARALLEL
+   !$OMP DO
    do j=2,ny
       do i=2,nx
+! test sur tout le tableau :
          if (flgzmx(i,j)) then 
-            debug_3D(i,j,11)=uxflgz(i,j)
+!            debug_3D(i,j,11)=uxflgz(i,j)
             uxbar(i,j)=uxflgz(i,j)+uxdef(i,j)                       ! on ajoute la deform. SIA
             ubx(i,j) = uxflgz(i,j)
-            do k=1,nz
+!            do k=1,nz
+!               ux(i,j,k)=ux(i,j,k)-ux(i,j,nz)+uxflgz(i,j)           ! remplace le glissement par uxflgz
+!            end do
+
+!            debug_3D(i,j,9)=uxdef(i,j)                   
+         endif
+
+         if (flgzmy(i,j)) then 
+!            debug_3D(i,j,12)=uyflgz(i,j)
+            uybar(i,j)=uyflgz(i,j)+uydef(i,j)                       ! on ajoute la deform. SIA
+            uby(i,j) = uyflgz(i,j)
+!            do k=1,nz
+!               uy(i,j,k)=uy(i,j,k)-uy(i,j,nz)+uyflgz(i,j)           ! remplace le glissement par uyflgz
+!            end do
+
+!            debug_3D(i,j,10)=uydef(i,j)              
+         endif
+      end do
+   end do
+        !$OMP END DO
+        !$OMP END PARALLEL
+
+! on ne recalcul ux que lorsque uxflgz est modifié soit tous les dtt
+        if (isynchro.eq.1) then
+!       !$OMP PARALLEL PRIVATE(rang,nb_taches)
+!       !$ rang=OMP_GET_THREAD_NUM()
+!       !$ nb_taches=OMP_GET_NUM_THREADS()
+!       !$OMP DO SCHEDULE(STATIC,NY/nb_taches)
+        !$OMP PARALLEL
+   !$OMP DO
+        do j=2,ny
+                do i=2,nx
+        if (flgzmx(i,j)) then 
+                        do k=1,nz
                ux(i,j,k)=ux(i,j,k)-ux(i,j,nz)+uxflgz(i,j)           ! remplace le glissement par uxflgz
             end do
-
-            debug_3D(i,j,9)=uxdef(i,j)                   
-         end if
-
+         endif
          if (flgzmy(i,j)) then 
-            debug_3D(i,j,12)=uyflgz(i,j)
-            uybar(i,j)=uyflgz(i,j)+uydef(i,j)                       ! on ajoute la deform. SIA
-            uby(i,j) = uyflgz(i,j)
-            do k=1,nz
+                        do k=1,nz
                uy(i,j,k)=uy(i,j,k)-uy(i,j,nz)+uyflgz(i,j)           ! remplace le glissement par uyflgz
             end do
-
-            debug_3D(i,j,10)=uydef(i,j)              
-         end if
-      end do
-   end do
-
-
+         endif
+       enddo
+   enddo
+        !$OMP END DO
+        !$OMP END PARALLEL 
+endif
 
 else ! SIA et L1 en 2 zones separe : pas d'action (uxbar est deja remis a uxflgz dans diffusiv)