FLEXPART CTBTO WO8
Main Page
Data Types List
Files
File List
File Members
All
Classes
Files
Functions
Variables
git
ctbto
flexpart_code
ew.f90
Go to the documentation of this file.
1
!**********************************************************************
2
! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010 *
3
! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, *
4
! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann *
5
! *
6
! This file is part of FLEXPART. *
7
! *
8
! FLEXPART is free software: you can redistribute it and/or modify *
9
! it under the terms of the GNU General Public License as published by*
10
! the Free Software Foundation, either version 3 of the License, or *
11
! (at your option) any later version. *
12
! *
13
! FLEXPART is distributed in the hope that it will be useful, *
14
! but WITHOUT ANY WARRANTY; without even the implied warranty of *
15
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
16
! GNU General Public License for more details. *
17
! *
18
! You should have received a copy of the GNU General Public License *
19
! along with FLEXPART. If not, see <http://www.gnu.org/licenses/>. *
20
!**********************************************************************
21
22
real
function
ew
(x)
23
24
!****************************************************************
25
!SAETTIGUNGSDAMPFDRUCK UEBER WASSER IN PA. X IN KELVIN.
26
!NACH DER GOFF-GRATCH-FORMEL.
27
!****************************************************************
28
29
implicit none
30
31
real
:: x, y, a, c, d
32
33
ew
=0.
34
if
(x.le.0.) stop
'sorry: t not in [k]'
35
y=373.16/x
36
a=-7.90298*(y-1.)
37
a=a+(5.02808*0.43429*alog(y))
38
c=(1.-(1./y))*11.344
39
c=-1.+(10.**c)
40
c=-1.3816*c/(10.**7)
41
d=(1.-y)*3.49149
42
d=-1.+(10.**d)
43
d=8.1328*d/(10.**3)
44
y=a+c+d
45
ew
=101324.6*(10.**y)
! Saettigungsdampfdruck in Pa
46
47
end function
ew
ew
real function ew(x)
Definition:
ew.f90:22
Generated on Sun Mar 27 2016 01:00:33 for FLEXPART CTBTO WO8 by
1.8.6