Changeset 2240 for branches/DEV_r2106_LOCEAN2010/NEMO/OFF_SRC/dtadyn.F90
- Timestamp:
- 2010-10-13T10:48:48+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2106_LOCEAN2010/NEMO/OFF_SRC/dtadyn.F90
r2082 r2240 568 568 DO jj = 2, jpjm1 569 569 DO ji = fs_2, fs_jpim1 ! vector opt. 570 #if defined key_zco571 zu = pu(ji ,jj ,jk) * umask(ji ,jj ,jk) * e2u(ji ,jj )572 zu1 = pu(ji-1,jj ,jk) * umask(ji-1,jj ,jk) * e2u(ji-1,jj )573 zv = pv(ji ,jj ,jk) * vmask(ji ,jj ,jk) * e1v(ji ,jj )574 zv1 = pv(ji ,jj-1,jk) * vmask(ji ,jj-1,jk) * e1v(ji ,jj-1)575 zet = 1. / ( e1t(ji,jj) * e2t(ji,jj) )576 #else577 570 zu = pu(ji ,jj ,jk) * umask(ji ,jj ,jk) * e2u(ji ,jj ) * fse3u(ji ,jj ,jk) 578 571 zu1 = pu(ji-1,jj ,jk) * umask(ji-1,jj ,jk) * e2u(ji-1,jj ) * fse3u(ji-1,jj ,jk) … … 580 573 zv1 = pv(ji ,jj-1,jk) * vmask(ji ,jj-1,jk) * e1v(ji ,jj-1) * fse3v(ji ,jj-1,jk) 581 574 zet = 1. / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) ) 582 #endif583 575 phdiv(ji,jj,jk) = ( zu - zu1 + zv - zv1 ) * zet 584 576 END DO
Note: See TracChangeset
for help on using the changeset viewer.