PageRenderTime 143ms CodeModel.GetById 91ms app.highlight 175ms RepoModel.GetById 37ms app.codeStats 1ms

/Utilities/otb6S/6SV1.0B/DUST.c

https://bitbucket.org/olahlou/otb
C | 721 lines | 596 code | 14 blank | 111 comment | 2 complexity | 3b36b2889328f52cb000464a266f38f2 MD5 | raw file
  1/* DUST.f -- translated by f2c (version 19970805).
  2   You must link the resulting object file with the libraries:
  3	-lf2c -lm   (in that order)
  4*/
  5
  6#ifdef __cplusplus
  7extern "C" {
  8#endif
  9/* OTB patches: replace "f2c.h" by "otb_6S.h" */
 10/*#include "f2c.h"*/
 11#include "otb_6S.h"
 12
 13/* Common Block Declarations */
 14
 15Extern struct {
 16    integer nquad;
 17} num_quad__;
 18
 19#define num_quad__1 num_quad__
 20
 21Extern struct {
 22    doublereal ph[20000]	/* was [20][1000] */, qh[20000]	/* was [20][
 23	    1000] */, uh[20000]	/* was [20][1000] */;
 24} sixs_aerbas__;
 25
 26#define sixs_aerbas__1 sixs_aerbas__
 27
 28Extern struct {
 29    doublereal ex[80]	/* was [4][20] */, sc[80]	/* was [4][20] */, 
 30	    asy[80]	/* was [4][20] */, vi[4];
 31} sixs_coef__;
 32
 33#define sixs_coef__1 sixs_coef__
 34
 35/*<       subroutine   dust >*/
 36/* Subroutine */ int dust_()
 37{
 38    /* Initialized data */
 39
 40    static doublereal asy_m__[20] = { .846,.838,.836,.832,.827,.826,.822,.817,
 41	    .812,.807,.802,.8,.792,.785,.778,.795,.808,.859,.892,.853 };
 42    static doublereal ex_m__[20] = { 8.230388,8.305723,8.328883,8.368396,
 43	    8.415968,8.42865,8.467864,8.519962,8.583242,8.632867,8.689027,
 44	    8.704309,8.803438,8.924636,9.300344,9.404915,9.368683,8.811917,
 45	    7.876924,7.002555 };
 46    static doublereal sc_m__[20] = { 5.5657,5.744132,5.786949,5.888061,
 47	    5.976042,6.028808,6.112733,6.215003,6.333235,6.450123,6.547941,
 48	    6.603794,6.771705,7.00242,7.719834,8.043072,8.08431,7.738148,
 49	    6.921192,6.187184 };
 50    static doublereal phr[1660]	/* was [20][83] */ = { .3243,.3504,.3625,
 51	    .3853,.4323,.4322,.4427,.496,.5164,.5181,.5654,.5712,.5994,.5925,
 52	    .3548,.236,.2183,.1136,.0644,.07657,.316,.341,.3534,.3726,.413,
 53	    .4125,.4235,.4724,.4898,.4938,.5352,.5403,.5672,.5612,.3322,.2178,
 54	    .1988,.1007,.05458,.07193,.3591,.3912,.4034,.4175,.4481,.4415,
 55	    .4548,.4868,.5048,.5082,.5364,.5328,.5564,.5433,.3156,.1964,.17,
 56	    .08035,.04514,.06322,.3476,.3864,.3995,.4135,.4391,.4331,.4487,
 57	    .4722,.492,.4987,.5185,.5126,.5333,.5175,.313,.194,.166,.08849,
 58	    .05591,.06634,.3226,.3531,.3637,.3787,.403,.4048,.4149,.4366,
 59	    .4535,.4599,.4779,.4783,.4909,.4768,.3146,.2077,.1823,.1064,
 60	    .05989,.07486,.2902,.3218,.3363,.3464,.3689,.3649,.3774,.3942,
 61	    .4158,.4226,.4371,.4312,.4489,.4387,.3242,.2326,.204,.1102,.05586,
 62	    .07879,.2596,.2983,.3122,.3207,.3346,.3299,.3463,.3518,.3775,
 63	    .3885,.3962,.3866,.4085,.4069,.3406,.2537,.2133,.1035,.05049,
 64	    .07773,.2419,.2747,.2775,.2931,.2956,.3059,.3185,.3219,.3353,
 65	    .3506,.3542,.3578,.3701,.3811,.3519,.2599,.2101,.09485,.04611,
 66	    .07451,.2232,.239,.2391,.255,.265,.2776,.2823,.2946,.3016,.3106,
 67	    .3207,.3308,.3386,.3563,.3468,.2539,.201,.08717,.04284,.07098,
 68	    .1846,.1973,.2059,.2139,.2299,.2307,.2372,.2565,.2648,.2717,.2835,
 69	    .287,.3018,.3207,.3281,.2416,.1911,.08152,.04062,.06784,.1428,
 70	    .1651,.1719,.1797,.1895,.1899,.2001,.2072,.2222,.2321,.2396,.2392,
 71	    .2571,.2746,.2985,.2274,.1825,.07739,.03894,.06531,.1168,.1363,
 72	    .1389,.1483,.1499,.1577,.1664,.1679,.1798,.1928,.1962,.2008,.2129,
 73	    .2299,.2611,.2136,.1742,.07472,.03775,.06342,.09709,.108,.1102,
 74	    .1176,.1248,.1289,.1339,.1415,.1499,.1568,.1648,.1685,.18,.1957,
 75	    .2251,.1981,.1662,.07332,.0369,.06207,.08103,.08892,.09256,.09732,
 76	    .1054,.1059,.1105,.1205,.1266,.1316,.1396,.1417,.1531,.1668,.1957,
 77	    .1817,.1576,.07256,.03647,.06114,.07004,.0784,.08045,.0853,.08953,
 78	    .09193,.09644,.1019,.1077,.1139,.119,.1217,.1307,.1423,.171,.1655,
 79	    .148,.07272,.03634,.06053,.0632,.07077,.0712,.07638,.07772,.08229,
 80	    .08606,.08837,.09357,.1003,.1033,.1072,.1135,.1235,.15,.1485,
 81	    .1385,.07329,.03638,.06011,.05814,.06358,.06493,.06854,.07041,
 82	    .07362,.07676,.07963,.08409,.08989,.09258,.09605,.1015,.1101,
 83	    .1331,.1334,.1282,.07381,.03668,.05981,.05326,.05933,.06051,
 84	    .06361,.06557,.06746,.07074,.073,.07747,.08219,.08497,.08694,
 85	    .09286,.1005,.1188,.1212,.1175,.07415,.03732,.0595,.05106,.0551,
 86	    .05595,.05873,.06245,.06315,.06528,.06925,.07286,.07531,.07979,
 87	    .0809,.08691,.09362,.1074,.1102,.1079,.07411,.03816,.05911,.0494,
 88	    .05225,.05397,.05571,.05991,.05925,.0613,.06632,.06935,.07122,
 89	    .07564,.07628,.08209,.08802,.09941,.1012,.09907,.07328,.03919,
 90	    .0585,.0487,.05194,.05247,.05508,.05656,.05816,.06028,.06339,
 91	    .06534,.06883,.07122,.07342,.07726,.08284,.0932,.09373,.09203,
 92	    .07137,.04024,.0577,.04778,.05138,.0516,.05421,.05516,.05746,
 93	    .05938,.06103,.06333,.06691,.06867,.07136,.07418,.07917,.0884,
 94	    .08759,.08621,.0694,.04118,.05674,.0482,.05099,.05119,.05365,
 95	    .05353,.05736,.05857,.05995,.06142,.06558,.06649,.07074,.07177,
 96	    .07662,.0847,.08329,.08117,.06678,.04189,.05577,.04788,.05019,
 97	    .05069,.05277,.05403,.05622,.0575,.06034,.06153,.06441,.06643,
 98	    .06956,.07148,.07603,.08216,.08024,.07705,.06417,.04234,.05485,
 99	    .04772,.05009,.05103,.05261,.05539,.05591,.05712,.05995,.06279,
100	    .06427,.0676,.06862,.07242,.07642,.08044,.0778,.07415,.06139,
101	    .04253,.05396,.04784,.04996,.05174,.05262,.05675,.05538,.05693,
102	    .06075,.06407,.06474,.06877,.06817,.0734,.07699,.07981,.07591,
103	    .07245,.0588,.04238,.05316,.04869,.05225,.05266,.05491,.05634,
104	    .05713,.05924,.06106,.06367,.06628,.06842,.06912,.07313,.07689,
105	    .07995,.07502,.07116,.05642,.04189,.05249,.04996,.05356,.05382,
106	    .05623,.05705,.0592,.06095,.06174,.06438,.06768,.06902,.07099,
107	    .07361,.07716,.0809,.075,.07028,.05471,.04135,.05211,.05202,
108	    .05519,.05511,.05787,.05789,.06115,.06269,.06349,.06546,.06929,
109	    .07022,.07334,.0749,.0785,.08223,.07553,.06998,.05351,.04079,.052,
110	    .0528,.05641,.05686,.05919,.06072,.06213,.06402,.06621,.06841,
111	    .07105,.07329,.07487,.07806,.08175,.08385,.07636,.07034,.05247,
112	    .04021,.0521,.05469,.05805,.05922,.06104,.06362,.06397,.06583,
113	    .06858,.07172,.07377,.07671,.07713,.08153,.0851,.08613,.07738,
114	    .07157,.05169,.03982,.05246,.05593,.06066,.06272,.06405,.06673,
115	    .06597,.06862,.07111,.07521,.07767,.08028,.07946,.08518,.08877,
116	    .08903,.079,.07347,.05156,.0396,.05313,.0591,.06323,.06472,.06664,
117	    .06976,.06955,.07187,.07485,.07843,.08059,.08373,.08371,.08879,
118	    .09248,.09256,.08176,.07534,.05198,.03966,.0542,.06304,.06633,
119	    .06745,.06981,.07261,.07405,.0756,.07902,.08167,.08416,.08727,
120	    .08901,.09263,.09669,.09657,.08511,.07742,.05276,.04003,.05562,
121	    .06599,.06987,.07095,.07354,.07714,.07774,.07968,.084,.08663,
122	    .08859,.09258,.09343,.09821,.1025,.1008,.08852,.08053,.05375,
123	    .04059,.05735,.06985,.07371,.0754,.07771,.08242,.082,.08408,
124	    .08898,.09257,.09403,.0988,.09861,.1046,.1088,.1064,.09255,.08427,
125	    .0552,.0414,.05946,.07474,.07894,.08057,.08325,.08687,.08777,
126	    .08997,.0945,.09782,.1004,.1043,.1053,.1104,.1148,.1129,.09769,
127	    .0881,.05741,.04262,.06211,.08012,.08533,.08557,.08967,.09191,
128	    .09513,.09726,.1006,.1034,.1069,.1104,.1133,.1168,.1216,.1194,
129	    .1034,.0926,.0601,.04421,.06525,.08582,.09124,.09202,.09601,
130	    .09852,.1018,.104,.1081,.1107,.1145,.118,.1214,.1248,.1299,.1269,
131	    .1099,.09819,.06304,.0461,.06887,.09162,.09803,.09913,.1033,.1069,
132	    .1088,.1117,.1161,.1199,.1234,.1276,.1294,.1347,.1397,.1367,.1173,
133	    .1043,.06663,.04841,.07309,.09846,.1067,.1079,.1124,.1152,.1176,
134	    .1213,.1238,.1294,.1341,.1375,.1387,.1448,.1498,.148,.1256,.1111,
135	    .07119,.05129,.07808,.1009,.1112,.1131,.1174,.1203,.1215,.1263,
136	    .1283,.1349,.14,.143,.1433,.1505,.1554,.1532,.1302,.1151,.07373,
137	    .05293,.08087,.105,.1168,.118,.1232,.1243,.1269,.1323,.1329,.1395,
138	    .146,.148,.1492,.1558,.161,.1587,.1349,.1194,.07642,.05471,.08386,
139	    .1138,.1253,.128,.1323,.1357,.1371,.1421,.1444,.1518,.1573,.1607,
140	    .1617,.1688,.1741,.17,.1455,.1291,.0822,.05867,.09045,.1246,.1344,
141	    .1375,.142,.1478,.148,.1526,.1591,.1649,.1691,.1747,.1753,.1836,
142	    .1896,.1837,.1573,.1399,.08908,.06332,.09808,.1362,.1455,.1498,
143	    .1538,.161,.1606,.1652,.1731,.1796,.1837,.19,.19,.1994,.2057,
144	    .2002,.1705,.1522,.09734,.06885,.107,.1473,.1616,.165,.1704,.1739,
145	    .176,.1823,.186,.1944,.2016,.2054,.206,.2154,.2219,.2178,.1857,
146	    .1664,.1067,.07527,.1172,.1611,.1759,.1798,.1854,.1898,.1924,
147	    .1986,.2023,.2118,.219,.2236,.2246,.2344,.2414,.2371,.2032,.1824,
148	    .1172,.08266,.129,.1732,.1898,.1974,.2008,.2107,.2063,.2143,.2216,
149	    .2338,.2387,.2462,.2423,.2573,.2644,.2589,.2229,.2005,.1298,
150	    .09136,.1427,.1897,.2075,.2154,.2193,.2298,.2261,.2342,.2414,
151	    .2546,.2603,.2678,.2649,.2798,.2876,.2826,.2451,.2217,.1444,.1016,
152	    .1586,.2101,.2249,.2317,.2374,.2505,.2484,.2549,.2662,.2767,.2813,
153	    .2913,.2913,.3044,.3134,.3082,.2707,.2458,.1611,.1135,.177,.2304,
154	    .2462,.2535,.2597,.274,.2719,.2788,.2908,.3023,.3073,.318,.3184,
155	    .3322,.342,.3389,.3002,.273,.1805,.1275,.1984,.254,.2702,.2763,
156	    .2847,.2986,.2982,.3057,.319,.3296,.336,.347,.3479,.3629,.3743,
157	    .3742,.3331,.3042,.2037,.1442,.2236,.2787,.2948,.3028,.3109,.3278,
158	    .3259,.3338,.3495,.3616,.3677,.3804,.3802,.3976,.4103,.4136,.3703,
159	    .3406,.2305,.1638,.2531,.306,.3259,.3322,.3431,.357,.3598,.3684,
160	    .3812,.3939,.4031,.4145,.4174,.4335,.4478,.4562,.4136,.3818,.2616,
161	    .1869,.2875,.3378,.3559,.3624,.3746,.3904,.3945,.4027,.4204,.4303,
162	    .4403,.4529,.459,.4739,.491,.5038,.4634,.4287,.2986,.2147,.3284,
163	    .3659,.3928,.4018,.4137,.43,.4314,.4431,.4577,.4737,.4856,.4978,
164	    .5009,.5205,.5389,.5594,.5194,.4834,.3423,.2479,.3765,.4035,.4306,
165	    .4391,.4531,.4704,.4746,.4859,.5039,.5181,.5309,.5447,.5516,.57,
166	    .5917,.6209,.5833,.5473,.3937,.2875,.4337,.4393,.4753,.4861,.5002,
167	    .5168,.5205,.535,.55,.5693,.5851,.5982,.6042,.6259,.65,.6897,
168	    .6567,.6205,.4552,.3358,.502,.4858,.5264,.5328,.5529,.5611,.5758,
169	    .5914,.6028,.6201,.6432,.6525,.6657,.6838,.712,.7698,.7427,.7037,
170	    .528,.3941,.5834,.536,.5803,.5849,.6089,.6133,.6349,.6516,.6593,
171	    .6798,.7082,.7156,.7324,.7504,.7826,.8596,.8407,.8011,.6154,.4649,
172	    .6815,.5962,.633,.633,.6637,.6722,.6999,.7137,.7294,.745,.7732,
173	    .7858,.8081,.8258,.8644,.958,.9513,.9173,.7208,.5523,.7996,.6485,
174	    .6854,.6967,.7219,.7469,.7585,.7757,.8073,.8255,.849,.8704,.8848,
175	    .9151,.9597,1.068,1.08,1.051,.8468,.6595,.9427,.7072,.7571,.7683,
176	    .7969,.819,.836,.8562,.8843,.906,.9353,.9557,.9761,1.005,1.056,
177	    1.197,1.232,1.205,1.,.7922,1.117,.7761,.8406,.8482,.8836,.8945,
178	    .9231,.948,.9682,.994,1.034,1.05,1.076,1.106,1.166,1.345,1.408,
179	    1.387,1.187,.9585,1.329,.8649,.9153,.9209,.9625,.985,1.018,1.039,
180	    1.073,1.096,1.134,1.16,1.191,1.225,1.294,1.516,1.61,1.603,1.415,
181	    1.167,1.591,.9362,1.007,1.022,1.062,1.095,1.115,1.145,1.179,1.217,
182	    1.258,1.288,1.311,1.361,1.439,1.708,1.848,1.858,1.694,1.43,1.913,
183	    1.036,1.116,1.127,1.176,1.201,1.236,1.269,1.307,1.341,1.392,1.423,
184	    1.459,1.508,1.603,1.928,2.126,2.161,2.039,1.764,2.312,1.154,1.225,
185	    1.239,1.292,1.334,1.371,1.401,1.456,1.493,1.543,1.587,1.627,1.687,
186	    1.799,2.19,2.459,2.522,2.466,2.193,2.81,1.261,1.362,1.39,1.441,
187	    1.497,1.516,1.56,1.621,1.679,1.733,1.787,1.813,1.901,2.028,2.507,
188	    2.857,2.956,3.001,2.746,3.431,1.424,1.534,1.552,1.623,1.665,1.714,
189	    1.762,1.821,1.881,1.957,2.008,2.055,2.145,2.3,2.887,3.333,3.489,
190	    3.671,3.466,4.213,1.594,1.73,1.763,1.836,1.902,1.937,1.999,2.071,
191	    2.155,2.236,2.305,2.346,2.467,2.649,3.357,3.923,4.144,4.518,4.407,
192	    5.199,1.832,1.992,2.026,2.118,2.191,2.24,2.313,2.399,2.498,2.599,
193	    2.682,2.734,2.878,3.099,3.955,4.668,4.967,5.602,5.65,6.449,2.152,
194	    2.347,2.388,2.502,2.582,2.658,2.745,2.844,2.962,3.098,3.19,3.27,
195	    3.433,3.705,4.749,5.632,6.022,7.003,7.306,8.048,2.617,2.856,2.909,
196	    3.052,3.169,3.255,3.363,3.503,3.656,3.816,3.948,4.039,4.257,4.596,
197	    5.869,6.929,7.419,8.848,9.537,10.1,3.335,3.662,3.74,3.928,4.088,
198	    4.189,4.342,4.528,4.736,4.95,5.119,5.228,5.517,5.941,7.494,8.742,
199	    9.362,11.34,12.59,12.75,4.608,5.08,5.173,5.455,5.645,5.828,6.041,
200	    6.274,6.551,6.868,7.074,7.25,7.603,8.148,10.,11.44,12.19,14.81,
201	    16.84,16.25,7.076,7.789,7.947,8.364,8.69,8.933,9.239,9.62,10.01,
202	    10.43,10.74,10.98,11.45,12.14,14.21,15.76,16.59,19.9,22.94,20.95,
203	    12.68,13.92,14.19,14.87,15.39,15.74,16.24,16.8,17.4,18.,18.42,
204	    18.69,19.33,20.12,21.96,23.21,24.03,27.93,32.14,27.4,28.42,30.64,
205	    31.13,32.27,33.19,33.7,34.42,35.29,36.08,36.73,37.25,37.49,38.06,
206	    38.51,38.05,37.66,38.01,41.85,47.11,36.86,89.63,93.19,93.92,95.21,
207	    96.27,96.5,96.88,97.25,97.21,96.55,96.09,95.41,93.86,90.71,77.68,
208	    70.09,68.38,69.31,74.26,50.85,546.3,522.3,516.2,499.2,484.4,474.2,
209	    459.7,441.3,421.1,399.6,382.1,370.6,342.6,304.7,205.8,159.3,147.5,
210	    130.1,125.1,65.26,5621.,4201.,3938.,3362.,2954.,2722.,2420.,2096.,
211	    1798.,1542.,1362.,1262.,1035.,791.2,361.8,235.3,207.6,165.9,149.9,
212	    69.8 };
213    static doublereal qhr[1660]	/* was [20][83] */ = { 1.268e-22,8.915e-22,
214	    -6.146e-22,-8.679e-22,5.549e-22,3.355e-22,2.263e-22,-1.172e-21,
215	    -5.743e-22,-4.904e-22,-1.071e-22,1.195e-21,-2.802e-22,7.227e-22,
216	    3.341e-22,4.569e-23,3.013e-23,-8.339e-23,-1.652e-22,-4.718e-23,
217	    .02476,.01047,.01969,.01142,.02641,.01874,.0132,.02621,.02658,
218	    .01744,.02586,.02337,.02409,.01962,1.668e-4,-.003494,-.006586,
219	    -.003001,-3.713e-5,-3.367e-4,.02255,.008794,.01817,.01278,.04434,
220	    .02124,.01774,.05035,.04799,.02926,.05239,.04146,.05438,.05306,
221	    .01905,.004438,-.006701,-.001501,.003969,9.818e-4,-.008459,
222	    -.01221,-.00623,-.007503,.01074,-2.99e-4,-1.102e-4,.01893,.02133,
223	    .01488,.02991,.02593,.03805,.04545,.03846,.01387,.003526,.009853,
224	    .00787,.005409,-.037,-.02791,-.03044,-.02461,-.03317,-.02018,
225	    -.01613,-.02101,-.01628,-.003591,-.005289,.004162,.007287,.02058,
226	    .0506,.01729,.01083,.01559,.004872,.007205,-.05391,-.04559,
227	    -.05209,-.04284,-.05365,-.03582,-.03241,-.03886,-.03471,-.01959,
228	    -.02169,-.008725,-.006781,.009166,.04795,.0162,.008649,.01043,
229	    6.642e-4,.004406,-.05697,-.05726,-.05398,-.05172,-.04567,-.04506,
230	    -.04148,-.03113,-.02654,-.01971,-.01274,-.008606,.002331,.01769,
231	    .03824,.01017,-.001759,.002577,-.00274,-4.729e-4,-.04863,-.04209,
232	    -.03391,-.03515,-.01946,-.03142,-.02489,-.0147,-8.06e-4,3.157e-4,
233	    .01198,.004025,.02494,.03513,.03434,-.00219,-.01359,-.003925,
234	    -.004686,-.005499,-.0216,-.009986,-.007471,-.003311,-.00409,
235	    -.001718,.006124,.00629,.01513,.02603,.02723,.02814,.03896,.04622,
236	    .03286,-.01421,-.02471,-.01012,-.006046,-.01012,.002527,.008919,
237	    .003891,.01311,.008491,.02246,.02474,.02117,.02537,.03606,.03617,
238	    .04608,.04654,.05277,.02928,-.02101,-.03343,-.01511,-.007327,
239	    -.01419,.007879,.01561,.01574,.02056,.01844,.02391,.02857,.02847,
240	    .03273,.0422,.04163,.04622,.05034,.05581,.0295,-.02372,-.0376,
241	    -.02013,-.008648,-.01783,.008199,.01519,.02093,.02022,.02527,
242	    .02167,.02612,.02779,.03695,.0412,.04372,.04163,.05053,.05464,
243	    .03271,-.02141,-.03836,-.02543,-.009896,-.02108,.01009,.01543,
244	    .01697,.01896,.02248,.0219,.02471,.02789,.03215,.03473,.03843,
245	    .03901,.04468,.04927,.03308,-.01531,-.0367,-.02928,-.01121,
246	    -.02385,.01381,.01286,.01185,.0154,.01768,.02166,.0214,.02736,
247	    .02589,.02781,.03191,.03714,.03791,.04313,.03264,-.01001,-.03092,
248	    -.03297,-.01258,-.02614,.0106,.01151,.01116,.01383,.01647,.0187,
249	    .01898,.02288,.02356,.02484,.02863,.03139,.03378,.03835,.03297,
250	    -.004689,-.02429,-.03565,-.01432,-.02769,.007513,.01267,.01147,
251	    .01481,.01332,.01663,.01882,.01786,.01992,.02373,.02415,.02669,
252	    .02862,.03267,.03119,.001363,-.01861,-.037,-.01624,-.02862,
253	    .006562,.009975,.01049,.01219,.01159,.01396,.01571,.01574,.01756,
254	    .02128,.02113,.02362,.02494,.02844,.02856,.005275,-.01251,-.03797,
255	    -.01813,-.02888,.004469,.007821,.009795,.01017,.01144,.01111,
256	    .01304,.01373,.01675,.01941,.01978,.0202,.02303,.026,.02664,
257	    .007626,-.008017,-.03668,-.01998,-.02848,.004342,.006247,.007055,
258	    .008202,.009839,.009925,.01121,.01353,.01476,.0165,.01789,.01838,
259	    .02103,.02379,.02484,.008819,-.004221,-.03487,-.02148,-.02766,
260	    .00387,.005102,.005709,.006948,.008875,.008866,.009889,.0123,
261	    .01358,.01463,.01657,.01688,.01951,.02197,.02333,.008749,-.001206,
262	    -.03197,-.02289,-.02647,.003488,.005426,.005082,.006972,.006834,
263	    .009426,.01011,.01004,.01168,.0138,.01451,.01596,.01736,.01969,
264	    .02166,.008743,2.898e-4,-.02826,-.0236,-.02509,.001873,.004434,
265	    .004872,.00619,.006969,.007567,.008898,.009056,.01147,.01319,
266	    .0141,.01429,.01671,.01876,.01996,.008694,.001564,-.02552,-.02407,
267	    -.02369,8.242e-4,.0047,.004787,.006439,.005712,.007128,.008921,
268	    .007886,.01027,.01287,.0128,.01335,.01536,.01747,.018,.008697,
269	    .001966,-.02208,-.02368,-.02232,1.707e-4,.002735,.003997,.004648,
270	    .006149,.00571,.007182,.007899,.01043,.01178,.01285,.01243,.01519,
271	    .01693,.01671,.008481,.001962,-.01931,-.0229,-.02109,-5.39e-4,
272	    .002602,.003696,.004435,.005588,.005228,.006796,.007281,.01003,
273	    .01128,.01251,.01168,.01491,.0167,.01536,.007835,.002216,-.01667,
274	    -.02176,-.02002,-3.612e-4,1.81e-4,.001473,.001918,.00611,.004137,
275	    .004625,.008963,.01019,.009016,.0129,.01164,.01537,.01727,.01453,
276	    .007277,.002135,-.01454,-.02038,-.01907,-2.495e-4,.001072,.001427,
277	    .00266,.004177,.004865,.005367,.008128,.008618,.009107,.01137,
278	    .01192,.01399,.01623,.01374,.006824,.001819,-.01253,-.01895,
279	    -.01834,-.001048,-1.217e-4,.001349,.001885,.004484,.003448,
280	    .004419,.008272,.008749,.009229,.01142,.01125,.01388,.0159,.0138,
281	    .006472,.001428,-.01133,-.01772,-.0178,-.001244,.001517,.001024,
282	    .00319,.002991,.004953,.006144,.006622,.00764,.00954,.01031,
283	    .01151,.01281,.0147,.01355,.006108,.001478,-.01019,-.01648,
284	    -.01737,-.002457,4.575e-4,.001414,.002531,.003738,.00401,.005411,
285	    .006194,.008194,.009797,.01065,.01126,.01291,.01446,.01344,
286	    .005573,.001519,-.009347,-.01539,-.01708,-.001549,.001173,
287	    4.724e-4,.003014,.002014,.005156,.006249,.006127,.006869,.009819,
288	    .009588,.01192,.01206,.01374,.01313,.005435,.001486,-.008573,
289	    -.01445,-.0169,-.004067,3.59e-4,.001628,.002541,.004061,.003211,
290	    .005222,.006022,.008858,.01013,.01143,.01064,.01379,.01538,.01261,
291	    .005348,.001357,-.008007,-.01367,-.01683,-.001796,-.001717,
292	    -9.738e-4,3.263e-4,.003613,.003659,.003798,.008311,.00832,.008269,
293	    .01127,.01222,.01385,.01575,.01256,.005386,.001147,-.007655,
294	    -.01312,-.01692,-.002753,8.142e-4,7.231e-4,.002842,.002326,
295	    .004971,.006125,.005975,.00766,.01005,.0105,.01246,.01311,.01494,
296	    .01244,.005309,.001426,-.007264,-.01278,-.01707,-.004492,
297	    -8.557e-4,.001648,.001856,.004274,.003029,.004896,.006272,.009486,
298	    .01096,.0121,.01146,.01437,.01548,.01323,.005225,.001829,-.006791,
299	    -.01244,-.01725,-.002752,-1.025e-4,-3.753e-4,.002172,.002807,
300	    .004759,.005991,.007114,.008386,.01057,.01144,.01251,.01403,
301	    .01542,.01401,.005663,.001888,-.006614,-.01223,-.01755,-.004822,
302	    .001248,.001837,.003738,.003406,.004766,.007145,.005577,.009288,
303	    .01216,.01216,.01208,.01477,.01628,.01399,.006047,.002242,
304	    -.006553,-.01218,-.01795,-.005504,-.001358,.001687,.001736,
305	    .005575,.002813,.005129,.00745,.01118,.01194,.01408,.01183,.01656,
306	    .01779,.01446,.006261,.00305,-.00638,-.01215,-.01836,-.003245,
307	    -7.803e-4,-4.635e-4,.001778,.004299,.004695,.00589,.008602,.01032,
308	    .01127,.01369,.01354,.01649,.01797,.01532,.007047,.003673,-.00626,
309	    -.01213,-.01885,-.003857,7.688e-4,.001536,.003486,.003577,.005471,
310	    .007235,.007773,.01004,.01264,.01331,.01443,.01624,.01813,.01538,
311	    .008086,.004362,-.006197,-.01218,-.01937,-.005823,-7.587e-4,
312	    .001999,.002549,.0059,.003873,.006232,.008895,.01187,.01296,.015,
313	    .0143,.01775,.01948,.01533,.008532,.005665,-.006101,-.01234,
314	    -.01993,-.002124,.001022,5.311e-4,.003672,.003044,.007237,.008208,
315	    .008417,.009991,.01291,.01354,.01616,.01668,.01875,.01614,.008645,
316	    .006338,-.006076,-.01244,-.02026,-.003315,-.002572,5.462e-5,
317	    7.203e-4,.006663,.004518,.005417,.01028,.01288,.01248,.01626,
318	    .01507,.01908,.02064,.01667,.009205,.006978,-.006109,-.01256,
319	    -.02061,-.00196,-4.315e-4,-1.677e-4,.002449,.004768,.006216,
320	    .007206,.01081,.01206,.01347,.016,.01574,.01936,.02149,.01881,
321	    .01082,.007626,-.005996,-.01278,-.02125,-.003327,.001019,.001806,
322	    .003997,.00532,.007103,.008669,.009788,.01274,.01478,.01651,
323	    .01641,.01989,.02214,.02027,.01264,.008406,-.005687,-.01301,
324	    -.02188,-.005509,.001875,.002745,.005103,.006569,.007084,.00976,
325	    .009921,.01348,.01571,.01701,.01743,.02018,.02241,.02104,.01413,
326	    .009918,-.005619,-.01334,-.02266,-.003118,-4.092e-4,.001526,
327	    .003269,.0072,.007509,.008739,.01222,.01401,.01539,.01762,.0195,
328	    .02084,.02331,.02245,.01593,.0114,-.005624,-.01371,-.02344,
329	    -.001082,-6.854e-4,-3.013e-4,.002811,.007851,.00896,.009252,
330	    .01378,.01469,.01537,.01858,.0209,.02196,.0245,.02493,.01813,
331	    .01239,-.005103,-.01394,-.02414,3.212e-4,.002913,-2.097e-4,
332	    .005948,.005473,.01228,.01279,.0138,.01326,.01705,.01768,.02259,
333	    .02164,.02494,.02767,.02008,.01399,-.0046,-.01419,-.02489,.001378,
334	    .005484,.001492,.008544,.004245,.01379,.01504,.014,.01307,.01936,
335	    .01778,.02357,.02215,.02608,.03004,.0216,.01624,-.004475,-.01459,
336	    -.02573,-4.333e-4,.00629,.005193,.009975,.006273,.01267,.01547,
337	    .01373,.01588,.02237,.02046,.02261,.02476,.02845,.03304,.0238,
338	    .01771,-.003889,-.01482,-.02651,-.002347,.007964,.007723,.01191,
339	    .007705,.0131,.01708,.01269,.0178,.02453,.02224,.02303,.02642,
340	    .02964,.03633,.02666,.01879,-.002786,-.01493,-.02712,-.001399,
341	    .008576,.006748,.0123,.006867,.01566,.01849,.0136,.01682,.02419,
342	    .02141,.0264,.02579,.02964,.03829,.02888,.02146,-.002153,-.01521,
343	    -.02791,-.002813,.0087,.007801,.01287,.007293,.01553,.01906,
344	    .01369,.01716,.0255,.02159,.02738,.02594,.03016,.03985,.03036,
345	    .02455,-.001486,-.01542,-.02868,-.004009,.009681,.008641,.01418,
346	    .006801,.01539,.0203,.01264,.01738,.02762,.02189,.0267,.02639,
347	    .03082,.0427,.03324,.02607,5.661e-6,-.01534,-.02913,-.006005,
348	    .007736,.008715,.01277,.008444,.01335,.01875,.01214,.01952,.0283,
349	    .02413,.02505,.02864,.03255,.04578,.03643,.02811,.001504,-.0154,
350	    -.02965,-.005521,.006657,.006408,.01147,.006439,.01344,.01797,
351	    .01224,.01837,.02694,.02357,.02563,.02866,.03317,.04765,.03835,
352	    .03165,.002807,-.0154,-.03011,-.005109,.001938,.002116,.006827,
353	    .00586,.01152,.01404,.01369,.01788,.0236,.02375,.02602,.02937,
354	    .0346,.04812,.04051,.0345,.004649,-.01501,-.03031,-.007876,
355	    5.279e-4,-3.111e-4,.005626,.003223,.01025,.01334,.01118,.01547,
356	    .02251,.02155,.02495,.02752,.03345,.04968,.04391,.03637,.006591,
357	    -.0147,-.03038,-.008446,-.005814,-.007032,-7.231e-4,3.031e-4,
358	    .008236,.008714,.01219,.01178,.01713,.01819,.02565,.02461,.03161,
359	    .05035,.04708,.03955,.009116,-.0142,-.03013,-.01587,-.01153,
360	    -.009027,-.005651,2.25e-4,.001398,.003267,.01015,.01121,.01383,
361	    .01756,.02134,.02407,.0314,.05016,.04807,.0436,.01215,-.01317,
362	    -.02967,-.02224,-.01495,-.01203,-.008878,-.003539,-.003515,
363	    -4.783e-4,.005442,.007955,.01067,.01458,.01729,.02125,.0285,
364	    .05161,.05076,.04602,.01426,-.01214,-.02892,-.02776,-.01527,
365	    -.01565,-.009875,-.01139,-.005077,-.001286,-.003573,.001397,
366	    .008326,.008183,.01336,.01523,.02306,.05113,.05418,.04828,.01796,
367	    -.01083,-.02755,-.02752,-.02602,-.02859,-.02113,-.0182,-.01056,
368	    -.01079,-.003473,-.005306,-.00201,.00276,.01032,.01101,.02084,
369	    .04874,.05633,.05033,.0228,-.008337,-.02562,-.03971,-.03755,
370	    -.0326,-.03108,-.01868,-.02413,-.02215,-.01,-.006195,-.007234,
371	    .001497,-2.555e-4,.009503,.01854,.04621,.0568,.05366,.02606,
372	    -.006306,-.02328,-.05038,-.04496,-.03945,-.03878,-.02639,-.03163,
373	    -.03003,-.01888,-.01378,-.01587,-.005853,-.006704,.002281,.01152,
374	    .04301,.05788,.05725,.03027,-.003129,-.02036,-.05438,-.05354,
375	    -.05263,-.04865,-.03951,-.03768,-.03895,-.02498,-.02609,-.02868,
376	    -.01693,-.01122,-.007148,.005739,.03629,.0568,.0602,.03632,
377	    .001217,-.01595,-.06782,-.07206,-.06436,-.06617,-.04185,-.0544,
378	    -.05641,-.03189,-.02939,-.04142,-.02007,-.02425,-.01019,.001794,
379	    .03227,.05861,.05813,.04196,.005563,-.01081,-.08576,-.07231,
380	    -.06736,-.0658,-.05879,-.06391,-.05817,-.05337,-.04444,-.04308,
381	    -.03546,-.03827,-.02554,-.01378,.02879,.05881,.05912,.04816,
382	    .01181,-.004895,-.09301,-.08516,-.08161,-.07857,-.07798,-.07361,
383	    -.07001,-.06653,-.0624,-.05402,-.05307,-.04682,-.04225,-.0281,
384	    .02171,.0561,.06185,.05429,.01875,.002784,-.1103,-.103,-.09714,
385	    -.09628,-.08801,-.09088,-.08741,-.07995,-.07218,-.06868,-.06244,
386	    -.06168,-.05133,-.03814,.02015,.05159,.06386,.06013,.02696,.01122,
387	    -.1294,-.1099,-.1062,-.1027,-.1085,-.1041,-.09479,-.1041,-.08915,
388	    -.0747,-.0791,-.0788,-.06732,-.05394,.01023,.04816,.06292,.06981,
389	    .03742,.02192,-.1358,-.1343,-.1375,-.1287,-.1288,-.1144,-.1152,
390	    -.1105,-.1097,-.1022,-.0968,-.0822,-.08226,-.06431,-.001314,
391	    .04755,.05653,.07773,.04885,.0338,-.1531,-.1605,-.1647,-.1558,
392	    -.1413,-.135,-.1404,-.1158,-.1203,-.127,-.1035,-.09563,-.08529,
393	    -.0636,-.004303,.04797,.05427,.08519,.06171,.04555,-.1867,-.1845,
394	    -.1723,-.1757,-.1492,-.1664,-.1631,-.1343,-.124,-.1337,-.1067,
395	    -.1193,-.08863,-.06949,-.01217,.04573,.05557,.09469,.07738,.05906,
396	    -.1947,-.211,-.212,-.2041,-.1788,-.1729,-.1828,-.1388,-.1513,
397	    -.1619,-.129,-.115,-.1055,-.07644,-.01671,.04337,.06068,.1001,
398	    .09339,.0714,-.2355,-.2412,-.2262,-.2289,-.185,-.2081,-.2092,
399	    -.1523,-.1526,-.1699,-.1278,-.1342,-.1014,-.0697,-.004811,.04462,
400	    .06457,.1114,.1113,.08277,-.2765,-.2572,-.2233,-.235,-.1963,
401	    -.2327,-.2178,-.1769,-.1536,-.1539,-.1274,-.1454,-.1001,-.07127,
402	    3.506e-4,.05837,.0675,.116,.1272,.09052,-.2621,-.2445,-.2653,
403	    -.2259,-.2498,-.1805,-.1823,-.1698,-.1977,-.1567,-.1616,-.08839,
404	    -.124,-.07447,.01281,.09104,.0557,.1219,.1371,.09304,-.2857,
405	    -.2903,-.2455,-.2563,-.1505,-.207,-.2122,-.1072,-.08945,-.1272,
406	    -.05051,-.07265,-.01593,.01157,.09823,.1244,.07216,.1272,.1355,
407	    .08736,-.2604,-.2884,-.09803,-.2237,.06993,-.2145,-.2078,.01727,
408	    .1337,-.008717,.1614,-.03975,.1754,.1504,.1908,.142,.1193,.0856,
409	    .07581,.05598,.08863,.07441,.1603,.09307,.2651,.1464,.1096,.1725,
410	    .2616,.1571,.2539,.1465,.2386,.1925,.1472,.08351,.07938,.03885,
411	    .02173,.01352,1.341e-18,-6.003e-18,-7.181e-18,-3.242e-18,
412	    -4.398e-18,1.28e-18,2.688e-18,5.777e-18,7.57e-19,2.719e-18,
413	    -2.801e-19,-9.238e-19,-8.006e-19,8.456e-19,-4.975e-19,2.65e-19,
414	    -1.204e-19,-3.065e-19,9.486e-20,-1.058e-19 };
415    static doublereal uhr[1660]	/* was [20][83] */ = { -.3243,-.3504,-.3625,
416	    -.3853,-.4323,-.4322,-.4427,-.496,-.5164,-.5181,-.5654,-.5712,
417	    -.5994,-.5925,-.3548,-.236,-.2183,-.1136,-.0644,-.07657,-.2275,
418	    -.2583,-.2725,-.2956,-.3391,-.3413,-.3553,-.4083,-.4297,-.4384,
419	    -.4828,-.4907,-.5223,-.5233,-.3133,-.2085,-.1922,-.09834,-.05294,
420	    -.07138,-.01795,-.04707,-.05735,-.07782,-.1088,-.1128,-.1299,
421	    -.1697,-.1933,-.2108,-.2465,-.2569,-.2911,-.3052,-.1695,-.1093,
422	    -.1014,-.04735,-.02004,-.0524,.04439,.02789,.02289,.01077,
423	    -.003171,-.007923,-.02048,-.04294,-.05802,-.07562,-.09581,-.1062,
424	    -.1285,-.1376,-.04147,-.009279,-.008594,.00231,.00555,-.02943,
425	    .03578,.03269,.03027,.02519,.01441,.01379,.008951,-.005037,
426	    -.01449,-.02196,-.03659,-.04171,-.05458,-.05303,.02811,.0544,
427	    .05487,.03021,.01412,-.01269,.01267,.0135,.007699,.01022,2.126e-4,
428	    .007758,.004957,-.004217,-.01248,-.01282,-.02283,-.01919,-.0296,
429	    -.02092,.05365,.08631,.08811,.038,.01288,-.004242,-.01034,-.01829,
430	    -.02313,-.01947,-.01919,-.01261,-.01808,-.01286,-.02362,-.0263,
431	    -.0262,-.01819,-.0265,-.01634,.05498,.09693,.09656,.03552,.00989,
432	    -.001328,-.03552,-.04433,-.04096,-.04326,-.0334,-.03633,-.04058,
433	    -.03018,-.03299,-.03815,-.03239,-.02964,-.03047,-.02247,.0491,
434	    .09204,.09036,.0308,.007319,-8.748e-4,-.05009,-.0421,-.03741,
435	    -.03956,-.03951,-.04292,-.03879,-.03998,-.03497,-.03195,-.03282,
436	    -.03479,-.02948,-.02243,.04089,.07919,.07855,.02686,.005413,
437	    -.001035,-.04099,-.03186,-.03494,-.03079,-.03483,-.02978,-.02715,
438	    -.03348,-.02895,-.02503,-.0255,-.02321,-.02136,-.01528,.02886,
439	    .06369,.06551,.024,.003973,-.001153,-.02721,-.03047,-.03292,
440	    -.02897,-.02856,-.02259,-.024,-.0187,-.02253,-.02079,-.01786,
441	    -.01096,-.01299,-.00669,.01992,.04884,.05402,.02242,.003066,
442	    -.001031,-.02594,-.02724,-.02654,-.02497,-.01972,-.02031,-.02033,
443	    -.01341,-.01395,-.01406,-.009588,-.007853,-.005201,-6.13e-4,
444	    .01745,.03684,.04436,.02143,.002644,-6.918e-4,-.02533,-.02107,
445	    -.01916,-.01821,-.01776,-.01749,-.01474,-.01375,-.01104,-.007127,
446	    -.007054,-.006274,-.002967,5.223e-4,.0156,.0287,.03628,.02109,
447	    .002585,-9.753e-5,-.02248,-.0183,-.01894,-.0161,-.01724,-.01337,
448	    -.01181,-.01335,-.01042,-.006479,-.00642,-.003348,-.002263,
449	    .001365,.01214,.02269,.03001,.02134,.002752,7.33e-4,-.0215,
450	    -.01933,-.0192,-.01711,-.01587,-.01375,-.01288,-.01116,-.009425,
451	    -.007255,-.005374,-.002809,-.001152,.002894,.01022,.0186,.02499,
452	    .02156,.00326,.001688,-.02233,-.02083,-.01882,-.01834,-.01416,
453	    -.01573,-.01461,-.009871,-.008193,-.007829,-.004197,-.004043,
454	    -9.091e-5,.003881,.01022,.01696,.02102,.02184,.004136,.002725,
455	    -.02209,-.02007,-.01973,-.01796,-.01489,-.01481,-.01396,-.01029,
456	    -.009074,-.008612,-.005056,-.00411,-9.361e-4,.003137,.009971,
457	    .01574,.01899,.02211,.005316,.003713,-.02106,-.01973,-.01992,
458	    -.01765,-.01631,-.01421,-.01361,-.011,-.01039,-.008582,-.006425,
459	    -.003913,-.002342,.001673,.01037,.01479,.01769,.02217,.006561,
460	    .00464,-.02203,-.01893,-.01825,-.01662,-.01659,-.01481,-.01308,
461	    -.01255,-.01045,-.007332,-.006812,-.005076,-.002902,9.269e-4,
462	    .01086,.0149,.01649,.02205,.007916,.005533,-.022,-.01882,-.0188,
463	    -.01674,-.01658,-.01419,-.01296,-.0128,-.01053,-.007874,-.006876,
464	    -.005094,-.002954,9.24e-4,.01077,.01502,.01642,.02163,.009371,
465	    .006384,-.02227,-.01979,-.01896,-.01769,-.01541,-.01492,-.01394,
466	    -.0121,-.009603,-.008479,-.006017,-.005268,-.002203,.001605,
467	    .01081,.01533,.01624,.02157,.01068,.007328,-.02192,-.01974,
468	    -.01868,-.0175,-.01529,-.01523,-.01404,-.01134,-.009431,-.008158,
469	    -.005717,-.005347,-.001818,.002073,.01098,.01562,.01643,.02118,
470	    .01203,.008258,-.02216,-.02,-.01895,-.01776,-.01377,-.01567,
471	    -.01414,-.01036,-.008129,-.008214,-.004357,-.005671,-4.951e-4,
472	    .00327,.01189,.01582,.01715,.0207,.01312,.009143,-.02116,-.01834,
473	    -.01813,-.01617,-.01455,-.01377,-.01242,-.01091,-.00865,-.006986,
474	    -.004872,-.004336,-9.19e-4,.00303,.01298,.01641,.01773,.02043,
475	    .01394,.01007,-.02043,-.01736,-.01723,-.01504,-.01482,-.01284,
476	    -.01119,-.01013,-.008813,-.005735,-.005091,-.003083,-.001151,
477	    .002866,.01421,.01715,.01847,.02035,.01455,.01101,-.01964,-.01565,
478	    -.01625,-.0134,-.01458,-.011,-.009432,-.009804,-.008393,-.00453,
479	    -.004568,-.001539,-4.745e-4,.003714,.01496,.01824,.01943,.02048,
480	    .01502,.0121,-.01853,-.01667,-.01593,-.01434,-.01217,-.01088,
481	    -.01012,-.007304,-.006066,-.004489,-.0022,4.284e-5,.001872,
482	    .006027,.01611,.01966,.02052,.02086,.0155,.0133,-.01754,-.01492,
483	    -.01448,-.01256,-.01104,-.009663,-.008492,-.005993,-.004523,
484	    -.002666,-4.823e-4,.001385,.003777,.008167,.01744,.02102,.02214,
485	    .02118,.01597,.01452,-.01771,-.0147,-.01257,-.01194,-.007639,
486	    -.009634,-.00797,-.004045,-.001066,-6.568e-4,.002991,.001904,
487	    .007191,.01134,.01963,.0226,.02411,.02152,.01638,.01583,-.01472,
488	    -.01251,-.01157,-.009849,-.00712,-.006501,-.005433,-.002252,
489	    -7.197e-5,.001308,.004298,.005102,.008798,.0132,.0222,.02504,
490	    .02589,.02245,.01691,.01733,-.01352,-.01047,-.009595,-.007649,
491	    -.005115,-.004269,-.002901,5.06e-5,.001953,.00363,.006428,.008036,
492	    .01108,.01583,.02459,.02775,.0277,.0236,.0175,.01902,-.009567,
493	    -.008317,-.009181,-.005715,-.003964,-2.647e-4,4.081e-6,.003228,
494	    .003462,.00545,.00834,.01239,.01324,.01807,.02753,.03076,.03023,
495	    .02479,.01823,.02087,-.009255,-.004739,-.004405,-.001671,
496	    -9.363e-4,.001572,.003413,.00512,.007174,.01033,.01206,.01458,
497	    .01705,.0218,.03145,.03411,.03331,.02609,.01908,.02289,-.006918,
498	    -.003179,-.002159,2.362e-4,.003409,.003747,.00573,.009411,.01161,
499	    .01366,.01669,.01793,.02177,.02645,.03648,.03785,.03648,.02761,
500	    .02012,.02513,-.002846,.002167,.001466,.0055,.004963,.009521,
501	    .01153,.01177,.01407,.01868,.01942,.02313,.02486,.02995,.04139,
502	    .04213,.04012,.02962,.02138,.02771,-.001439,.005453,.00707,
503	    .009494,.01041,.01195,.01514,.016,.01987,.02418,.02527,.02711,
504	    .03082,.03603,.04626,.04695,.04457,.03204,.02297,.03066,.004171,
505	    .006826,.008585,.01093,.01626,.0161,.01751,.02377,.02597,.02711,
506	    .03205,.03347,.038,.0434,.05233,.05213,.04944,.03453,.02479,
507	    .03391,.008762,.01329,.01377,.01738,.02074,.02243,.02428,.02858,
508	    .03146,.03375,.03803,.03973,.0445,.0504,.05903,.05795,.05497,
509	    .03741,.02682,.03753,.0109,.01912,.02235,.02426,.02805,.02677,
510	    .03088,.03357,.03938,.04277,.04582,.04503,.05224,.05794,.06685,
511	    .06486,.06108,.04108,.02919,.04172,.0184,.0237,.02656,.02912,
512	    .03467,.03418,.03701,.04242,.04644,.04937,.0534,.05426,.06014,
513	    .06595,.07655,.07288,.06756,.04526,.03198,.04649,.02884,.0312,
514	    .03128,.03616,.03989,.045,.04562,.05172,.05292,.05631,.0608,
515	    .06591,.06843,.0753,.08663,.08128,.07513,.04974,.03513,.05182,
516	    .02866,.03354,.03672,.03928,.04751,.04617,.0483,.05684,.06022,
517	    .06125,.06801,.06919,.07551,.08223,.09116,.08575,.07953,.05223,
518	    .03685,.05475,.03381,.04061,.04147,.04609,.04928,.05283,.05523,
519	    .06051,.06336,.0671,.07166,.07529,.07977,.08715,.09689,.09042,
520	    .08416,.055,.03872,.05789,.0396,.04866,.05297,.05536,.06204,
521	    .06081,.06469,.06992,.07667,.0797,.08498,.08497,.09305,.1,.109,
522	    .1016,.09372,.06123,.04291,.06486,.05157,.05869,.0606,.06555,
523	    .0721,.07332,.07658,.08332,.0878,.09125,.09705,.09891,.1059,.1134,
524	    .1238,.1145,.104,.06833,.04782,.07282,.06634,.06984,.07029,.07693,
525	    .08275,.08783,.08945,.09825,.1,.1041,.1104,.1154,.1202,.1289,
526	    .1401,.128,.1166,.07631,.05345,.08192,.07588,.08732,.0862,.09465,
527	    .09633,.1031,.1072,.1103,.1155,.1221,.1264,.1306,.1368,.1456,
528	    .1576,.1429,.1312,.08562,.05995,.09229,.08726,.1039,.1049,.1124,
529	    .1117,.1186,.1249,.1248,.133,.1425,.1444,.1477,.1554,.1646,.1777,
530	    .1605,.1468,.09653,.06758,.1043,.09928,.1181,.1238,.1284,.1331,
531	    .1329,.1413,.1436,.1551,.1631,.1667,.166,.1779,.187,.2002,.1804,
532	    .1645,.109,.07645,.1182,.1166,.1346,.1412,.1458,.1533,.1523,.1601,
533	    .166,.1766,.1837,.1891,.1895,.2012,.2113,.2237,.2024,.1859,.1234,
534	    .08673,.1341,.1371,.1507,.1578,.1627,.1769,.1735,.1794,.1921,
535	    .2015,.2047,.2154,.2151,.2285,.2396,.2495,.2278,.2099,.1403,
536	    .09874,.1527,.1568,.1733,.1821,.1869,.2014,.1964,.2045,.216,.2285,
537	    .2335,.2433,.2408,.2573,.2686,.2811,.2574,.2362,.1598,.1129,.1743,
538	    .182,.1985,.2055,.213,.226,.2247,.2329,.2439,.2559,.2629,.2723,
539	    .2716,.2878,.3003,.3166,.2903,.2671,.1824,.1295,.1994,.2087,.2222,
540	    .2309,.238,.2561,.2529,.2599,.2774,.2885,.2932,.3065,.3057,.3235,
541	    .3377,.355,.3266,.3036,.2092,.149,.2288,.2357,.2526,.2616,.2697,
542	    .2881,.2857,.2936,.3101,.3234,.3292,.343,.3426,.3616,.3771,.3976,
543	    .3693,.3446,.2404,.1722,.2633,.2673,.2853,.2932,.3039,.3213,.3218,
544	    .3306,.3485,.3596,.3681,.3811,.3843,.4016,.4195,.446,.4191,.3907,
545	    .2769,.1999,.304,.297,.3221,.3324,.3427,.3611,.3602,.3712,.3871,
546	    .4029,.4129,.426,.4279,.448,.4673,.5007,.475,.4449,.3203,.2329,
547	    .352,.3342,.3598,.3703,.3819,.4039,.4026,.4133,.4347,.4494,.4581,
548	    .475,.478,.4997,.5221,.5613,.538,.5082,.3718,.2725,.4092,.3698,
549	    .4059,.4198,.4307,.451,.4483,.4634,.4801,.5017,.5142,.5291,.5297,
550	    .5559,.5802,.6306,.6111,.5807,.4326,.3206,.4772,.418,.4596,.4665,
551	    .4853,.4938,.5066,.5224,.5322,.5508,.5726,.5819,.5929,.6121,.6406,
552	    .7101,.6965,.6637,.5048,.3785,.5583,.4708,.5129,.5173,.5405,.5458,
553	    .5661,.5817,.5914,.61,.6361,.6447,.661,.6787,.7117,.7984,.7935,
554	    .7606,.5921,.4492,.6563,.5307,.5654,.5667,.5953,.6068,.6303,.6434,
555	    .6611,.6772,.7018,.7165,.7352,.7554,.7941,.8966,.9035,.8755,.697,
556	    .5363,.7741,.584,.6197,.6307,.6551,.6804,.6903,.7071,.7383,.7563,
557	    .778,.7996,.8123,.8431,.8877,1.006,1.032,1.008,.8221,.643,.9167,
558	    .6446,.6905,.7018,.7291,.7523,.7684,.7865,.8166,.8359,.8627,.884,
559	    .9047,.9325,.9832,1.133,1.183,1.162,.9752,.7755,1.091,.7125,.774,
560	    .7826,.8157,.8294,.8541,.8777,.9,.9253,.9614,.9791,1.002,1.034,
561	    1.093,1.28,1.358,1.342,1.161,.9413,1.303,.8015,.8514,.8565,.8973,
562	    .9173,.9506,.9707,1.002,1.024,1.063,1.086,1.117,1.149,1.217,1.451,
563	    1.558,1.557,1.388,1.149,1.564,.876,.9416,.9551,.9948,1.025,1.049,
564	    1.076,1.109,1.142,1.183,1.211,1.239,1.282,1.36,1.64,1.794,1.812,
565	    1.667,1.411,1.886,.9746,1.049,1.059,1.107,1.132,1.168,1.197,1.238,
566	    1.267,1.315,1.346,1.384,1.43,1.524,1.857,2.071,2.113,2.011,1.745,
567	    2.284,1.09,1.159,1.172,1.224,1.265,1.3,1.329,1.383,1.419,1.466,
568	    1.51,1.548,1.606,1.717,2.117,2.403,2.472,2.437,2.174,2.782,1.199,
569	    1.294,1.32,1.371,1.425,1.446,1.486,1.546,1.601,1.652,1.705,1.734,
570	    1.816,1.942,2.431,2.797,2.904,2.972,2.726,3.402,1.36,1.463,1.479,
571	    1.549,1.592,1.639,1.684,1.745,1.801,1.872,1.925,1.97,2.058,2.212,
572	    2.808,3.271,3.435,3.639,3.444,4.184,1.526,1.657,1.691,1.76,1.825,
573	    1.858,1.917,1.989,2.069,2.149,2.215,2.254,2.372,2.553,3.274,3.859,
574	    4.088,4.486,4.385,5.17,1.764,1.917,1.946,2.038,2.104,2.16,2.228,
575	    2.31,2.401,2.502,2.579,2.638,2.772,2.991,3.864,4.599,4.909,5.569,
576	    5.628,6.419,2.077,2.262,2.3,2.412,2.491,2.567,2.649,2.749,2.858,
577	    2.988,3.08,3.163,3.318,3.59,4.649,5.559,5.96,6.969,7.283,8.02,
578	    2.534,2.764,2.815,2.954,3.067,3.154,3.257,3.395,3.541,3.698,3.826,
579	    3.919,4.13,4.467,5.758,6.851,7.354,8.813,9.514,10.07,3.247,3.565,
580	    3.638,3.824,3.976,4.083,4.231,4.409,4.609,4.824,4.986,5.103,5.378,
581	    5.8,7.369,8.653,9.292,11.3,12.56,12.73,4.514,4.972,5.059,5.34,
582	    5.526,5.714,5.92,6.153,6.42,6.732,6.937,7.12,7.463,8.007,9.871,
583	    11.35,12.12,14.77,16.81,16.22,6.968,7.674,7.829,8.243,8.562,8.808,
584	    9.113,9.488,9.875,10.3,10.6,10.84,11.31,12.,14.08,15.66,16.51,
585	    19.87,22.92,20.93,12.56,13.79,14.06,14.74,15.25,15.61,16.11,16.67,
586	    17.26,17.86,18.28,18.55,19.19,19.99,21.85,23.12,23.96,27.9,32.12,
587	    27.39,28.29,30.5,30.97,32.12,33.04,33.57,34.27,35.16,35.93,36.58,
588	    37.12,37.37,37.93,38.39,37.97,37.6,37.95,41.82,47.1,36.85,89.49,
589	    93.05,93.77,95.07,96.14,96.37,96.75,97.14,97.09,96.43,95.98,95.31,
590	    93.76,90.63,77.64,70.06,68.35,69.3,74.26,50.85,546.2,522.2,516.1,
591	    499.1,484.4,474.1,459.6,441.2,421.,399.6,382.1,370.6,342.5,304.7,
592	    205.7,159.3,147.5,130.1,125.1,65.26,5621.,4201.,3938.,3362.,2954.,
593	    2722.,2420.,2096.,1798.,1542.,1362.,1262.,1035.,791.2,361.8,235.3,
594	    207.6,165.9,149.9,69.8 };
595
596    /* System generated locals */
597    integer i__1;
598
599    /* Local variables */
600    integer i__, j;
601    doublereal vi_m__;
602
603/* - to vary the number of quadratures */
604/*<       include "paramdef.inc" >*/
605/*<       integer nquad >*/
606/*<       parameter (nt_p_max=100,nqmax_p=1000,nqdef_p=83) ! do not change >*/
607/* Attention */
608/* mu2_p has to be equal to (mu_p-1)*2 */
609/*<       common /num_quad/ nquad >*/
610/*<       real ph,qh,uh >*/
611/*<       common /sixs_aerbas/ ph(20,nqmax_p),qh(20,nqmax_p),uh(20,nqmax_p) >*/
612/*<       real phr(20,nqdef_p),qhr(20,nqdef_p),uhr(20,nqdef_p) >*/
613/* - to vary the number of quadratures */
614/*<       real ex,sc,asy,vi >*/
615/*<       common /sixs_coef/ ex(4,20),sc(4,20),asy(4,20),vi(4) >*/
616/*<       real ex_m(20),sc_m(20),asy_m(20),vi_m >*/
617/*<       integer i,j >*/
618/* Dust model */
619/* Parameters and refractive indices - G.A. d'Almeida, 'Atmospheric */
620/* aerosols. Global climatology and radiative characteristics', 1991. */
621/* Extinction coefficients are calculated in km-1 */
622/* *************** MIE (asy) ****************** */
623/*<    >*/
624/* *************** MIE (ext&sca) ****************** */
625/*<    >*/
626/* ************** MIE (phase fun ph) *************** */
627/*<    >*/
628/*<    >*/
629/*<    >*/
630/*<    >*/
631/*<    >*/
632/*<    >*/
633/*<    >*/
634/*<    >*/
635/*<    >*/
636/*<    >*/
637/*<    >*/
638/*<    >*/
639/*<    >*/
640/*<    >*/
641/*<    >*/
642/*<    >*/
643/*<    >*/
644/*<    >*/
645/*<    >*/
646/*<    >*/
647/* ************** MIE (phase fun qh) *************** */
648/*<    >*/
649/*<    >*/
650/*<    >*/
651/*<    >*/
652/*<    >*/
653/*<    >*/
654/*<    >*/
655/*<    >*/
656/*<    >*/
657/*<    >*/
658/*<    >*/
659/*<    >*/
660/*<    >*/
661/*<    >*/
662/*<    >*/
663/*<    >*/
664/*<    >*/
665/*<    >*/
666/*<    >*/
667/*<    >*/
668/* ************** MIE (phase fun uh) *************** */
669/*<    >*/
670/*<    >*/
671/*<    >*/
672/*<    >*/
673/*<    >*/
674/*<    >*/
675/*<    >*/
676/*<    >*/
677/*<    >*/
678/*<    >*/
679/*<    >*/
680/*<    >*/
681/*<    >*/
682/*<    >*/
683/*<    >*/
684/*<    >*/
685/*<    >*/
686/*<    >*/
687/*<    >*/
688/*<    >*/
689/*<       vi_m=18.2859898 >*/
690    vi_m__ = 18.2859898;
691/*<       vi(1)=vi_m >*/
692    sixs_coef__1.vi[0] = vi_m__;
693/*<       do 1 i=1,20 >*/
694    for (i__ = 1; i__ <= 20; ++i__) {
695/*<        asy(1,i)=asy_m(i) >*/
696	sixs_coef__1.asy[(i__ << 2) - 4] = asy_m__[i__ - 1];
697/*<        ex(1,i)=ex_m(i) >*/
698	sixs_coef__1.ex[(i__ << 2) - 4] = ex_m__[i__ - 1];
699/*<        sc(1,i)=sc_m(i) >*/
700	sixs_coef__1.sc[(i__ << 2) - 4] = sc_m__[i__ - 1];
701/*<       do 1 j=1,nquad >*/
702	i__1 = num_quad__1.nquad;
703	for (j = 1; j <= i__1; ++j) {
704/*<       ph(i,j)=phr(i,j) >*/
705	    sixs_aerbas__1.ph[i__ + j * 20 - 21] = phr[i__ + j * 20 - 21];
706/*<       qh(i,j)=qhr(i,j) >*/
707	    sixs_aerbas__1.qh[i__ + j * 20 - 21] = qhr[i__ + j * 20 - 21];
708/*<       uh(i,j)=uhr(i,j) >*/
709	    sixs_aerbas__1.uh[i__ + j * 20 - 21] = uhr[i__ + j * 20 - 21];
710/*<     1 continue >*/
711/* L1: */
712	}
713    }
714/*<       return >*/
715    return 0;
716/*<       end >*/
717} /* dust_ */
718
719#ifdef __cplusplus
720	}
721#endif