calcsigmaW <- function(ustar,zL,bw){ phiW <- ifelse(zL<0, (1 - 3*zL)^(1/3), rep(1,length(zL)) ) return(bw*ustar*phiW) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.