Changeset 15145 for NEMO/trunk/src/OCE/DYN/divhor.F90
- Timestamp:
- 2021-07-26T18:16:45+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DYN/divhor.F90
r15058 r15145 79 79 ENDIF 80 80 ! 81 DO_3D_OVR( nn_hls-1, nn_hls , nn_hls-1, nn_hls, 1, jpkm1 )!== Horizontal divergence ==!81 DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpkm1 ) !== Horizontal divergence ==! 82 82 ! round brackets added to fix the order of floating point operations 83 83 ! needed to ensure halo 1 - halo 2 compatibility 84 hdiv(ji,jj,jk) = ( ( e2u(ji ,jj) * e3u(ji ,jj,jk,Kmm) * uu(ji ,jj,jk,Kmm) & 85 & - e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm) * uu(ji-1,jj,jk,Kmm) & 84 hdiv(ji,jj,jk) = ( ( e2u(ji ,jj) * e3u(ji ,jj,jk,Kmm) * uu(ji ,jj,jk,Kmm) & ! Warning: with qco, e3u uses r3u that 85 & - e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm) * uu(ji-1,jj,jk,Kmm) & ! is not defined in jpi 86 86 & ) & ! bracket for halo 1 - halo 2 compatibility 87 & + ( e1v(ji,jj ) * e3v(ji,jj ,jk,Kmm) * vv(ji,jj ,jk,Kmm) & 88 & - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vv(ji,jj-1,jk,Kmm) & 87 & + ( e1v(ji,jj ) * e3v(ji,jj ,jk,Kmm) * vv(ji,jj ,jk,Kmm) & ! Warning: with qco, e3v uses r3v that 88 & - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vv(ji,jj-1,jk,Kmm) & ! is not defined in jpj 89 89 & ) & ! bracket for halo 1 - halo 2 compatibility 90 90 & ) * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm)
Note: See TracChangeset
for help on using the changeset viewer.