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

Timestamp:

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

Author:

dubos

Message:

OpenMP fixes for DCMIP

File:

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

codes/icosagcm/trunk/src/etat0_dcmip2.f90

@@ -1 +1 @@
 MODULE etat0_dcmip2_mod
-
 ! test cases DCMIP 2012, category 2 : Orographic gravity waves
-
   USE icosa
+  IMPLICIT NONE
   PRIVATE
 
@@ -33 +32 @@
   SUBROUTINE compute_etat0(ngrid,lon,lat, phis, ps, Temp, ulon, ulat)
     USE disvert_mod
-    IMPLICIT NONE
+    USE omp_para
     INTEGER, INTENT(IN)    :: ngrid
     REAL(rstd),INTENT(IN)  :: lon(ngrid)
@@ -46 +45 @@
 
     ! Hexagons : ps,phis,temp
-    DO l=1,llm
+    DO l=ll_begin,ll_end
        ! The surface pressure is not set yet so we provide the hybrid coefficients
        hyam = .5*(ap(l)+ap(l+1))/preff