- Timestamp:
- 2015-12-16T10:25:22+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_merge_2015/NEMOGCM/NEMO/OPA_SRC/DYN/divhor.F90
r5836 r6060 36 36 37 37 !! * Substitutions 38 # include "domzgr_substitute.h90"39 38 # include "vectopt_loop_substitute.h90" 40 39 !!---------------------------------------------------------------------- … … 74 73 DO jj = 2, jpjm1 75 74 DO ji = fs_2, fs_jpim1 ! vector opt. 76 hdivn(ji,jj,jk) = ( e2u(ji ,jj) * fse3u_n(ji ,jj,jk) * un(ji ,jj,jk) &77 & - e2u(ji-1,jj) * fse3u_n(ji-1,jj,jk) * un(ji-1,jj,jk) &78 & + e1v(ji,jj ) * fse3v_n(ji,jj ,jk) * vn(ji,jj ,jk) &79 & - e1v(ji,jj-1) * fse3v_n(ji,jj-1,jk) * vn(ji,jj-1,jk) ) &80 & / ( e1e2t(ji,jj) * fse3t_n(ji,jj,jk) )75 hdivn(ji,jj,jk) = ( e2u(ji ,jj) * e3u_n(ji ,jj,jk) * un(ji ,jj,jk) & 76 & - e2u(ji-1,jj) * e3u_n(ji-1,jj,jk) * un(ji-1,jj,jk) & 77 & + e1v(ji,jj ) * e3v_n(ji,jj ,jk) * vn(ji,jj ,jk) & 78 & - e1v(ji,jj-1) * e3v_n(ji,jj-1,jk) * vn(ji,jj-1,jk) ) & 79 & / ( e1e2t(ji,jj) * e3t_n(ji,jj,jk) ) 81 80 END DO 82 81 END DO
Note: See TracChangeset
for help on using the changeset viewer.