Changeset 74 for trunk/SOURCES/litho-0.4.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/litho-0.4.f90 (modified) (1 diff)

trunk/SOURCES/litho-0.4.f90

@@ -37 +37 @@
 
  
-      subroutine litho
+subroutine litho
+  !$ USE OMP_LIB
+  USE module3D_phy
+  USE param_phy_mod
+  USE ISO_DECLAR ! module de declaration des variables specifiques a l'isostasie
 
-      USE module3D_phy
-      USE param_phy_mod
-      USE ISO_DECLAR ! module de declaration des variables specifiques a l'isostasie
+  implicit none
 
-      implicit none
-
-      INTEGER :: II,SOM1,SOM2
-      REAL, dimension(:,:), allocatable :: WLOC
+  INTEGER :: II,SOM1,SOM2
+  REAL, dimension(:,:), allocatable :: WLOC
 
 
-!----- allocation de WLOC  et de croix -----------
+  !----- allocation de WLOC  et de croix -----------
 
-      if (.not.allocated(WLOC)) THEN
-         allocate(WLOC(-LBLOC:LBLOC,-LBLOC:LBLOC),stat=err)
-         if (err/=0) then
-            print *,"Erreur a l'allocation du tableau WLOC",err
-            stop 4
-         end if
-      end if
+  if (.not.allocated(WLOC)) THEN
+     allocate(WLOC(-LBLOC:LBLOC,-LBLOC:LBLOC),stat=err)
+     if (err/=0) then
+        print *,"Erreur a l'allocation du tableau WLOC",err
+        stop 4
+     end if
+  end if
 
-!----- fin de l'allocation --------------
+  !----- fin de l'allocation --------------
 
-!     calcul de la deflexion par sommation des voisins appartenant
-!     au bloc de taille 2*LBLOC
-       som1=0.
-       som2=0.
+  !     calcul de la deflexion par sommation des voisins appartenant
+  !     au bloc de taille 2*LBLOC
+  som1=0.
+  som2=0.
 
-!      On somme aussi les contributions des points exterieurs au domaine
-!      lorsque la charge est due a l'ocean. On suppose alors  que
-!      ces points ont la meme charge que les limites
-
-boucleij : do J=1,NY
-           do I=1,NX
+  !      On somme aussi les contributions des points exterieurs au domaine
+  !      lorsque la charge est due a l'ocean. On suppose alors  que
+  !      ces points ont la meme charge que les limites
+  !$OMP PARALLEL
+  !$OMP DO PRIVATE(Wloc)
+  boucleij : do J=1,NY
+     do I=1,NX
 
 
-              W1(I,J)=0. 
-              ii=0
+        W1(I,J)=0. 
+ !   ii=0
 
-! Wloc : enfoncement centre sur ij du a chaque charges autour
-      Wloc(:,:) = We(:,:) * charge(i-lbloc:i+lbloc,j-lbloc:j+lbloc)
+ ! Wloc : enfoncement centre sur ij du a chaque charges autour
+        Wloc(:,:) = We(:,:) * charge(i-lbloc:i+lbloc,j-lbloc:j+lbloc)
 
 
-      W1(i,j) = sum (Wloc(:,:))   ! effet en ij de toutes les charges autour (distance < lbloc)
+        W1(i,j) = sum (Wloc(:,:))   ! effet en ij de toutes les charges autour (distance < lbloc)
 
-!    ------------------------ FIN DE L'INTEGRATION SUR LE PAVE LBLOC
-              som1=som1+W1(I,J)
-              som2=som2-charge(I,J)/ROMG
+ !    ------------------------ FIN DE L'INTEGRATION SUR LE PAVE LBLOC
+ !   som1=som1+W1(I,J)
+ !   som2=som2-charge(I,J)/ROMG
 
-      end do
-      end do  boucleij
-       
-!     write(6,*)'enfoncement total avec rigidite de la lithosphere:',
-!    &  som1
-!     write(6,*)'enfoncement total avec isostasie locale :',
-!    &  som2
+     end do
+  end do  boucleij
+  !$OMP END DO
+  !$OMP END PARALLEL
+
+  !     write(6,*)'enfoncement total avec rigidite de la lithosphere:',
+  !    &  som1
+  !     write(6,*)'enfoncement total avec isostasie locale :',
+  !    &  som2
 
 
-    end subroutine litho
-      
+end subroutine litho
+