Changeset 5845 for branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
- Timestamp:
- 2015-10-31T08:40:45+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
r5836 r5845 30 30 31 31 !! * Substitutions 32 # include "domzgr_substitute.h90"33 32 # include "vectopt_loop_substitute.h90" 34 33 !!---------------------------------------------------------------------- … … 77 76 DO jk = 1, jpkm1 ! ====================== ! 78 77 ! ! horizontal volume fluxes 79 zfu(:,:,jk) = 0.25 * e2u(:,:) * fse3u(:,:,jk) * un(:,:,jk)80 zfv(:,:,jk) = 0.25 * e1v(:,:) * fse3v(:,:,jk) * vn(:,:,jk)78 zfu(:,:,jk) = 0.25 * e2u(:,:) * e3u_n(:,:,jk) * un(:,:,jk) 79 zfv(:,:,jk) = 0.25 * e1v(:,:) * e3v_n(:,:,jk) * vn(:,:,jk) 81 80 ! 82 81 DO jj = 1, jpjm1 ! horizontal momentum fluxes at T- and F-point … … 90 89 DO jj = 2, jpjm1 ! divergence of horizontal momentum fluxes 91 90 DO ji = fs_2, fs_jpim1 ! vector opt. 92 zbu = e1e2u(ji,jj) * fse3u(ji,jj,jk)93 zbv = e1e2v(ji,jj) * fse3v(ji,jj,jk)91 zbu = e1e2u(ji,jj) * e3u_n(ji,jj,jk) 92 zbv = e1e2v(ji,jj) * e3v_n(ji,jj,jk) 94 93 ! 95 94 ua(ji,jj,jk) = ua(ji,jj,jk) - ( zfu_t(ji+1,jj ,jk) - zfu_t(ji ,jj ,jk) & … … 144 143 DO ji = fs_2, fs_jpim1 ! vector opt. 145 144 ua(ji,jj,jk) = ua(ji,jj,jk) - ( zfu_uw(ji,jj,jk) - zfu_uw(ji,jj,jk+1) ) & 146 & / ( e1e2u(ji,jj) * fse3u(ji,jj,jk) )145 & / ( e1e2u(ji,jj) * e3u_n(ji,jj,jk) ) 147 146 va(ji,jj,jk) = va(ji,jj,jk) - ( zfv_vw(ji,jj,jk) - zfv_vw(ji,jj,jk+1) ) & 148 & / ( e1e2v(ji,jj) * fse3v(ji,jj,jk) )147 & / ( e1e2v(ji,jj) * e3v_n(ji,jj,jk) ) 149 148 END DO 150 149 END DO
Note: See TracChangeset
for help on using the changeset viewer.