- Timestamp:
- 2013-04-09T18:34:38+02:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3858_NOC_ZTC/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilap.F90
r3294 r3865 110 110 DO jj = 1, jpjm1 111 111 DO ji = 1, fs_jpim1 ! vector opt. 112 zeeu(ji,jj) = e2u(ji,jj) * fse3u(ji,jj,jk) / e1u(ji,jj) * umask(ji,jj,jk)113 zeev(ji,jj) = e1v(ji,jj) * fse3v(ji,jj,jk) / e2v(ji,jj) * vmask(ji,jj,jk)112 zeeu(ji,jj) = re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk) * umask(ji,jj,jk) 113 zeev(ji,jj) = re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk) * vmask(ji,jj,jk) 114 114 END DO 115 115 END DO … … 133 133 DO jj = 2, jpjm1 ! Second derivative (divergence) time the eddy diffusivity coefficient 134 134 DO ji = fs_2, fs_jpim1 ! vector opt. 135 zbtr = 1.0 / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )135 zbtr = 1.0 / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 136 136 zlt(ji,jj) = fsahtt(ji,jj,jk) * zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) & 137 137 & + ztv(ji,jj,jk) - ztv(ji,jj-1,jk) ) … … 151 151 DO ji = fs_2, fs_jpim1 ! vector opt. 152 152 ! horizontal diffusive trends 153 zbtr = 1.0 / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )153 zbtr = 1.0 / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 154 154 ztra = zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji,jj-1,jk) ) 155 155 ! add it to the general tracer trends
Note: See TracChangeset
for help on using the changeset viewer.