Directory: | ./ |
---|---|
File: | rad/suaerh.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 13 | 13 | 100.0% |
Branches: | 2 | 4 | 50.0% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | 1 | SUBROUTINE SUAERH | |
2 | |||
3 | !**** *SUAERH* - SETS UP THE AEROSOL HORIZONTAL DISTRIBUTION | ||
4 | |||
5 | ! PURPOSE. | ||
6 | ! -------- | ||
7 | |||
8 | ! THIS ROUTINE FETCHES VALUES OF A T10 SPECTRAL DISTRIBUTION | ||
9 | ! FOR FOUR AEROSOL TYPES OF DIFFERENT ORIGINS (SEA,LAND,URBAN AREAS | ||
10 | ! AND DESERTS). THE VALUES TO BE OBTAINED ARE BETWEEN ZERO AND ONE. | ||
11 | |||
12 | !** INTERFACE. | ||
13 | ! ---------- | ||
14 | |||
15 | ! THERE ARE EIGHT DUMMY ARGUMENTS: *PAESC*, *PAESS*, *PAELC*, | ||
16 | ! *PAELS*, *PAEUC*, *PAEUS*, *PAEDC* AND *PAEDS*ARE ARRAYS FOR THE | ||
17 | ! T10 DISTRIBUTIONS (*S FOR SEA, *L FOR LAND, *U FOR URBAN AND *D | ||
18 | ! FOR DESERT, *C FOR COSINE AND *S FOR SINE). | ||
19 | |||
20 | ! METHOD. | ||
21 | ! ------- | ||
22 | |||
23 | ! NONE. | ||
24 | |||
25 | ! EXTERNALS. | ||
26 | ! ---------- | ||
27 | |||
28 | ! NONE. | ||
29 | |||
30 | ! REFERENCE. | ||
31 | ! ---------- | ||
32 | |||
33 | ! NONE. | ||
34 | |||
35 | ! AUTHOR | ||
36 | ! ------ | ||
37 | ! J.-F. GELEYN E.C.M.W.F. 29/09/82. | ||
38 | |||
39 | ! MODIFICATIONS | ||
40 | ! ------------- | ||
41 | ! J.-J. MORCRETTE E.C.M.W.F. 91/07/14 ADAPTATION TO I.F.S. | ||
42 | ! M.Hamrud 01-Oct-2003 CY28 Cleaning | ||
43 | |||
44 | ! ------------------------------------------------------------------ | ||
45 | |||
46 | USE PARKIND1 ,ONLY : JPIM ,JPRB | ||
47 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK | ||
48 | |||
49 | USE YOEAERD , ONLY : RCAEROS ,RAESC ,RAESS ,RAELC ,& | ||
50 | & RAELS ,RAEUC ,RAEUS ,RAEDC ,RAEDS | ||
51 | |||
52 | IMPLICIT NONE | ||
53 | |||
54 | ! ------------------------------------------------------------------ | ||
55 | |||
56 | !* DATA STATEMENTS. | ||
57 | ! ---- ----------- | ||
58 | |||
59 | ! *RAE S/L/U/D C/S* (SEE ABOVE). | ||
60 | |||
61 | REAL(KIND=JPRB) :: ZHOOK_HANDLE | ||
62 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | IF (LHOOK) CALL DR_HOOK('SUAERH',0,ZHOOK_HANDLE) |
63 | RAESC = (/& | ||
64 | & +.6688E+00_JPRB,-.1172E+00_JPRB,-.1013E+00_JPRB,+.1636E-01_JPRB,-.3699E-01_JPRB,+.1775E-01_JPRB,& | ||
65 | & -.9635E-02_JPRB,+.1290E-02_JPRB,+.4681E-04_JPRB,-.9106E-04_JPRB,+.9355E-04_JPRB,& | ||
66 | & -.7076E-01_JPRB,-.1782E-01_JPRB,+.1856E-01_JPRB,+.1372E-01_JPRB,+.8210E-04_JPRB,& | ||
67 | & +.2149E-02_JPRB,+.4856E-03_JPRB,+.2231E-03_JPRB,+.1824E-03_JPRB,+.1960E-05_JPRB,& | ||
68 | & +.2057E-01_JPRB,+.2703E-01_JPRB,+.2424E-01_JPRB,+.9716E-02_JPRB,& | ||
69 | & +.1312E-02_JPRB,-.8846E-03_JPRB,-.3347E-03_JPRB,+.6231E-04_JPRB,+.6397E-04_JPRB,& | ||
70 | & -.3341E-02_JPRB,-.1295E-01_JPRB,-.4598E-02_JPRB,& | ||
71 | & +.3242E-03_JPRB,+.8122E-03_JPRB,-.2975E-03_JPRB,-.7757E-04_JPRB,+.7793E-04_JPRB,& | ||
72 | & +.4455E-02_JPRB,-.1584E-01_JPRB,& | ||
73 | & -.2551E-02_JPRB,+.1174E-02_JPRB,+.1335E-04_JPRB,+.5112E-04_JPRB,+.5605E-04_JPRB,& | ||
74 | & +.7412E-04_JPRB,& | ||
75 | & +.1857E-02_JPRB,-.1917E-03_JPRB,+.4460E-03_JPRB,+.1767E-04_JPRB,-.5281E-04_JPRB,& | ||
76 | & -.5043E-03_JPRB,+.2467E-03_JPRB,-.2497E-03_JPRB,-.2377E-04_JPRB,-.3954E-04_JPRB,& | ||
77 | & +.2666E-03_JPRB,-.8186E-03_JPRB,-.1441E-03_JPRB,-.1904E-04_JPRB,& | ||
78 | & +.3337E-03_JPRB,-.1696E-03_JPRB,-.2503E-04_JPRB,& | ||
79 | & +.1239E-03_JPRB,-.9983E-04_JPRB,& | ||
80 | & -.5283E-04_JPRB & | ||
81 | 1 | & /) | |
82 | RAESS = (/& | ||
83 | & -.3374E-01_JPRB,-.3247E-01_JPRB,-.1012E-01_JPRB,+.6002E-02_JPRB,+.5190E-02_JPRB,& | ||
84 | & +.7784E-03_JPRB,-.1090E-02_JPRB,+.3294E-03_JPRB,+.1719E-03_JPRB,-.5866E-05_JPRB,& | ||
85 | & -.4124E-03_JPRB,-.3742E-01_JPRB,-.5054E-02_JPRB,+.3430E-02_JPRB,& | ||
86 | & +.5513E-03_JPRB,-.6235E-03_JPRB,+.2892E-03_JPRB,-.9730E-04_JPRB,+.7078E-04_JPRB,& | ||
87 | & -.3300E-01_JPRB,+.5104E-03_JPRB,-.2156E-02_JPRB,& | ||
88 | & -.3194E-02_JPRB,-.5079E-03_JPRB,-.5517E-03_JPRB,+.4632E-04_JPRB,+.5369E-04_JPRB,& | ||
89 | & -.2731E-01_JPRB,+.5126E-02_JPRB,& | ||
90 | & +.2241E-02_JPRB,-.5789E-03_JPRB,-.3048E-03_JPRB,-.1774E-03_JPRB,+.1946E-05_JPRB,& | ||
91 | & -.8247E-02_JPRB,& | ||
92 | & +.2338E-02_JPRB,+.1021E-02_JPRB,+.1575E-04_JPRB,+.2612E-05_JPRB,+.1995E-04_JPRB,& | ||
93 | & -.1319E-02_JPRB,+.1384E-02_JPRB,-.4159E-03_JPRB,-.2337E-03_JPRB,+.5764E-04_JPRB,& | ||
94 | & +.1495E-02_JPRB,-.3727E-03_JPRB,+.6075E-04_JPRB,-.4642E-04_JPRB,& | ||
95 | & +.5368E-03_JPRB,-.7619E-04_JPRB,+.3774E-04_JPRB,& | ||
96 | & +.1206E-03_JPRB,-.4104E-06_JPRB,& | ||
97 | & +.2158E-04_JPRB & | ||
98 | 1 | & /) | |
99 | RAELC = (/& | ||
100 | & +.1542E+00_JPRB,+.8245E-01_JPRB,-.1879E-03_JPRB,+.4864E-02_JPRB,-.5527E-02_JPRB,-.7966E-02_JPRB,& | ||
101 | & -.2683E-02_JPRB,-.2011E-02_JPRB,-.8889E-03_JPRB,-.1058E-03_JPRB,-.1614E-04_JPRB,& | ||
102 | & +.4206E-01_JPRB,+.1912E-01_JPRB,-.9476E-02_JPRB,-.6780E-02_JPRB,+.1767E-03_JPRB,& | ||
103 | & -.5422E-03_JPRB,-.7753E-03_JPRB,-.2106E-03_JPRB,-.9870E-04_JPRB,-.1721E-04_JPRB,& | ||
104 | & -.9536E-02_JPRB,-.9580E-02_JPRB,-.1050E-01_JPRB,-.5747E-02_JPRB,& | ||
105 | & -.1282E-02_JPRB,+.2248E-03_JPRB,+.1694E-03_JPRB,-.4782E-04_JPRB,-.2441E-04_JPRB,& | ||
106 | & +.5781E-03_JPRB,+.6212E-02_JPRB,+.1921E-02_JPRB,& | ||
107 | & -.1102E-02_JPRB,-.8145E-03_JPRB,+.2497E-03_JPRB,+.1539E-03_JPRB,-.2538E-04_JPRB,& | ||
108 | & -.3993E-02_JPRB,+.9777E-02_JPRB,& | ||
109 | & +.4837E-03_JPRB,-.1304E-02_JPRB,+.2417E-04_JPRB,-.1370E-04_JPRB,-.3731E-05_JPRB,& | ||
110 | & +.1922E-02_JPRB,& | ||
111 | & -.5167E-03_JPRB,+.4295E-03_JPRB,-.1888E-03_JPRB,+.2427E-04_JPRB,+.4012E-04_JPRB,& | ||
112 | & +.1529E-02_JPRB,-.2120E-03_JPRB,+.8166E-04_JPRB,+.2579E-04_JPRB,+.3488E-04_JPRB,& | ||
113 | & +.2140E-03_JPRB,+.2274E-03_JPRB,-.3447E-05_JPRB,-.1075E-04_JPRB,& | ||
114 | & -.1018E-03_JPRB,+.2864E-04_JPRB,+.3442E-04_JPRB,& | ||
115 | & -.1002E-03_JPRB,+.7117E-04_JPRB,& | ||
116 | & +.2045E-04_JPRB & | ||
117 | 1 | & /) | |
118 | RAELS = (/& | ||
119 | & +.1637E-01_JPRB,+.1935E-01_JPRB,+.1080E-01_JPRB,+.2784E-02_JPRB,+.1606E-03_JPRB,& | ||
120 | & +.1860E-02_JPRB,+.1263E-02_JPRB,-.2707E-03_JPRB,-.2290E-03_JPRB,-.9761E-05_JPRB,& | ||
121 | & -.7317E-02_JPRB,+.2465E-01_JPRB,+.6799E-02_JPRB,-.1913E-02_JPRB,& | ||
122 | & +.1382E-02_JPRB,+.6691E-03_JPRB,+.1414E-03_JPRB,+.3527E-04_JPRB,-.5210E-04_JPRB,& | ||
123 | & +.1873E-01_JPRB,+.2977E-02_JPRB,+.4650E-02_JPRB,& | ||
124 | & +.2509E-02_JPRB,+.3680E-03_JPRB,+.1481E-03_JPRB,-.6594E-04_JPRB,-.5634E-04_JPRB,& | ||
125 | & +.1592E-01_JPRB,-.1875E-02_JPRB,& | ||
126 | & -.1093E-02_JPRB,+.3022E-03_JPRB,+.2625E-03_JPRB,+.3252E-04_JPRB,-.3803E-04_JPRB,& | ||
127 | & +.4218E-02_JPRB,& | ||
128 | & -.1843E-02_JPRB,-.1351E-02_JPRB,-.2952E-03_JPRB,-.8171E-05_JPRB,-.1473E-04_JPRB,& | ||
129 | & +.9076E-03_JPRB,-.1057E-02_JPRB,+.2676E-03_JPRB,+.1307E-03_JPRB,-.3628E-04_JPRB,& | ||
130 | & -.9158E-03_JPRB,+.4335E-03_JPRB,+.2927E-04_JPRB,+.6602E-04_JPRB,& | ||
131 | & -.3570E-03_JPRB,+.5760E-04_JPRB,-.3465E-04_JPRB,& | ||
132 | & -.8535E-04_JPRB,-.2011E-04_JPRB,& | ||
133 | & +.6612E-06_JPRB & | ||
134 | 1 | & /) | |
135 | RAEUC = (/& | ||
136 | & +.8005E-01_JPRB,+.7095E-01_JPRB,+.2014E-01_JPRB,-.1412E-01_JPRB,-.2425E-01_JPRB,-.1332E-01_JPRB,& | ||
137 | & -.2904E-02_JPRB,+.5068E-03_JPRB,+.9369E-03_JPRB,+.4114E-03_JPRB,+.7549E-04_JPRB,& | ||
138 | & +.1922E-01_JPRB,+.2534E-01_JPRB,+.2088E-01_JPRB,+.1064E-01_JPRB,+.1063E-02_JPRB,& | ||
139 | & -.2526E-02_JPRB,-.2091E-02_JPRB,-.9660E-03_JPRB,-.2030E-03_JPRB,+.3865E-04_JPRB,& | ||
140 | & -.9900E-02_JPRB,-.5964E-02_JPRB,+.2223E-02_JPRB,+.4941E-02_JPRB,& | ||
141 | & +.3277E-02_JPRB,+.1038E-02_JPRB,-.1480E-03_JPRB,-.2844E-03_JPRB,-.1208E-03_JPRB,& | ||
142 | & +.3999E-02_JPRB,+.6282E-02_JPRB,+.2813E-02_JPRB,& | ||
143 | & +.1475E-02_JPRB,+.4571E-03_JPRB,-.1349E-03_JPRB,-.9011E-04_JPRB,-.1936E-04_JPRB,& | ||
144 | & +.1994E-02_JPRB,+.3540E-02_JPRB,& | ||
145 | & +.8837E-03_JPRB,+.1992E-03_JPRB,+.3092E-04_JPRB,-.7979E-04_JPRB,-.2664E-04_JPRB,& | ||
146 | & -.5006E-04_JPRB,& | ||
147 | & +.6447E-03_JPRB,+.5550E-03_JPRB,+.1197E-03_JPRB,+.6657E-04_JPRB,+.1488E-04_JPRB,& | ||
148 | & -.9141E-04_JPRB,-.2896E-03_JPRB,-.1561E-03_JPRB,-.6524E-04_JPRB,-.1559E-04_JPRB,& | ||
149 | & -.1082E-03_JPRB,-.4126E-03_JPRB,-.1732E-03_JPRB,-.8286E-04_JPRB,& | ||
150 | & -.1993E-04_JPRB,+.3850E-04_JPRB,+.2870E-04_JPRB,& | ||
151 | & +.4493E-04_JPRB,+.4721E-04_JPRB,& | ||
152 | & +.1338E-04_JPRB & | ||
153 | 1 | & /) | |
154 | RAEUS = (/& | ||
155 | & +.6646E-02_JPRB,+.8373E-02_JPRB,+.5463E-02_JPRB,+.4554E-02_JPRB,+.3301E-02_JPRB,& | ||
156 | & +.5725E-03_JPRB,-.7482E-03_JPRB,-.6222E-03_JPRB,-.2603E-03_JPRB,-.5127E-04_JPRB,& | ||
157 | & -.3849E-04_JPRB,+.9741E-02_JPRB,+.8190E-02_JPRB,+.5712E-02_JPRB,& | ||
158 | & +.3039E-02_JPRB,+.5290E-03_JPRB,-.2044E-03_JPRB,-.2309E-03_JPRB,-.1160E-03_JPRB,& | ||
159 | & +.9160E-02_JPRB,+.1286E-01_JPRB,+.1170E-01_JPRB,& | ||
160 | & +.5491E-02_JPRB,+.1393E-02_JPRB,-.6288E-04_JPRB,-.2715E-03_JPRB,-.1047E-03_JPRB,& | ||
161 | & +.4873E-02_JPRB,+.3545E-02_JPRB,& | ||
162 | & +.3069E-02_JPRB,+.1819E-02_JPRB,+.6947E-03_JPRB,+.1416E-03_JPRB,-.1538E-04_JPRB,& | ||
163 | & -.4351E-03_JPRB,& | ||
164 | & -.1907E-02_JPRB,-.5774E-03_JPRB,-.2247E-03_JPRB,+.5345E-04_JPRB,+.9052E-04_JPRB,& | ||
165 | & -.3972E-04_JPRB,-.9665E-04_JPRB,+.7912E-04_JPRB,-.1094E-04_JPRB,-.6776E-05_JPRB,& | ||
166 | & +.2724E-03_JPRB,+.1973E-03_JPRB,+.6837E-04_JPRB,+.4313E-04_JPRB,& | ||
167 | & -.7174E-05_JPRB,+.8527E-05_JPRB,-.2160E-05_JPRB,& | ||
168 | & -.7852E-04_JPRB,+.3453E-06_JPRB,& | ||
169 | & -.2402E-05_JPRB & | ||
170 | 1 | & /) | |
171 | RAEDC = (/& | ||
172 | & +.2840E-01_JPRB,+.1775E-01_JPRB,-.1069E-01_JPRB,-.1553E-01_JPRB,-.3299E-02_JPRB,+.3583E-02_JPRB,& | ||
173 | & +.2274E-02_JPRB,+.5767E-04_JPRB,-.3678E-03_JPRB,-.1050E-03_JPRB,+.2133E-04_JPRB,& | ||
174 | & +.2326E-01_JPRB,+.1566E-01_JPRB,-.3130E-02_JPRB,-.8253E-02_JPRB,-.2615E-02_JPRB,& | ||
175 | & +.1247E-02_JPRB,+.1059E-02_JPRB,+.1196E-03_JPRB,-.1303E-03_JPRB,-.5094E-04_JPRB,& | ||
176 | & +.1185E-01_JPRB,+.7238E-02_JPRB,-.1562E-02_JPRB,-.3665E-02_JPRB,& | ||
177 | & -.1182E-02_JPRB,+.4678E-03_JPRB,+.4448E-03_JPRB,+.8307E-04_JPRB,-.3468E-04_JPRB,& | ||
178 | & +.5273E-02_JPRB,+.3037E-02_JPRB,-.4014E-03_JPRB,& | ||
179 | & -.1202E-02_JPRB,-.4647E-03_JPRB,+.5148E-04_JPRB,+.1014E-03_JPRB,+.2996E-04_JPRB,& | ||
180 | & +.2505E-02_JPRB,+.1495E-02_JPRB,& | ||
181 | & +.2438E-03_JPRB,-.1223E-03_JPRB,-.7669E-04_JPRB,-.1638E-04_JPRB,+.1869E-05_JPRB,& | ||
182 | & +.1094E-02_JPRB,& | ||
183 | & +.6131E-03_JPRB,+.1508E-03_JPRB,+.1765E-04_JPRB,+.1360E-05_JPRB,-.7998E-06_JPRB,& | ||
184 | & +.4475E-03_JPRB,+.2737E-03_JPRB,+.6430E-04_JPRB,-.6759E-05_JPRB,-.6761E-05_JPRB,& | ||
185 | & +.1992E-03_JPRB,+.1531E-03_JPRB,+.4828E-04_JPRB,+.5103E-06_JPRB,& | ||
186 | & +.7454E-04_JPRB,+.5917E-04_JPRB,+.2152E-04_JPRB,& | ||
187 | & +.9300E-05_JPRB,+.9790E-05_JPRB,& | ||
188 | & -.8853E-05_JPRB & | ||
189 | 1 | & /) | |
190 | RAEDS = (/& | ||
191 | & +.9815E-02_JPRB,+.8436E-02_JPRB,+.1087E-02_JPRB,-.2717E-02_JPRB,-.1755E-02_JPRB,& | ||
192 | & -.1559E-03_JPRB,+.2367E-03_JPRB,+.8808E-04_JPRB,+.2001E-05_JPRB,-.1244E-05_JPRB,& | ||
193 | & +.1041E-01_JPRB,+.8039E-02_JPRB,+.1005E-02_JPRB,-.1981E-02_JPRB,& | ||
194 | & -.1090E-02_JPRB,+.1595E-05_JPRB,+.1787E-03_JPRB,+.4644E-04_JPRB,-.1052E-04_JPRB,& | ||
195 | & +.6593E-02_JPRB,+.3983E-02_JPRB,-.1527E-03_JPRB,& | ||
196 | & -.1235E-02_JPRB,-.5078E-03_JPRB,+.3649E-04_JPRB,+.1005E-03_JPRB,+.3182E-04_JPRB,& | ||
197 | & +.3225E-02_JPRB,+.1672E-02_JPRB,& | ||
198 | & -.7752E-04_JPRB,-.4312E-03_JPRB,-.1872E-03_JPRB,-.1666E-04_JPRB,+.1872E-04_JPRB,& | ||
199 | & +.1133E-02_JPRB,& | ||
200 | & +.5643E-03_JPRB,+.7747E-04_JPRB,-.2980E-04_JPRB,-.2092E-04_JPRB,-.8590E-05_JPRB,& | ||
201 | & +.2988E-03_JPRB,+.6714E-04_JPRB,-.6249E-05_JPRB,+.1052E-04_JPRB,+.8790E-05_JPRB,& | ||
202 | & +.1569E-03_JPRB,-.1175E-04_JPRB,-.3033E-04_JPRB,-.9777E-06_JPRB,& | ||
203 | & +.1101E-03_JPRB,+.6827E-05_JPRB,-.1023E-04_JPRB,& | ||
204 | & +.4231E-04_JPRB,+.4905E-05_JPRB,& | ||
205 | & +.6229E-05_JPRB & | ||
206 | 1 | & /) | |
207 | |||
208 | ! ------------------------------------------------------------------ | ||
209 | 1 | RCAEROS= 0.1462E-16_JPRB | |
210 | ! ------------------------------------------------------------------ | ||
211 | |||
212 |
1/2✓ Branch 0 taken 1 times.
✗ Branch 1 not taken.
|
1 | IF (LHOOK) CALL DR_HOOK('SUAERH',1,ZHOOK_HANDLE) |
213 | 1 | END SUBROUTINE SUAERH | |
214 |