source: flexpart.git/src/dynamic_viscosity.f90 @ 3481cc1

10.4.1_peseiGFS_025bugfixes+enhancementsdevrelease-10release-10.4.1scaling-bug
Last change on this file since 3481cc1 was 3481cc1, checked in by Ignacio Pisso <ip@…>, 4 years ago

move license from headers to a different file

  • Property mode set to 100644
File size: 307 bytes
Line 
1! Function calculates dynamic viscosity of air (kg/m/s) as function of
2! temperature (K) using Sutherland's formula
3
4real function viscosity(t)
5
6  implicit none
7
8  real :: t
9  real,parameter :: c=120.,t_0=291.15,eta_0=1.827e-5
10
11  viscosity=eta_0*(t_0+c)/(t+c)*(t/t_0)**1.5
12
13  return
14
15end function viscosity
Note: See TracBrowser for help on using the repository browser.
hosted by ZAMG