dynamic_viscosity.f90 307 Bytes
Newer Older
Matthias Langer's avatar
 
Matthias Langer committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
! Function calculates dynamic viscosity of air (kg/m/s) as function of
! temperature (K) using Sutherland's formula

real function viscosity(t)

  implicit none

  real :: t
  real,parameter :: c=120.,t_0=291.15,eta_0=1.827e-5

  viscosity=eta_0*(t_0+c)/(t+c)*(t/t_0)**1.5

  return

end function viscosity