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 |
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 |
IF (LHOOK) CALL DR_HOOK('SUAERH',1,ZHOOK_HANDLE) |
213 |
|
1 |
END SUBROUTINE SUAERH |