Changeset 186 for codes/icosagcm/trunk/src/theta_rhodz.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/theta_rhodz.f90 (modified) (5 diffs)

codes/icosagcm/trunk/src/theta_rhodz.f90

@@ -16 +16 @@
 
     DO ind=1,ndomain
+      IF (.NOT. assigned_domain(ind)) CYCLE
       CALL swap_dimensions(ind)
       CALL swap_geometry(ind)
@@ -39 +40 @@
 
     DO ind=1,ndomain
+      IF (.NOT. assigned_domain(ind)) CYCLE
       CALL swap_dimensions(ind)
       CALL swap_geometry(ind)
@@ -62 +64 @@
 
     DO ind=1,ndomain
+      IF (.NOT. assigned_domain(ind)) CYCLE
       CALL swap_dimensions(ind)
       CALL swap_geometry(ind)
@@ -82 +85 @@
     INTEGER :: i,j,ij,l
     REAL(rstd),ALLOCATABLE,SAVE :: p(:,:)
+!$OMP THREADPRIVATE(p)
 
     ALLOCATE( p(iim*jjm,llm+1))
@@ -109 +113 @@
     INTEGER :: i,j,ij,l
     REAL(rstd),SAVE,ALLOCATABLE :: p(:,:)
+!$OMP THREADPRIVATE(p)
 
     ALLOCATE( p(iim*jjm,llm+1))