Changeset 5845 for branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.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/TRA/traadv_qck.F90
r5836 r5845 39 39 40 40 !! * Substitutions 41 # include "domzgr_substitute.h90"42 41 # include "vectopt_loop_substitute.h90" 43 42 !!---------------------------------------------------------------------- … … 171 170 DO ji = fs_2, fs_jpim1 ! vector opt. 172 171 zdir = 0.5 + SIGN( 0.5, pun(ji,jj,jk) ) ! if pun > 0 : zdir = 1 otherwise zdir = 0 173 zdx = ( zdir * e1t(ji,jj) + ( 1. - zdir ) * e1t(ji+1,jj) ) * e2u(ji,jj) * fse3u(ji,jj,jk)172 zdx = ( zdir * e1t(ji,jj) + ( 1. - zdir ) * e1t(ji+1,jj) ) * e2u(ji,jj) * e3u_n(ji,jj,jk) 174 173 zwx(ji,jj,jk) = ABS( pun(ji,jj,jk) ) * zdt / zdx ! (0<zc_cfl<1 : Courant number on x-direction) 175 174 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 … … 217 216 DO jj = 2, jpjm1 218 217 DO ji = fs_2, fs_jpim1 ! vector opt. 219 zbtr = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) )218 zbtr = 1. / ( e1e2t(ji,jj) * e3t_n(ji,jj,jk) ) 220 219 ! horizontal advective trends 221 220 ztra = - zbtr * ( zwx(ji,jj,jk) - zwx(ji-1,jj,jk) ) … … 294 293 DO ji = fs_2, fs_jpim1 ! vector opt. 295 294 zdir = 0.5 + SIGN( 0.5, pvn(ji,jj,jk) ) ! if pun > 0 : zdir = 1 otherwise zdir = 0 296 zdx = ( zdir * e2t(ji,jj) + ( 1. - zdir ) * e2t(ji,jj+1) ) * e1v(ji,jj) * fse3v(ji,jj,jk)295 zdx = ( zdir * e2t(ji,jj) + ( 1. - zdir ) * e2t(ji,jj+1) ) * e1v(ji,jj) * e3v_n(ji,jj,jk) 297 296 zwy(ji,jj,jk) = ABS( pvn(ji,jj,jk) ) * zdt / zdx ! (0<zc_cfl<1 : Courant number on x-direction) 298 297 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 … … 341 340 DO jj = 2, jpjm1 342 341 DO ji = fs_2, fs_jpim1 ! vector opt. 343 zbtr = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) )342 zbtr = 1. / ( e1e2t(ji,jj) * e3t_n(ji,jj,jk) ) 344 343 ! horizontal advective trends 345 344 ztra = - zbtr * ( zwy(ji,jj,jk) - zwy(ji,jj-1,jk) ) … … 413 412 DO ji = fs_2, fs_jpim1 ! vector opt. 414 413 pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) - ( zwz(ji,jj,jk) - zwz(ji,jj,jk+1) ) & 415 & / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) )414 & / ( e1e2t(ji,jj) * e3t_n(ji,jj,jk) ) 416 415 END DO 417 416 END DO
Note: See TracChangeset
for help on using the changeset viewer.