GCC Code Coverage Report


Directory: ./
File: rad/suaer15.f90
Date: 2022-01-11 19:19:34
Exec Total Coverage
Lines: 0 21 0.0%
Branches: 0 30 0.0%

Line Branch Exec Source
1 SUBROUTINE SUAER15
2
3 !**** *SUAER15* - INITIALIZE COMMON YOMAER15
4
5 ! PURPOSE.
6 ! --------
7 ! INITIALIZE YOMAER15, 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 YOMAER15
21
22 ! METHOD.
23 ! -------
24 ! SEE DOCUMENTATION
25
26 ! EXTERNALS.
27 ! ----------
28
29 ! REFERENCE.
30 ! ----------
31
32 ! 04-06: F. Bouyssel. Meteo-France
33 ! 05-09: A. Alias - PAESOD (black_carbon) is added to PAELAN (oraganic)
34 ! - the sulfate is put in a separate type
35 ! (see Hu Ron Ming work) - P.Marquet
36
37 !=======================================================================
38 !-- The (old) five aerosol types were respectively:
39
40 ! 1/ continental average (+desert) 2/ maritime
41 ! 3/ urban 4/ volcanic active
42 ! 5/ stratospheric background
43
44 !=======================================================================
45
46 !-- The (new) six aerosol types are respectively:
47
48 ! 1/ continental average 2/ maritime
49 ! 3/ desert 4/ volcanic active
50 ! 5/ stratospheric background 6/ sulfate
51
52 ! ------------------------------------------------------------------
53
54 USE PARKIND1 ,ONLY : JPIM ,JPRB
55 USE YOMHOOK ,ONLY : LHOOK, DR_HOOK
56
57 USE YOMRAD15 , ONLY : LNEWAER15
58 USE YOMAER15 , ONLY : RAER15 ,TAUA15 ,RPIZA15 ,RCGA15
59
60 IMPLICIT NONE
61
62 ! ------------------------------------------------------------------
63
64 REAL(KIND=JPRB) :: ZHOOK_HANDLE
65 IF (LHOOK) CALL DR_HOOK('SUAER15',0,ZHOOK_HANDLE)
66 IF (.NOT.LNEWAER15) THEN
67
68 RAER15=RESHAPE((/&
69 & .038520_JPRB, .037196_JPRB, .040532_JPRB, .054934_JPRB, .038520_JPRB ,&
70 & .12613_JPRB , .18313_JPRB , .10357_JPRB , .064106_JPRB, .126130_JPRB ,&
71 & .012579_JPRB, .013649_JPRB, .018652_JPRB, .025181_JPRB, .012579_JPRB ,&
72 & .011890_JPRB, .016142_JPRB, .021105_JPRB, .028908_JPRB, .011890_JPRB ,&
73 & .013792_JPRB, .026810_JPRB, .052203_JPRB, .066338_JPRB, .013792_JPRB ,&
74 & .012579_JPRB, .013649_JPRB, .018652_JPRB, .025181_JPRB, .012579_JPRB /)&
75 & ,SHAPE=(/5,6/))
76
77 TAUA15=RESHAPE((/&
78 & .730719_JPRB,.730719_JPRB,&
79 & .912819_JPRB,.912819_JPRB,&
80 & .725059_JPRB,.725059_JPRB,&
81 & .745405_JPRB,.745405_JPRB,&
82 & .682188_JPRB,.682188_JPRB,&
83 & 1.35059_JPRB,.725059_JPRB /)&
84 & ,SHAPE=(/2,6/))
85
86 RPIZA15=RESHAPE((/&
87 & .872212_JPRB,.872212_JPRB,&
88 & .982545_JPRB,.982545_JPRB,&
89 & .623143_JPRB,.623143_JPRB,&
90 & .944887_JPRB,.944887_JPRB,&
91 & .997975_JPRB,.997975_JPRB,&
92 & .999999_JPRB,.999999_JPRB /)&
93 & ,SHAPE=(/2,6/))
94
95 RCGA15=RESHAPE((/&
96 & .647596_JPRB,.647596_JPRB,&
97 & .739002_JPRB,.739002_JPRB,&
98 & .580845_JPRB,.580845_JPRB,&
99 & .662657_JPRB,.662657_JPRB,&
100 & .624246_JPRB,.624246_JPRB,&
101 & .680845_JPRB,.680845_JPRB /)&
102 & ,SHAPE=(/2,6/))
103
104 ELSE
105
106 RAER15( :, 1)= (/&
107 & .036271_JPRB, .030153_JPRB, .017343_JPRB, .015002_JPRB, .008806_JPRB /)
108 RAER15( :, 2)= (/&
109 & .026561_JPRB, .032657_JPRB, .017977_JPRB, .014210_JPRB, .016775_JPRB /)
110 RAER15( :, 3)= (/&
111 & .014897_JPRB, .016359_JPRB, .019789_JPRB, .030777_JPRB, .013341_JPRB /)
112 RAER15( :, 4)= (/&
113 & .011890_JPRB, .016142_JPRB, .021105_JPRB, .028908_JPRB, .011890_JPRB /)
114 RAER15( :, 5)= (/&
115 & .013792_JPRB, .026810_JPRB, .052203_JPRB, .066338_JPRB, .013792_JPRB /)
116 RAER15( :, 6)= (/&
117 & .012579_JPRB, .013649_JPRB, .018652_JPRB, .025181_JPRB, .012579_JPRB /)
118
119 TAUA15(1, :)= (/&
120 & 1.69446_JPRB, 1.11855_JPRB, 1.09212_JPRB, 1.03858_JPRB, 1.12044_JPRB, 1.35059_JPRB /)
121 TAUA15(2, :)= (/&
122 & 0.40174_JPRB, 0.89383_JPRB, 0.89546_JPRB, 0.51143_JPRB, 0.32646_JPRB, .725059_JPRB /)
123
124 RPIZA15(1, :)= (/&
125 & .9148907_JPRB, .9956173_JPRB, .7504584_JPRB, .9401905_JPRB, .9999999_JPRB, .999999_JPRB/)
126 RPIZA15(2, :)= (/&
127 & .8814597_JPRB, .9920407_JPRB, .9239428_JPRB, .9515548_JPRB, .9938563_JPRB, .999999_JPRB/)
128
129 RCGA15(1, :)= (/&
130 & 0.729019_JPRB, 0.803129_JPRB, 0.784592_JPRB, .7008249_JPRB, .7270548_JPRB, .680845_JPRB/)
131 RCGA15(2, :)= (/&
132 & 0.663224_JPRB, 0.793746_JPRB, 0.696315_JPRB, .6608509_JPRB, .6318786_JPRB, .680845_JPRB/)
133
134 ENDIF
135
136 ! ----------------------------------------------------------------
137
138 IF (LHOOK) CALL DR_HOOK('SUAER15',1,ZHOOK_HANDLE)
139 END SUBROUTINE SUAER15
140