libf/phylmd/FCTTRE.f90

module FCTTRE

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

  ! This COMDECK includes the thermodynamical functions for the cycle
  ! 39 ECMWF Physics package. Consistent with YOMCST basic physics
  ! 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.

  implicit none

  LOGICAL, PARAMETER:: thermcep= .TRUE.

contains

  REAL function FOEEW(T, DEL)

    use yoethf, only: R3LES, R3IES, R4LES, R4IES
    use YOMCST, only: rtt

    REAL, intent(in):: T, DEL

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

    FOEEW = EXP((R3LES * (1. - DEL) + R3IES * DEL) * (T - RTT) &
         / (T - (R4LES * (1. - DEL) + R4IES * DEL)))

  end function FOEEW

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

  REAL function FOEDE(T, DEL, P5ARG, QS, PCOARG)

    use yoethf, only: R4LES, R4IES

    REAL, intent(in):: T, DEL
    real, intent(in):: P5ARG, QS, PCOARG

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

    FOEDE = QS*PCOARG*P5ARG / (T-(R4LES*(1.-DEL)+R4IES*DEL))**2

  end function FOEDE

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

  REAL function qsats(t)

    REAL, intent(in):: T

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

    qsats = 100.0 * 0.622 &
         * 10.**(2.07023 - 0.00320991 * t - 2484.896 / t + 3.56654 * LOG10(t))

  end function qsats

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

  REAL function qsatl(t)

    REAL, intent(in):: T

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

    qsatl = 100.0 * 0.622 * 10.**(23.8319 - 2948.964 / t &
         - 5.028 * LOG10(t) - 29810.16 * EXP(- 0.0699382 * t) &
         + 25.21935 * EXP(- 2999.924 / t))

  end function qsatl

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

  REAL function dqsats(t, qs)

    use YOMCST, only: RLVTT, rcpd

    REAL, intent(in):: T, qs

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

    dqsats = RLVTT / RCPD * qs * (3.56654/t &
         +2484.896*LOG(10.)/t**2 - 0.00320991*LOG(10.))

  end function dqsats

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

  REAL function dqsatl(t, qs)

    use YOMCST, only: RLVTT, rcpd

    REAL, intent(in):: T, qs

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

    dqsatl = RLVTT / RCPD * qs * LOG(10.) &
         * (2948.964/t**2-5.028/LOG(10.)/t &
         +25.21935*2999.924/t**2*EXP(-2999.924/t) &
         +29810.16*0.0699382*EXP(-0.0699382*t))

  end function dqsatl

end module FCTTRE
1	module FCTTRE
2
3	! From phylmd/FCTTRE.inc, version 1.2 2004/06/22 11:45:32
4
5	! This COMDECK includes the thermodynamical functions for the cycle
6	! 39 ECMWF Physics package. Consistent with YOMCST basic physics
7	! constants, assuming the partial pressure of water vapour is given
8	! by a first order Taylor expansion of Qs(T) with respect to
9	! temperature, using constants in yoethf.
10
11	implicit none
12
13	LOGICAL, PARAMETER:: thermcep= .TRUE.
14
15	contains
16
17	REAL function FOEEW(T, DEL)
18
19	use yoethf, only: R3LES, R3IES, R4LES, R4IES
20	use YOMCST, only: rtt
21
22	REAL, intent(in):: T, DEL
23
24	!-----------------------
25
26	FOEEW = EXP((R3LES * (1. - DEL) + R3IES * DEL) * (T - RTT) &
27	/ (T - (R4LES * (1. - DEL) + R4IES * DEL)))
28
29	end function FOEEW
30
31	!******************************************
32
33	REAL function FOEDE(T, DEL, P5ARG, QS, PCOARG)
34
35	use yoethf, only: R4LES, R4IES
36
37	REAL, intent(in):: T, DEL
38	real, intent(in):: P5ARG, QS, PCOARG
39
40	!-----------------------
41
42	FOEDE = QSPCOARGP5ARG / (T-(R4LES(1.-DEL)+R4IESDEL))**2
43
44	end function FOEDE
45
46	!******************************************
47
48	REAL function qsats(t)
49
50	REAL, intent(in):: T
51
52	!-----------------------
53
54	qsats = 100.0 * 0.622 &
55	* 10.*(2.07023 - 0.00320991 t - 2484.896 / t + 3.56654 * LOG10(t))
56
57	end function qsats
58
59	!******************************************
60
61	REAL function qsatl(t)
62
63	REAL, intent(in):: T
64
65	!-----------------------
66
67	qsatl = 100.0 * 0.622 * 10.**(23.8319 - 2948.964 / t &
68	- 5.028 * LOG10(t) - 29810.16 * EXP(- 0.0699382 * t) &
69	+ 25.21935 * EXP(- 2999.924 / t))
70
71	end function qsatl
72
73	!******************************************
74
75	REAL function dqsats(t, qs)
76
77	use YOMCST, only: RLVTT, rcpd
78
79	REAL, intent(in):: T, qs
80
81	!-----------------------
82
83	dqsats = RLVTT / RCPD * qs * (3.56654/t &
84	+2484.896LOG(10.)/t2 - 0.00320991LOG(10.))
85
86	end function dqsats
87
88	!******************************************
89
90	REAL function dqsatl(t, qs)
91
92	use YOMCST, only: RLVTT, rcpd
93
94	REAL, intent(in):: T, qs
95
96	!-----------------------
97
98	dqsatl = RLVTT / RCPD * qs * LOG(10.) &
99	* (2948.964/t**2-5.028/LOG(10.)/t &
100	+25.219352999.924/t2EXP(-2999.924/t) &
101	+29810.160.0699382EXP(-0.0699382*t))
102
103	end function dqsatl
104
105	end module FCTTRE