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/PISCES/P4Z/p4zint.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/PISCES/P4Z/p4zint.F90

    r7698 r7753  
    2121   !!---------------------------------------------------------------------- 
    2222   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    23    !! $Id$ 
     23   !! $Id$  
    2424   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    2525   !!---------------------------------------------------------------------- 
     
    3636      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
    3737      ! 
    38       INTEGER  :: ji, jj, jk             ! dummy loop indices 
     38      INTEGER  :: ji, jj                 ! dummy loop indices 
    3939      REAL(wp) :: zvar                   ! local variable 
    4040      !!--------------------------------------------------------------------- 
     
    4444      ! Computation of phyto and zoo metabolic rate 
    4545      ! ------------------------------------------- 
    46 !$OMP PARALLEL 
    47 !$OMP DO schedule(static) private(jk,jj,ji) 
    48       DO jk = 1, jpk 
    49          DO jj = 1, jpj 
    50             DO ji = 1, jpi 
    51                tgfunc (ji,jj,jk) = EXP( 0.063913 * tsn(ji,jj,jk,jp_tem) ) 
    52                tgfunc2(ji,jj,jk) = EXP( 0.07608  * tsn(ji,jj,jk,jp_tem) ) 
    53             END DO 
    54          END DO 
    55       END DO 
     46      tgfunc (:,:,:) = EXP( 0.063913 * tsn(:,:,:,jp_tem) ) 
     47      tgfunc2(:,:,:) = EXP( 0.07608  * tsn(:,:,:,jp_tem) ) 
    5648 
    5749      ! Computation of the silicon dependant half saturation  constant for silica uptake 
    5850      ! --------------------------------------------------- 
    59 !$OMP DO schedule(static) private(jj,ji,zvar) 
    6051      DO ji = 1, jpi 
    6152         DO jj = 1, jpj 
     
    6657      ! 
    6758      IF( nday_year == nyear_len(1) ) THEN 
    68 !$OMP DO schedule(static) private(jj,ji) 
    69          DO jj = 1, jpj 
    70             DO ji = 1, jpi 
    71                xksi   (ji,jj) = xksimax(ji,jj) 
    72                xksimax(ji,jj) = 0._wp 
    73             END DO 
    74          END DO 
     59         xksi   (:,:) = xksimax(:,:) 
     60         xksimax(:,:) = 0._wp 
    7561      ENDIF 
    76 !$OMP END PARALLEL 
    7762      ! 
    7863      IF( nn_timing == 1 )  CALL timing_stop('p4z_int') 
Note: See TracChangeset for help on using the changeset viewer.