- 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/TOP_SRC/TRP/trcadv.F90
r5836 r6060 62 62 63 63 !! * Substitutions 64 # include "domzgr_substitute.h90"65 64 # include "vectopt_loop_substitute.h90" 66 65 !!---------------------------------------------------------------------- … … 109 108 ! !== effective transport ==! 110 109 DO jk = 1, jpkm1 111 zun(:,:,jk) = e2u (:,:) * fse3u(:,:,jk) * un(:,:,jk) ! eulerian transport112 zvn(:,:,jk) = e1v (:,:) * fse3v(:,:,jk) * vn(:,:,jk)110 zun(:,:,jk) = e2u (:,:) * e3u_n(:,:,jk) * un(:,:,jk) ! eulerian transport 111 zvn(:,:,jk) = e1v (:,:) * e3v_n(:,:,jk) * vn(:,:,jk) 113 112 zwn(:,:,jk) = e1e2t(:,:) * wn(:,:,jk) 114 113 END DO … … 231 230 CALL ctl_stop( 'trc_adv_init: force 2nd order FCT scheme, 4th order does not exist with sub-timestepping' ) 232 231 ENDIF 233 IF( lk_vvl) THEN232 IF( .NOT.ln_linssh ) THEN 234 233 CALL ctl_stop( 'trc_adv_init: vertical sub-timestepping not allow in non-linear free surface' ) 235 234 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.