- Timestamp:
- 02/13/12 08:23:51 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
vendor/nemo/current/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90
r4 r19 42 42 !!---------------------------------------------------------------------- 43 43 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 44 !! $Id: traadv_qck.F90 330 0 2012-02-08 09:54:42Z cbricaud $44 !! $Id: traadv_qck.F90 3301 2012-02-08 16:56:27Z cbricaud $ 45 45 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 46 46 !!---------------------------------------------------------------------- … … 177 177 DO ji = fs_2, fs_jpim1 ! vector opt. 178 178 zdir = 0.5 + SIGN( 0.5, pun(ji,jj,jk) ) ! if pun > 0 : zdir = 1 otherwise zdir = 0 179 zdx = zdir * e1t(ji,jj) + ( 1. - zdir ) * e1t(ji+1,jj)179 zdx = ( zdir * e1t(ji,jj) + ( 1. - zdir ) * e1t(ji+1,jj) ) * e2u(ji,jj) * fse3u(ji,jj,jk) 180 180 zwx(ji,jj,jk) = ABS( pun(ji,jj,jk) ) * zdt / zdx ! (0<zc_cfl<1 : Courant number on x-direction) 181 181 zfc(ji,jj,jk) = zdir * ptb(ji ,jj,jk,jn) + ( 1. - zdir ) * ptb(ji+1,jj,jk,jn) ! FC in the x-direction for T … … 302 302 DO ji = fs_2, fs_jpim1 ! vector opt. 303 303 zdir = 0.5 + SIGN( 0.5, pvn(ji,jj,jk) ) ! if pun > 0 : zdir = 1 otherwise zdir = 0 304 zdx = zdir * e2t(ji,jj) + ( 1. - zdir ) * e2t(ji,jj+1)304 zdx = ( zdir * e2t(ji,jj) + ( 1. - zdir ) * e2t(ji,jj+1) ) * e1v(ji,jj) * fse3v(ji,jj,jk) 305 305 zwy(ji,jj,jk) = ABS( pvn(ji,jj,jk) ) * zdt / zdx ! (0<zc_cfl<1 : Courant number on x-direction) 306 306 zfc(ji,jj,jk) = zdir * ptb(ji,jj ,jk,jn) + ( 1. - zdir ) * ptb(ji,jj+1,jk,jn) ! FC in the x-direction for T
Note: See TracChangeset
for help on using the changeset viewer.