Changeset 71 for trunk/SOURCES/moy_mxmy.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/moy_mxmy.f90 (modified) (2 diffs)

trunk/SOURCES/moy_mxmy.f90

@@ -18 +18 @@
 subroutine moy_mxmy(n1,n2,X2D,X_mx,X_my)
 ! fait la moyenne ponderee d'un tableau X2D sur les mailles staggered
-
+ !$ USE OMP_LIB
 implicit none
 integer, intent(in) :: n1,n2   !< dimension des tableaux
@@ -29 +29 @@
 
 
+  !$OMP PARALLEL PRIVATE(i_moins1,j_moins1,i_plus1,j_plus1)
+  !$OMP DO
+        do j=1,n2
+                do i=1,n1
+                        i_moins1=max(1,i-1)
+                        j_moins1=max(1,j-1)
+                        i_plus1=min(n1,i+1)
+                        j_plus1=min(n2,j+1)
 
-do j=1,n2
-   do i=1,n1
-      i_moins1=max(1,i-1)
-      j_moins1=max(1,j-1)
-      i_plus1=min(n1,i+1)
-      j_plus1=min(n2,j+1)
-
-      X_mx(i,j)=0.25*(X2D(i,j)+X2D(i_moins1,j)) &
+                        X_mx(i,j)=0.25*(X2D(i,j)+X2D(i_moins1,j)) &
            + 0.125*((X2D(i_moins1,j_plus1)+X2D(i,j_plus1))   &
            +       (X2D(i_moins1,j_moins1)+X2D(i,j_moins1)))
 
-      X_my(i,j)=0.25*(X2D(i,j)+X2D(i,j_moins1)) &
+                        X_my(i,j)=0.25*(X2D(i,j)+X2D(i,j_moins1)) &
            + 0.125*((X2D(i_plus1,j_moins1)+X2D(i_plus1,j))   &
            +       (X2D(i_moins1,j_moins1)+X2D(i_moins1,j)))
-
-
-   end do
-end do
+                end do
+        end do
+  !$OMP END DO
+  !$OMP END PARALLEL
 
 return