source: LMDZ5/trunk/libf/phymar/suaerh.F90 @ 2381

Last change on this file since 2381 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR

Integration of MAR physics
File size: 12.2 KB
Line 
1SUBROUTINE SUAERH
2
3
4!**** *SUAERH* - SETS UP THE AEROSOL HORIZONTAL DISTRIBUTION
5
6!     PURPOSE.
7!     --------
8
9!          THIS ROUTINE FETCHES VALUES OF A T10 SPECTRAL DISTRIBUTION
10!     FOR FOUR AEROSOL TYPES OF DIFFERENT ORIGINS (SEA,LAND,URBAN AREAS
11!     AND DESERTS). THE VALUES TO BE OBTAINED ARE BETWEEN ZERO AND ONE.
12
13!**   INTERFACE.
14!     ----------
15
16!          THERE ARE EIGHT DUMMY ARGUMENTS: *PAESC*, *PAESS*, *PAELC*,
17!     *PAELS*, *PAEUC*, *PAEUS*, *PAEDC* AND *PAEDS*ARE ARRAYS FOR THE
18!     T10 DISTRIBUTIONS (*S FOR SEA, *L FOR LAND, *U FOR URBAN AND *D
19!     FOR DESERT, *C FOR COSINE AND *S FOR SINE).
20
21!     METHOD.
22!     -------
23
24!          NONE.
25
26!     EXTERNALS.
27!     ----------
28
29!          NONE.
30
31!     REFERENCE.
32!     ----------
33
34!          NONE.
35
36!     AUTHOR
37!     ------
38!     J.-F. GELEYN      E.C.M.W.F.     29/09/82.
39
40!     MODIFICATIONS
41!     -------------
42!     J.-J. MORCRETTE   E.C.M.W.F.     91/07/14   ADAPTATION TO I.F.S.
43
44!     ------------------------------------------------------------------
45
46
47#include "tsmbkind.h"
48
49USE YOEAERD  , ONLY : RCAEROS  ,RAESC    ,RAESS    ,RAELC    ,&
50            &RAELS    ,RAEUC    ,RAEUS    ,RAEDC    ,RAEDS
51
52
53IMPLICIT NONE
54
55
56!     ------------------------------------------------------------------
57
58!*    DATA STATEMENTS.
59!     ---- -----------
60
61!          *RAE S/L/U/D C/S*  (SEE ABOVE).
62
63RAESC = (/&
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         &/)
82RAESS = (/&
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         &/)
99RAELC = (/&
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         &/)
118RAELS = (/&
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         &/)
135RAEUC = (/&
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         &/)
154RAEUS = (/&
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         &/)
171RAEDC = (/&
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         &/)
190RAEDS = (/&
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         &/)
207
208!     ------------------------------------------------------------------
209RCAEROS= 0.1462E-16_JPRB
210!     ------------------------------------------------------------------
211
212RETURN
213END SUBROUTINE SUAERH
Note: See TracBrowser for help on using the repository browser.