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

Language C Lines 722
MD5 Hash 3b36b2889328f52cb000464a266f38f2 Estimated Cost $15,691 (why?)
Repository https://bitbucket.org/olahlou/otb View Raw File View Project SPDX
  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
/* DUST.f -- translated by f2c (version 19970805).
   You must link the resulting object file with the libraries:
	-lf2c -lm   (in that order)
*/

#ifdef __cplusplus
extern "C" {
#endif
/* OTB patches: replace "f2c.h" by "otb_6S.h" */
/*#include "f2c.h"*/
#include "otb_6S.h"

/* Common Block Declarations */

Extern struct {
    integer nquad;
} num_quad__;

#define num_quad__1 num_quad__

Extern struct {
    doublereal ph[20000]	/* was [20][1000] */, qh[20000]	/* was [20][
	    1000] */, uh[20000]	/* was [20][1000] */;
} sixs_aerbas__;

#define sixs_aerbas__1 sixs_aerbas__

Extern struct {
    doublereal ex[80]	/* was [4][20] */, sc[80]	/* was [4][20] */, 
	    asy[80]	/* was [4][20] */, vi[4];
} sixs_coef__;

#define sixs_coef__1 sixs_coef__

/*<       subroutine   dust >*/
/* Subroutine */ int dust_()
{
    /* Initialized data */

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

    /* System generated locals */
    integer i__1;

    /* Local variables */
    integer i__, j;
    doublereal vi_m__;

/* - to vary the number of quadratures */
/*<       include "paramdef.inc" >*/
/*<       integer nquad >*/
/*<       parameter (nt_p_max=100,nqmax_p=1000,nqdef_p=83) ! do not change >*/
/* Attention */
/* mu2_p has to be equal to (mu_p-1)*2 */
/*<       common /num_quad/ nquad >*/
/*<       real ph,qh,uh >*/
/*<       common /sixs_aerbas/ ph(20,nqmax_p),qh(20,nqmax_p),uh(20,nqmax_p) >*/
/*<       real phr(20,nqdef_p),qhr(20,nqdef_p),uhr(20,nqdef_p) >*/
/* - to vary the number of quadratures */
/*<       real ex,sc,asy,vi >*/
/*<       common /sixs_coef/ ex(4,20),sc(4,20),asy(4,20),vi(4) >*/
/*<       real ex_m(20),sc_m(20),asy_m(20),vi_m >*/
/*<       integer i,j >*/
/* Dust model */
/* Parameters and refractive indices - G.A. d'Almeida, 'Atmospheric */
/* aerosols. Global climatology and radiative characteristics', 1991. */
/* Extinction coefficients are calculated in km-1 */
/* *************** MIE (asy) ****************** */
/*<    >*/
/* *************** MIE (ext&sca) ****************** */
/*<    >*/
/* ************** MIE (phase fun ph) *************** */
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/* ************** MIE (phase fun qh) *************** */
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/* ************** MIE (phase fun uh) *************** */
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<    >*/
/*<       vi_m=18.2859898 >*/
    vi_m__ = 18.2859898;
/*<       vi(1)=vi_m >*/
    sixs_coef__1.vi[0] = vi_m__;
/*<       do 1 i=1,20 >*/
    for (i__ = 1; i__ <= 20; ++i__) {
/*<        asy(1,i)=asy_m(i) >*/
	sixs_coef__1.asy[(i__ << 2) - 4] = asy_m__[i__ - 1];
/*<        ex(1,i)=ex_m(i) >*/
	sixs_coef__1.ex[(i__ << 2) - 4] = ex_m__[i__ - 1];
/*<        sc(1,i)=sc_m(i) >*/
	sixs_coef__1.sc[(i__ << 2) - 4] = sc_m__[i__ - 1];
/*<       do 1 j=1,nquad >*/
	i__1 = num_quad__1.nquad;
	for (j = 1; j <= i__1; ++j) {
/*<       ph(i,j)=phr(i,j) >*/
	    sixs_aerbas__1.ph[i__ + j * 20 - 21] = phr[i__ + j * 20 - 21];
/*<       qh(i,j)=qhr(i,j) >*/
	    sixs_aerbas__1.qh[i__ + j * 20 - 21] = qhr[i__ + j * 20 - 21];
/*<       uh(i,j)=uhr(i,j) >*/
	    sixs_aerbas__1.uh[i__ + j * 20 - 21] = uhr[i__ + j * 20 - 21];
/*<     1 continue >*/
/* L1: */
	}
    }
/*<       return >*/
    return 0;
/*<       end >*/
} /* dust_ */

#ifdef __cplusplus
	}
#endif
Back to Top