Changeset 353 for codes/icosagcm/trunk/src/etat0_dcmip5.f90

Timestamp:

08/06/15 10:48:45 (9 years ago)

Author:

dubos

Message:

OpenMP fixes for DCMIP

File:

• codes/icosagcm/trunk/src/etat0_dcmip5.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/etat0_dcmip5.f90

@@ -22 +22 @@
 
   INTEGER, SAVE :: dcmip5_testcase
+!$OMP THREADPRIVATE(dcmip5_testcase)
 
   PUBLIC getin_etat0, compute_etat0
@@ -41 +42 @@
   SUBROUTINE compute_etat0(ngrid,lon,lat, phis, ps, Temp, ulon, ulat, q)
     USE disvert_mod
+    USE omp_para
     INTEGER, INTENT(IN)    :: ngrid
     REAL(rstd),INTENT(IN)  :: lon(ngrid)
@@ -61 +63 @@
     END DO
 
-    DO l=1,llm
+    DO l=ll_begin,ll_end
        aa=.5*(ap(l)+ap(l+1))
        bb=.5*(bp(l)+bp(l+1))