New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7753 for trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trczdf.F90 – NEMO

Ignore:
Timestamp:
2017-03-03T12:46:59+01:00 (7 years ago)
Author:
mocavero
Message:

Reverting trunk to remove OpenMP

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trczdf.F90

    r7698 r7753  
    4040   !!---------------------------------------------------------------------- 
    4141   !! NEMO/TOP 3.7 , NEMO Consortium (2015) 
    42    !! $Id$ 
     42   !! $Id$  
    4343   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4444   !!---------------------------------------------------------------------- 
     
    5353      INTEGER, INTENT( in ) ::  kt      ! ocean time-step index 
    5454      ! 
    55       INTEGER               ::  jk, jn, jj, ji 
     55      INTEGER               ::  jk, jn 
    5656      CHARACTER (len=22)    :: charout 
    5757      REAL(wp), POINTER, DIMENSION(:,:,:,:) ::   ztrtrd   ! 4D workspace 
     
    6262      IF( l_trdtrc )  THEN 
    6363         CALL wrk_alloc( jpi, jpj, jpk, jptra, ztrtrd ) 
    64 !$OMP PARALLEL DO schedule(static) private(jn,jk,jj,ji) 
    65          DO jn = 1, jptra 
    66             DO jk = 1, jpk 
    67                DO jj = 1, jpj 
    68                   DO ji = 1, jpi 
    69                      ztrtrd(ji,jj,jk,jn)  = tra(ji,jj,jk,jn) 
    70                   END DO 
    71                END DO 
    72             END DO 
    73          END DO 
     64         ztrtrd(:,:,:,:)  = tra(:,:,:,:) 
    7465      ENDIF 
    7566 
     
    8172      IF( l_trdtrc )   THEN                      ! save the vertical diffusive trends for further diagnostics 
    8273         DO jn = 1, jptra 
    83 !$OMP PARALLEL DO schedule(static) private(jk,jj,ji) 
    8474            DO jk = 1, jpkm1 
    85                DO jj = 1, jpj 
    86                   DO ji = 1, jpi 
    87                      ztrtrd(ji,jj,jk,jn) = ( ( tra(ji,jj,jk,jn) - trb(ji,jj,jk,jn) ) / r2dttrc ) - ztrtrd(ji,jj,jk,jn) 
    88                   END DO 
    89                END DO 
     75               ztrtrd(:,:,jk,jn) = ( ( tra(:,:,jk,jn) - trb(:,:,jk,jn) ) / r2dttrc ) - ztrtrd(:,:,jk,jn) 
    9076            END DO 
    9177            CALL trd_tra( kt, 'TRC', jn, jptra_zdf, ztrtrd(:,:,:,jn) ) 
Note: See TracChangeset for help on using the changeset viewer.