swh:1:snp:bbc0b94b1c0463480b615c86f1ae2df3fcc700c0
Tip revision: e1145574dc02ebe9e0d96419755e4702aef85172 authored by Karline Soetaert on 16 September 2010, 00:00:00 UTC
version 1.3
version 1.3
Tip revision: e114557
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)
}