Changeset 4032 for branches/2013/dev_r4017_CMCC_MPI_OpenMP/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

Timestamp:

2013-09-21T16:15:51+02:00 (11 years ago)

Author:

epico

Message:

The hybrid parallel version OpenMP+MPI has been introduced for the GYRE configuration. A further refinement is needed to remove dependencies in those loops not parallelized yet (whenever possible)

File:

• branches/2013/dev_r4017_CMCC_MPI_OpenMP/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90 (modified) (2 diffs)

branches/2013/dev_r4017_CMCC_MPI_OpenMP/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

@@ -136 +136 @@
       IF( neuler == 0 .AND. kt == nit000 ) THEN       ! Euler time-stepping at first time-step (only swap)
          DO jn = 1, jpts
+!$OMP PARALLEL DO schedule(static) private(jk)
             DO jk = 1, jpkm1
                tsn(:,:,jk,jn) = tsa(:,:,jk,jn)    
@@ -154 +155 @@
       ! trends computation
       IF( l_trdtra ) THEN      ! trend of the Asselin filter (tb filtered - tb)/dt     
+!$OMP PARALLEL DO schedule(static) private(jk, zfact)
          DO jk = 1, jpkm1
             zfact = 1.e0 / r2dtra(jk)