-
Notifications
You must be signed in to change notification settings - Fork 146
/
radiation_cloud_overlap.F90
116 lines (96 loc) · 4.05 KB
/
radiation_cloud_overlap.F90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
module module_radiation_cloud_overlap
use physparam, only : kind_phys
implicit none
public :: cmp_dcorr_lgth
interface cmp_dcorr_lgth
module procedure cmp_dcorr_lgth_hogan
module procedure cmp_dcorr_lgth_oreopoulos
end interface
contains
! ######################################################################################
! Hogan et al. (2010)
! "Effect of improving representation of horizontal and vertical cloud structure on the
! Earth's global radiation budget. Part I: Review and parametrization"
! https://rmets.onlinelibrary.wiley.com/doi/full/10.1002/qj.647
! ######################################################################################
subroutine cmp_dcorr_lgth_hogan(nCol, lat, con_pi, dcorr_lgth)
! Inputs
integer, intent(in) :: &
nCol ! Number of horizontal grid-points
real(kind_phys), intent(in) :: &
con_pi ! Physical constant: Pi
real(kind_phys), dimension(nCol), intent(in) :: &
lat ! Latitude
! Outputs
real(kind_phys), dimension(nCol),intent(out) :: &
dcorr_lgth ! Decorrelation length
! Local variables
integer :: iCol
! Parameters
real(kind_phys),parameter :: min_dcorr = 0.6 ! (see section 2.3)
do iCol =1,nCol
dcorr_lgth(iCol) = max(min_dcorr, 2.78-4.6*abs(lat(iCol)/con_pi)) ! (eq. 13)
enddo
end subroutine cmp_dcorr_lgth_hogan
! ######################################################################################
! Oreopoulos et al. (2012)
! "Radiative impacts of cloud heterogeneity and overlap in an
! atmospheric General Circulation Model"
! 10.5194/acp-12-9097-2012
! ######################################################################################
subroutine cmp_dcorr_lgth_oreopoulos(nCol, lat, juldat, yearlength, dcorr_lgth)
! Inputs
integer, intent(in) :: &
nCol, & ! Number of horizontal grid-points
yearlength ! Number of days in year
real(kind_phys), intent(in) :: &
juldat ! Julian date
real(kind_phys), dimension(nCol), intent(in) :: &
lat ! Latitude
! Outputs
real(kind_phys), dimension(nCol),intent(out) :: &
dcorr_lgth ! Decorrelation length (km)
! Parameters for the Gaussian fits per Eqs. (10) and (11) (See Table 1)
real(kind_phys), parameter :: & !
am1 = 1.4315_kind_phys, & !
am2 = 2.1219_kind_phys, & !
am4 = -25.584_kind_phys, & !
amr = 7.0_kind_phys !
! Local variables
integer :: iCol
real(kind_phys) :: am3
do iCol = 1, nCol
if (juldat .gt. 181._kind_phys) then
am3 = -4._kind_phys * amr * (juldat - 272._kind_phys) / yearlength ! (eq. 11a)
else
am3 = 4._kind_phys * amr * (juldat - 91._kind_phys) / yearlength ! (eq. 11b)
endif
dcorr_lgth(iCol) = am1 + am2 * exp( -(lat(iCol) - am3)**2 / am4**2) ! (eq. 10)
enddo
end subroutine cmp_dcorr_lgth_oreopoulos
! ######################################################################################
!
! ######################################################################################
subroutine get_alpha_exp(nCol, nLay, dzlay, dcorr_lgth, alpha)
! Inputs
integer, intent(in) :: &
nCol, & ! Number of horizontal grid points
nLay ! Number of vertical grid points
real(kind_phys), dimension(nCol), intent(in) :: &
dcorr_lgth ! Decorrelation length (km)
real(kind_phys), dimension(nCol,nLay), intent(in) :: &
dzlay !
! Outputs
real(kind_phys), dimension(nCol,nLay) :: &
alpha ! Cloud overlap parameter
! Local variables
integer :: iCol,iLay
do iCol = 1, nCol
alpha(iCol,1) = 0.0d0
do iLay = 2, nLay
alpha(iCol,iLay) = exp( -(dzlay(iCol,iLay)) / dcorr_lgth(iCol))
enddo
enddo
return
end subroutine get_alpha_exp
end module module_radiation_cloud_overlap