PageRenderTime 59ms CodeModel.GetById 18ms RepoModel.GetById 0ms app.codeStats 0ms

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