- Timestamp:
- 2020-11-27T17:26:33+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/tickets_icb_1900
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/tickets_icb_1900
- Property svn:externals
-
NEMO/branches/2020/tickets_icb_1900/src/OCE/TRD/trdglo.F90
r13237 r13899 86 86 ! 87 87 CASE( 'TRA' ) !== Tracers (T & S) ==! 88 DO_3D _11_11( 1, jpkm1)88 DO_3D( 1, 1, 1, 1, 1, jpkm1 ) ! global sum of mask volume trend and trend*T (including interior mask) 89 89 zvm = e1e2t(ji,jj) * e3t(ji,jj,jk,Kmm) * tmask(ji,jj,jk) * tmask_i(ji,jj) 90 90 zvt = ptrdx(ji,jj,jk) * zvm … … 115 115 ! 116 116 CASE( 'DYN' ) !== Momentum and KE ==! 117 DO_3D _10_10(1, jpkm1 )117 DO_3D( 1, 0, 1, 0, 1, jpkm1 ) 118 118 zvt = ptrdx(ji,jj,jk) * tmask_i(ji+1,jj) * tmask_i(ji,jj) * umask(ji,jj,jk) & 119 119 & * e1e2u (ji,jj) * e3u(ji,jj,jk,Kmm) … … 127 127 IF( ktrd == jpdyn_zdf ) THEN ! zdf trend: compute separately the surface forcing trend 128 128 z1_2rho0 = 0.5_wp / rho0 129 DO_2D _10_10129 DO_2D( 1, 0, 1, 0 ) 130 130 zvt = ( utau_b(ji,jj) + utau(ji,jj) ) * tmask_i(ji+1,jj) * tmask_i(ji,jj) * umask(ji,jj,jk) & 131 131 & * z1_2rho0 * e1e2u(ji,jj) … … 211 211 212 212 zcof = 0.5_wp / rho0 ! Density flux at u and v-points 213 DO_3D _10_10(1, jpkm1 )213 DO_3D( 1, 0, 1, 0, 1, jpkm1 ) 214 214 zkx(ji,jj,jk) = zcof * e2u(ji,jj) * e3u(ji,jj,jk,Kmm) & 215 215 & * uu(ji,jj,jk,Kmm) * ( rhop(ji,jj,jk) + rhop(ji+1,jj,jk) ) … … 218 218 END_3D 219 219 220 DO_3D _00_00( 1, jpkm1 )220 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) ! Density flux divergence at t-point 221 221 zkepe(ji,jj,jk) = - ( zkz(ji,jj,jk) - zkz(ji ,jj ,jk+1) & 222 222 & + zkx(ji,jj,jk) - zkx(ji-1,jj ,jk ) & … … 527 527 tvolv = 0._wp 528 528 529 DO_3D _00_00(1, jpk )529 DO_3D( 0, 0, 0, 0, 1, jpk ) 530 530 tvolu = tvolu + e1u(ji,jj) * e2u(ji,jj) * e3u(ji,jj,jk,Kmm) & 531 531 & * tmask_i(ji+1,jj ) * tmask_i(ji,jj) * umask(ji,jj,jk)
Note: See TracChangeset
for help on using the changeset viewer.