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

Timestamp:

01/09/14 09:56:11 (11 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_dcmip1.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/etat0_dcmip1.f90

@@ -4 +4 @@
 
   REAL(rstd), SAVE  :: h0=1.
+!$OMP THREADPRIVATE(h0)
   REAL(rstd), SAVE  :: lon0=3*pi/2
+!$OMP THREADPRIVATE(lon0)
   REAL(rstd), SAVE  :: lat0=0.0
+!$OMP THREADPRIVATE(lat0)
   REAL(rstd), SAVE  :: alpha=0.0
+!$OMP THREADPRIVATE(alpha)
   REAL(rstd), SAVE  :: R0 
+!$OMP THREADPRIVATE(R0)
   REAL(rstd), SAVE  :: lat1=0.
+!$OMP THREADPRIVATE(lat1)
   REAL(rstd), SAVE  :: lat2=0.
+!$OMP THREADPRIVATE(lat2)
   REAL(rstd), SAVE  :: lon1=pi/6
+!$OMP THREADPRIVATE(lon1)
   REAL(rstd), SAVE  :: lon2=-pi/6
+!$OMP THREADPRIVATE(lon2)
   REAL(rstd), SAVE  :: latc1=0.
+!$OMP THREADPRIVATE(latc1)
   REAL(rstd), SAVE  :: latc2=0.
+!$OMP THREADPRIVATE(latc2)
   REAL(rstd), SAVE  :: lonc1=5*pi/6
+!$OMP THREADPRIVATE(lonc1)
   REAL(rstd), SAVE  :: lonc2=7*pi/6
+!$OMP THREADPRIVATE(lonc2)
   REAL(rstd), SAVE  :: zt=1000.0
+!$OMP THREADPRIVATE(zt)
   REAL(rstd), SAVE  :: rt
+!$OMP THREADPRIVATE(rt)
   REAL(rstd), SAVE  :: zc=5000.0
+!$OMP THREADPRIVATE(zc)
 
   PUBLIC etat0
@@ -47 +63 @@
     
     DO ind=1,ndomain
+      IF (.NOT. assigned_domain(ind)) CYCLE
       CALL swap_dimensions(ind)
       CALL swap_geometry(ind)