Changeset 74 for trunk/SOURCES/taubed-0.3.f90

Timestamp:

06/24/16 09:19:59 (8 years ago)

Author:

dumas

Message:

OpenMP parallelization in remplimat-shelves-tabTu, litho and taubed

File:

• trunk/SOURCES/taubed-0.3.f90 (modified) (5 diffs)

trunk/SOURCES/taubed-0.3.f90

@@ -43 +43 @@
 subroutine taubed()
 
+  !$USE OMP_LIB 
   USE module3D_phy
   USE param_phy_mod
@@ -52 +53 @@
   ! NLITH est defini dans isostasie et permet le choix du modele d'isostasie
 
+  
   if (NLITH.eq.1) then
      !       avec rigidite de la lithosphere
+     !$OMP PARALLEL
+     !$OMP DO
      do J=1,NY 
         do I=1,NX
@@ -65 +69 @@
         end do
      end do
-
+     !$OMP END DO
 
      ! il faut remplir CHARGE dans les parties a l'exterieur de la grille :
      ! a l'exterieur de la grille CHARGE est egale a la valeur sur le bord
-
+     !$OMP DO
      do J=1,NY
         CHARGE(1-LBLOC:0,J)=CHARGE(1,J)      ! bord W
         CHARGE(NX+1:NX+LBLOC,J)=CHARGE(NX,J) ! bord E
      end do
+     !$OMP END DO
+     !$OMP DO
      do I=1,NX
         CHARGE(I,1-LBLOC:0)=CHARGE(I,1)      ! bord S
         CHARGE(I,NY+1:NY+LBLOC)=CHARGE(I,NY) ! bord N
      end do
+     !$OMP END DO
 
      ! valeur dans les quatres coins exterieurs au domaine       
+          !$OMP WORKSHARE
      CHARGE(1-LBLOC:0,1-LBLOC:0)=CHARGE(1,1)           ! coin SW
      CHARGE(1-LBLOC:0,NY+1:NY+LBLOC)=CHARGE(1,NY)      ! coin NW
      CHARGE(NX+1:NX+LBLOC,1-LBLOC:0)=CHARGE(NX,1)      ! coin SE
      CHARGE(NX+1:NX+LBLOC,NY+1:NY+LBLOC)=CHARGE(NX,NY) ! coin NE
-
+     !$OMP END WORKSHARE
+     !$OMP END PARALLEL 
      call litho
 
   else
      !     enfoncement local
+     !$OMP PARALLEL
+     !$OMP DO
      do J=1,NY
         do I=1,NX
@@ -100 +111 @@
         end do
      end do
+     !$OMP END DO
+     !$OMP END PARALLEL
   endif
 
   !     decroissance exponentielle de l'enfoncement
+  !$OMP PARALLEL
+  !$OMP DO
   do J=1,NY
      do I=1,NX
@@ -109 +124 @@
      end do
   end do
+  !$OMP END DO
+  !$OMP END PARALLEL
 
 end subroutine taubed