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 457 for trunk/NEMO/OPA_SRC/TRA/trabbl_adv.h90 – NEMO

Ignore:
Timestamp:
2006-05-10T19:01:19+02:00 (18 years ago)
Author:
opalod
Message:

nemo_v1_update_049:RB: reorganization of tracers part, remove traadv_cen2_atsk.h90 traldf_iso_zps.F90 trazdf_iso.F90 trazdf_iso_vopt.F90, change atsk routines to jki

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/TRA/trabbl_adv.h90

    r409 r457  
    7373         zgdrho, zbtr, zta, zsa            !    "         "  
    7474      REAL(wp), DIMENSION(jpi,jpj) ::   & 
    75          ztnb, zsnb, zdep, ztbb, zsbb,  &  ! temporary workspace arrays 
     75         ztnb, zsnb, zdep, ztbb, zsbb,  &  !    "                  " 
    7676         zahu, zahv                        !    "                  " 
    7777      REAL(wp), DIMENSION(jpi,jpj) ::   &  ! temporary workspace arrays 
     
    120120      ! mbathy= number of w-level, minimum value=1 (cf dommsk.F) 
    121121 
    122 #if defined key_vectopt_loop   &&   ! defined key_autotasking 
     122#if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
    123123      jj = 1 
    124124      DO ji = 1, jpij   ! vector opt. (forced unrolling) 
     
    133133            zsbb(ji,jj) = sb(ji,jj,ik) * tmask(ji,jj,1)    ! masked before S at the ocean bottom 
    134134            zdep(ji,jj) = fsdept(ji,jj,ik)                 ! depth of the ocean bottom T-level 
    135 #if ! defined key_vectopt_loop   ||   defined key_autotasking 
    136          END DO 
    137 #endif 
    138       END DO 
    139 #if defined key_vectopt_loop   &&   ! defined key_autotasking  
     135#if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
     136         END DO 
     137#endif 
     138      END DO 
     139#if defined key_vectopt_loop   &&   ! defined key_mpp_omp  
    140140      jj = 1 
    141141      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    260260          u_bbl(:,:,:) = 0.e0 
    261261          v_bbl(:,:,:) = 0.e0 
    262 # if defined key_vectopt_loop   &&   ! defined key_autotasking 
     262# if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
    263263      jj = 1 
    264264      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    273273               v_bbl(ji,jj,ikv) = zalphay(ji,jj) * vn(ji,jj,ikv) * vmask(ji,jj,1) 
    274274            ENDIF 
    275 # if ! defined key_vectopt_loop   ||   defined key_autotasking 
     275# if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    276276         END DO 
    277277# endif 
     
    285285      ! ... Second order centered tracer flux at u and v-points 
    286286 
    287 # if defined key_vectopt_loop   &&   ! defined key_autotasking 
     287# if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
    288288      jj = 1 
    289289      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    310310            zwz(ji,jj) = ( ( zfvj + ABS( zfvj ) ) * zsbb(ji  ,jj  )   & 
    311311               &          +( zfvj - ABS( zfvj ) ) * zsbb(ji  ,jj+1) ) * 0.5 
    312 #if ! defined key_vectopt_loop   ||   defined key_autotasking 
     312#if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    313313         END DO 
    314314#endif 
    315315        END DO 
    316 # if defined key_vectopt_loop   &&   ! defined key_autotasking 
     316# if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
    317317      jj = 1 
    318318      DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    332332            ta(ji,jj,ik) = ta(ji,jj,ik) + zta 
    333333            sa(ji,jj,ik) = sa(ji,jj,ik) + zsa 
    334 #if ! defined key_vectopt_loop   ||   defined key_autotasking 
     334#if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    335335         END DO 
    336336#endif 
     
    340340      ! BBL lateral advection tracers trends 
    341341      IF( l_trdtra )   THEN 
    342 #  if defined key_vectopt_loop   &&   ! defined key_autotasking 
     342#  if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
    343343         jj = 1 
    344344         DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    350350            tladbbl(ji,jj) = ta(ji,jj,ik) - ztdta(ji,jj,ik) 
    351351            sladbbl(ji,jj) = sa(ji,jj,ik) - ztdsa(ji,jj,ik) 
    352 #  if ! defined key_vectopt_loop   ||   defined key_autotasking 
     352#  if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    353353            END DO 
    354354#  endif 
     
    386386 
    387387      ! ... horizontal bottom divergence 
    388 # if defined key_vectopt_loop   &&   ! defined key_autotasking 
     388# if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
    389389      jj = 1 
    390390      DO ji = 1, jpij-jpi   ! vector opt. (forced unrolling) 
     
    397397            zwu(ji,jj) = zalphax(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,iku)  
    398398            zwv(ji,jj) = zalphay(ji,jj) * e1v(ji,jj) * fse3v(ji,jj,ikv)  
    399 #if ! defined key_vectopt_loop   ||   defined key_autotasking 
     399#if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    400400         END DO 
    401401#endif 
    402402        END DO 
    403403 
    404 # if defined key_vectopt_loop   &&   ! defined key_autotasking 
     404# if defined key_vectopt_loop   &&   ! defined key_mpp_omp 
    405405      jj = 1 
    406406      DO ji = jpi+2, jpij-jpi-1   ! vector opt. (forced unrolling) 
     
    418418               &   ) / zbt 
    419419 
    420 # if ! defined key_vectopt_loop   ||   defined key_autotasking 
     420# if ! defined key_vectopt_loop   ||   defined key_mpp_omp 
    421421         END DO 
    422422# endif 
Note: See TracChangeset for help on using the changeset viewer.