Changeset 71 for trunk/SOURCES/Temperature-routines/Qprod_icetemp.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/Temperature-routines/Qprod_icetemp.f90 (modified) (10 diffs)

trunk/SOURCES/Temperature-routines/Qprod_icetemp.f90

@@ -17 +17 @@
 
   use Icetemp_declar
+!$ USE OMP_LIB
+
 
   implicit none
@@ -41 +43 @@
   real, dimension(Nx,Ny)      :: Tox           !< Contraintes Sur Maille Mx
   real, dimension(Nx,Ny)      :: Toy           !< Contraintes Sur Maille Mx
+  real      :: Chalk_1       !< Utilise Pour Le Calcul De Chalk : Glace Posée
+  real      :: Chalk_2       !< Utilise Pour Le Calcul De Chalk : Ice Streams Et Ice
 
 
@@ -52 +56 @@
      ! Q_prod_demi(T_m%T,T_m%Tpmp,Deform_m%Tobmx,Deform_m%Tobmy,Geom_m%H,Geom_m%Hmx,Geom_m%Hmy,mask_flot_m%Flot,mask_flot_m%Flotmx,mask_flot_m%Flotmy,mask_flot_m%Flgzmx ,mask_flot_m%Flgzmy,Geom_m%Slop_x,Geom_m%Slop_y,Ice_flow_m%Ddx,Ice_flow_m%Ddy,&
      !      Ice_flow_m%Ddbx,Ice_flow_m%Ddby,Deform_m%Epsxx,Deform_m%Epsyy,Deform_m%Epsxy,mask_flot_m%Gzmx,mask_flot_m%Gzmy,Deform_m%Btt,Ice_flow_m%Uxbar,Ice_flow_m%Uybar,therm_var_m%Phid,Deform_m%Glen,Deform_m%Visc)
-
-     do K=2,Nz
-        do J=2,Ny-1
-           do I=2,Nx-1
+     
+
+
+    !$OMP PARALLEL PRIVATE(chalk_1,chalk_2)
+    !$OMP DO COLLAPSE(2)
+    do K=2,Nz
+    do J=2,Ny-1
+            do I=2,Nx-1
 
               ! Calcul De La Chaleur De Deformation Selon Xx Yy Zz Et Xy
@@ -80 +88 @@
            end do
         end do
-     end do
+!        !$ i_min=min(i_min,j)
+!        !$ i_max=max(i_max,j)
+     end do
+        !$OMP END DO NOWAIT
+!       !$ print*,'nb_taches = ', nb_taches
+!   !$ print*,"Rang : ",rang," i_min : ",i_min," i_max : ",i_max
 
      !  Partie Sia   Calcul De La Chaleur Produite Sur Chaque Demi Maille
-     do L=1,size(Btt,4) !N1poly,N2poly
-        do J=2,Ny
+                !$OMP DO COLLAPSE(2)
+                do L=1,size(Btt,4) !N1poly,N2poly
+      do J=2,Ny
            do I=2,Nx
               !          Ffx A 3 Dimensions !
@@ -92 +106 @@
         end do
      end do
-
-     do L=1,size(Btt,4)!N1poly,N2poly
-        do K=2,Nz
-           do J=2,Ny
-              do I=2,Nx
-                 if ((.not.Flotmx(I,J)).and.(.not.fleuvemx(I,J))) then ! grounded et slowssa
-                    Chalx(I,J,K,L)=(Btt(I-1,J,K,L)+Btt(I,J,K,L))*Ffx(I,J,L) !&
+     !$OMP END DO
+     
+     !$OMP DO COLLAPSE(2)
+          do K=2,Nz
+       do J=2,Ny
+              do I=2,Nx
+                     do L=1,size(Btt,4)!N1poly,N2poly
+                    if ((.not.Flotmx(I,J)).and.(.not.fleuvemx(I,J))) then ! grounded et slowssa
+                        Chalx(I,J,K,L)=(Btt(I-1,J,K,L)+Btt(I,J,K,L))*Ffx(I,J,L) !&
                     !                      *Ro*G*Ee(K)**(Glen(L)+1)/Cp(I,J,K)
 
-                 else if (fleuvemx(I,J)) then ! Ice Streams
-                    Chalx(I,J,K,L)=0.
-                 else                     ! Ice Shelves  
-                    Chalx(I,J,K,L)=0.
-                 endif
-                 if ((.not.Flotmy(I,J)).and.(.not.fleuvemy(I,J))) then ! grounded et slowssa
-                    Chaly(I,J,K,L)=(Btt(I,J-1,K,L)+Btt(I,J,K,L))*Ffy(I,J,L) !&
+                    else if (fleuvemx(I,J)) then ! Ice Streams
+                        Chalx(I,J,K,L)=0.
+                    else                     ! Ice Shelves  
+                        Chalx(I,J,K,L)=0.
+                    endif
+                    if ((.not.Flotmy(I,J)).and.(.not.fleuvemy(I,J))) then ! grounded et slowssa
+                        Chaly(I,J,K,L)=(Btt(I,J-1,K,L)+Btt(I,J,K,L))*Ffy(I,J,L) !&
                     !                      *Ro*G*Ee(K)**(Glen(L)+1)/Cp(I,J,K)
-                 else if (fleuvemy(I,J)) then  ! Ice Streams
-                    Chaly(I,J,K,L)=0.
-                 else                      ! Ice Shelves
-                    Chaly(I,J,K,L)=0.
-                 endif
-
-              end do
-           end do
-        end do
-     end do
+                    else if (fleuvemy(I,J)) then  ! Ice Streams
+                        Chaly(I,J,K,L)=0.
+                    else                      ! Ice Shelves
+                        Chaly(I,J,K,L)=0.
+                    endif
+
+                   end do
+            end do
+      end do
+     end do
+     !$OMP END DO
+     
 
      !         Nouvelle Formulation De Chaldef_maj(I,J,K), Le 4 Vient Des Moyennes
@@ -125 +143 @@
      !         Ancienne Formulation Chal=(Ro*G*H(I,J))**4*(Sx2+Sy2)*(Sx*Sx+Sy*Sy) 
 
-     do K=2,Nz
-        do J=1,Ny-1
-           do I=1,Nx-1
-
+    !$OMP DO COLLAPSE(2)
+    do K=2,Nz
+      do J=1,Ny-1
+            do I=1,Nx-1
               ! Modif Christophe Mars 2000 : Chalx Et Chaly Sont A 4 Dim
               Chaldef_maj(I,J,K)= 0.
@@ -152 +170 @@
         end do
      end do
+     !$OMP END DO NOWAIT
+          
      ! Chaleur Produite A La Base Par Le Glissement
+     !$OMP DO
      do J=2,Ny
         do I=2,Nx
@@ -168 +189 @@
            endif
 
-        end do
-     end do
-
+        end do   
+     end do
+     !$OMP END DO
+          
      !  Boundary Condition Ice-Rock Interface 
      K=Nz
 
+     !$OMP DO
      do J=2,Ny-1
         do I=2,Nx-1
@@ -227 +250 @@
         end do
      end do
-
+     !$OMP END DO
+     !$OMP END PARALLEL
+     
 
   case (2)             ! Q_prod_centre : Calcul Avec La Somme Des Carres
@@ -904 +929 @@
 
   end select
+
 end subroutine Qprod_ice