MAIC-2  Revision 19
 All Classes Files Functions Variables
p_sat.f90
Go to the documentation of this file.
1 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2 !
3 ! Function : p _ s a t
4 !
5 !> @file
6 !!
7 !! Computation of the water-vapour saturation pressure.
8 !!
9 !! @section Copyright
10 !!
11 !! Copyright 2010-2013 Ralf Greve, Bjoern Grieger, Oliver J. Stenzel
12 !!
13 !! @section License
14 !!
15 !! This file is part of MAIC-2.
16 !!
17 !! MAIC-2 is free software: you can redistribute it and/or modify
18 !! it under the terms of the GNU General Public License as published by
19 !! the Free Software Foundation, either version 3 of the License, or
20 !! (at your option) any later version.
21 !!
22 !! MAIC-2 is distributed in the hope that it will be useful,
23 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
24 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
25 !! GNU General Public License for more details.
26 !!
27 !! You should have received a copy of the GNU General Public License
28 !! along with MAIC-2. If not, see <http://www.gnu.org/licenses/>.
29 !<
30 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
31 
32 !-------------------------------------------------------------------------------
33 !> Computation of the water-vapour saturation pressure.
34 !<------------------------------------------------------------------------------
35 elemental function p_sat(temp_surf) result(psat)
36 
37 use maic2_types
38 
39 implicit none
40 
41 real(dp), intent(in) :: temp_surf
42 
43 real(dp) :: psat
44 
45 ! psat = exp( -5504.4088_dp * temp_surf**(-1) &
46 ! - 3.5704628_dp &
47 ! - 1.7337458e-2_dp * temp_surf &
48 ! + 6.5204209e-6_dp * temp_surf**2 &
49 ! + 6.1295027_dp * log(temp_surf) )
50 
51  psat = 610.66_dp * exp( 21.875_dp*(temp_surf-273.16_dp)/ &
52  (temp_surf- 7.65_dp) )
53 
54 end function p_sat
55 !