GCC Code Coverage Report
Directory: ./ Exec Total Coverage
File: phylmd/rrtm/suaerl.F90 Lines: 10 10 100.0 %
Date: 2023-06-30 12:56:34 Branches: 14 16 87.5 %

Line Branch Exec Source
1
1
SUBROUTINE SUAERL
2
3
!**** *SUAERL*   - INITIALIZE COMMON YOEAER
4
5
!     PURPOSE.
6
!     --------
7
!           INITIALIZE YOEAER, THE COMMON THAT CONTAINS THE
8
!           RADIATIVE CHARACTERISTICS OF THE AEROSOLS
9
10
!**   INTERFACE.
11
!     ----------
12
!              -----        -----
13
14
!        EXPLICIT ARGUMENTS :
15
!        --------------------
16
!        NONE
17
18
!        IMPLICIT ARGUMENTS :
19
!        --------------------
20
!        COMMON YOEAER
21
22
!     METHOD.
23
!     -------
24
!        SEE DOCUMENTATION
25
26
!     EXTERNALS.
27
!     ----------
28
29
!     REFERENCE.
30
!     ----------
31
!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE "IFS MODEL"
32
33
!     AUTHOR.
34
!     -------
35
!        JEAN-JACQUES MORCRETTE *ECMWF*
36
37
!     MODIFICATIONS.
38
!     --------------
39
!        ORIGINAL : 88-02-15
40
!        99-05-25   JJMorcrette   Revised Aerosol LW properties (OPAC)
41
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
42
43
!     ------------------------------------------------------------------
44
45
USE PARKIND1  ,ONLY : JPIM     ,JPRB
46
USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
47
48
USE YOESW    , ONLY : RAER
49
50
!      ----------------------------------------------------------------
51
52
!*       1.    SHORTWAVE COEFFICIENTS
53
!              ----------------------
54
55
!      Initialised in SUAERSn depending on NSW number of SW spectral
56
!      intervals
57
58
!      ----------------------------------------------------------------
59
60
!*       2.    LONGWAVE COEFFICIENTS
61
!              ---------------------
62
63
!=======================================================================
64
!-- The (old) five aerosol types were respectively:
65
66
!  1/ continental average (+desert)       2/ maritime
67
!  3/ urban                               4/ volcanic active
68
!  5/ stratospheric background
69
70
!      RAER  = RESHAPE((/
71
!     &          .038520, .037196, .040532, .054934, .038520
72
!     &        , .126130, .18313 , .10357 , .064106, .126130
73
!     &        , .012579, .013649, .018652, .025181, .012579
74
!     &        , .011890, .016142, .021105, .028908, .011890
75
!     &        , .013792, .026810, .052203, .066338, .013792 /)
76
!     & ,SHAPE=(/5,5/))
77
78
!=======================================================================
79
80
!-- The six aerosol types are respectively:
81
82
!  1/ continental average                 2/ maritime
83
!  3/ desert                              4/ urban
84
!  5/ volcanic active                     6/ stratospheric background
85
86
! The quantities given are:
87
! TAU : ratio of average optical thickness in interval to that at 0.55
88
!       micron
89
! PIZA: average single scattering albedo
90
! CGA : average asymmetry factor
91
92
! computed from Hess and Koepke (con, mar, des, urb)
93
!          from Bonnel et al.   (vol, str)
94
95
!-- data are entered for the 6 spectral intervals of the LW scheme (line)
96
!   and the different types of aerosols (column)
97
98
IMPLICIT NONE
99
REAL(KIND=JPRB) :: ZHOOK_HANDLE
100
1
IF (LHOOK) CALL DR_HOOK('SUAERL',0,ZHOOK_HANDLE)
101
RAER( :, 1)= (/&
102
7
 & .036271_JPRB, .030153_JPRB, .017343_JPRB, .015002_JPRB, .008806_JPRB, .006865_JPRB /)
103
RAER( :, 2)= (/&
104
7
 & .026561_JPRB, .032657_JPRB, .017977_JPRB, .014210_JPRB, .016775_JPRB, .022123_JPRB /)
105
RAER( :, 3)= (/&
106
7
 & .014897_JPRB, .016359_JPRB, .019789_JPRB, .030777_JPRB, .013341_JPRB, .014321_JPRB /)
107
RAER( :, 4)= (/&
108
7
 & .001863_JPRB, .002816_JPRB, .002355_JPRB, .002557_JPRB, .001774_JPRB, .001780_JPRB /)
109
RAER( :, 5)= (/&
110
7
 & .011890_JPRB, .016142_JPRB, .021105_JPRB, .028908_JPRB, .011890_JPRB, .011890_JPRB /)
111
RAER( :, 6)= (/&
112
7
 & .013792_JPRB, .026810_JPRB, .052203_JPRB, .066338_JPRB, .013792_JPRB, .013792_JPRB /)
113
114
!      ----------------------------------------------------------------
115
116
1
IF (LHOOK) CALL DR_HOOK('SUAERL',1,ZHOOK_HANDLE)
117
1
END SUBROUTINE SUAERL