- Timestamp:
- 2018-09-12T15:59:13+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90
r7806 r10115 296 296 ! 297 297 IF ( iom_use("eken") ) THEN 298 rke(:,:,j k) = 0._wp ! kinetic energy298 rke(:,:,jpk) = 0._wp ! kinetic energy 299 299 DO jk = 1, jpkm1 300 300 DO jj = 2, jpjm1 301 301 DO ji = fs_2, fs_jpim1 ! vector opt. 302 zztmp = 1._wp / ( e1e2t(ji,jj) * fse3t(ji,jj,jk))303 zztmpx = 0.5 * ( un(ji-1,jj,jk) * un(ji-1,jj,jk) * e 2u(ji-1,jj) * fse3u(ji-1,jj,jk) &304 & + un(ji ,jj,jk) * un(ji ,jj,jk) * e 2u(ji ,jj) * fse3u(ji ,jj,jk) ) &302 zztmp = 1 / (e1e2t(ji,jj) * fse3t(ji,jj,jk)) 303 zztmpx = 0.5 * ( un(ji-1,jj,jk) * un(ji-1,jj,jk) * e1u(ji-1,jj) * e2u(ji-1,jj) * fse3u(ji-1,jj,jk) & 304 & + un(ji ,jj,jk) * un(ji ,jj,jk) * e1u(ji, jj) * e2u(ji ,jj) * fse3u(ji ,jj,jk) ) & 305 305 & * zztmp 306 306 ! 307 zztmpy = 0.5 * ( vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1v(ji,jj-1) * fse3v(ji,jj-1,jk) &308 & + vn(ji,jj ,jk) * vn(ji,jj ,jk) * e1v(ji,jj ) * fse3v(ji,jj ,jk) ) &307 zztmpy = 0.5 * ( vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1v(ji,jj-1) * e2v(ji,jj-1) * fse3v(ji,jj-1,jk) & 308 & + vn(ji,jj ,jk) * vn(ji,jj ,jk) * e1v(ji,jj ) * e2v(ji,jj ) * fse3v(ji,jj ,jk) ) & 309 309 & * zztmp 310 310 !
Note: See TracChangeset
for help on using the changeset viewer.