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 2122 – NEMO

Changeset 2122


Ignore:
Timestamp:
2010-09-24T17:38:54+02:00 (14 years ago)
Author:
cetlod
Message:

correction of minor bug in trabbl.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/trabbl.F90

    r2106 r2122  
    102102      ENDIF 
    103103 
    104       IF( l_bbl )   CALL bbl( kt, 'TRA' )       !* bbl coef. and transport (only if not already done trcbbl) 
     104      IF( l_bbl )   CALL bbl( kt, 'TRA' )       !* bbl coef. and transport (only if not already done in trcbbl) 
    105105 
    106106 
     
    166166      INTEGER  ::   ji, jj, jn   ! dummy loop indices 
    167167      INTEGER  ::   ik           ! local integers 
    168       REAL(wp) ::   zbtr, ztra   ! local scalars 
     168      REAL(wp) ::   zbtr         ! local scalars 
    169169      !!---------------------------------------------------------------------- 
    170170      ! 
     
    216216      !! 
    217217      INTEGER  ::   ji, jj, jk, jn           ! dummy loop indices 
    218       INTEGER  ::   ik                       ! local integers 
    219       INTEGER  ::   iis , iid , ijs , ijd    !   -       - 
     218      INTEGER  ::   iis , iid , ijs , ijd    ! local integers 
    220219      INTEGER  ::   ikus, ikud, ikvs, ikvd   !   -       - 
    221220      REAL(wp) ::   zbtr, ztra               ! local scalars 
     
    228227# if defined key_vectopt_loop 
    229228         DO jj = 1, 1 
    230             DO ji = jpi+1, jpij-jpi-1   ! vector opt. (forced unrolling) 
     229            DO ji = 1, jpij-jpi-1   ! vector opt. (forced unrolling) 
    231230# else 
    232          DO jj = 2, jpjm1 
     231         DO jj = 1, jpjm1 
    233232            DO ji = 1, jpim1            ! CAUTION start from i=1 to update i=2 when cyclic east-west 
    234233# endif 
     
    418417         ! 
    419418         CASE( 1 )                                   != use of upper velocity 
    420             DO jj = 2, jpjm1                                 ! criteria: grad(rho).grad(h)<0  and grad(rho).grad(h)<0 
     419            DO jj = 1, jpjm1                                 ! criteria: grad(rho).grad(h)<0  and grad(rho).grad(h)<0 
    421420               DO ji = 1, fs_jpim1   ! vector opt. 
    422421                  !                                               ! i-direction 
     
    451450         CASE( 2 )                                 != bbl velocity = F( delta rho ) 
    452451            zgbbl = grav * rn_gambbl 
    453             DO jj = 2, jpjm1                            ! criteria: rho_up > rho_down 
     452            DO jj = 1, jpjm1                            ! criteria: rho_up > rho_down 
    454453               DO ji = 1, fs_jpim1   ! vector opt. 
    455454                  !                                         ! i-direction 
     
    555554      zmbk(:,:) = FLOAT( mbkv_d(:,:) )   ;   CALL lbc_lnk(zmbk,'V',1.)   ;   mbkv_d(:,:) = MAX( INT( zmbk(:,:) ), 1 ) 
    556555 
    557       DO jj = 1, jpj                !* sign of grad(H) at u- and v-points 
    558          DO ji = 1, jpi 
     556                                        !* sign of grad(H) at u- and v-points 
     557      mgrhu(jpi,:) = 0.    ;    mgrhu(:,jpj) = 0.   ;    mgrhv(jpi,:) = 0.    ;    mgrhv(:,jpj) = 0. 
     558      DO jj = 1, jpjm1                 
     559         DO ji = 1, jpim1 
    559560            mgrhu(ji,jj) = INT(  SIGN( 1.e0, fsdept_0(ji+1,jj,mbkt(ji+1,jj)) - fsdept_0(ji,jj,mbkt(ji,jj)) )  ) 
    560561            mgrhv(ji,jj) = INT(  SIGN( 1.e0, fsdept_0(ji,jj+1,mbkt(ji,jj+1)) - fsdept_0(ji,jj,mbkt(ji,jj)) )  ) 
Note: See TracChangeset for help on using the changeset viewer.