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

Timestamp:

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

Author:

dubos

Message:

OpenMP fixes for DCMIP

File:

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

codes/icosagcm/trunk/src/etat0_dcmip3.f90

@@ -11 +11 @@
     USE dcmip_initial_conditions_test_1_2_3
     USE disvert_mod
+    USE omp_para
     INTEGER, INTENT(IN) :: ngrid
     REAL(rstd), INTENT(IN) :: lon(ngrid)
@@ -28 +29 @@
             dummy,dummy,dummy,dummy,phis(ij),ps(ij),dummy,dummy)
     END DO
-    DO l=1,llm
+    DO l=ll_begin,ll_end
        DO ij=1,ngrid
           pp = .5*(ap(l)+ap(l+1)) + .5*(bp(l)+bp(l+1))*ps(ij) ! full-layer pressure
@@ -34 +35 @@
                ulon(ij,l),ulat(ij,l),dummy,Temp(ij,l),dummy,dummy,dummy,dummy)
        END DO
-    END DO    
-    q=0.
+    END DO
+    q(:,:,:)=0.
   END SUBROUTINE compute_etat0