- Timestamp:
- 2014-04-06T17:28:25+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90
r3294 r4616 47 47 !! 48 48 !! ** Method : The now vertical advection of momentum is given by: 49 !! w dz(u) = ua + 1/(e1 u*e2u*e3u) mk+1[ mi(e1t*e2t*wn) dk(un) ]50 !! w dz(v) = va + 1/(e1 v*e2v*e3v) mk+1[ mj(e1t*e2t*wn) dk(vn) ]49 !! w dz(u) = ua + 1/(e1e2u*e3u) mk+1[ mi(e1e2t*wn) dk(un) ] 50 !! w dz(v) = va + 1/(e1e2v*e3v) mk+1[ mj(e1e2t*wn) dk(vn) ] 51 51 !! Add this trend to the general trend (ua,va): 52 52 !! (ua,va) = (ua,va) + w dz(u,v) … … 83 83 DO jj = 2, jpj ! vertical fluxes 84 84 DO ji = fs_2, jpi ! vector opt. 85 zww(ji,jj) = 0.25 * e1 t(ji,jj) *e2t(ji,jj) * wn(ji,jj,jk)85 zww(ji,jj) = 0.25 * e1e2t(ji,jj) * wn(ji,jj,jk) 86 86 END DO 87 87 END DO … … 106 106 DO ji = fs_2, fs_jpim1 ! vector opt. 107 107 ! ! vertical momentum advective trends 108 zua = - ( zwuw(ji,jj,jk) + zwuw(ji,jj,jk+1) ) / ( e1 u(ji,jj) *e2u(ji,jj) * fse3u(ji,jj,jk) )109 zva = - ( zwvw(ji,jj,jk) + zwvw(ji,jj,jk+1) ) / ( e1 v(ji,jj) *e2v(ji,jj) * fse3v(ji,jj,jk) )108 zua = - ( zwuw(ji,jj,jk) + zwuw(ji,jj,jk+1) ) / ( e1e2u(ji,jj) * fse3u(ji,jj,jk) ) 109 zva = - ( zwvw(ji,jj,jk) + zwvw(ji,jj,jk+1) ) / ( e1e2v(ji,jj) * fse3v(ji,jj,jk) ) 110 110 ! ! add the trends to the general momentum trends 111 111 ua(ji,jj,jk) = ua(ji,jj,jk) + zua
Note: See TracChangeset
for help on using the changeset viewer.