trunk/phylmd/FCTTRE.f

module FCTTRE

  ! From phylmd/FCTTRE.inc, version 1.2, 2004/06/22 11:45:32

  ! This module includes the thermodynamical functions for the cycle
  ! 39 ECMWF physics package. Consistent with "SUPHEC_M" basic
  ! physical constants, assuming the partial pressure of water vapour
  ! is given by a first order Taylor expansion of "Qs(T)" with respect
  ! to temperature, using constants in "yoethf_m".

  ! Probably from Buck, 1981, Journal of Applied Meteorology, volume
  ! 20, number 12, page 1527.

  implicit none

contains

  elemental REAL function FOEEW(T, ICE)

    use yoethf_m, only: R3LES, R3IES, R4LES, R4IES
    use SUPHEC_M, only: rtt

    REAL, intent(in):: T
    logical, intent(in):: ICE ! else liquid

    !-----------------------

    FOEEW = exp(merge(R3IES / (T - R4IES), R3lES / (T - R4lES), ice) &
         * (T - RTT))

  end function FOEEW

  !******************************************

  pure REAL function FOEDE(T, ICE, P5ARG, QS, PCOARG)

    use yoethf_m, only: R4LES, R4IES

    REAL, intent(in):: T
    logical, intent(in):: ICE ! else liquid
    real, intent(in):: P5ARG, QS, PCOARG

    !-----------------------

    FOEDE = QS * PCOARG * P5ARG / (T - merge(R4IES, R4lES, ice))**2

  end function FOEDE

end module FCTTRE
1	guez	3	module FCTTRE
2
3	guez	78	! From phylmd/FCTTRE.inc, version 1.2, 2004/06/22 11:45:32
4	guez	3
5	guez	43	! This module includes the thermodynamical functions for the cycle
6	guez	78	! 39 ECMWF physics package. Consistent with "SUPHEC_M" basic
7	guez	43	! physical constants, assuming the partial pressure of water vapour
8			! is given by a first order Taylor expansion of "Qs(T)" with respect
9			! to temperature, using constants in "yoethf_m".
10	guez	3
11	guez	70	! Probably from Buck, 1981, Journal of Applied Meteorology, volume
12			! 20, number 12, page 1527.
13	guez	43
14	guez	3	implicit none
15
16			contains
17
18	guez	109	elemental REAL function FOEEW(T, ICE)
19	guez	3
20	guez	38	use yoethf_m, only: R3LES, R3IES, R4LES, R4IES
21			use SUPHEC_M, only: rtt
22	guez	3
23	guez	70	REAL, intent(in):: T
24	guez	109	logical, intent(in):: ICE ! else liquid
25	guez	3
26			!-----------------------
27
28	guez	109	FOEEW = exp(merge(R3IES / (T - R4IES), R3lES / (T - R4lES), ice) &
29	guez	103	* (T - RTT))
30	guez	3
31			end function FOEEW
32
33			!******************************************
34
35	guez	324	pure REAL function FOEDE(T, ICE, P5ARG, QS, PCOARG)
36	guez	3
37	guez	38	use yoethf_m, only: R4LES, R4IES
38	guez	3
39	guez	103	REAL, intent(in):: T
40	guez	109	logical, intent(in):: ICE ! else liquid
41	guez	37	real, intent(in):: P5ARG, QS, PCOARG
42	guez	3
43			!-----------------------
44
45	guez	109	FOEDE = QS * PCOARG * P5ARG / (T - merge(R4IES, R4lES, ice))**2
46	guez	3
47			end function FOEDE
48
49			end module FCTTRE