source: trunk/src/ew.f90 @ 28

Last change on this file since 28 was 4, checked in by mlanger, 11 years ago
File size: 2.0 KB
Line 
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
22real 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
47end function ew
Note: See TracBrowser for help on using the repository browser.
hosted by ZAMG