/Utilities/otb6S/6SV1.0B/STM.c
C | 707 lines | 581 code | 14 blank | 112 comment | 2 complexity | 8cd2693e49052b582be84d915cd3e46f 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
- /* STM.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 stm >*/
- /* Subroutine */ int stm_()
- {
- /* Initialized data */
- static doublereal asy_m__[20] = { .787,.795,.796,.8,.802,.804,.806,.808,
- .809,.81,.81,.809,.807,.802,.77,.736,.723,.685,.647,.469 };
- static doublereal ex_m__[20] = { 1.323426,1.367378,1.375664,1.389866,
- 1.397082,1.397282,1.394358,1.384032,1.363314,1.332585,1.301248,
- 1.278673,1.211386,1.101411,.720679,.5124845,.4508164,.3256089,
- .2397076,.06640586 };
- static doublereal sc_m__[20] = { 1.323426,1.367378,1.375664,1.389866,
- 1.397082,1.397282,1.394358,1.384032,1.363314,1.332585,1.301248,
- 1.278673,1.211386,1.101411,.720679,.5124845,.4508164,.3256089,
- .2397076,.06640586 };
- static doublereal phr[1660] /* was [20][83] */ = { .2975,.2458,.2347,
- .2072,.1879,.1755,.1638,.1489,.1415,.1337,.1284,.1253,.1225,.12,
- .128,.142,.1485,.1709,.1989,.4039,.2696,.2255,.2165,.1924,.1756,
- .165,.1551,.1422,.1359,.1293,.1248,.1223,.1201,.1185,.1274,.1416,
- .1482,.1707,.1987,.4037,.2051,.1731,.1684,.1526,.1427,.1366,.1306,
- .1236,.1198,.1163,.114,.113,.1128,.1132,.1252,.14,.1468,.1696,
- .1978,.4031,.1882,.1522,.1476,.1348,.1286,.1235,.1178,.1141,.1103,
- .1075,.1066,.1062,.1067,.1081,.122,.1375,.1445,.1677,.1962,.402,
- .2239,.1759,.1687,.1519,.1445,.1363,.1276,.1226,.1156,.1102,.1085,
- .1073,.1062,.1062,.1188,.1346,.1417,.1653,.1942,.4004,.2685,.2136,
- .2047,.1832,.1724,.1615,.1496,.141,.1302,.1214,.1173,.1148,.1107,
- .1076,.1161,.1315,.1387,.1626,.1917,.3985,.2867,.2382,.2294,.2073,
- .1934,.1828,.1696,.1576,.1448,.1339,.1275,.1239,.1169,.1106,.114,
- .1285,.1356,.1596,.1889,.3961,.2805,.2436,.2353,.2165,.2016,.193,
- .1804,.167,.154,.1427,.1349,.1308,.122,.1135,.1121,.1255,.1325,
- .1564,.1859,.3934,.2628,.2355,.2287,.2133,.2005,.1932,.182,.1696,
- .1574,.1464,.1385,.1342,.1249,.1151,.1102,.1226,.1294,.1532,.1828,
- .3903,.2453,.222,.2174,.2042,.1948,.1873,.1777,.1674,.1565,.1459,
- .1387,.1344,.1252,.1151,.1082,.1197,.1264,.15,.1796,.387,.2311,
- .2096,.2058,.1942,.1868,.1796,.1711,.1624,.1524,.1427,.1362,.1322,
- .1236,.1137,.1059,.1168,.1234,.1468,.1764,.3834,.2201,.2001,.1961,
- .1853,.1778,.1718,.1641,.1559,.1467,.138,.132,.1284,.1204,.1112,
- .1034,.114,.1204,.1437,.1732,.3797,.2119,.1918,.1878,.1771,.1693,
- .1638,.1567,.1489,.1404,.1325,.127,.1238,.1164,.108,.1008,.1111,
- .1175,.1406,.1701,.3758,.2031,.1835,.1798,.169,.1614,.1558,.1491,
- .1417,.1341,.1267,.1216,.1187,.1121,.1044,.09812,.1084,.1147,
- .1377,.1671,.3718,.1933,.1749,.1713,.1608,.1539,.148,.1416,.1347,
- .1279,.1208,.1162,.1133,.1076,.1006,.09541,.1057,.112,.1349,.1643,
- .3678,.1849,.1656,.1619,.1523,.1462,.1402,.1342,.1278,.1215,.115,
- .1109,.1082,.1029,.09679,.09278,.1032,.1095,.1323,.1616,.3638,
- .1764,.1561,.1522,.1436,.1377,.133,.1271,.1213,.115,.1093,.1057,
- .1035,.09837,.09304,.09029,.1009,.1072,.13,.1592,.3598,.1653,
- .1466,.1431,.135,.1296,.1258,.1203,.1149,.109,.1039,.1007,.0988,
- .09419,.08947,.08797,.09875,.1051,.1278,.1569,.3559,.152,.138,
- .1349,.1273,.1221,.1183,.1139,.1088,.1037,.09893,.09599,.09428,
- .09039,.08618,.08585,.09687,.1033,.126,.155,.3522,.1401,.1298,
- .1269,.1204,.1149,.1115,.1078,.103,.09871,.09453,.09168,.0901,
- .08685,.08326,.08397,.09528,.1018,.1244,.1534,.3487,.1305,.1211,
- .1188,.1131,.1084,.1055,.102,.09774,.0941,.0905,.08792,.08647,
- .08359,.0807,.08238,.09398,.1005,.1231,.152,.3454,.1216,.1129,
- .1111,.1063,.1026,.1002,.09688,.0932,.0899,.08683,.08467,.08337,
- .08068,.07847,.0811,.09298,.09963,.1222,.1511,.3425,.1134,.1061,
- .1045,.1006,.0975,.09529,.09251,.08922,.08621,.08361,.08184,
- .08071,.07822,.07653,.08013,.09231,.09906,.1217,.1506,.34,.1062,
- .1004,.0992,.09577,.09308,.09099,.08873,.08576,.08313,.08086,
- .07938,.07843,.07626,.0749,.07948,.09197,.09884,.1215,.1504,.3379,
- .1001,.09533,.09465,.09133,.08933,.08735,.08542,.08295,.08063,
- .07857,.07732,.07649,.07474,.07362,.07912,.09199,.09899,.1218,
- .1508,.3364,.09488,.09133,.09069,.08774,.08627,.08439,.0827,
- .08069,.07866,.07673,.07569,.07493,.07361,.07273,.07908,.0924,
- .09953,.1226,.1516,.3355,.0912,.08829,.08751,.08502,.08377,.08211,
- .08061,.07886,.07713,.07539,.07448,.0738,.0728,.07225,.07939,
- .09323,.1005,.1238,.153,.3352,.0889,.08568,.08517,.08293,.08174,
- .08043,.079,.07756,.07601,.07453,.07368,.07311,.07231,.07216,
- .08008,.09449,.1019,.1256,.1549,.3358,.08686,.0837,.08347,.08154,
- .08043,.07921,.07789,.07675,.07535,.07411,.07333,.07282,.0722,
- .07243,.08118,.09622,.1038,.1279,.1576,.3372,.08502,.08269,.08239,
- .08091,.07987,.07857,.07737,.07632,.07517,.07409,.07341,.07295,
- .07251,.07305,.08271,.09843,.1062,.1309,.1609,.3396,.08471,.08238,
- .08178,.08072,.07963,.07864,.07742,.07636,.07538,.07445,.07386,
- .07353,.07322,.07399,.08468,.1011,.1092,.1345,.165,.343,.08559,
- .08235,.08146,.08075,.07973,.07928,.07789,.07698,.07595,.07518,
- .07469,.07455,.07434,.0753,.0871,.1044,.1129,.1389,.1699,.3476,
- .08624,.08286,.08207,.08128,.08058,.08015,.07881,.07807,.07705,
- .07637,.076,.07594,.07591,.07704,.09,.1083,.1171,.1441,.1758,
- .3536,.0863,.08428,.08392,.08266,.08203,.08126,.08032,.07955,
- .07869,.07808,.07779,.07772,.07795,.07925,.09343,.1129,.1221,
- .1502,.1828,.361,.0877,.0865,.08605,.08473,.08377,.08312,.08242,
- .08152,.08075,.08021,.08003,.08,.08044,.08194,.09745,.1182,.1279,
- .1573,.1909,.37,.09053,.08908,.08827,.08716,.08592,.08568,.08499,
- .08411,.08326,.08279,.08273,.08283,.08339,.08516,.1022,.1245,
- .1347,.1656,.2003,.3808,.0936,.09184,.09132,.09004,.0889,.08859,
- .08796,.08729,.08638,.08594,.08597,.0862,.08688,.08896,.1076,
- .1316,.1425,.175,.2111,.3934,.09716,.09539,.09505,.09369,.09265,
- .09215,.09154,.091,.09012,.08972,.08983,.09014,.09097,.0934,.114,
- .1399,.1515,.1859,.2235,.4082,.1019,.0999,.09934,.09813,.09692,
- .0966,.09587,.0953,.0945,.09418,.09435,.09471,.09573,.09855,.1213,
- .1493,.1618,.1983,.2377,.4253,.1071,.1049,.1045,.1032,.1021,.1016,
- .1008,.1003,.09966,.09941,.09965,.09998,.1012,.1045,.1297,.1601,
- .1736,.2125,.2539,.4449,.1127,.1105,.1105,.1089,.1086,.1073,.1066,
- .1062,.1057,.1055,.1058,.1061,.1076,.1113,.1393,.1725,.1871,.2286,
- .2724,.4672,.1161,.1138,.1138,.1122,.112,.1107,.1098,.1094,.1091,
- .1088,.1092,.1095,.1111,.1151,.1447,.1793,.1945,.2375,.2826,.4795,
- .1202,.1174,.1174,.1158,.1155,.1143,.1134,.1128,.1126,.1124,.1128,
- .1131,.1149,.1192,.1504,.1866,.2025,.247,.2934,.4925,.1294,.126,
- .1255,.124,.1233,.1226,.1214,.1205,.1203,.1205,.1208,.1213,.1232,
- .1282,.1631,.2028,.2201,.2679,.3172,.5209,.1389,.1356,.1348,.1334,
- .1325,.1318,.1305,.1295,.1293,.1297,.1301,.1306,.1328,.1386,.1776,
- .2213,.2402,.2917,.3441,.5529,.1504,.1461,.1454,.1439,.1431,.1421,
- .1408,.1401,.1398,.1401,.1407,.1413,.1439,.1505,.1942,.2424,.2631,
- .3187,.3745,.5886,.1641,.1585,.1579,.156,.1552,.1542,.1527,.1522,
- .1518,.1521,.1529,.1536,.1566,.1642,.2133,.2666,.2892,.3495,.4089,
- .6283,.1786,.1734,.1725,.1704,.1693,.1679,.1665,.166,.1657,.1659,
- .167,.1678,.1714,.18,.2352,.2944,.319,.3844,.4478,.6724,.1954,
- .1899,.1893,.1866,.1854,.1836,.1824,.182,.1816,.1818,.1832,.1842,
- .1884,.1981,.2605,.3262,.3531,.4241,.4916,.7211,.2161,.2092,.2083,
- .2052,.2036,.2022,.2009,.2005,.1998,.2002,.2018,.2032,.208,.2191,
- .2898,.3627,.3922,.4692,.541,.7748,.2396,.2325,.2303,.2273,.2246,
- .224,.2225,.2215,.2207,.2215,.2233,.2251,.2308,.2435,.3237,.4046,
- .4369,.5204,.5967,.8339,.2653,.2584,.256,.2523,.2495,.2486,.2471,
- .2459,.2453,.2463,.2485,.2506,.2573,.2718,.363,.4529,.4882,.5785,
- .6595,.8986,.2963,.2868,.2851,.2805,.2789,.2769,.2752,.2745,.2739,
- .275,.2779,.2802,.2882,.305,.4088,.5084,.5471,.6445,.7301,.9693,
- .3324,.3206,.3183,.314,.3123,.3103,.3081,.3075,.307,.3086,.3121,
- .3148,.3243,.3439,.4621,.5725,.6147,.7194,.8095,1.046,.3735,.3611,
- .3579,.3537,.3508,.3492,.3467,.3457,.3457,.3481,.3521,.3554,.3665,
- .3898,.5243,.6464,.6923,.8044,.8987,1.13,.4204,.407,.4044,.3995,
- .3963,.3941,.3916,.3904,.3913,.3945,.3993,.4031,.4163,.4439,.5969,
- .7316,.7814,.9008,.9988,1.221,.475,.4599,.4576,.4524,.4492,.4465,
- .4438,.4433,.4449,.4491,.4548,.4594,.4752,.508,.6818,.8301,.8837,
- 1.01,1.111,1.319,.5394,.5231,.5197,.5147,.5101,.5083,.5057,.5058,
- .508,.5136,.5204,.5261,.5453,.584,.7814,.9437,1.001,1.134,1.237,
- 1.424,.6129,.5974,.5936,.5876,.5824,.5808,.5788,.5794,.5831,.5902,
- .5986,.6054,.6291,.6745,.8982,1.075,1.136,1.274,1.377,1.537,.6964,
- .683,.6804,.6726,.6689,.6659,.6651,.667,.6728,.6816,.6923,.7005,
- .7296,.7827,1.036,1.227,1.291,1.433,1.534,1.659,.7964,.7828,.7808,
- .7726,.7705,.7676,.7677,.7718,.7794,.7908,.8047,.8151,.8502,.9124,
- 1.197,1.402,1.469,1.611,1.708,1.788,.9181,.9016,.8981,.892,.89,
- .8897,.8912,.8971,.9069,.9226,.9404,.9538,.9954,1.069,1.388,1.604,
- 1.672,1.813,1.902,1.925,1.057,1.043,1.039,1.035,1.034,1.036,1.04,
- 1.047,1.061,1.083,1.105,1.122,1.171,1.258,1.612,1.837,1.905,2.039,
- 2.117,2.069,1.217,1.209,1.207,1.206,1.207,1.211,1.217,1.229,1.249,
- 1.278,1.306,1.326,1.385,1.487,1.875,2.105,2.171,2.293,2.354,2.222,
- 1.407,1.404,1.407,1.409,1.417,1.421,1.432,1.451,1.479,1.516,1.551,
- 1.574,1.646,1.764,2.185,2.412,2.474,2.577,2.615,2.382,1.633,1.639,
- 1.645,1.655,1.669,1.677,1.694,1.722,1.761,1.807,1.851,1.88,1.965,
- 2.101,2.549,2.765,2.818,2.894,2.901,2.548,1.9,1.927,1.934,1.955,
- 1.975,1.991,2.017,2.056,2.108,2.167,2.22,2.255,2.357,2.511,2.974,
- 3.167,3.207,3.244,3.212,2.721,2.216,2.274,2.287,2.322,2.351,2.376,
- 2.416,2.47,2.537,2.611,2.676,2.72,2.838,3.008,3.472,3.624,3.645,
- 3.631,3.55,2.898,2.599,2.692,2.716,2.769,2.817,2.854,2.911,2.984,
- 3.07,3.163,3.242,3.294,3.428,3.613,4.05,4.14,4.136,4.054,3.914,
- 3.079,3.062,3.206,3.242,3.324,3.4,3.452,3.53,3.627,3.736,3.85,
- 3.945,4.004,4.153,4.345,4.717,4.72,4.681,4.515,4.303,3.263,3.634,
- 3.85,3.899,4.025,4.133,4.205,4.308,4.433,4.569,4.707,4.815,4.881,
- 5.039,5.229,5.484,5.365,5.282,5.012,4.714,3.447,4.352,4.669,4.737,
- 4.921,5.062,5.16,5.293,5.448,5.615,5.776,5.892,5.962,6.121,6.29,
- 6.355,6.076,5.937,5.542,5.146,3.631,5.27,5.724,5.823,6.07,6.251,
- 6.379,6.547,6.735,6.928,7.103,7.221,7.29,7.433,7.55,7.333,6.85,
- 6.644,6.1,5.591,3.81,6.481,7.113,7.245,7.562,7.793,7.949,8.151,
- 8.367,8.573,8.746,8.853,8.91,9.008,9.033,8.417,7.679,7.393,6.678,
- 6.045,3.984,8.163,8.98,9.145,9.53,9.812,9.985,10.21,10.43,10.62,
- 10.77,10.84,10.87,10.88,10.75,9.594,8.551,8.172,7.267,6.499,4.15,
- 10.59,11.55,11.73,12.16,12.46,12.64,12.84,13.03,13.16,13.22,13.22,
- 13.19,13.05,12.7,10.84,9.446,8.964,7.852,6.942,4.304,14.17,15.15,
- 15.32,15.71,15.94,16.07,16.19,16.26,16.24,16.15,16.01,15.9,15.52,
- 14.85,12.13,10.34,9.745,8.416,7.364,4.444,19.56,20.22,20.29,20.46,
- 20.46,20.47,20.38,20.19,19.9,19.54,19.18,18.94,18.24,17.15,13.41,
- 11.19,10.49,8.942,7.749,4.567,27.63,27.22,27.05,26.63,26.18,25.9,
- 25.44,24.8,24.07,23.3,22.63,22.21,21.07,19.45,14.6,11.97,11.15,
- 9.406,8.086,4.671,39.22,36.37,35.76,34.19,32.98,32.2,31.12,29.84,
- 28.48,27.15,26.09,25.44,23.79,21.6,15.63,12.62,11.71,9.788,8.359,
- 4.753,54.07,46.93,45.59,42.35,40.09,38.6,36.73,34.64,32.56,30.59,
- 29.13,28.24,26.09,23.35,16.42,13.11,12.13,10.07,8.557,4.811,68.07,
- 55.89,53.77,48.85,45.6,43.45,40.85,38.07,35.4,32.92,31.15,30.07,
- 27.56,24.44,16.89,13.4,12.37,10.23,8.671,4.844,72.5,58.54,56.16,
- 50.7,47.14,44.78,41.96,38.99,36.14,33.52,31.67,30.54,27.93,24.71,
- 17.,13.46,12.42,10.27,8.698,4.852 };
- static doublereal qhr[1660] /* was [20][83] */ = { -5.308e-22,6.787e-22,
- 5.459e-22,-2.16e-22,-2.354e-22,-1.055e-22,1.958e-23,-1.308e-22,
- -8.726e-24,-4.912e-23,9.602e-23,2.191e-23,-1.487e-22,-1.464e-22,
- -3.694e-22,-3.902e-22,-5.161e-23,1.407e-22,2.306e-22,-1.974e-22,
- -.003847,-.001662,-8.407e-4,-5.478e-4,7.911e-5,-4.079e-4,4.175e-5,
- 5.053e-4,6.411e-4,6.492e-4,7.231e-4,6.181e-4,7.015e-4,6.29e-4,
- 3.515e-4,2.073e-4,1.652e-4,7.84e-5,1.529e-5,-1.772e-4,-.007671,
- -.003169,-4.804e-4,-3.328e-4,.001869,2.357e-4,.001454,.002848,
- .003394,.003305,.003548,.00319,.003409,.003055,.00174,.001033,
- 8.25e-4,3.879e-4,6.618e-5,-9.349e-4,-.001357,9.185e-4,.003452,
- .002842,.005575,.004274,.004916,.006212,.007113,.006903,.007227,
- .006894,.007016,.0064,.003816,.002299,.001839,8.488e-4,1.025e-4,
- -.002303,.004828,.003963,.003932,.00443,.005783,.006725,.006743,
- .007468,.008528,.008964,.009379,.009547,.009635,.009252,.006001,
- .003696,.002962,.001322,4.286e-5,-.004288,-.00222,-.003095,
- -.00457,-.001803,-.001751,.001356,.002535,.003675,.00534,.007371,
- .008188,.009106,.009834,.01043,.007727,.004902,.003931,.001656,
- -2.063e-4,-.006899,-.02526,-.02183,-.0215,-.01715,-.01525,-.01153,
- -.008462,-.005366,-.00199,.001685,.003559,.005059,.007324,.009451,
- .008609,.005663,.00453,.001711,-7.331e-4,-.01014,-.05228,-.04467,
- -.04196,-.03619,-.03073,-.0271,-.02263,-.01737,-.01168,-.006733,
- -.003404,-.001485,.002669,.006562,.008516,.00583,.004625,.001388,
- -.001609,-.01403,-.07486,-.06291,-.06035,-.05231,-.04557,-.04117,
- -.03546,-.02903,-.0219,-.01568,-.01127,-.008772,-.003157,.002449,
- .007521,.005367,.004163,6.31e-4,-.00288,-.01856,-.0878,-.07432,
- -.07299,-.06302,-.05729,-.05124,-.04487,-.03802,-.0307,-.02337,
- -.01855,-.01539,-.009005,-.002067,.005816,.004316,.003161,
- -5.766e-4,-.00457,-.02375,-.09331,-.07984,-.07823,-.06842,-.06318,
- -.05667,-.05053,-.04379,-.03643,-.02887,-.02402,-.02067,-.01392,
- -.006299,.003628,.00277,.001677,-.00222,-.006683,-.02957,-.09431,
- -.08006,-.07741,-.06906,-.06339,-.05834,-.05253,-.04637,-.03897,
- -.03203,-.02736,-.02437,-.01754,-.00987,.001173,8.345e-4,-2.09e-4,
- -.004263,-.009205,-.03605,-.08953,-.07549,-.07284,-.06578,-.06055,
- -.05661,-.05135,-.04593,-.03928,-.03319,-.02894,-.02639,-.01995,
- -.01269,-.001373,-.001384,-.002412,-.006662,-.01212,-.04316,
- -.07919,-.06806,-.06617,-.06007,-.05588,-.05225,-.04799,-.04323,
- -.03808,-.03286,-.02921,-.0269,-.02134,-.0148,-.003888,-.003793,
- -.004854,-.009368,-.01538,-.05089,-.06584,-.05945,-.05821,-.05313,
- -.0499,-.04638,-.04344,-.03939,-.03571,-.03147,-.02851,-.02638,
- -.02185,-.01624,-.006304,-.006318,-.007469,-.01233,-.01898,
- -.05925,-.05399,-.04954,-.04899,-.04509,-.04325,-.03989,-.03806,
- -.03532,-.03252,-.02924,-.02708,-.02531,-.02166,-.0171,-.008584,
- -.008901,-.0102,-.01552,-.02288,-.0682,-.04485,-.03957,-.03932,
- -.03694,-.03617,-.03397,-.03257,-.03117,-.02891,-.02652,-.02515,
- -.02393,-.02101,-.01748,-.0107,-.0115,-.01301,-.01889,-.02706,
- -.07773,-.03392,-.03153,-.03126,-.03008,-.02967,-.02861,-.02774,
- -.02684,-.02539,-.02388,-.02303,-.02224,-.0201,-.01758,-.01265,
- -.01409,-.01587,-.02241,-.03148,-.08783,-.02352,-.02509,-.02512,
- -.02451,-.0241,-.02347,-.02352,-.02267,-.02219,-.0215,-.02094,
- -.02042,-.01907,-.01751,-.01444,-.01667,-.01875,-.02606,-.03613,
- -.09849,-.01671,-.01935,-.01961,-.01964,-.01923,-.01903,-.01961,
- -.01905,-.01929,-.01926,-.019,-.0187,-.01799,-.01733,-.01609,
- -.01921,-.02165,-.02983,-.04099,-.1097,-.01125,-.01391,-.01464,
- -.01505,-.01552,-.01524,-.016,-.01609,-.0168,-.01712,-.01729,
- -.01712,-.01694,-.01705,-.01762,-.02172,-.02455,-.0337,-.04606,
- -.1214,-.006411,-.00939,-.01062,-.01111,-.01253,-.01208,-.01298,
- -.01366,-.0146,-.01519,-.01575,-.01569,-.01593,-.0167,-.01906,
- -.0242,-.02745,-.03767,-.0513,-.1335,-.003637,-.006723,-.007453,
- -.008535,-.009657,-.009792,-.01077,-.01151,-.01258,-.01354,
- -.01428,-.01445,-.01499,-.01632,-.02041,-.02666,-.03036,-.04174,
- -.05673,-.1462,-.002257,-.004789,-.005056,-.006538,-.007314,
- -.007959,-.008974,-.009679,-.01088,-.01214,-.01298,-.01335,
- -.01418,-.01596,-.02168,-.0291,-.03327,-.04592,-.06234,-.1592,
- -4.036e-4,-.002548,-.003402,-.004446,-.005819,-.006207,-.007248,
- -.008249,-.009605,-.01093,-.01191,-.01231,-.0135,-.01563,-.02288,
- -.03154,-.03619,-.0502,-.06812,-.1727,.001497,-.001083,-.002244,
- -.003014,-.004774,-.004666,-.005826,-.007079,-.008601,-.009904,
- -.011,-.01135,-.01293,-.01535,-.02404,-.03399,-.03913,-.05459,
- -.07408,-.1866,.002087,-4.758e-4,-.001182,-.002147,-.003636,
- -.00367,-.004808,-.006049,-.007678,-.009012,-.01015,-.01056,
- -.0124,-.0151,-.02518,-.03647,-.04211,-.0591,-.08021,-.2009,
- .001657,4.041e-4,-2.04e-4,-.001212,-.002503,-.003007,-.00389,
- -.005228,-.006826,-.008237,-.009345,-.009909,-.01191,-.01489,
- -.02632,-.03897,-.04513,-.06375,-.08653,-.2154,.002143,.00128,
- 4.527e-4,-5.118e-4,-.001788,-.002176,-.003038,-.004549,-.006152,
- -.007603,-.008677,-.009308,-.0115,-.01471,-.02746,-.04152,-.04821,
- -.06852,-.09304,-.2303,.003304,.001548,8.72e-4,-2.022e-4,-.00133,
- -.001423,-.00243,-.003857,-.005612,-.007065,-.008121,-.008767,
- -.01115,-.01456,-.02863,-.04411,-.05136,-.07344,-.09976,-.2455,
- .003609,.001781,.001457,2.086e-4,-7.149e-4,-.001069,-.001931,
- -.003234,-.005069,-.006523,-.007609,-.00835,-.01085,-.01442,
- -.02982,-.04676,-.0546,-.07852,-.1067,-.2608,.003288,.002567,
- .002287,9.618e-4,-6.828e-5,-5.962e-4,-.00132,-.002838,-.00453,
- -.005954,-.007149,-.008019,-.01058,-.0143,-.03103,-.04947,-.05792,
- -.08376,-.1138,-.2764,.003797,.003456,.002885,.001733,4.573e-4,
- 9.352e-5,-7.169e-4,-.002501,-.004033,-.005442,-.006763,-.007705,
- -.01035,-.0142,-.03228,-.05226,-.06136,-.08918,-.1212,-.2921,
- .004872,.003663,.003056,.002072,.001058,6.808e-4,-3.549e-4,
- -.001975,-.003557,-.005021,-.00642,-.007394,-.01016,-.01416,
- -.03357,-.05515,-.06492,-.0948,-.1288,-.3079,.006129,.00386,
- .003389,.002448,.001663,.001268,-1.914e-5,-.001434,-.003079,
- -.004591,-.006104,-.007096,-.009979,-.01414,-.03492,-.05815,
- -.06864,-.1006,-.1366,-.3238,.006769,.004835,.004199,.003363,
- .002294,.001913,5.437e-4,-.001045,-.002562,-.004111,-.005782,
- -.006827,-.009791,-.01415,-.03635,-.06129,-.07252,-.1066,-.1448,
- -.3397,.007335,.005663,.004951,.004169,.003129,.002525,.001106,
- -5.853e-4,-.001973,-.003656,-.005421,-.006575,-.009579,-.01417,
- -.03789,-.06459,-.07659,-.1129,-.1531,-.3555,.008643,.006322,
- .005611,.004819,.003931,.003139,.001661,1.863e-5,-.001383,
- -.003222,-.005057,-.006302,-.009358,-.01422,-.03954,-.06805,
- -.08087,-.1194,-.1617,-.3713,.009557,.007667,.006817,.005963,
- .004575,.003855,.002412,6.278e-4,-8.448e-4,-.002738,-.004709,
- -.005988,-.009143,-.01429,-.04131,-.07171,-.08539,-.1261,-.1706,
- -.3868,.01068,.008904,.008115,.007034,.005319,.004686,.003195,
- .001332,-2.873e-4,-.00228,-.004337,-.005596,-.008919,-.01439,
- -.04323,-.07558,-.09015,-.133,-.1797,-.4021,.01298,.009905,
- .008826,.0078,.006089,.005611,.004025,.002225,2.886e-4,-.001886,
- -.003936,-.005117,-.008682,-.01452,-.04528,-.07966,-.09515,-.1402,
- -.1891,-.417,.01395,.01067,.009268,.008361,.006411,.006124,
- .004532,.002699,5.708e-4,-.001658,-.003722,-.004872,-.008559,
- -.01459,-.04635,-.08179,-.09775,-.1439,-.1939,-.4243,.01447,
- .01162,.009985,.009085,.006748,.006635,.005086,.003165,8.697e-4,
- -.001385,-.003488,-.004635,-.008428,-.01467,-.04747,-.08398,
- -.1004,-.1477,-.1987,-.4315,.01504,.01323,.01195,.01044,.007799,
- .007467,.006131,.004089,.001622,-7.561e-4,-.002916,-.004161,
- -.008117,-.01481,-.0498,-.08854,-.1059,-.1555,-.2086,-.4456,
- .01661,.01408,.01341,.01124,.009324,.008339,.007058,.005098,
- .002563,-1.249e-4,-.002227,-.003613,-.007749,-.01495,-.05226,
- -.09335,-.1117,-.1635,-.2186,-.459,.01898,.01571,.01459,.01252,
- .01057,.009799,.008325,.006137,.003436,6.209e-4,-.001541,-.002982,
- -.007388,-.01511,-.05486,-.09842,-.1177,-.1718,-.2288,-.4717,
- .02061,.01827,.01651,.01454,.01163,.01147,.009937,.00714,.004277,
- .001598,-8.513e-4,-.002308,-.007034,-.01524,-.0576,-.1037,-.1239,
- -.1803,-.2391,-.4835,.02232,.01991,.01849,.01614,.01332,.01303,
- .01133,.008324,.005371,.002669,-2.911e-5,-.001572,-.006626,
- -.01535,-.0605,-.1093,-.1304,-.1891,-.2494,-.4945,.02514,.02125,
- .01987,.0176,.01535,.01479,.01269,.009743,.006689,.003765,
- 9.395e-4,-7.213e-4,-.006153,-.01543,-.06359,-.1151,-.1372,-.1981,
- -.2598,-.5044,.02842,.02362,.02171,.01983,.01722,.01684,.01443,
- .01118,.008146,.005048,.00203,2.08e-4,-.005615,-.01545,-.06688,
- -.1212,-.1442,-.2072,-.2702,-.5131,.03015,.02592,.0246,.02229,
- .01963,.01866,.01622,.01284,.009911,.006578,.003312,.001176,
- -.004972,-.01541,-.0704,-.1275,-.1514,-.2164,-.2805,-.5205,.03158,
- .02787,.0274,.02464,.02237,.02041,.01805,.01499,.01188,.008207,
- .004719,.002305,-.004239,-.01534,-.07414,-.1341,-.1589,-.2257,
- -.2906,-.5264,.03669,.03089,.02975,.02746,.0245,.02324,.02049,
- .01729,.01382,.009949,.006144,.003727,-.003435,-.01528,-.07812,
- -.1408,-.1666,-.235,-.3004,-.5308,.04103,.03509,.03303,.03088,
- .02689,.02671,.0236,.0196,.01599,.01199,.00777,.005314,-.002463,
- -.01518,-.0823,-.1478,-.1743,-.2441,-.3097,-.5335,.04226,.03916,
- .03747,.03451,.03038,.02937,.02675,.02248,.01869,.01435,.009783,
- .00698,-.001231,-.01506,-.08667,-.1549,-.1822,-.2531,-.3186,
- -.5344,.04606,.04221,.04136,.03771,.03432,.03233,.02983,.02611,
- .02173,.01684,.0121,.008974,2.256e-4,-.01491,-.09116,-.1621,-.19,
- -.2617,-.3268,-.5333,.05303,.04604,.04479,.04113,.0383,.03682,
- .03365,.0299,.02496,.01957,.01466,.01142,.00185,-.01471,-.09575,
- -.1693,-.1977,-.2698,-.3341,-.5302,.05746,.05197,.04979,.04632,
- .04288,.04179,.03846,.0337,.02862,.02287,.01759,.01417,.003656,
- -.01436,-.1004,-.1763,-.2052,-.2774,-.3404,-.5249,.061,.05751,
- .05591,.05222,.04833,.04665,.04341,.03818,.0328,.02673,.02099,
- .01723,.005654,-.01383,-.105,-.1831,-.2123,-.2841,-.3456,-.5173,
- .06848,.06183,.06082,.05728,.05407,.05236,.04857,.04368,.03737,
- .03088,.02479,.02076,.00785,-.01309,-.1095,-.1895,-.2188,-.2898,
- -.3493,-.5073,.07725,.06828,.06577,.06361,.06,.0591,.05476,.0496,
- .04247,.03556,.02896,.02469,.01038,-.01211,-.1139,-.1954,-.2247,
- -.2944,-.3514,-.4949,.08308,.07651,.07395,.07184,.06706,.06618,
- .06171,.05564,.0486,.04115,.03368,.02885,.01347,-.01083,-.118,
- -.2005,-.2296,-.2975,-.3517,-.4801,.08673,.08355,.08361,.07958,
- .07582,.07311,.06882,.06262,.05585,.04731,.03906,.03336,.01719,
- -.009279,-.1217,-.2047,-.2334,-.2989,-.35,-.4629,.09333,.09061,
- .0912,.08707,.08512,.08086,.0767,.07113,.06356,.05366,.04492,
- .03862,.02135,-.007482,-.1247,-.2078,-.2358,-.2984,-.346,-.4432,
- .1045,.1001,.09799,.09636,.09317,.09079,.08611,.0806,.07117,
- .06057,.0511,.04488,.0258,-.005331,-.1269,-.2094,-.2365,-.2958,
- -.3395,-.4213,.1155,.1107,.1079,.1069,.1017,.102,.09663,.09015,
- .07944,.06847,.058,.05183,.03074,-.002583,-.128,-.2093,-.2353,
- -.2909,-.3305,-.3971,.1211,.1202,.1194,.1174,.1138,.1128,.1074,.1,
- .08932,.07729,.06612,.05905,.03639,9.063e-4,-.1279,-.2072,-.2319,
- -.2833,-.3187,-.3709,.1233,.1292,.1303,.128,.1278,.123,.1184,
- .1111,.1002,.08679,.07506,.06651,.04269,.005089,-.1263,-.2029,
- -.2261,-.2731,-.3041,-.3429,.1321,.1388,.1403,.1395,.1402,.1346,
- .1299,.1231,.1111,.09672,.08401,.07458,.04938,.009833,-.1231,
- -.1963,-.2178,-.26,-.2867,-.3134,.1445,.1497,.1506,.1521,.1511,
- .148,.1421,.1351,.1219,.1067,.09276,.08323,.05638,.01506,-.1179,
- -.1871,-.2068,-.2442,-.2666,-.2827,.1527,.1612,.1618,.1646,.1629,
- .1612,.1546,.1468,.1328,.1166,.1016,.09185,.06358,.0207,-.1108,
- -.1753,-.1931,-.2256,-.244,-.2512,.1531,.1729,.1742,.1762,.1762,
- .172,.1665,.1578,.1436,.1258,.1104,.09967,.07066,.02665,-.1018,
- -.161,-.1768,-.2046,-.2192,-.2193,.1504,.1808,.1867,.1857,.1899,
- .1799,.176,.1674,.1534,.1338,.1183,.1061,.077,.03257,-.09103,
- -.1444,-.1582,-.1814,-.1928,-.1877,.1585,.1857,.1955,.1931,.2005,
- .1877,.1827,.1744,.1602,.1394,.1237,.1108,.08167,.03777,-.07876,
- -.126,-.1377,-.1566,-.1651,-.1568,.1756,.1927,.1989,.1996,.2039,
- .1954,.1872,.1776,.1618,.1413,.1249,.1131,.08353,.04143,-.0655,
- -.1062,-.1159,-.131,-.1371,-.1271,.1892,.2037,.2008,.205,.2001,
- .1984,.1885,.1755,.1574,.1384,.1214,.1112,.08191,.04292,-.0519,
- -.0858,-.09356,-.1052,-.1094,-.09931,.1897,.2113,.2088,.2055,
- .1944,.1912,.1819,.1647,.1485,.1301,.1133,.1032,.07671,.04191,
- -.03872,-.06571,-.07167,-.08025,-.08302,-.07396,.1712,.2034,.2211,
- .1947,.188,.168,.1625,.1434,.1349,.1157,.1004,.08785,.06772,
- .03801,-.02677,-.04688,-.05119,-.05716,-.05887,-.05161,.1322,
- .1811,.2187,.1714,.1713,.1281,.1309,.1119,.1131,.09404,.0814,
- .06592,.05425,.03086,-.0167,-.03025,-.0331,-.03691,-.03788,
- -.03277,.08504,.1435,.1758,.1326,.133,.08054,.09003,.07395,.07987,
- .06447,.05582,.04075,.0365,.02082,-.008919,-.0167,-.01832,-.02041,
- -.0209,-.01789,.03936,.08311,.09568,.07461,.0731,.03755,.04509,
- .03589,.04026,.03161,.02757,.01806,.01764,.01004,-.003585,
- -.006904,-.007594,-.00846,-.008646,-.007351,.008367,.01985,.02146,
- .01751,.01672,.00796,.009761,.007629,.008715,.006695,.005893,
- .003567,.003711,.002105,-6.734e-4,-.001321,-.001456,-.001622,
- -.001656,-.001403,1.901e-20,-6.852e-20,-2.92e-22,6.497e-20,
- 4.95e-20,9.892e-20,1.818e-20,-9.169e-22,-3.362e-21,6.922e-20,
- 2.684e-20,2.626e-20,-3.55e-20,2.511e-20,-5.888e-21,4.135e-20,
- -1.341e-21,-1.894e-20,1.482e-20,-2.286e-21 };
- static doublereal uhr[1660] /* was [20][83] */ = { -.2975,-.2458,-.2347,
- -.2072,-.1879,-.1755,-.1638,-.1489,-.1415,-.1337,-.1284,-.1253,
- -.1225,-.12,-.128,-.142,-.1485,-.1709,-.1989,-.4039,-.2668,-.2241,
- -.2152,-.1914,-.1748,-.1643,-.1546,-.1418,-.1356,-.1291,-.1246,
- -.1221,-.12,-.1184,-.1274,-.1416,-.1482,-.1707,-.1987,-.4037,
- -.1549,-.1447,-.1427,-.1329,-.1262,-.1222,-.1193,-.1143,-.1128,
- -.111,-.1098,-.1093,-.1102,-.1117,-.1248,-.1399,-.1467,-.1695,
- -.1977,-.4031,-.01893,-.04556,-.04866,-.05698,-.06112,-.06482,
- -.07024,-.07374,-.0794,-.08412,-.08708,-.08905,-.09449,-.1007,
- -.1202,-.1368,-.144,-.1675,-.1961,-.402,.08761,.03753,.03153,
- .01086,-6.59e-5,-.009411,-.02173,-.03142,-.04385,-.05464,-.06128,
- -.06548,-.07561,-.0867,-.1139,-.1325,-.1402,-.1645,-.1937,-.4004,
- .149,.09297,.08391,.05998,.0454,.03384,.01788,.00398,-.01264,
- -.02736,-.03663,-.04236,-.05627,-.07146,-.1061,-.127,-.1353,
- -.1608,-.1907,-.3983,.169,.1209,.1109,.08847,.07199,.06082,.0448,
- .02856,.01078,-.005128,-.01602,-.02251,-.03862,-.05647,-.09731,
- -.1207,-.1296,-.1563,-.187,-.3957,.1621,.1259,.1176,.09828,.08237,
- .07282,.05866,.04264,.02581,.01067,-6.693e-4,-.007273,-.02391,
- -.04287,-.08809,-.1137,-.1232,-.1512,-.1828,-.3926,.1425,.1161,
- .1105,.09476,.08222,.074,.0622,.04865,.03391,.0202,.009705,
- .003367,-.01248,-.03132,-.07884,-.1062,-.1164,-.1457,-.1781,
- -.3891,.1202,.09999,.0965,.08454,.07657,.06903,.05955,.04937,
- .03702,.02497,.01606,.01027,-.004072,-.02191,-.06988,-.09867,
- -.1093,-.1397,-.173,-.3851,.1006,.08441,.08184,.07332,.06805,
- .06219,.05464,.04702,.03695,.02688,.01949,.01454,.001893,-.01438,
- -.06141,-.09113,-.1022,-.1336,-.1676,-.3806,.08386,.07229,.06955,
- .06387,.05898,.05552,.04967,.04327,.03532,.0274,.02112,.01709,
- .00611,-.008345,-.05355,-.08376,-.09513,-.1273,-.162,-.3757,.0708,
- .06217,.05952,.05578,.0514,.04945,.04499,.03953,.03338,.02725,
- .02193,.01863,.009222,-.003461,-.04633,-.07664,-.08821,-.121,
- -.1563,-.3703,.05982,.05317,.05163,.04867,.04579,.04406,.04066,
- .03665,.03176,.02679,.02246,.01961,.01167,5.524e-4,-.03973,
- -.06985,-.08152,-.1147,-.1504,-.3645,.05096,.04587,.04572,.04299,
- .04188,.03959,.03709,.03458,.03059,.02633,.0229,.02037,.01368,
- .003912,-.0337,-.06339,-.07508,-.1084,-.1444,-.3582,.04534,.04108,
- .04123,.0392,.03887,.03676,.03475,.03295,.0297,.02608,.02328,
- .02115,.01539,.006786,-.02821,-.05728,-.06891,-.1023,-.1385,
- -.3514,.04095,.03883,.03829,.03724,.03623,.03534,.0336,.03169,
- .02902,.02615,.02365,.02194,.0169,.009315,-.02319,-.05149,-.06302,
- -.09627,-.1324,-.3442,.03683,.03724,.03675,.03592,.03462,.03409,
- .03288,.03099,.02881,.02647,.02419,.02271,.01836,.01161,-.01858,
- -.04602,-.05738,-.09034,-.1264,-.3364,.03467,.03524,.03552,.03453,
- .03409,.03301,.03214,.03081,.02907,.02689,.02493,.02351,.0198,
- .01374,-.01432,-.04082,-.05197,-.08453,-.1204,-.3282,.03363,
- .03361,.03436,.03352,.03395,.03255,.03166,.03096,.02951,.02739,
- .02577,.02441,.02119,.01573,-.01035,-.03587,-.04678,-.0788,-.1143,
- -.3194,.03287,.03266,.03347,.03301,.03388,.0325,.0316,.03132,
- .02992,.02801,.02666,.02544,.02249,.01761,-.00661,-.03112,-.04175,
- -.07314,-.1081,-.31,.03248,.03251,.03299,.03301,.03375,.03281,
- .032,.03176,.03035,.02876,.02756,.02654,.0238,.0194,-.00305,
- -.02653,-.03685,-.06751,-.1019,-.3,.03227,.03315,.03325,.03357,
- .03376,.03336,.03272,.03214,.03096,.02966,.02852,.02766,.02518,
- .02116,3.643e-4,-.02207,-.03205,-.06189,-.09558,-.2894,.03202,
- .03373,.03383,.03408,.03412,.03393,.03344,.03268,.03179,.03063,
- .02956,.0288,.02665,.0229,.003671,-.01769,-.02731,-.05623,-.0891,
- -.278,.03221,.03388,.03404,.03426,.03465,.03451,.03408,.03362,
- .03274,.03161,.03068,.03001,.02817,.02468,.006915,-.01333,-.02258,
- -.05048,-.08243,-.2659,.03315,.03429,.03451,.03476,.03528,.0352,
- .03483,.03472,.0338,.03272,.03189,.03133,.02972,.02653,.01015,
- -.008958,-.01781,-.04461,-.07552,-.253,.03341,.03542,.03568,
- .03596,.03607,.03611,.03588,.03576,.03498,.03407,.03322,.03275,
- .03133,.0285,.01342,-.004508,-.01296,-.03855,-.06831,-.2391,
- .03416,.0368,.03694,.03742,.03709,.03734,.03718,.03689,.03632,
- .03558,.03473,.03432,.03304,.03059,.0168,5.975e-5,-.007946,
- -.03226,-.06073,-.2243,.03596,.03789,.03806,.03875,.03849,.03875,
- .03852,.03828,.03788,.03721,.03646,.03607,.03493,.0328,.02031,
- .00479,-.00273,-.02567,-.05272,-.2084,.0375,.03897,.03948,.04012,
- .04035,.04016,.03995,.03996,.03968,.03899,.0384,.03801,.03701,
- .03514,.02401,.009731,.002748,-.01872,-.04419,-.1914,.03891,
- .04054,.04132,.04176,.04245,.04187,.04174,.04184,.04165,.04098,
- .04053,.04016,.0393,.03764,.02791,.01494,.008549,-.01132,-.03507,
- -.1732,.04129,.04285,.04339,.04385,.04439,.04421,.044,.04396,
- .04372,.04322,.04283,.04258,.0418,.04035,.03206,.02049,.01474,
- -.003382,-.02525,-.1536,.04417,.04571,.0458,.04635,.04649,.04679,
- .04662,.04639,.04606,.04572,.04539,.04526,.04458,.04334,.03651,
- .02646,.02139,.005188,-.01464,-.1325,.04684,.04849,.04869,.04896,
- .04928,.04939,.04938,.04919,.04885,.0485,.04833,.04821,.04769,
- .04666,.04133,.03294,.02859,.0145,-.003121,-.1098,.05013,.05137,
- .05183,.05179,.05249,.05231,.05239,.05242,.05201,.05159,.05162,
- .0515,.05114,.05035,.0466,.04003,.03646,.02466,.009448,-.08535,
- .05405,.05498,.05537,.05527,.05587,.05573,.05583,.05607,.0555,
- .05508,.05523,.05519,.05494,.05445,.05241,.04783,.0451,.03581,
- .02321,-.05904,.05841,.05913,.05952,.05936,.05974,.05966,.05974,
- .06009,.05942,.05907,.05925,.05931,.05917,.05902,.05887,.05646,
- .05467,.04809,.03833,-.03069,.06347,.06371,.064,.06393,.06419,
- .06418,.06417,.06449,.06385,.0636,.06377,.06391,.06391,.06414,
- .06608,.06605,.0653,.06166,.05501,-1.341e-4,.06904,.06905,.06899,
- .06915,.06921,.06938,.06925,.06936,.06886,.06873,.06888,.06905,
- .06922,.06987,.07416,.07674,.07716,.07669,.07343,.03282,.07469,
- .07498,.07493,.07503,.07508,.07512,.0749,.07485,.07458,.07452,
- .07468,.07481,.07522,.07631,.08324,.0887,.09044,.09342,.09383,
- .06835,.08139,.08126,.08148,.0814,.0818,.08147,.08112,.08116,
- .08105,.08102,.08123,.08135,.08197,.08355,.09346,.1021,.1053,
- .1121,.1165,.1067,.08556,.08467,.08488,.08483,.08536,.08505,
- .08453,.08465,.08454,.08456,.0848,.08494,.08565,.08752,.09906,
- .1095,.1135,.1222,.1287,.1269,.09018,.08842,.08848,.08854,.08906,
- .08894,.08822,.08836,.08823,.08833,.08857,.08878,.08956,.09173,
- .105,.1173,.1221,.1329,.1416,.148,.09939,.09722,.09678,.09708,
- .09724,.09748,.09652,.09644,.09634,.09665,.09687,.09717,.09817,
- .101,.1181,.1344,.1409,.1563,.1696,.1925,.1083,.1073,.1068,.1069,
- .1069,.1067,.1059,.1056,.1057,.1061,.1063,.1067,.108,.1116,.1329,
- .1537,.1622,.1826,.2008,.2405,.1199,.1181,.118,.1177,.1176,.1172,
- .1164,.1162,.1163,.1168,.1171,.1175,.1192,.1236,.1497,.1757,.1863,
- .2122,.2355,.2923,.134,.1307,.1305,.13,.1296,.1294,.1285,.1284,
- .1284,.1288,.1294,.1299,.132,.1373,.1689,.2007,.2136,.2455,.2742,
- .348,.1482,.1455,.145,.1444,.1436,.1432,.1424,.1423,.1422,.1426,
- .1435,.1441,.1467,.153,.191,.2293,.2446,.2831,.3175,.408,.1645,
- .162,.1618,.1605,.1597,.1588,.1582,.1583,.1581,.1584,.1596,.1605,
- .1637,.171,.2164,.2619,.2799,.3255,.3658,.4725,.1852,.1811,.1805,
- .1788,.1778,.1774,.1766,.1767,.1761,.1766,.1781,.1793,.1831,.1918,
- .2457,.2992,.3202,.3733,.4199,.5419,.2084,.2039,.2023,.2005,.1987,
- .1988,.198,.1976,.1969,.1977,.1994,.201,.2057,.2159,.2796,.342,
- .3663,.4273,.4803,.6164,.2329,.2296,.2281,.2254,.2236,.2231,.2224,
- .2216,.2213,.2223,.2243,.2261,.2319,.2439,.319,.3911,.419,.4884,
- .5479,.6963,.2629,.2581,.2569,.2536,.2526,.2511,.2503,.2497,.2496,
- .2507,.2533,.2554,.2624,.2767,.3648,.4476,.4793,.5574,.6235,.782,
- .2995,.2912,.2894,.2864,.2855,.2842,.2827,.2824,.2822,.2838,.2871,
- .2896,.298,.3151,.4182,.5127,.5485,.6355,.7079,.8737,.3403,.3305,
- .3283,.3254,.3236,.3227,.3207,.3202,.3204,.3228,.3266,.3297,.3396,
- .3604,.4804,.5877,.6277,.7238,.8023,.9717,.3861,.3762,.3744,.3709,
- .3687,.367,.3649,.3643,.3655,.3686,.3731,.3767,.3887,.414,.553,
- .6742,.7186,.8237,.9077,1.076,.4398,.4288,.4271,.4233,.421,.4185,
- .4166,.4165,.4184,.4225,.4279,.4321,.4469,.4774,.6381,.774,.8229,
- .9366,1.025,1.188,.5032,.4909,.4887,.4847,.4812,.4795,.4776,.4782,
- .4808,.4862,.4927,.4979,.5162,.5526,.7378,.8892,.9424,1.064,1.156,
- 1.306,.5757,.5646,.5619,.5569,.5526,.5512,.5499,.5508,.555,.5619,
- .5699,.5763,.5991,.6423,.8549,1.022,1.08,1.208,1.302,1.432,.6589,
- .6497,.6477,.641,.638,.6354,.6353,.6374,.6435,.6521,.6625,.6704,
- .6986,.7496,.9926,1.175,1.237,1.371,1.465,1.565,.758,.748,.7467,
- .7396,.7385,.7359,.7367,.7413,.7489,.7601,.7738,.7839,.8181,.8784,
- 1.155,1.352,1.417,1.553,1.645,1.705,.8782,.8653,.8629,.8577,.8571,
- .8567,.8588,.8654,.8752,.8905,.9083,.9213,.9622,1.034,1.346,1.557,
- 1.623,1.759,1.844,1.852,1.017,1.006,1.002,.9997,.9993,1.002,1.006,
- 1.014,1.028,1.049,1.072,1.088,1.137,1.222,1.571,1.792,1.859,1.99,
- 2.065,2.007,1.175,1.17,1.169,1.169,1.171,1.175,1.182,1.194,1.214,
- 1.243,1.271,1.291,1.349,1.45,1.835,2.063,2.128,2.248,2.307,2.169,
- 1.362,1.364,1.368,1.371,1.379,1.383,1.395,1.414,1.443,1.479,1.515,
- 1.538,1.609,1.727,2.147,2.373,2.435,2.537,2.574,2.337,1.588,1.597,
- 1.604,1.614,1.63,1.637,1.655,1.684,1.723,1.77,1.813,1.842,1.928,
- 2.064,2.512,2.729,2.782,2.857,2.865,2.511,1.853,1.882,1.89,1.912,
- 1.933,1.949,1.976,2.016,2.068,2.127,2.181,2.217,2.319,2.473,2.94,
- 3.134,3.175,3.212,3.181,2.69,2.167,2.227,2.24,2.276,2.306,2.332,
- 2.373,2.428,2.496,2.571,2.636,2.68,2.799,2.971,3.439,3.594,3.616,
- 3.603,3.524,2.873,2.547,2.643,2.667,2.722,2.771,2.808,2.866,2.94,
- 3.027,3.121,3.202,3.253,3.39,3.576,4.02,4.114,4.11,4.031,3.892,
- 3.06,3.006,3.153,3.191,3.274,3.351,3.403,3.483,3.581,3.692,3.808,
- 3.904,3.963,4.114,4.31,4.691,4.697,4.659,4.496,4.285,3.248,3.574,
- 3.794,3.844,3.972,4.081,4.155,4.26,4.386,4.524,4.664,4.774,4.84,
- 5.002,5.196,5.46,5.346,5.264,4.996,4.7,3.436,4.289,4.609,4.679,
- 4.864,5.008,5.107,5.243,5.401,5.57,5.733,5.852,5.923,6.085,6.258,
- 6.334,6.06,5.923,5.53,5.135,3.622,5.202,5.661,5.76,6.012,6.196,
- 6.326,6.496,6.688,6.884,7.062,7.183,7.253,7.399,7.522,7.316,6.837,
- 6.632,6.091,5.583,3.805,6.408,7.047,7.181,7.502,7.737,7.895,8.1,
- 8.32,8.53,8.707,8.817,8.876,8.978,9.008,8.403,7.669,7.384,6.672,
- 6.04,3.98,8.085,8.911,9.079,9.47,9.756,9.932,10.16,10.39,10.58,
- 10.73,10.81,10.84,10.85,10.73,9.584,8.544,8.166,7.262,6.495,4.147,
- 10.5,11.48,11.67,12.1,12.41,12.58,12.79,12.99,13.12,13.19,13.19,
- 13.17,13.03,12.68,10.84,9.442,8.96,7.849,6.94,4.303,14.08,15.08,
- 15.25,15.65,15.89,16.03,16.15,16.22,16.21,16.12,15.98,15.88,15.51,
- 14.84,12.13,10.34,9.743,8.415,7.362,4.443,19.47,20.15,20.22,20.4,
- 20.41,20.43,20.34,20.16,19.88,19.52,19.16,18.93,18.23,17.14,13.4,
- 11.19,10.49,8.941,7.749,4.567,27.55,27.16,27.,26.59,26.14,25.87,
- 25.41,24.78,24.05,23.28,22.62,22.2,21.06,19.45,14.6,11.97,11.15,
- 9.406,8.085,4.671,39.18,36.33,35.73,34.17,32.96,32.18,31.11,29.83,
- 28.47,27.14,26.09,25.43,23.79,21.6,15.63,12.62,11.71,9.788,8.359,
- 4.753,54.06,46.92,45.58,42.34,40.09,38.6,36.73,34.64,32.56,30.59,
- 29.13,28.24,26.09,23.35,16.42,13.11,12.13,10.07,8.557,4.811,68.07,
- 55.89,53.77,48.85,45.6,43.45,40.85,38.07,35.4,32.92,31.15,30.07,
- 27.56,24.44,16.89,13.4,12.37,10.23,8.671,4.844,72.5,58.54,56.16,
- 50.7,47.14,44.78,41.96,38.99,36.14,33.52,31.67,30.54,27.93,24.71,
- 17.,13.46,12.42,10.27,8.698,4.852 };
- /* System generated locals */
- integer i__1;
- /* Local variables */
- integer i__, j;
- /* - 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) >*/
- /*< integer i,j >*/
- /* Aerosol stratospheric model */
- /* Parameters (STM-params.dat) are calculated with PINATUBO.f */
- /* for one-year-old aerosol. PINATUBO.f is written based on */
- /* Russel et al., 'Global to microscale evolution of the Pinatubo */
- /* volcanic aerosol derived from diverse measurements and analyses', */
- /* Journal of Geophysical Research, 101(D13), pp. 18745-18763, 1996. */
- /* Extinction coefficients are calculated in km-1 */
- /* *************** MIE (asy) ****************** */
- /*< >*/
- /* *************** MIE (ext&sca) ****************** */
- /*< >*/
- /* ************** MIE (phase fun ph) *************** */
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /* ************** MIE (phase fun qh) *************** */
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /* ************** MIE (phase fun uh) *************** */
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< >*/
- /*< 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 >*/
- } /* stm_ */
- #ifdef __cplusplus
- }
- #endif