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/OPA_SRC/DYN/dynzdf.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/OPA_SRC/DYN/dynzdf.F90

    r7698 r7753  
    5353      !! 
    5454      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
    55       INTEGER  ::   ji, jj, jk   ! dummy loop indices 
    5655      ! 
    5756      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv 
     
    6766      IF( l_trddyn )   THEN                      ! temporary save of ta and sa trends 
    6867         CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv )  
    69 !$OMP PARALLEL DO schedule(static) private(jk,jj,ji) 
    70          DO jk = 1, jpk 
    71             DO jj = 1, jpj 
    72                DO ji = 1, jpi 
    73                   ztrdu(ji,jj,jk) = ua(ji,jj,jk) 
    74                   ztrdv(ji,jj,jk) = va(ji,jj,jk) 
    75                END DO 
    76             END DO 
    77          END DO 
     68         ztrdu(:,:,:) = ua(:,:,:) 
     69         ztrdv(:,:,:) = va(:,:,:) 
    7870      ENDIF 
    7971 
     
    8678 
    8779      IF( l_trddyn )   THEN                      ! save the vertical diffusive trends for further diagnostics 
    88 !$OMP PARALLEL DO schedule(static) private(jk,jj,ji) 
    89          DO jk = 1, jpk 
    90             DO jj = 1, jpj 
    91                DO ji = 1, jpi 
    92                   ztrdu(ji,jj,jk) = ( ua(ji,jj,jk) - ub(ji,jj,jk) ) / r2dt - ztrdu(ji,jj,jk) 
    93                   ztrdv(ji,jj,jk) = ( va(ji,jj,jk) - vb(ji,jj,jk) ) / r2dt - ztrdv(ji,jj,jk) 
    94                END DO 
    95             END DO 
    96          END DO 
     80         ztrdu(:,:,:) = ( ua(:,:,:) - ub(:,:,:) ) / r2dt - ztrdu(:,:,:) 
     81         ztrdv(:,:,:) = ( va(:,:,:) - vb(:,:,:) ) / r2dt - ztrdv(:,:,:) 
    9782         CALL trd_dyn( ztrdu, ztrdv, jpdyn_zdf, kt ) 
    9883         CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv )  
Note: See TracChangeset for help on using the changeset viewer.