- Timestamp:
- 2015-09-13T09:42:41+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilap.F90
r5147 r5737 112 112 DO jj = 1, jpjm1 113 113 DO ji = 1, fs_jpim1 ! vector opt. 114 zeeu(ji,jj) = re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk) * umask(ji,jj,jk)115 zeev(ji,jj) = re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk) * vmask(ji,jj,jk)114 zeeu(ji,jj) = e2_e1u(ji,jj) * fse3u_n(ji,jj,jk) * umask(ji,jj,jk) 115 zeev(ji,jj) = e1_e2v(ji,jj) * fse3v_n(ji,jj,jk) * vmask(ji,jj,jk) 116 116 END DO 117 117 END DO … … 145 145 DO jj = 2, jpjm1 ! Second derivative (divergence) time the eddy diffusivity coefficient 146 146 DO ji = fs_2, fs_jpim1 ! vector opt. 147 zbtr = 1.0 / ( e1 2t(ji,jj) * fse3t_n(ji,jj,jk) )147 zbtr = 1.0 / ( e1e2t(ji,jj) * fse3t_n(ji,jj,jk) ) 148 148 zlt(ji,jj) = fsahtt(ji,jj,jk) * zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) & 149 149 & + ztv(ji,jj,jk) - ztv(ji,jj-1,jk) ) … … 163 163 DO ji = fs_2, fs_jpim1 ! vector opt. 164 164 ! horizontal diffusive trends 165 zbtr = 1.0 / ( e1 2t(ji,jj) * fse3t_n(ji,jj,jk) )165 zbtr = 1.0 / ( e1e2t(ji,jj) * fse3t_n(ji,jj,jk) ) 166 166 ztra = zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji,jj-1,jk) ) 167 167 ! add it to the general tracer trends
Note: See TracChangeset
for help on using the changeset viewer.