viscosity {marelac} | R Documentation |
Calculates the shear viscosity of water, in centipoise.
Valid for 0<T<30 and 0<S<36.
Based on the code "CANDI" by B.P. Boudreau
viscosity(S=35, T=25, P=0)
S |
Salinity, - |
T |
Temperature, degrees C |
P |
Pressure, atm |
The details given in the original code by B. Boudreau are repeated here:
Uses the equation given by Kukulka et al. (1987).
shear visocisity in centipoise
Karline Soetaert <k.soetaert@nioo.knaw.nl>
Based on the FORTRAN implementation of the diagenetic model "CANDI" of B.P. Boudreau:
Boudreau BP,
A method-of-lines code for carbon and nutrient diagenesis in aquatic sediments.
COMPUTERS & GEOSCIENCES 22 (5): 479-496 JUN 1996
plot(0:30,viscosity(S=35,T=0:30,P=1),xlab="temperature",ylab="centipoise", main="shear viscosity of water",type="l") lines(0:30,viscosity(S=0,T=0:30,P=1),col="red") lines(0:30,viscosity(S=35,T=0:30,P=100),col="blue") legend("topright",col=c("black","red","blue"),lty=1, legend=c("S=35,P=1","S=0,P=1","S=35,P=100"))