- Timestamp:
- 2010-10-13T10:48:48+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynldf_lap.F90
r1954 r2240 50 50 !! difu = 1/e1u di[ahmt hdivb] - 1/(e2u*e3u) dj-1[e3f ahmf rotb] 51 51 !! difv = 1/e2v dj[ahmt hdivb] + 1/(e1v*e3v) di-1[e3f ahmf rotb] 52 !! If lk_zco=T, e3f=e3u=e3v, the vertical scale factor are simplified53 52 !! in the rotational part of the diffusion. 54 53 !! Add this before trend to the general trend (ua,va): … … 87 86 DO jj = 2, jpjm1 88 87 DO ji = fs_2, fs_jpim1 ! vector opt. 89 #if defined key_zco90 ! horizontal diffusive trends91 ze2u = rotb (ji,jj,jk)*fsahmf(ji,jj,jk)92 ze1v = hdivb(ji,jj,jk)*fsahmt(ji,jj,jk)93 zua = - ( ze2u - rotb (ji,jj-1,jk)*fsahmf(ji,jj-1,jk) ) / e2u(ji,jj) &94 + ( hdivb(ji+1,jj,jk)*fsahmt(ji+1,jj,jk) - ze1v ) / e1u(ji,jj)95 96 zva = + ( ze2u - rotb (ji-1,jj,jk)*fsahmf(ji-1,jj,jk) ) / e1v(ji,jj) &97 + ( hdivb(ji,jj+1,jk)*fsahmt(ji,jj+1,jk) - ze1v ) / e2v(ji,jj)98 #else99 88 ze2u = rotb (ji,jj,jk)*fsahmf(ji,jj,jk)*fse3f(ji,jj,jk) 100 89 ze1v = hdivb(ji,jj,jk)*fsahmt(ji,jj,jk) … … 105 94 zva = + ( ze2u - rotb (ji-1,jj,jk)*fsahmf(ji-1,jj,jk)*fse3f(ji-1,jj,jk) ) / ( e1v(ji,jj) * fse3v(ji,jj,jk) ) & 106 95 + ( hdivb(ji,jj+1,jk)*fsahmt(ji,jj+1,jk) - ze1v ) / e2v(ji,jj) 107 #endif108 96 109 97 ! add it to the general momentum trends
Note: See TracChangeset
for help on using the changeset viewer.