Changeset 76 for trunk/SOURCES/Temperature-routines/prop_th_icetemp.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/Temperature-routines/prop_th_icetemp.f90 (modified) (4 diffs)

trunk/SOURCES/Temperature-routines/prop_th_icetemp.f90

@@ -15 +15 @@
 
 Subroutine Thermal_prop_icetemp
+!$ USE OMP_LIB
   Use Icetemp_declar 
 
@@ -21 +22 @@
   If (Itracebug.Eq.1) Write(Num_tracebug,*)' Entree Dans Routine Thermal_prop_icetemp'
 
+!$OMP PARALLEL
+!$OMP DO COLLAPSE(2)
   Do K=1,Nz  
      Do J=1,Ny
@@ -29 +32 @@
            !  Attention Pour La Conductivite C'Est La Formule De Yin-Chao Yen
 
-           Ct(I,J,K)=6.62E7
-           Cp(I,J,K)=2009. 
+!           Ct(I,J,K)=6.62E7
+!           Cp(I,J,K)=2009. 
            Tpmp(I,J,K)=-0.00087*(K-1)*1./(Nz-1)*H(I,J)!De=1/(Nz_m-1)
            Cp(I,J,K)=(2115.3+7.79293*T(I,J,K)) 
@@ -38 +41 @@
      End Do
   End Do
-
+!$OMP END DO
+!$OMP END PARALLEL
 End Subroutine Thermal_prop_icetemp