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

codes/icosagcm/trunk/src/etat0_dcmip4.f90

@@ -12 +12 @@
   REAL(rstd),PARAMETER :: up0=1
   REAL(rstd) :: latw=2*Pi/9
+!$OMP THREADPRIVATE(latw)
   REAL(rstd) :: pw=34000
+!$OMP THREADPRIVATE(pw)
   REAL(rstd) :: q0=0.021
+!$OMP THREADPRIVATE(q0)
   REAL(rstd) :: lonc
+!$OMP THREADPRIVATE(lonc)
   REAL(rstd) :: latc
+!$OMP THREADPRIVATE(latc)
+
+  INTEGER,SAVE :: testcase
+!$OMP THREADPRIVATE(testcase)  
+
   PUBLIC  etat0
 CONTAINS
@@ -36 +45 @@
     REAL(rstd),POINTER :: q(:,:,:)
     INTEGER :: ind
+
+    testcase=1
+    CALL getin("dcmip4_testcase",testcase)
     
     DO ind=1,ndomain
+      IF (.NOT. assigned_domain(ind)) CYCLE
       CALL swap_dimensions(ind)
       CALL swap_geometry(ind)
@@ -81 +94 @@
   REAL(rstd) :: lonx,latx
   REAL(rstd) :: dthetaodeta_ave, dthetaodeta, dthetaodlat, duodeta, K, r
-  INTEGER :: testcase
   
     lonc=Pi/9
     latc=2*Pi/9 
-  
-    testcase=1
-    CALL getin("dcmip4_testcase",testcase)
-  
   
     DO l=1,llm