- Timestamp:
- 2016-11-21T09:55:07+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90
r6140 r7277 298 298 ! 299 299 INTEGER :: ji, jj, jk ! dummy loop indices 300 REAL(wp) :: zaht, zah t_min, z1_f20 ! local scalar300 REAL(wp) :: zaht, zahf, zaht_min, z1_f20 ! local scalar 301 301 !!---------------------------------------------------------------------- 302 302 ! … … 327 327 DO jj = 1, jpj 328 328 DO ji = 1, jpi 329 zaht = ( 1._wp - MIN( 1._wp , ABS( ff(ji,jj) * z1_f20 ) ) ) * ( rn_aht_0 - zaht_min ) 329 !!gm CAUTION : here we assume lat/lon grid in 20deg N/S band (like all ORCA cfg) 330 !! ==>>> The Coriolis value is identical for t- & u_points, and for v- and f-points 331 zaht = ( 1._wp - MIN( 1._wp , ABS( ff_t(ji,jj) * z1_f20 ) ) ) * ( rn_aht_0 - zaht_min ) 332 zahf = ( 1._wp - MIN( 1._wp , ABS( ff_f(ji,jj) * z1_f20 ) ) ) * ( rn_aht_0 - zaht_min ) 330 333 ahtu(ji,jj,1) = ( MAX( zaht_min, ahtu(ji,jj,1) ) + zaht ) * umask(ji,jj,1) ! min value zaht_min 331 ahtv(ji,jj,1) = ( MAX( zaht_min, ahtv(ji,jj,1) ) + zah t) * vmask(ji,jj,1) ! increase within 20S-20N334 ahtv(ji,jj,1) = ( MAX( zaht_min, ahtv(ji,jj,1) ) + zahf ) * vmask(ji,jj,1) ! increase within 20S-20N 332 335 END DO 333 336 END DO … … 555 558 END DO 556 559 557 !!gm IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA R2558 !!gm DO jj = 2, jpjm1559 !!gm DO ji = fs_2, fs_jpim1 ! vector opt.560 !!gm ! Take the minimum between aeiw and 1000 m2/s over shelves (depth shallower than 650 m)561 !!gm IF( mbkt(ji,jj) <= 20 ) zaeiw(ji,jj) = MIN( zaeiw(ji,jj), 1000. )562 !!gm END DO563 !!gm END DO564 !!gm ENDIF565 566 560 ! !== Bound on eiv coeff. ==! 567 561 z1_f20 = 1._wp / ( 2._wp * omega * sin( rad * 20._wp ) ) 568 562 DO jj = 2, jpjm1 569 563 DO ji = fs_2, fs_jpim1 ! vector opt. 570 zzaei = MIN( 1._wp, ABS( ff (ji,jj) * z1_f20 ) ) * zaeiw(ji,jj) ! tropical decrease564 zzaei = MIN( 1._wp, ABS( ff_t(ji,jj) * z1_f20 ) ) * zaeiw(ji,jj) ! tropical decrease 571 565 zaeiw(ji,jj) = MIN( zzaei , paei0 ) ! Max value = paei0 572 566 END DO
Note: See TracChangeset
for help on using the changeset viewer.