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