10.4.1_peseiGFS_025bugfixes+enhancementsdevrelease-10release-10.4.1scaling-bug
|
Last change
on this file was
92fab65,
checked in by Ignacio Pisso <ip@…>, 4 years ago
|
|
add SPDX-License-Identifier to all .f90 files
|
-
Property mode set to
100644
|
|
File size:
423 bytes
|
| Line | |
|---|
| 1 | ! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt |
|---|
| 2 | ! SPDX-License-Identifier: GPL-3.0-or-later |
|---|
| 3 | |
|---|
| 4 | ! Function calculates dynamic viscosity of air (kg/m/s) as function of |
|---|
| 5 | ! temperature (K) using Sutherland's formula |
|---|
| 6 | |
|---|
| 7 | real function viscosity(t) |
|---|
| 8 | |
|---|
| 9 | implicit none |
|---|
| 10 | |
|---|
| 11 | real :: t |
|---|
| 12 | real,parameter :: c=120.,t_0=291.15,eta_0=1.827e-5 |
|---|
| 13 | |
|---|
| 14 | viscosity=eta_0*(t_0+c)/(t+c)*(t/t_0)**1.5 |
|---|
| 15 | |
|---|
| 16 | return |
|---|
| 17 | |
|---|
| 18 | end function viscosity |
|---|
Note: See
TracBrowser
for help on using the repository browser.
