Changeset 77 for branches/iLoveclim/SOURCES/diffusiv-polyn-0.6.f90

Timestamp:

07/05/16 11:20:52 (8 years ago)

Author:

dumas

Message:

Merge branche iLOVECLIM sur rev 76

Location:

branches/iLoveclim

Files:

• . (modified) (1 prop)
• SOURCES/diffusiv-polyn-0.6.f90 (modified) (12 diffs)

branches/iLoveclim/SOURCES/diffusiv-polyn-0.6.f90

@@ -34 +34 @@
 
 !     ===============================================================
-
+!$ USE OMP_LIB
 USE module3D_phy
 
@@ -44 +44 @@
 REAL :: INV_4DX ! inverse de dx pour eviter les divisions =1/(4*dx)
 REAL :: INV_4DY ! inverse de dy pour eviter les divisions =1/(4*dy)
-INTEGER :: ii
-integer :: ll
 
 if (itracebug.eq.1)  call tracebug(' Entree dans routine diffusiv')
@@ -60 +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 :
@@ -81 +82 @@
 ! Calcul de la vitesse de glissement, qu'elle soit calculee par diagno ou par sliding
 !------------------------------------------------------------------------------------
-! 
+!$OMP PARALLEL
+!$OMP WORKSHARE
 where (flgzmx(:,:))
     ubx(:,:) = uxflgz(:,:)
@@ -93 +95 @@
     uby(:,:) = ddby(:,:)* (-sdy(:,:))
 endwhere
-
-
+!$OMP END WORKSHARE
+!$OMP END PARALLEL
 if (itracebug.eq.1) call tracebug(' diffusiv  apres calcul glissement')
 
@@ -101 +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)) &
@@ -135 +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 
@@ -166 +169 @@
     end do
    end do ydef
-
+!$OMP END DO
+
+!$OMP DO
 xdef: do j=1,ny
        do i=2,nx 
@@ -194 +199 @@
     end do
    end do xdef
+!$OMP END DO
 
 if (itracebug.eq.1) call tracebug(' diffusiv  avant limit')
@@ -201 +207 @@
 
 ! la limitation selon x
+!$OMP WORKSHARE
 where(.not.flgzmx(:,:))
    uxbar(:,:)=min(uxbar(:,:),ultot)
@@ -211 +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
@@ -223 +232 @@
    end do
 end do
-
+!$OMP END DO
+!$OMP END PARALLEL
 if (itracebug.eq.1)  call tracebug(' fin diffusiv ')