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

Timestamp:

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

Author:

dubos

Message:

OpenMP fixes for DCMIP

File:

• codes/icosagcm/trunk/src/etat0_dcmip1.f90 (modified) (5 diffs)

codes/icosagcm/trunk/src/etat0_dcmip1.f90

@@ -105 +105 @@
 
   SUBROUTINE compute_etat0_ncar(icase,ngrid,lon,lat, q)
-  USE icosa
   USE disvert_mod
-  IMPLICIT NONE  
+  USE omp_para
   INTEGER, INTENT(IN) :: icase, ngrid
   REAL(rstd),INTENT(IN) :: lon(ngrid),lat(ngrid)  
@@ -158 +157 @@
 
     SUBROUTINE cosine_bell_1(hx)
-    IMPLICIT NONE 
     REAL(rstd) :: hx(ngrid,llm)
     REAL(rstd) :: rr1,rr2   
     INTEGER :: n,l
-    DO l=1,llm 
+    DO l=ll_begin,ll_end 
        DO n=1,ngrid
           CALL dist_lonlat(lon0,lat0,lon(n),lat(n),rr1)  ! GC distance from center 
@@ -181 +179 @@
     REAL(rstd) :: rr1,rr2,dd1,dd2,dd1t1,dd1t2,dd2t1
     INTEGER :: n,l
-    DO l=1,llm
+    DO l=ll_begin,ll_end
        DO n=1,ngrid
           CALL dist_lonlat(lonc1,latc1,lon(n),lat(n),rr1)  ! GC distance from center 
@@ -208 +206 @@
     REAL(rstd) :: rr1,rr2,dd1,dd2,dd1t1,dd1t2,dd2t1
     INTEGER :: n,l
-    DO l=1,llm
+    DO l=ll_begin,ll_end
        DO n=1,ngrid
           CALL dist_lonlat(lonc1,latc1,lon(n),lat(n),rr1)  ! GC distance from center 
@@ -245 +243 @@
       INTEGER :: n,l
       
-      DO l=1,llm
+      DO l=ll_begin,ll_end
          IF ( ( zz1 .LT. zrl(l) ) .and. ( zrl(l) .LT. zz2 ) )  THEN 
             hx(:,l) = 0.5*(1. + cos(2*pi*(zrl(l)-zz0)/(zz2-zz1)))