Changeset 7525 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

Timestamp:

2017-01-04T17:47:47+01:00 (7 years ago)

Author:

mocavero

Message:

changes after review

File:

• branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90 (modified) (2 diffs)

branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

@@ -108 +108 @@
       IF( l_trdtrc )  THEN             ! trends: store now fields before the Asselin filter application
          CALL wrk_alloc( jpi, jpj, jpk, jptra, ztrdt )
+!$OMP PARALLEL DO schedule(static) private(jn,jk,jj,ji)
          DO jn = 1, jptra
-!$OMP PARALLEL DO schedule(static) private(jk,jj,ji)
             DO jk = 1, jpk
                DO jj = 1, jpj
@@ -121 +121 @@
       !                                ! Leap-Frog + Asselin filter time stepping
       IF( neuler == 0 .AND. kt == nittrc000 ) THEN    ! Euler time-stepping at first time-step (only swap)
+!$OMP PARALLEL DO schedule(static) private(jn,jk,jj,ji)
          DO jn = 1, jptra
-!$OMP PARALLEL DO schedule(static) private(jk)
             DO jk = 1, jpkm1
-               trn(:,:,jk,jn) = tra(:,:,jk,jn)
+               DO jj = 1, jpj
+                  DO ji = 1, jpi
+                     trn(ji,jj,jk,jn) = tra(ji,jj,jk,jn)
+                  END DO
+               END DO
             END DO
          END DO