Line |
Branch |
Exec |
Source |
1 |
|
|
module sens_heat_rain_m |
2 |
|
|
|
3 |
|
|
implicit none |
4 |
|
|
|
5 |
|
|
contains |
6 |
|
|
|
7 |
|
467812 |
real function sens_heat_rain(rain, t, q, rhoa, xlv, t_int, p) |
8 |
|
|
|
9 |
|
|
! Computes heat flux due to rainfall, in W m-2, positive |
10 |
|
|
! upward. |
11 |
|
|
|
12 |
|
|
! Do not declare this function as elemental because it may include |
13 |
|
|
! YOMCST.h, which contains an OpenMP directive, and that only |
14 |
|
|
! works in OpenMP 5.0. |
15 |
|
|
|
16 |
|
|
use const, only: cpa, cpw, rgas |
17 |
|
|
use esat_m, only: esat |
18 |
|
|
|
19 |
|
|
include "YOMCST.h" |
20 |
|
|
! for eps_w |
21 |
|
|
|
22 |
|
|
real, intent(in):: rain ! rain mass flux, in kg m-2 s-1 |
23 |
|
|
real, intent(in):: t ! air temperature at 10 m, in K |
24 |
|
|
real, intent(in):: q ! specific humidity at 10 m |
25 |
|
|
real, intent(in):: rhoa ! density of moist air (kg / m3) |
26 |
|
|
real, intent(in):: xlv ! latent heat of evaporation (J / kg) |
27 |
|
|
real, intent(in):: t_int ! interface temperature, in K |
28 |
|
|
real, intent(in):: p ! surface pressure, in Pa |
29 |
|
|
|
30 |
|
|
! Local: |
31 |
|
|
|
32 |
|
|
real es ! saturation pressure of wator vapor, in Pa |
33 |
|
|
real alfac ! wet bulb factor |
34 |
|
|
real dwat ! water vapour diffusivity |
35 |
|
|
real dtmp ! heat diffusivity |
36 |
|
|
real q_int ! specific (saturation) humidity at ocean interface |
37 |
|
|
real t_celsius ! air temperature at 10 m, in Celsius degrees |
38 |
|
|
|
39 |
|
|
real wetc |
40 |
|
|
! derivative of saturated mass fraction of water vapor with |
41 |
|
|
! respect to temperature, at constant total pressure |
42 |
|
|
|
43 |
|
|
!--------------------------------------------------------------------- |
44 |
|
|
|
45 |
|
467812 |
es = esat(t_int, p) * 0.98 ! reduced for salinity, Kraus 1972 page 46 |
46 |
|
467812 |
q_int = eps_w * (es / (p - (1. - eps_w) * es)) |
47 |
|
467812 |
wetc = eps_w * xlv * q_int / (rgas * t_int**2) |
48 |
|
467812 |
dwat = 2.11e-5 * (t / 273.15)**1.94 |
49 |
|
467812 |
t_celsius = t - 273.15 |
50 |
|
|
dtmp = (1. + 3.309e-3 * t_celsius - 1.44e-6 * t_celsius**2) * 0.02411 & |
51 |
|
467812 |
/ (rhoa * cpa) |
52 |
|
|
|
53 |
|
|
! Gosnell 1995 k0991, equation (11): |
54 |
|
467812 |
alfac = 1. / (1. + (wetc * xlv * dwat) / (cpa * dtmp)) |
55 |
|
|
|
56 |
|
|
! Gosnell 1995 k0991, equation (12): |
57 |
|
467812 |
sens_heat_rain = rain * alfac * cpw * (t_int - t + (q_int - q) * xlv / cpa) |
58 |
|
|
|
59 |
|
467812 |
end function sens_heat_rain |
60 |
|
|
|
61 |
|
|
end module sens_heat_rain_m |
62 |
|
|
|