Changeset 76 for trunk/SOURCES/diffusiv-polyn-0.6.f90

Timestamp:

06/29/16 16:21:13 (8 years ago)

Author:

dumas

Message:

OpenMP parallelization in conserv-mass-adv-diff_sept2009_mod.f90, diffusiv-polyn-0.6.f90, dragging_neff_slope_mod.f90, eaubasale-0.5_mod.f90 and relaxation_water_diffusion.f90

File:

• trunk/SOURCES/diffusiv-polyn-0.6.f90 (modified) (11 diffs)

trunk/SOURCES/diffusiv-polyn-0.6.f90

@@ -34 +34 @@
 
 !     ===============================================================
-
+!$ USE OMP_LIB
 USE module3D_phy
 
@@ -58 +58 @@
 
 ! initialisation de la diffusion
+!$OMP PARALLEL
+!$OMP WORKSHARE
 diffmx(:,:)=0.
 diffmy(:,:)=0.
-
+!$OMP END WORKSHARE
+!$OMP END PARALLEL
 
 ! calcul de SDX, SDY et de la pente au carre en mx et en my :
@@ -79 +82 @@
 ! Calcul de la vitesse de glissement, qu'elle soit calculee par diagno ou par sliding
 !------------------------------------------------------------------------------------
-! 
+!$OMP PARALLEL
+!$OMP WORKSHARE
 where (flgzmx(:,:))
     ubx(:,:) = uxflgz(:,:)
@@ -91 +95 @@
     uby(:,:) = ddby(:,:)* (-sdy(:,:))
 endwhere
-
-
+!$OMP END WORKSHARE
+!$OMP END PARALLEL
 if (itracebug.eq.1) call tracebug(' diffusiv  apres calcul glissement')
 
@@ -99 +103 @@
 ! Deformation SIA : Calcul de ddy  et ddx pour chaque valeur de iglen
 !------------------------------------------------------------------
-
 do iglen=n1poly,n2poly
    glenexp=max(0.,(glen(iglen)-1.)/2.)
-
-
+  !$OMP PARALLEL 
+  !$OMP DO
    do j=1,ny
       do i=1,nx
-
          if (.not.flotmy(i,j)) then !  On calcule quand la deformation même pour les ice streams
             ddy(i,j,iglen)=((slope2my(i,j)**  &  ! slope2my calcule en debut de routine
                  glenexp)*(rog)**glen(iglen)) &
                  *hmy(i,j)**(glen(iglen)+1)
-
-
          endif
       end do
    end do
+   !$OMP END DO
+   !$OMP END PARALLEL
 end do
+
+
 
 do iglen=n1poly,n2poly
    glenexp=max(0.,(glen(iglen)-1.)/2.)
-
-
+   !$OMP PARALLEL
+   !$OMP DO
    do j=1,ny  
       do i=1,nx
-
          if (.not.flotmx(i,j)) then      ! bug y->x corrige le 5 avril 2008 (
-                                    
             ddx(i,j,iglen)=((slope2mx(i,j)**  & ! slope2mx calcule en debut de routine
                  glenexp)*(rog)**glen(iglen)) &
@@ -133 +135 @@
       end do
    end do
+   !$OMP END DO
+   !$OMP END PARALLEL
 end do
 
+
 ! vitesse due a la déformation-----------------------------------------------
-
-
+!$OMP PARALLEL
+!$OMP DO
 ydef: do j=2,ny
        do i=1,nx 
@@ -164 +169 @@
     end do
    end do ydef
-
+!$OMP END DO
+
+!$OMP DO
 xdef: do j=1,ny
        do i=2,nx 
@@ -192 +199 @@
     end do
    end do xdef
+!$OMP END DO
 
 if (itracebug.eq.1) call tracebug(' diffusiv  avant limit')
@@ -199 +207 @@
 
 ! la limitation selon x
+!$OMP WORKSHARE
 where(.not.flgzmx(:,:))
    uxbar(:,:)=min(uxbar(:,:),ultot)
@@ -209 +218 @@
    uybar(:,:)=max(uybar(:,:),-ultot)
 end where
+!$OMP END WORKSHARE
 
 !    cas particulier des sommets ou il ne reste plus de neige.
+!$OMP DO
 do j=2,ny-1
    do i=2,nx-1
@@ -221 +232 @@
    end do
 end do
-
+!$OMP END DO
+!$OMP END PARALLEL
 if (itracebug.eq.1)  call tracebug(' fin diffusiv ')