dynamic_viscosity.f90 423 Bytes
Newer Older
1 2
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
3

Matthias Langer's avatar
 
Matthias Langer committed
4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
! 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