Changeset 1119 for trunk/NEMO/TOP_SRC/TRP/trcadv_tvd.F90
- Timestamp:
- 2008-06-20T17:17:41+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/TRP/trcadv_tvd.F90
r941 r1119 14 14 !!---------------------------------------------------------------------- 15 15 USE oce_trc ! ocean dynamics and active tracers variables 16 USE tr c ! ocean passive tracers variables16 USE trp_trc ! ocean passive tracers variables 17 17 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 18 18 USE trcbbl ! advective passive tracers in the BBL … … 342 342 INTEGER :: ikm1 343 343 REAL(wp), DIMENSION (jpi,jpj,jpk) :: zbetup, zbetdo 344 REAL(wp) :: zpos, zneg, zbt, za, zb, zc, zbig, z rtrn, z2dtt344 REAL(wp) :: zpos, zneg, zbt, za, zb, zc, zbig, z2dtt 345 345 !!---------------------------------------------------------------------- 346 346 347 347 zbig = 1.e+40 348 zrtrn = 1.e-15349 348 zbetup(:,:,:) = 0.e0 ; zbetdo(:,:,:) = 0.e0 350 349 … … 410 409 ! up & down beta terms 411 410 zbt = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) / z2dtt 412 zbetup(ji,jj,jk) = ( zbetup(ji,jj,jk) - paft(ji,jj,jk) ) / (zpos+ zrtrn) * zbt413 zbetdo(ji,jj,jk) = ( paft(ji,jj,jk) - zbetdo(ji,jj,jk) ) / (zneg+ zrtrn) * zbt411 zbetup(ji,jj,jk) = ( zbetup(ji,jj,jk) - paft(ji,jj,jk) ) / (zpos+rtrn) * zbt 412 zbetdo(ji,jj,jk) = ( paft(ji,jj,jk) - zbetdo(ji,jj,jk) ) / (zneg+rtrn) * zbt 414 413 END DO 415 414 END DO
Note: See TracChangeset
for help on using the changeset viewer.