- Timestamp:
- 2014-04-06T17:28:25+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90
r4499 r4616 191 191 z0u = SIGN( 0.5, pun(ji,jj,jk) ) 192 192 zalpha = 0.5 - z0u 193 zu = z0u - 0.5 * pun(ji,jj,jk) * zdt / ( e1 u(ji,jj) *e2u(ji,jj) * fse3u(ji,jj,jk) )193 zu = z0u - 0.5 * pun(ji,jj,jk) * zdt / ( e1e2u(ji,jj) * fse3u(ji,jj,jk) ) 194 194 zzwx = ptb(ji+1,jj,jk,jn) + xind(ji,jj,jk) * (zu * zslpx(ji+1,jj,jk)) 195 195 zzwy = ptb(ji ,jj,jk,jn) + xind(ji,jj,jk) * (zu * zslpx(ji ,jj,jk)) … … 198 198 z0v = SIGN( 0.5, pvn(ji,jj,jk) ) 199 199 zalpha = 0.5 - z0v 200 zv = z0v - 0.5 * pvn(ji,jj,jk) * zdt / ( e1 v(ji,jj) *e2v(ji,jj) * fse3v(ji,jj,jk) )200 zv = z0v - 0.5 * pvn(ji,jj,jk) * zdt / ( e1e2v(ji,jj) * fse3v(ji,jj,jk) ) 201 201 zzwx = ptb(ji,jj+1,jk,jn) + xind(ji,jj,jk) * (zv * zslpy(ji,jj+1,jk)) 202 202 zzwy = ptb(ji,jj ,jk,jn) + xind(ji,jj,jk) * (zv * zslpy(ji,jj ,jk)) … … 212 212 DO jj = 2, jpjm1 213 213 DO ji = fs_2, fs_jpim1 ! vector opt. 214 zbtr = 1. / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) )214 zbtr = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 215 215 ! horizontal advective trends 216 216 ztra = - zbtr * ( zwx(ji,jj,jk) - zwx(ji-1,jj ,jk ) & … … 270 270 DO jj = 2, jpjm1 271 271 DO ji = fs_2, fs_jpim1 ! vector opt. 272 zbtr = 1. / ( e1 t(ji,jj) *e2t(ji,jj) * fse3w(ji,jj,jk+1) )272 zbtr = 1. / ( e1e2t(ji,jj) * fse3w(ji,jj,jk+1) ) 273 273 z0w = SIGN( 0.5, pwn(ji,jj,jk+1) ) 274 274 zalpha = 0.5 + z0w … … 285 285 DO jj = 2, jpjm1 286 286 DO ji = fs_2, fs_jpim1 ! vector opt. 287 zbtr = 1. / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) )287 zbtr = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 288 288 ! vertical advective trends 289 289 ztra = - zbtr * ( zwx(ji,jj,jk) - zwx(ji,jj,jk+1) )
Note: See TracChangeset
for help on using the changeset viewer.