Changeset 8464
- Timestamp:
- 2017-08-25T10:31:40+02:00 (7 years ago)
- Location:
- branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90
r7494 r8464 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 = r1_e1e2t(ji,jj) / fse3t(ji,jj,jk) 303 zztmpx = 0.5 * ( un(ji-1,jj,jk) * un(ji-1,jj,jk) * e1e2u(ji-1,jj) * fse3u(ji-1,jj,jk) & 304 & + un(ji ,jj,jk) * un(ji ,jj,jk) * e1e2u(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) * e1 v(ji,jj-1) * fse3v(ji,jj-1,jk) &308 & + vn(ji,jj ,jk) * vn(ji,jj ,jk) * e1 v(ji,jj ) * fse3v(ji,jj ,jk) ) &307 zztmpy = 0.5 * ( vn(ji,jj-1,jk) * vn(ji,jj-1,jk) * e1e2v(ji,jj-1) * fse3v(ji,jj-1,jk) & 308 & + vn(ji,jj ,jk) * vn(ji,jj ,jk) * e1e2v(ji,jj ) * fse3v(ji,jj ,jk) ) & 309 309 & * zztmp 310 310 ! -
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/step.F90
r7494 r8464 234 234 IF( lk_diaar5 ) CALL dia_ar5( kstp ) ! ar5 diag 235 235 IF( lk_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis 236 236 IF( lk_diaprod ) CALL dia_prod( kstp ) ! ocean model: product diagnostics 237 237 CALL dia_wri( kstp ) ! ocean model: outputs 238 238 !
Note: See TracChangeset
for help on using the changeset viewer.