https://github.com/cran/ReacTran
Tip revision: 4945515d3d9eb03bb3f931bb41d2f124904f90bc authored by Karline Soetaert on 15 August 2017, 21:13:03 UTC
version 1.4.3.1
version 1.4.3.1
Tip revision: 4945515
fiadeiro.R
##==============================================================================
## The Fiadeiro and Veronis scheme: advective finite difference weights (AFDW)
##==============================================================================
fiadeiro <- function(v, D, dx.aux=NULL, grid=list(dx.aux=dx.aux)) {
Pe <- abs(v)*grid$dx.aux/D
sigma <- (1+(1/tanh(Pe)-1/Pe))/2
return(sigma)
}