Philippe 23/02/2016: lfi2cdf: modif for JPHEXT/=1
[MNH-git_open_source-lfs.git] / LIBTOOLS / tools / diachro / src / POS / esat.f90
1 !-----------------------------------------------------------------
2 !--------------- special set of characters for SCCS information
3 !-----------------------------------------------------------------
4 !      @(#) Lib:/opt/local/MESONH/sources/post/s.esat.f90, Version:1.2, Date:98/10/02, Last modified:98/06/04
5 !-----------------------------------------------------------------
6 !     ######spl
7       FUNCTION ESAT(PT)
8 !     #################
9 !
10 !!****  *ESAT* - Computes the saturation water vapor pressure
11 !!            
12 !!
13 !!    PURPOSE
14 !!    -------
15 !       Computes the saturation water vapor pressure at a given temperature,
16 !      used in the emagram routine of TRACE
17 !
18 !!**  METHOD
19 !!    ------ 
20 !!      Analytical formula of Tetens (1930)
21 !!
22 !!    EXTERNAL
23 !!    --------
24 !!      NONE 
25 !!
26 !!    IMPLICIT ARGUMENTS
27 !!    ------------------
28 !!      NONE
29 !!
30 !!    REFERENCE
31 !!    ---------
32 !!
33 !!      Among many others, see for instance:
34 !!       Bluestein H. B., 1992, "Synoptic-Dynamic Meteorology in mid-latitudes"
35 !!       Volume 1, Priciples of Kinematics and Dynamics, Section 4.3, p. 195,
36 !!       Oxford University Press.
37 !!
38 !!    AUTHOR
39 !!    ------
40 !!      - Initial version Peridot TRACE Program, P.Bougeault *Meteo-France*,
41 !!      modified by R. Benoit (mc2, april 91) for the PYREX Oracle data base.
42 !!      - Present version J. Duron    * Laboratoire d'Aerologie *
43 !!
44 !!    MODIFICATIONS
45 !!    -------------
46 !!      Original       06/06/94
47 !!      Updated   PM   10/01/95
48 !-------------------------------------------------------------------------------
49 !
50 !*       0.    DECLARATIONS
51 !              ------------
52 !
53 IMPLICIT NONE
54 !
55 !*       0.1   Declaration of argument and result
56 !
57 REAL,INTENT(IN)                :: PT
58 REAL                           :: ESAT
59 !
60 !*       0.2   Declaration of local variables
61 !
62 !
63 REAL        :: ZABZ, ZTC
64 !
65 !-------------------------------------------------------------------------------
66 !
67 !*       1.    CALCULATION OF ESAT
68 !              -------------------
69 !
70 ! ESAT (MILLIBARS), PT (KELVIN)
71 !
72 ZABZ=273.16
73 ZTC = PT-ZABZ
74 ESAT = 6.1078*EXP((17.2693882*ZTC)/(ZTC+237.3))
75 !
76 !------------------------------------------------------------------------------
77 !
78 !*       2.     EXIT
79 !               ----
80 !
81 RETURN
82 END FUNCTION ESAT