Changeset 14072 for NEMO/trunk/src/OCE/TRA/trabbl.F90
- Timestamp:
- 2020-12-04T08:48:38+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/TRA/trabbl.F90
r13982 r14072 31 31 USE trdtra ! trends: active tracers 32 32 ! 33 USE iom ! IOM library 33 USE iom ! IOM library 34 34 USE in_out_manager ! I/O manager 35 35 USE lbclnk ! ocean lateral boundary conditions 36 36 USE prtctl ! Print control 37 37 USE timing ! Timing 38 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 38 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 39 39 40 40 IMPLICIT NONE … … 200 200 zptb(ji,jj) = pt(ji,jj,ik,jn) ! bottom before T and S 201 201 END_2D 202 ! 202 ! 203 203 DO_2D( 0, 0, 0, 0 ) ! Compute the trend 204 204 ik = mbkt(ji,jj) ! bottom T-level index … … 399 399 za = zab(ji+1,jj,jp_tem) + zab(ji,jj,jp_tem) ! 2*(alpha,beta) at u-point 400 400 zb = zab(ji+1,jj,jp_sal) + zab(ji,jj,jp_sal) 401 ! ! 2*masked bottom density gradient 401 ! ! 2*masked bottom density gradient 402 402 zgdrho = ( za * ( zts(ji+1,jj,jp_tem) - zts(ji,jj,jp_tem) ) & 403 403 - zb * ( zts(ji+1,jj,jp_sal) - zts(ji,jj,jp_sal) ) ) * umask(ji,jj,1) … … 523 523 END_2D 524 524 ! converte into REAL to use lbc_lnk ; impose a min value of 1 as a zero can be set in lbclnk 525 zmbku(:,:) = REAL( mbku_d(:,:), wp ) ; zmbkv(:,:) = REAL( mbkv_d(:,:), wp ) 526 CALL lbc_lnk_multi( 'trabbl', zmbku,'U',1.0_wp, zmbkv,'V',1.0_wp) 525 zmbku(:,:) = REAL( mbku_d(:,:), wp ) ; zmbkv(:,:) = REAL( mbkv_d(:,:), wp ) 526 CALL lbc_lnk_multi( 'trabbl', zmbku,'U',1.0_wp, zmbkv,'V',1.0_wp) 527 527 mbku_d(:,:) = MAX( INT( zmbku(:,:) ), 1 ) ; mbkv_d(:,:) = MAX( NINT( zmbkv(:,:) ), 1 ) 528 528 !
Note: See TracChangeset
for help on using the changeset viewer.