- 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/dynadv_ubs.F90
r4153 r4616 183 183 DO jj = 2, jpjm1 ! divergence of horizontal momentum fluxes 184 184 DO ji = fs_2, fs_jpim1 ! vector opt. 185 zbu = e1 u(ji,jj) *e2u(ji,jj) * fse3u(ji,jj,jk)186 zbv = e1 v(ji,jj) *e2v(ji,jj) * fse3v(ji,jj,jk)185 zbu = e1e2u(ji,jj) * fse3u(ji,jj,jk) 186 zbv = e1e2v(ji,jj) * fse3v(ji,jj,jk) 187 187 ! 188 188 ua(ji,jj,jk) = ua(ji,jj,jk) - ( zfu_t(ji+1,jj ,jk) - zfu_t(ji ,jj ,jk) & … … 205 205 DO jk = 1, jpkm1 ! ==================== ! 206 206 ! ! Vertical volume fluxesÊ 207 zfw(:,:,jk) = 0.25 * e1 t(:,:) *e2t(:,:) * wn(:,:,jk)207 zfw(:,:,jk) = 0.25 * e1e2t(:,:) * wn(:,:,jk) 208 208 ! 209 209 IF( jk == 1 ) THEN ! surface/bottom advective fluxes … … 235 235 DO ji = fs_2, fs_jpim1 ! vector opt. 236 236 ua(ji,jj,jk) = ua(ji,jj,jk) - ( zfu_uw(ji,jj,jk) - zfu_uw(ji,jj,jk+1) ) & 237 & / ( e1 u(ji,jj) *e2u(ji,jj) * fse3u(ji,jj,jk) )237 & / ( e1e2u(ji,jj) * fse3u(ji,jj,jk) ) 238 238 va(ji,jj,jk) = va(ji,jj,jk) - ( zfv_vw(ji,jj,jk) - zfv_vw(ji,jj,jk+1) ) & 239 & / ( e1 v(ji,jj) *e2v(ji,jj) * fse3v(ji,jj,jk) )239 & / ( e1e2v(ji,jj) * fse3v(ji,jj,jk) ) 240 240 END DO 241 241 END DO
Note: See TracChangeset
for help on using the changeset viewer.