Changeset 186 for codes/icosagcm/trunk/src/etat0_academic.f90

Timestamp:

01/09/14 09:56:11 (10 years ago)

Author:

ymipsl

Message:

Add new openMP parallelism based on distribution of domains on threads. There is no more limitation of number of threads by MPI process.

YM

File:

• codes/icosagcm/trunk/src/etat0_academic.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/etat0_academic.f90

@@ -10 +10 @@
   USE kinetic_mod
   IMPLICIT NONE
-    TYPE(t_field),POINTER :: f_ps(:)
-    TYPE(t_field),POINTER :: f_phis(:)
-    TYPE(t_field),POINTER :: f_theta_rhodz(:)
-    TYPE(t_field),POINTER :: f_u(:)
-    TYPE(t_field),POINTER :: f_q(:)
-    TYPE(t_field),POINTER :: f_Ki(:)
-    TYPE(t_field),POINTER :: f_temp(:)
+    TYPE(t_field),POINTER,SAVE :: f_ps(:)
+    TYPE(t_field),POINTER,SAVE :: f_phis(:)
+    TYPE(t_field),POINTER,SAVE :: f_theta_rhodz(:)
+    TYPE(t_field),POINTER,SAVE :: f_u(:)
+    TYPE(t_field),POINTER,SAVE :: f_q(:)
+    TYPE(t_field),POINTER,SAVE :: f_Ki(:)
+    TYPE(t_field),POINTER,SAVE :: f_temp(:)
   
     REAL(rstd),POINTER :: Ki(:,:)
@@ -58 +58 @@
     
     DO ind=1,ndomain
+      IF (.NOT. assigned_domain(ind)) CYCLE
       CALL swap_dimensions(ind)
       CALL swap_geometry(ind)