Changeset 2122
- Timestamp:
- 2010-09-24T17:38:54+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/TRA/trabbl.F90
r2106 r2122 102 102 ENDIF 103 103 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) 105 105 106 106 … … 166 166 INTEGER :: ji, jj, jn ! dummy loop indices 167 167 INTEGER :: ik ! local integers 168 REAL(wp) :: zbtr , ztra! local scalars168 REAL(wp) :: zbtr ! local scalars 169 169 !!---------------------------------------------------------------------- 170 170 ! … … 216 216 !! 217 217 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 220 219 INTEGER :: ikus, ikud, ikvs, ikvd ! - - 221 220 REAL(wp) :: zbtr, ztra ! local scalars … … 228 227 # if defined key_vectopt_loop 229 228 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) 231 230 # else 232 DO jj = 2, jpjm1231 DO jj = 1, jpjm1 233 232 DO ji = 1, jpim1 ! CAUTION start from i=1 to update i=2 when cyclic east-west 234 233 # endif … … 418 417 ! 419 418 CASE( 1 ) != use of upper velocity 420 DO jj = 2, jpjm1 ! criteria: grad(rho).grad(h)<0 and grad(rho).grad(h)<0419 DO jj = 1, jpjm1 ! criteria: grad(rho).grad(h)<0 and grad(rho).grad(h)<0 421 420 DO ji = 1, fs_jpim1 ! vector opt. 422 421 ! ! i-direction … … 451 450 CASE( 2 ) != bbl velocity = F( delta rho ) 452 451 zgbbl = grav * rn_gambbl 453 DO jj = 2, jpjm1 ! criteria: rho_up > rho_down452 DO jj = 1, jpjm1 ! criteria: rho_up > rho_down 454 453 DO ji = 1, fs_jpim1 ! vector opt. 455 454 ! ! i-direction … … 555 554 zmbk(:,:) = FLOAT( mbkv_d(:,:) ) ; CALL lbc_lnk(zmbk,'V',1.) ; mbkv_d(:,:) = MAX( INT( zmbk(:,:) ), 1 ) 556 555 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 559 560 mgrhu(ji,jj) = INT( SIGN( 1.e0, fsdept_0(ji+1,jj,mbkt(ji+1,jj)) - fsdept_0(ji,jj,mbkt(ji,jj)) ) ) 560 561 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.