SkyMapMod 0.1.16__py3-none-any.whl → 0.1.18__py3-none-any.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (10) hide show
  1. SkyMapMod/__init__.py +3 -1
  2. SkyMapMod/airglow.py +0 -1
  3. SkyMapMod/albedo_of_planets.py +23 -0
  4. SkyMapMod/planets.py +131 -0
  5. SkyMapMod/star_catalogues.py +2 -0
  6. SkyMapMod/zodiac.py +12 -55
  7. {skymapmod-0.1.16.dist-info → skymapmod-0.1.18.dist-info}/METADATA +1 -1
  8. {skymapmod-0.1.16.dist-info → skymapmod-0.1.18.dist-info}/RECORD +10 -8
  9. {skymapmod-0.1.16.dist-info → skymapmod-0.1.18.dist-info}/WHEEL +0 -0
  10. {skymapmod-0.1.16.dist-info → skymapmod-0.1.18.dist-info}/licenses/LICENSE +0 -0
SkyMapMod/__init__.py CHANGED
@@ -6,4 +6,6 @@ from .modtran_default_kp_transparency import *
6
6
  from .solar_spectrum import *
7
7
  from .band_V_data import *
8
8
  from .solar_radio_flux import *
9
- from .star_catalogues import *
9
+ from .star_catalogues import *
10
+ from .albedo_of_planets import *
11
+ from .planets import *
SkyMapMod/airglow.py CHANGED
@@ -1,6 +1,5 @@
1
1
  #для работы с собственным свечением атмосферы. возможно, сюда можно досыпать поправку на радиопоток
2
2
  import numpy as np
3
- import math
4
3
  import matplotlib.pyplot as plt
5
4
 
6
5
  from .solar_radio_flux import fluxdate
SkyMapMod/albedo_of_planets.py ADDED
@@ -0,0 +1,23 @@
1
+ import numpy as np
2
+
3
+ #все длины волн в нм
4
+ #все коэффициенты отражения в долях единицы
5
+
6
+ venus_alb_wl = np.array([300.0, 300.5, 301.0, 301.5, 302.0, 302.5, 303.0, 303.5, 304.0, 304.5, 305.0, 305.5, 306.0, 306.5, 307.0, 307.5, 308.0, 308.5, 309.0, 309.5, 310.0, 310.5, 311.0, 311.5, 312.0, 312.5, 313.0, 313.5, 314.0, 314.5, 315.0, 315.5, 316.0, 316.5, 317.0, 317.5, 318.0, 318.5, 319.0, 319.5, 320.0, 320.5, 321.0, 321.5, 322.0, 322.5, 323.0, 323.5, 324.0, 324.5, 325.0, 325.5, 326.0, 326.5, 327.0, 327.5, 328.0, 328.5, 329.0, 329.5, 330.0, 330.5, 331.0, 331.5, 332.0, 332.5, 333.0, 333.5, 334.0, 334.5, 335.0, 335.5, 336.0, 336.5, 337.0, 337.5, 338.0, 338.5, 339.0, 339.5, 340.0, 340.5, 341.0, 341.5, 342.0, 342.5, 343.0, 343.5, 344.0, 344.5, 345.0, 345.5, 346.0, 346.5, 347.0, 347.5, 348.0, 348.5, 349.0, 349.5, 350.0, 350.5, 351.0, 351.5, 352.0, 352.5, 353.0, 353.5, 354.0, 354.5, 355.0, 355.5, 356.0, 356.5, 357.0, 357.5, 358.0, 358.5, 359.0, 359.5, 360.0, 360.5, 361.0, 361.5, 362.0, 362.5, 363.0, 363.5, 364.0, 364.5, 365.0, 365.5, 366.0, 366.5, 367.0, 367.5, 368.0, 368.5, 369.0, 369.5, 370.0, 370.5, 371.0, 371.5, 372.0, 372.5, 373.0, 373.5, 374.0, 374.5, 375.0, 375.5, 376.0, 376.5, 377.0, 377.5, 378.0, 378.5, 379.0, 379.5, 380.0, 381.0, 382.0, 383.5, 385.0, 385.5, 386.5, 388.0, 388.5, 390.0, 390.0, 391.0, 392.0, 394.0, 395.0, 397.0, 398.0, 399.5, 400.5, 402.0, 403.0, 404.0, 406.0, 411.5, 417.5, 419.5, 420.0, 421.5, 424.5, 427.0, 428.0, 430.0, 431.5, 432.0, 434.0, 436.5, 437.5, 439.0, 439.0, 441.0, 442.0, 443.5, 445.0, 445.0, 445.5, 447.0, 448.5, 449.5, 451.0, 454.0, 455.0, 458.0, 460.5, 462.5, 464.0, 466.5, 470.0, 473.5, 474.5, 477.0, 478.0, 479.5, 481.5, 483.0, 484.0, 486.0, 487.5, 489.5, 491.0, 492.0, 494.0, 494.5, 496.0, 497.5, 498.5, 500.0, 502.5, 503.5, 505.5, 506.5, 508.0, 511.0, 512.0, 513.0, 514.0, 516.5, 518.0, 519.5, 522.5, 523.5, 524.0, 526.5, 529.0, 535.0, 537.0, 541.5, 542.5, 544.5, 545.5, 547.5, 550.0, 552.0, 554.0, 557.5, 561.0, 562.5, 567.0, 568.0, 571.0, 573.0, 574.5, 576.5, 580.5, 585.0, 589.0, 592.0, 593.5, 597.5, 599.5, 604.0, 604.0, 605.0, 607.0, 610.0, 612.0, 613.5, 616.5, 619.5, 623.5, 627.0, 630.0, 633.5, 635.5, 638.0, 639.5, 640.0, 641.0, 645.5, 647.0, 648.5, 650.5, 651.0, 653.0, 655.0, 656.0, 657.0, 658.5, 660.5, 662.5, 665.0, 668.0, 670.0, 672.0, 675.0, 676.5, 678.0, 679.5, 680.5, 681.5, 683.0, 686.0, 689.0, 690.5, 692.0, 693.0, 696.0, 698.0, 702.5, 703.5, 707.0, 709.0, 711.0, 712.5, 714.5, 715.0, 716.5, 718.0, 720.0, 721.5, 723.5, 726.5, 730.0, 731.0, 733.0, 734.0, 736.0, 737.5, 739.0, 740.0, 741.0, 742.0, 743.0, 744.0, 745.0, 746.0, 747.0, 748.0, 749.0, 750.0, 751.0, 752.0, 753.0, 754.0, 755.0, 756.0, 757.0, 758.0, 759.0, 760.0, 761.0, 762.0, 763.0, 764.0, 765.0, 766.0, 767.0, 768.0, 769.0, 770.0, 771.0, 772.0, 773.0, 774.0, 775.0, 776.0, 777.0, 778.0, 779.0, 780.0, 781.0, 782.0, 783.0, 784.0, 785.0, 786.0, 787.0, 788.0, 789.0, 790.0, 791.0, 792.0, 793.0, 794.0, 795.0, 796.0, 797.0, 798.0, 799.0, 800.0, 801.0, 802.0, 803.0, 804.0, 805.0, 806.0, 807.0, 808.0, 809.0, 810.0, 811.0, 812.0, 813.0, 814.0, 815.0, 816.0, 817.0, 818.0, 819.0, 820.0, 821.0, 822.0, 823.0, 824.0, 825.0, 826.0, 827.0, 828.0, 829.0, 830.0, 831.0, 832.0, 833.0, 834.0, 835.0, 836.0, 837.0, 838.0, 839.0, 840.0, 841.0, 842.0, 843.0, 844.0, 845.0, 846.0, 847.0, 848.0, 849.0, 850.0, 851.0, 852.0, 853.0, 854.0, 855.0, 856.0, 857.0, 858.0, 859.0, 860.0, 861.0, 862.0, 863.0, 864.0, 865.0, 866.0, 867.0, 868.0, 869.0, 870.0, 871.0, 872.0, 873.0, 874.0, 875.0, 876.0, 877.0, 878.0, 879.0, 880.0, 881.0, 882.0, 883.0, 884.0, 885.0, 886.0, 887.0, 888.0, 889.0, 890.0, 891.0, 892.0, 893.0, 894.0, 895.0, 896.0, 897.0, 898.0, 899.0, 900.0, 901.0, 902.0, 903.0, 904.0, 905.0, 906.0, 907.0, 908.0, 909.0, 910.0, 911.0, 912.0, 913.0, 914.0, 915.0, 916.0, 917.0, 918.0, 919.0, 920.0, 921.0, 922.0, 923.0, 924.0, 925.0, 926.0, 927.0, 928.0, 929.0, 930.0, 931.0, 932.0, 933.0, 934.0, 935.0, 936.0, 937.0, 938.0, 939.0, 940.0, 941.0, 942.0, 943.0, 944.0, 945.0, 946.0, 947.0, 948.0, 949.0, 950.0, 951.0, 952.0, 953.0, 954.0, 955.0, 956.0, 957.0, 958.0, 959.0, 960.0, 961.0, 962.0, 963.0, 964.0, 965.0, 966.0, 967.0, 968.0, 969.0, 970.0, 971.0, 972.0, 973.0, 974.0, 975.0, 976.0, 977.0, 978.0, 979.0, 980.0, 981.0, 982.0, 983.0, 984.0, 985.0, 986.0, 987.0, 988.0, 989.0, 990.0, 991.0, 992.0, 993.0, 994.0, 995.0, 996.0, 997.0, 998.0, 999.0])
7
+
8
+ venus_alb_rf = np.array([0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.09263657957244642, 0.11402518106174298, 0.1520372178424143, 0.14729755329845673, 0.0808001690395348, 0.104560424347522, 0.12119478891916702, 0.09982804599041128, 0.15921411188668533, 0.16397563499118373, 0.19010390102443786, 0.1853605933870569, 0.29701047753668597, 0.2732720807892396, 0.2970323360972268, 0.29705055156434423, 0.261428384069481, 0.2566887195255234, 0.29470075630619463, 0.285210497938009, 0.32559783163079425, 0.33748524547163483, 0.36600338079069694, 0.4064271454177171, 0.43022747475336254, 0.44686912551185454, 0.4349999271381315, 0.43976145024263, 0.47541276248488107, 0.477809917957536, 0.5134466578260933, 0.5467153869693834, 0.5728545822829081, 0.551484196260729, 0.5443728778980808, 0.5657687655742244, 0.5847784270579836, 0.5942941870801335, 0.6156718592891596, 0.6275629162234235, 0.5990666394649025, 0.6157046471299709, 0.6370968917126911, 0.6370968917126911, 0.6252240502455444, 0.6466126517348412, 0.6466235810151115, 0.6276321349984699, 0.6490207364877665, 0.6799214549057897, 0.663305305801262, 0.682329539658715, 0.6823477551258323, 0.6681141891202658, 0.6895027906095623, 0.7109023213791295, 0.7228043075936639, 0.7204581554289378, 0.7323419261763549, 0.7323601416434724, 0.7204945863631727, 0.7371325940282414, 0.7418977602261632, 0.7727875493639158, 0.7965514477653264, 0.8060672077874764, 0.7775745741223787, 0.7395843959022486, 0.7514754528365126, 0.7704851143202716, 0.7633737959576236, 0.7823798143479591, 0.7728931990731971, 0.7871559098261516, 0.8132914620462528, 0.8323047666234353, 0.803823062238608, 0.8157068329860251, 0.8038449207991489, 0.7967263162496538, 0.8109926700960319, 0.7943874502717747, 0.8110218148434198, 0.8086574472115763, 0.8371682963437914, 0.849062996371479, 0.8371974410911793, 0.799207262871049, 0.8277326843769581, 0.8229893767395771, 0.8491212858662548, 0.8705208166358218, 0.827783687684887, 0.8610852046689885, 0.8777304985209041, 0.8611362079769174, 0.8468953557845038, 0.8469099281581978, 0.8659195896419568, 0.8825648834938724, 0.8707066144004196, 0.8588447022135436, 0.873114699153345, 0.8873919822799936, 0.8707940486425834, 0.8613037902743979, 0.8684661119449747, 0.8851004765166198, 0.8851259781705842, 0.8542616906867959, 0.8804008860003206, 0.8875449922037801, 0.8566952770936858, 0.8828599740611749, 0.8923939495504423, 0.854411057517159, 0.871052708275651, 0.8734571499351529, 0.9091048190839806, 0.911512903836906, 0.911512903836906, 0.8972684085510689, 0.8925360301939583, 0.9186861547877534, 0.9139537764306428, 0.8925870335018872, 0.9187371580956822, 0.8736320184194803, 0.906918963029888, 0.9306974338049925, 0.8998440756014747, 0.9378743278492634, 0.9307666525800388, 0.9070318989260161, 0.9070428282062866, 0.8951736298325634, 0.9236844789647786, 0.9237209098990135, 0.8881023854975737, 0.9047403931626423, 0.8928784809757662, 0.9071339055418737, 0.9261544963059033, 0.9760503038339915, 0.9499293239875843, 0.9499402532678548, 0.9071958381300731, 0.911961004327995, 0.9594815149439693, 0.9595033735045102, 0.928646372207569, 0.9286645876746864, 0.9453062384331784, 0.9405811462629148, 0.902587324949361, 0.8835958789327194, 0.8313502761464815, 0.831361205426752, 0.8242426008772569, 0.8670088745755796, 0.9786733310989028, 0.9288139545050493, 0.8860695393672675, 0.9003358932136456, 0.9264714454337467, 0.9478709762033137, 0.9526397854946592, 0.9407960887749005, 0.9099281581976888, 0.9170795505879953, 0.9170977660551127, 0.8149745712079041, 0.8411137665214287, 0.8435036358072366, 0.8316344374335136, 0.8744007111318363, 0.8530339682030806, 0.8245449776314064, 0.8198089561808724, 0.8720800606210746, 0.8816067499234951, 0.8816322515774595, 0.8626371624673943, 0.8579047841102837, 0.8745391486819287, 0.8888055025283068, 0.914948340935255, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254, 0.9149592702155254])
9
+
10
+
11
+ mars_alb_wl = np.array([300.0, 300.5, 301.0, 301.5, 302.0, 302.5, 303.0, 303.5, 304.0, 304.5, 305.0, 305.5, 306.0, 306.5, 307.0, 307.5, 308.0, 308.5, 309.0, 309.5, 310.0, 310.5, 311.0, 311.5, 312.0, 312.5, 313.0, 313.5, 314.0, 314.5, 315.0, 315.5, 316.0, 316.5, 317.0, 317.5, 318.0, 318.5, 319.0, 319.5, 320.0, 320.5, 321.0, 321.5, 322.0, 322.5, 323.0, 323.5, 324.0, 324.5, 325.0, 325.5, 326.0, 326.5, 327.0, 327.5, 328.0, 328.5, 329.0, 329.5, 330.0, 330.5, 331.0, 331.5, 332.0, 332.5, 333.0, 333.5, 334.0, 334.5, 335.0, 335.5, 336.0, 336.5, 337.0, 337.5, 338.0, 338.5, 339.0, 339.5, 340.0, 340.5, 341.0, 341.5, 342.0, 342.5, 343.0, 343.5, 344.0, 344.5, 345.0, 345.5, 346.0, 346.5, 347.0, 347.5, 348.0, 348.5, 349.0, 349.5, 350.0, 350.5, 351.0, 351.5, 352.0, 352.5, 353.0, 353.5, 354.0, 354.5, 355.0, 355.5, 356.0, 356.5, 357.0, 357.5, 358.0, 358.5, 359.0, 359.5, 360.0, 360.5, 361.0, 361.5, 362.0, 362.5, 363.0, 363.5, 364.0, 364.5, 365.0, 365.5, 366.0, 366.5, 367.0, 367.5, 368.0, 368.5, 369.0, 369.5, 370.0, 370.5, 371.0, 371.5, 372.0, 372.5, 373.0, 373.5, 374.0, 374.5, 375.0, 375.5, 376.0, 376.5, 377.0, 377.5, 378.0, 378.5, 379.0, 379.5, 380.0, 380.5, 381.0, 381.5, 382.0, 382.5, 383.0, 383.5, 384.0, 384.5, 385.0, 385.5, 386.0, 386.5, 387.0, 387.5, 388.0, 388.5, 389.0, 389.5, 390.0, 390.5, 391.0, 391.5, 392.0, 392.5, 393.0, 393.5, 394.0, 394.5, 395.0, 395.5, 396.0, 401.0, 406.5, 411.5, 417.0, 422.0, 427.5, 433.0, 438.5, 444.0, 449.5, 455.0, 460.5, 466.0, 471.5, 477.0, 483.0, 488.5, 494.5, 500.0, 506.0, 511.5, 517.5, 523.5, 529.5, 535.0, 541.0, 547.0, 553.0, 559.0, 565.5, 571.5, 577.5, 583.5, 590.0, 596.0, 602.5, 608.5, 615.0, 621.0, 627.5, 640.5, 643.0, 647.0, 653.5, 660.0, 666.5, 673.0, 679.5, 682.0, 693.0, 699.5, 704.0, 716.5, 729.5, 742.0, 754.5, 767.0, 779.5, 792.0, 805.0, 803.0, 829.0, 856.5, 883.5, 910.0, 937.5, 965.5, 992.0, 993.0, 994.0, 995.0, 996.0, 997.0, 998.0, 999.0])
12
+
13
+ mars_alb_rf = np.array([0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0370416, 0.0367995, 0.0390063, 0.0439349, 0.0416433, 0.0429654, 0.045467, 0.0498745, 0.0542356, 0.0566804, 0.058356, 0.0603925, 0.06358, 0.0683027, 0.0700244, 0.0720945, 0.0778272, 0.0799044, 0.084825, 0.0878545, 0.0924617, 0.0965058, 0.102426, 0.1064, 0.111605, 0.117631, 0.126016, 0.132353, 0.137337, 0.149238, 0.158461, 0.170818, 0.178154, 0.191217, 0.202846, 0.208467, 0.224776, 0.23087, 0.24077, 0.24184, 0.246895, 0.252237, 0.251927, 0.256119, 0.260163, 0.264544, 0.264351, 0.270838, 0.272577, 0.277526, 0.274634, 0.28391, 0.291076, 0.298498, 0.306176, 0.306296, 0.310309, 0.305621, 0.306522, 0.304822, 0.306767, 0.306008, 0.30278, 0.30234, 0.299554, 0.302506, 0.303779, 0.30435, 0.305426, 0.305426, 0.305426, 0.305426, 0.305426, 0.305426, 0.305426, 0.305426])
14
+
15
+
16
+ jupiter_alb_wl = np.array([300.4, 300.8, 301.20000000000005, 301.59999999999997, 302.0, 302.4, 302.8, 303.20000000000005, 303.59999999999997, 304.0, 304.4, 304.8, 305.20000000000005, 305.59999999999997, 306.0, 306.4, 306.8, 307.2, 307.59999999999997, 308.0, 308.40000000000003, 308.8, 309.2, 309.59999999999997, 310.0, 310.40000000000003, 310.8, 311.2, 311.59999999999997, 312.0, 312.40000000000003, 312.8, 313.2, 313.59999999999997, 314.0, 314.40000000000003, 314.8, 315.2, 315.59999999999997, 316.0, 316.40000000000003, 316.8, 317.2, 317.59999999999997, 318.0, 318.40000000000003, 318.79999999999995, 319.2, 319.6, 320.0, 320.40000000000003, 320.79999999999995, 321.2, 321.6, 322.0, 322.40000000000003, 322.79999999999995, 323.2, 323.6, 324.0, 324.40000000000003, 324.79999999999995, 325.2, 325.6, 326.0, 326.40000000000003, 326.79999999999995, 327.2, 327.6, 328.0, 328.40000000000003, 328.79999999999995, 329.2, 329.6, 330.0, 330.40000000000003, 330.79999999999995, 331.2, 331.6, 332.0, 332.4, 332.8, 333.2, 333.6, 334.0, 334.4, 334.8, 335.2, 335.6, 336.0, 336.4, 336.8, 337.2, 337.6, 338.0, 338.4, 338.8, 339.2, 339.6, 340.0, 340.4, 340.8, 341.2, 341.6, 342.0, 342.4, 342.8, 343.2, 343.6, 344.0, 344.4, 344.8, 345.2, 345.6, 346.0, 346.4, 346.8, 347.2, 347.6, 348.0, 348.4, 348.8, 349.2, 349.6, 350.0, 350.4, 350.8, 351.2, 351.6, 352.0, 352.4, 352.8, 353.2, 353.6, 354.0, 354.4, 354.8, 355.2, 355.6, 356.0, 356.4, 356.8, 357.20000000000005, 357.59999999999997, 358.0, 358.4, 358.8, 359.20000000000005, 359.59999999999997, 360.0, 360.4, 360.8, 361.20000000000005, 361.59999999999997, 362.0, 362.4, 362.8, 363.20000000000005, 363.59999999999997, 364.0, 364.4, 364.8, 365.20000000000005, 365.59999999999997, 366.0, 366.4, 366.8, 367.20000000000005, 367.59999999999997, 368.0, 368.4, 368.8, 369.2, 369.59999999999997, 370.0, 370.40000000000003, 370.8, 371.2, 371.59999999999997, 372.0, 372.40000000000003, 372.8, 373.2, 373.59999999999997, 374.0, 374.40000000000003, 374.8, 375.2, 375.59999999999997, 376.0, 376.40000000000003, 376.8, 377.2, 377.59999999999997, 378.0, 378.40000000000003, 378.8, 379.2, 379.59999999999997, 380.0, 380.40000000000003, 380.8, 381.2, 381.6, 382.0, 382.40000000000003, 382.79999999999995, 383.2, 383.6, 384.0, 384.40000000000003, 384.79999999999995, 385.2, 385.6, 386.0, 386.40000000000003, 386.79999999999995, 387.2, 387.6, 388.0, 388.40000000000003, 388.79999999999995, 389.2, 389.6, 390.0, 390.40000000000003, 390.79999999999995, 391.2, 391.6, 392.0, 392.40000000000003, 392.79999999999995, 393.2, 393.6, 394.0, 394.4, 394.8, 395.2, 395.6, 396.0, 396.4, 396.8, 397.2, 397.6, 398.0, 398.4, 398.8, 399.2, 399.6, 400.0, 400.4, 400.8, 401.2, 401.6, 402.0, 402.4, 402.8, 403.2, 403.6, 404.0, 404.4, 404.8, 405.2, 405.6, 406.0, 406.4, 406.8, 407.2, 407.6, 408.0, 408.4, 408.8, 409.2, 409.6, 410.0, 410.4, 410.8, 411.2, 411.6, 412.0, 412.4, 412.8, 413.2, 413.6, 414.0, 414.4, 414.8, 415.2, 415.6, 416.0, 416.4, 416.8, 417.2, 417.6, 418.0, 418.4, 418.8, 419.20000000000005, 419.59999999999997, 420.0, 420.4, 420.8, 421.20000000000005, 421.59999999999997, 422.0, 422.4, 422.8, 423.20000000000005, 423.59999999999997, 424.0, 424.4, 424.8, 425.20000000000005, 425.59999999999997, 426.0, 426.4, 426.8, 427.20000000000005, 427.59999999999997, 428.0, 428.4, 428.8, 429.20000000000005, 429.59999999999997, 430.0, 430.4, 430.8, 431.20000000000005, 431.59999999999997, 432.0, 432.40000000000003, 432.8, 433.2, 433.59999999999997, 434.0, 434.40000000000003, 434.8, 435.2, 435.59999999999997, 436.0, 436.40000000000003, 436.8, 437.2, 437.59999999999997, 438.0, 438.40000000000003, 438.8, 439.2, 439.59999999999997, 440.0, 440.40000000000003, 440.8, 441.2, 441.59999999999997, 442.0, 442.40000000000003, 442.8, 443.2, 443.6, 444.0, 444.40000000000003, 444.79999999999995, 445.2, 445.6, 446.0, 446.40000000000003, 446.79999999999995, 447.2, 447.6, 448.0, 448.40000000000003, 448.79999999999995, 449.2, 449.6, 450.0, 450.40000000000003, 450.79999999999995, 451.2, 451.6, 452.0, 452.40000000000003, 452.79999999999995, 453.2, 453.6, 454.0, 454.40000000000003, 454.79999999999995, 455.2, 455.6, 456.0, 456.4, 456.8, 457.2, 457.6, 458.0, 458.4, 458.8, 459.2, 459.6, 460.0, 460.4, 460.8, 461.2, 461.6, 462.0, 462.4, 462.8, 463.2, 463.6, 464.0, 464.4, 464.8, 465.2, 465.6, 466.0, 466.4, 466.8, 467.2, 467.6, 468.0, 468.4, 468.8, 469.2, 469.6, 470.0, 470.4, 470.8, 471.2, 471.6, 472.0, 472.4, 472.8, 473.2, 473.6, 474.0, 474.4, 474.8, 475.2, 475.6, 476.0, 476.4, 476.8, 477.2, 477.6, 478.0, 478.4, 478.8, 479.2, 479.6, 480.0, 480.4, 480.8, 481.2, 481.59999999999997, 482.0, 482.4, 482.8, 483.20000000000005, 483.59999999999997, 484.0, 484.4, 484.8, 485.20000000000005, 485.59999999999997, 486.0, 486.4, 486.8, 487.20000000000005, 487.59999999999997, 488.0, 488.4, 488.8, 489.20000000000005, 489.59999999999997, 490.0, 490.4, 490.8, 491.20000000000005, 491.59999999999997, 492.0, 492.4, 492.8, 493.20000000000005, 493.59999999999997, 494.0, 494.40000000000003, 494.8, 495.2, 495.59999999999997, 496.0, 496.40000000000003, 496.8, 497.2, 497.59999999999997, 498.0, 498.40000000000003, 498.8, 499.2, 499.59999999999997, 500.0, 500.4, 500.8, 501.2, 501.6, 502.0, 502.4, 502.8, 503.2, 503.6, 504.0, 504.4, 504.8, 505.2, 505.6, 506.0, 506.4, 506.8, 507.2, 507.6, 508.0, 508.4, 508.8, 509.2, 509.6000000000001, 510.0, 510.4, 510.8, 511.2, 511.6000000000001, 512.0, 512.4, 512.8000000000001, 513.2, 513.5999999999999, 514.0, 514.4, 514.8000000000001, 515.2, 515.5999999999999, 516.0, 516.4, 516.8000000000001, 517.2, 517.5999999999999, 518.0, 518.4, 518.8000000000001, 519.2, 519.5999999999999, 520.0, 520.4, 520.8000000000001, 521.2, 521.5999999999999, 522.0, 522.4, 522.8000000000001, 523.2, 523.5999999999999, 524.0, 524.4, 524.8000000000001, 525.2, 525.5999999999999, 526.0, 526.4, 526.8000000000001, 527.2, 527.5999999999999, 528.0, 528.4, 528.8000000000001, 529.2, 529.5999999999999, 530.0, 530.4, 530.8000000000001, 531.2, 531.5999999999999, 532.0, 532.4, 532.8000000000001, 533.2, 533.5999999999999, 534.0, 534.4, 534.8000000000001, 535.2, 535.5999999999999, 536.0, 536.4, 536.8000000000001, 537.2, 537.6, 538.0, 538.4, 538.8, 539.2, 539.6, 540.0, 540.4, 540.8, 541.2, 541.6, 542.0, 542.4, 542.8, 543.2, 543.6, 544.0, 544.4, 544.8, 545.2, 545.6, 546.0, 546.4, 546.8, 547.2, 547.6, 548.0, 548.4, 548.8, 549.2, 549.6, 550.0, 550.4, 550.8, 551.2, 551.6, 552.0, 552.4, 552.8, 553.2, 553.6, 554.0, 554.4, 554.8, 555.2, 555.6, 556.0, 556.4, 556.8, 557.2, 557.6, 558.0, 558.4, 558.8, 559.2, 559.6, 560.0, 560.4, 560.8, 561.2, 561.6, 562.0, 562.4, 562.8, 563.2, 563.6, 564.0, 564.4, 564.8, 565.2, 565.6, 566.0, 566.4, 566.8, 567.2, 567.6, 568.0, 568.4, 568.8, 569.2, 569.6, 570.0, 570.4, 570.8, 571.2, 571.6, 572.0, 572.4, 572.8, 573.2, 573.6, 574.0, 574.4, 574.8, 575.2, 575.6, 576.0, 576.4, 576.8, 577.2, 577.6, 578.0, 578.4, 578.8, 579.2, 579.6, 580.0, 580.4, 580.8, 581.2, 581.6, 582.0, 582.4, 582.8, 583.2, 583.6, 584.0, 584.4, 584.8, 585.2, 585.6, 586.0, 586.4, 586.8, 587.2, 587.6, 588.0, 588.4000000000001, 588.8, 589.1999999999999, 589.6, 590.0, 590.4000000000001, 590.8, 591.1999999999999, 591.6, 592.0, 592.4000000000001, 592.8, 593.1999999999999, 593.6, 594.0, 594.4000000000001, 594.8, 595.1999999999999, 595.6, 596.0, 596.4000000000001, 596.8, 597.1999999999999, 597.6, 598.0, 598.4000000000001, 598.8, 599.1999999999999, 599.6, 600.0, 600.4000000000001, 600.8, 601.1999999999999, 601.6, 602.0, 602.4000000000001, 602.8, 603.1999999999999, 603.6, 604.0, 604.4000000000001, 604.8, 605.1999999999999, 605.6, 606.0, 606.4000000000001, 606.8, 607.1999999999999, 607.6, 608.0, 608.4000000000001, 608.8, 609.1999999999999, 609.6, 610.0, 610.4000000000001, 610.8, 611.1999999999999, 611.6, 612.0, 612.4000000000001, 612.8, 613.1999999999999, 613.6, 614.0, 614.4, 614.8000000000001, 615.1999999999999, 615.6, 616.0, 616.4, 616.8000000000001, 617.1999999999999, 617.6, 618.0, 618.4, 618.8000000000001, 619.1999999999999, 619.6, 620.0, 620.4, 620.8000000000001, 621.1999999999999, 621.6, 622.0, 622.4, 622.8000000000001, 623.1999999999999, 623.6, 624.0, 624.4, 624.8000000000001, 625.1999999999999, 625.6, 626.0, 626.4, 626.8000000000001, 627.1999999999999, 627.6, 628.0, 628.4, 628.8000000000001, 629.1999999999999, 629.6, 630.0, 630.4, 630.8000000000001, 631.1999999999999, 631.6, 632.0, 632.4, 632.8000000000001, 633.1999999999999, 633.6, 634.0, 634.4, 634.8000000000001, 635.1999999999999, 635.6, 636.0, 636.4, 636.8000000000001, 637.2, 637.5999999999999, 638.0, 638.4, 638.8000000000001, 639.2, 639.5999999999999, 640.0, 640.4, 640.8000000000001, 641.2, 641.5999999999999, 642.0, 642.4, 642.8000000000001, 643.2, 643.5999999999999, 644.0, 644.4, 644.8000000000001, 645.2, 645.5999999999999, 646.0, 646.4, 646.8000000000001, 647.2, 647.5999999999999, 648.0, 648.4, 648.8000000000001, 649.2, 649.5999999999999, 650.0, 650.4, 650.8000000000001, 651.2, 651.5999999999999, 652.0, 652.4, 652.8000000000001, 653.2, 653.5999999999999, 654.0, 654.4, 654.8000000000001, 655.2, 655.5999999999999, 656.0, 656.4, 656.8000000000001, 657.2, 657.5999999999999, 658.0, 658.4, 658.8000000000001, 659.2, 659.5999999999999, 660.0, 660.4, 660.8000000000001, 661.2, 661.5999999999999, 662.0, 662.4, 662.8, 663.2, 663.6, 664.0, 664.4, 664.8, 665.2, 665.6, 666.0, 666.4, 666.8, 667.2, 667.6, 668.0, 668.4, 668.8, 669.2, 669.6, 670.0, 670.4, 670.8, 671.2, 671.6, 672.0, 672.4, 672.8, 673.2, 673.6, 674.0, 674.4, 674.8, 675.2, 675.6, 676.0, 676.4, 676.8, 677.2, 677.6, 678.0, 678.4, 678.8, 679.2, 679.6, 680.0, 680.4, 680.8, 681.2, 681.6, 682.0, 682.4, 682.8, 683.2, 683.6, 684.0, 684.4, 684.8, 685.2, 685.6, 686.0, 686.4, 686.8, 687.2, 687.6, 688.0, 688.4, 688.8, 689.2, 689.6, 690.0, 690.4, 690.8, 691.2, 691.6, 692.0, 692.4, 692.8, 693.2, 693.6, 694.0, 694.4, 694.8, 695.2, 695.6, 696.0, 696.4, 696.8, 697.2, 697.6, 698.0, 698.4, 698.8, 699.2, 699.6, 700.0, 700.4, 700.8, 701.2, 701.6, 702.0, 702.4, 702.8, 703.2, 703.6, 704.0, 704.4, 704.8, 705.2, 705.6, 706.0, 706.4, 706.8, 707.2, 707.6, 708.0, 708.4, 708.8, 709.2, 709.6, 710.0, 710.4, 710.8, 711.2, 711.6, 712.0, 712.4, 712.8, 713.1999999999999, 713.6, 714.0, 714.4000000000001, 714.8, 715.1999999999999, 715.6, 716.0, 716.4000000000001, 716.8, 717.1999999999999, 717.6, 718.0, 718.4000000000001, 718.8, 719.1999999999999, 719.6, 720.0, 720.4000000000001, 720.8, 721.1999999999999, 721.6, 722.0, 722.4000000000001, 722.8, 723.1999999999999, 723.6, 724.0, 724.4000000000001, 724.8, 725.1999999999999, 725.6, 726.0, 726.4000000000001, 726.8, 727.1999999999999, 727.6, 728.0, 728.4000000000001, 728.8, 729.1999999999999, 729.6, 730.0, 730.4000000000001, 730.8, 731.1999999999999, 731.6, 732.0, 732.4000000000001, 732.8, 733.1999999999999, 733.6, 734.0, 734.4000000000001, 734.8, 735.1999999999999, 735.6, 736.0, 736.4000000000001, 736.8, 737.1999999999999, 737.6, 738.0, 738.4, 738.8000000000001, 739.1999999999999, 739.6, 740.0, 740.4, 740.8000000000001, 741.1999999999999, 741.6, 742.0, 742.4, 742.8000000000001, 743.1999999999999, 743.6, 744.0, 744.4, 744.8000000000001, 745.1999999999999, 745.6, 746.0, 746.4, 746.8000000000001, 747.1999999999999, 747.6, 748.0, 748.4, 748.8000000000001, 749.1999999999999, 749.6, 750.0, 750.4, 750.8000000000001, 751.1999999999999, 751.6, 752.0, 752.4, 752.8000000000001, 753.1999999999999, 753.6, 754.0, 754.4, 754.8000000000001, 755.1999999999999, 755.6, 756.0, 756.4, 756.8000000000001, 757.1999999999999, 757.6, 758.0, 758.4, 758.8000000000001, 759.1999999999999, 759.6, 760.0, 760.4, 760.8000000000001, 761.1999999999999, 761.6, 762.0, 762.4, 762.8000000000001, 763.2, 763.5999999999999, 764.0, 764.4, 764.8000000000001, 765.2, 765.5999999999999, 766.0, 766.4, 766.8000000000001, 767.2, 767.5999999999999, 768.0, 768.4, 768.8000000000001, 769.2, 769.5999999999999, 770.0, 770.4, 770.8000000000001, 771.2, 771.5999999999999, 772.0, 772.4, 772.8000000000001, 773.2, 773.5999999999999, 774.0, 774.4, 774.8000000000001, 775.2, 775.5999999999999, 776.0, 776.4, 776.8000000000001, 777.2, 777.5999999999999, 778.0, 778.4, 778.8000000000001, 779.2, 779.5999999999999, 780.0, 780.4, 780.8000000000001, 781.2, 781.5999999999999, 782.0, 782.4, 782.8000000000001, 783.2, 783.5999999999999, 784.0, 784.4, 784.8000000000001, 785.2, 785.5999999999999, 786.0, 786.4, 786.8000000000001, 787.2, 787.6, 788.0, 788.4, 788.8, 789.2, 789.6, 790.0, 790.4, 790.8, 791.2, 791.6, 792.0, 792.4, 792.8, 793.2, 793.6, 794.0, 794.4, 794.8, 795.2, 795.6, 796.0, 796.4, 796.8, 797.2, 797.6, 798.0, 798.4, 798.8, 799.2, 799.6, 800.0, 800.4, 800.8, 801.2, 801.6, 802.0, 802.4, 802.8, 803.2, 803.6, 804.0, 804.4, 804.8, 805.2, 805.6, 806.0, 806.4, 806.8, 807.2, 807.6, 808.0, 808.4, 808.8, 809.2, 809.6, 810.0, 810.4, 810.8, 811.2, 811.6, 812.0, 812.4, 812.8, 813.2, 813.6, 814.0, 814.4, 814.8, 815.2, 815.6, 816.0, 816.4, 816.8, 817.2, 817.6, 818.0, 818.4, 818.8, 819.2, 819.6, 820.0, 820.4, 820.8, 821.2, 821.6, 822.0, 822.4, 822.8, 823.2, 823.6, 824.0, 824.4, 824.8, 825.2, 825.6, 826.0, 826.4, 826.8, 827.2, 827.6, 828.0, 828.4, 828.8, 829.2, 829.6, 830.0, 830.4, 830.8, 831.2, 831.6, 832.0, 832.4, 832.8, 833.2, 833.6, 834.0, 834.4, 834.8, 835.2, 835.6, 836.0, 836.4, 836.8, 837.2, 837.6, 838.0, 838.4000000000001, 838.8, 839.1999999999999, 839.6, 840.0, 840.4000000000001, 840.8, 841.1999999999999, 841.6, 842.0, 842.4000000000001, 842.8, 843.1999999999999, 843.6, 844.0, 844.4000000000001, 844.8, 845.1999999999999, 845.6, 846.0, 846.4000000000001, 846.8, 847.1999999999999, 847.6, 848.0, 848.4000000000001, 848.8, 849.1999999999999, 849.6, 850.0, 850.4000000000001, 850.8, 851.1999999999999, 851.6, 852.0, 852.4000000000001, 852.8, 853.1999999999999, 853.6, 854.0, 854.4000000000001, 854.8, 855.1999999999999, 855.6, 856.0, 856.4000000000001, 856.8, 857.1999999999999, 857.6, 858.0, 858.4000000000001, 858.8, 859.1999999999999, 859.6, 860.0, 860.4000000000001, 860.8, 861.1999999999999, 861.6, 862.0, 862.4000000000001, 862.8, 863.1999999999999, 863.6, 864.0, 864.4, 864.8000000000001, 865.1999999999999, 865.6, 866.0, 866.4, 866.8000000000001, 867.1999999999999, 867.6, 868.0, 868.4, 868.8000000000001, 869.1999999999999, 869.6, 870.0, 870.4, 870.8000000000001, 871.1999999999999, 871.6, 872.0, 872.4, 872.8000000000001, 873.1999999999999, 873.6, 874.0, 874.4, 874.8000000000001, 875.1999999999999, 875.6, 876.0, 876.4, 876.8000000000001, 877.1999999999999, 877.6, 878.0, 878.4, 878.8000000000001, 879.1999999999999, 879.6, 880.0, 880.4, 880.8000000000001, 881.1999999999999, 881.6, 882.0, 882.4, 882.8000000000001, 883.1999999999999, 883.6, 884.0, 884.4, 884.8000000000001, 885.1999999999999, 885.6, 886.0, 886.4, 886.8000000000001, 887.2, 887.5999999999999, 888.0, 888.4, 888.8000000000001, 889.2, 889.5999999999999, 890.0, 890.4, 890.8000000000001, 891.2, 891.5999999999999, 892.0, 892.4, 892.8000000000001, 893.2, 893.5999999999999, 894.0, 894.4, 894.8000000000001, 895.2, 895.5999999999999, 896.0, 896.4, 896.8000000000001, 897.2, 897.5999999999999, 898.0, 898.4, 898.8000000000001, 899.2, 899.5999999999999, 900.0, 900.4, 900.8000000000001, 901.2, 901.5999999999999, 902.0, 902.4, 902.8000000000001, 903.2, 903.5999999999999, 904.0, 904.4, 904.8000000000001, 905.2, 905.5999999999999, 906.0, 906.4, 906.8000000000001, 907.2, 907.5999999999999, 908.0, 908.4, 908.8000000000001, 909.2, 909.5999999999999, 910.0, 910.4, 910.8000000000001, 911.2, 911.5999999999999, 912.0, 912.4, 912.8, 913.2, 913.6, 914.0, 914.4, 914.8, 915.2, 915.6, 916.0, 916.4, 916.8, 917.2, 917.6, 918.0, 918.4, 918.8, 919.2, 919.6, 920.0, 920.4, 920.8, 921.2, 921.6, 922.0, 922.4, 922.8, 923.2, 923.6, 924.0, 924.4, 924.8, 925.2, 925.6, 926.0, 926.4, 926.8, 927.2, 927.6, 928.0, 928.4, 928.8, 929.2, 929.6, 930.0, 930.4, 930.8, 931.2, 931.6, 932.0, 932.4, 932.8, 933.2, 933.6, 934.0, 934.4, 934.8, 935.2, 935.6, 936.0, 936.4, 936.8, 937.2, 937.6, 938.0, 938.4, 938.8, 939.2, 939.6, 940.0, 940.4, 940.8, 941.2, 941.6, 942.0, 942.4, 942.8, 943.2, 943.6, 944.0, 944.4, 944.8, 945.2, 945.6, 946.0, 946.4, 946.8, 947.2, 947.6, 948.0, 948.4, 948.8, 949.2, 949.6, 950.0, 950.4, 950.8, 951.2, 951.6, 952.0, 952.4, 952.8, 953.2, 953.6, 954.0, 954.4, 954.8, 955.2, 955.6, 956.0, 956.4, 956.8, 957.2, 957.6, 958.0, 958.4, 958.8, 959.2, 959.6, 960.0, 960.4, 960.8, 961.2, 961.6, 962.0, 962.4, 962.8, 963.1999999999999, 963.6, 964.0, 964.4000000000001, 964.8, 965.1999999999999, 965.6, 966.0, 966.4000000000001, 966.8, 967.1999999999999, 967.6, 968.0, 968.4000000000001, 968.8, 969.1999999999999, 969.6, 970.0, 970.4000000000001, 970.8, 971.1999999999999, 971.6, 972.0, 972.4000000000001, 972.8, 973.1999999999999, 973.6, 974.0, 974.4000000000001, 974.8, 975.1999999999999, 975.6, 976.0, 976.4000000000001, 976.8, 977.1999999999999, 977.6, 978.0, 978.4000000000001, 978.8, 979.1999999999999, 979.6, 980.0, 980.4000000000001, 980.8, 981.1999999999999, 981.6, 982.0, 982.4000000000001, 982.8, 983.1999999999999, 983.6, 984.0, 984.4000000000001, 984.8, 985.1999999999999, 985.6, 986.0, 986.4000000000001, 986.8, 987.1999999999999, 987.6, 988.0, 988.4, 988.8000000000001, 989.1999999999999, 989.6, 990.0, 990.4, 990.8000000000001, 991.1999999999999, 991.6, 992.0, 992.4, 992.8000000000001, 993.1999999999999, 993.6, 994.0, 994.4, 994.8000000000001, 995.1999999999999, 995.6, 996.0, 996.4, 996.8000000000001, 997.1999999999999, 997.6, 998.0, 998.4, 998.8000000000001, 999.1999999999999, 999.6, 1000.0, 1000.4, 1000.8, 1001.2, 1001.6, 1002.0, 1002.4, 1002.8, 1003.2, 1003.6, 1004.0, 1004.4, 1004.8, 1005.2, 1005.6, 1006.0, 1006.4, 1006.8, 1007.2, 1007.6, 1008.0, 1008.4, 1008.8, 1009.2, 1009.6, 1010.0, 1010.4, 1010.8, 1011.2, 1011.6, 1012.0, 1012.4, 1012.8, 1013.2, 1013.6, 1014.0, 1014.4, 1014.8, 1015.2, 1015.6, 1016.0, 1016.4, 1016.8, 1017.2000000000002, 1017.6, 1018.0, 1018.4, 1018.8, 1019.2000000000002, 1019.6, 1020.0, 1020.4, 1020.8, 1021.2000000000002, 1021.6, 1022.0, 1022.4, 1022.8, 1023.2000000000002, 1023.6, 1024.0, 1024.4, 1024.8, 1025.1999999999998, 1025.6000000000001, 1026.0, 1026.4, 1026.8, 1027.1999999999998, 1027.6000000000001, 1028.0, 1028.4, 1028.8, 1029.1999999999998, 1029.6000000000001, 1030.0, 1030.4, 1030.8, 1031.1999999999998, 1031.6000000000001, 1032.0, 1032.4, 1032.8, 1033.1999999999998, 1033.6000000000001, 1034.0, 1034.4, 1034.8, 1035.1999999999998, 1035.6000000000001, 1036.0, 1036.4, 1036.8, 1037.1999999999998, 1037.6000000000001, 1038.0, 1038.4, 1038.8, 1039.1999999999998, 1039.6000000000001, 1040.0, 1040.4, 1040.8, 1041.1999999999998, 1041.6000000000001, 1042.0, 1042.4, 1042.8, 1043.1999999999998, 1043.6000000000001, 1044.0, 1044.4, 1044.8, 1045.1999999999998, 1045.6000000000001, 1046.0, 1046.4, 1046.8, 1047.1999999999998, 1047.6000000000001, 1048.0, 1048.4, 1048.8, 1049.1999999999998, 1049.6000000000001, 1050.0])
17
+
18
+ jupiter_alb_rf = np.array([0.2128, 0.2084, 0.2215, 0.2265, 0.2281, 0.2454, 0.2585, 0.2667, 0.2695, 0.2721, 0.2709, 0.2681, 0.2691, 0.2705, 0.2722, 0.2733, 0.2724, 0.2718, 0.2713, 0.2703, 0.2698, 0.2691, 0.2682, 0.2706, 0.272, 0.273, 0.2729, 0.2717, 0.2708, 0.2704, 0.2697, 0.2697, 0.2706, 0.2706, 0.2704, 0.2716, 0.2718, 0.2707, 0.2707, 0.2731, 0.2743, 0.2731, 0.2717, 0.2706, 0.2717, 0.2739, 0.2735, 0.2738, 0.2749, 0.2737, 0.2727, 0.2727, 0.2742, 0.2746, 0.2738, 0.2744, 0.2765, 0.2768, 0.2762, 0.2759, 0.2763, 0.2782, 0.2796, 0.2793, 0.2781, 0.2785, 0.2787, 0.2787, 0.2796, 0.281, 0.2814, 0.282, 0.282, 0.2822, 0.2821, 0.2829, 0.2846, 0.2841, 0.2843, 0.2852, 0.2864, 0.2866, 0.2867, 0.2877, 0.2891, 0.2874, 0.2873, 0.2897, 0.288, 0.2903, 0.2921, 0.2901, 0.2934, 0.2922, 0.2921, 0.2932, 0.2913, 0.293, 0.2952, 0.2933, 0.2929, 0.2946, 0.2949, 0.2967, 0.2961, 0.2959, 0.2959, 0.2961, 0.2964, 0.298, 0.3021, 0.3005, 0.2994, 0.299, 0.3018, 0.3038, 0.3034, 0.3031, 0.3059, 0.3054, 0.3047, 0.3043, 0.3055, 0.3079, 0.307, 0.3047, 0.3051, 0.3067, 0.3078, 0.3054, 0.3101, 0.3111, 0.3085, 0.3094, 0.3095, 0.3097, 0.3089, 0.3091, 0.3113, 0.3115, 0.3114, 0.3161, 0.3198, 0.3173, 0.3176, 0.3204, 0.3214, 0.319, 0.3173, 0.316, 0.3179, 0.3225, 0.3229, 0.3187, 0.3226, 0.3233, 0.3178, 0.3207, 0.3231, 0.3202, 0.3211, 0.325, 0.3245, 0.3222, 0.3213, 0.3215, 0.3237, 0.3247, 0.3245, 0.3271, 0.3268, 0.3294, 0.328, 0.3287, 0.3283, 0.3292, 0.3345, 0.3309, 0.3277, 0.3357, 0.3379, 0.3337, 0.335, 0.3447, 0.3415, 0.3385, 0.3474, 0.3434, 0.3366, 0.3415, 0.3418, 0.3413, 0.3395, 0.3385, 0.3402, 0.339, 0.3423, 0.3435, 0.3439, 0.3462, 0.342, 0.343, 0.3445, 0.3487, 0.3508, 0.353, 0.3557, 0.355, 0.3537, 0.356, 0.3525, 0.3497, 0.3509, 0.3515, 0.358, 0.354, 0.3512, 0.3579, 0.3552, 0.3554, 0.3571, 0.361, 0.3595, 0.359, 0.3605, 0.361, 0.361, 0.3583, 0.3587, 0.3629, 0.3637, 0.3641, 0.3786, 0.384, 0.3717, 0.3711, 0.3703, 0.3682, 0.3696, 0.3698, 0.3737, 0.386, 0.3853, 0.3749, 0.3708, 0.3718, 0.3731, 0.3709, 0.3731, 0.3741, 0.3744, 0.3771, 0.3744, 0.3745, 0.3741, 0.3741, 0.3738, 0.3785, 0.382, 0.3781, 0.383, 0.387, 0.381, 0.3838, 0.383, 0.3858, 0.387, 0.3872, 0.3908, 0.3898, 0.3888, 0.389, 0.388, 0.3902, 0.394, 0.3955, 0.3918, 0.394, 0.393, 0.3948, 0.3962, 0.3959, 0.3995, 0.4006, 0.3968, 0.4026, 0.4025, 0.4011, 0.4036, 0.4017, 0.4003, 0.4024, 0.4057, 0.4076, 0.4061, 0.4069, 0.41, 0.4088, 0.4085, 0.4137, 0.4113, 0.4087, 0.409, 0.4132, 0.4123, 0.4137, 0.4238, 0.4149, 0.417, 0.4188, 0.4152, 0.4184, 0.4221, 0.4186, 0.4231, 0.4215, 0.4179, 0.4296, 0.4271, 0.4212, 0.4241, 0.4259, 0.427, 0.4235, 0.4318, 0.434, 0.4376, 0.4415, 0.4282, 0.4225, 0.4355, 0.4357, 0.425, 0.4259, 0.4361, 0.4349, 0.4269, 0.4331, 0.435, 0.4326, 0.4322, 0.4326, 0.4352, 0.4372, 0.4335, 0.4456, 0.4447, 0.4381, 0.4387, 0.4375, 0.444, 0.4455, 0.4364, 0.441, 0.4411, 0.4411, 0.444, 0.4445, 0.4469, 0.4442, 0.4462, 0.4455, 0.4436, 0.4487, 0.4529, 0.4495, 0.4496, 0.4506, 0.4467, 0.4523, 0.454, 0.4503, 0.4538, 0.4532, 0.4534, 0.4515, 0.4505, 0.4534, 0.4553, 0.454, 0.4541, 0.4602, 0.4614, 0.459, 0.4559, 0.4564, 0.4581, 0.4583, 0.4594, 0.4602, 0.4587, 0.4592, 0.4602, 0.4619, 0.4627, 0.4624, 0.4642, 0.4637, 0.4617, 0.4619, 0.462, 0.4655, 0.47, 0.4692, 0.4684, 0.4688, 0.4691, 0.4699, 0.47, 0.4703, 0.471, 0.4716, 0.4718, 0.472, 0.473, 0.472, 0.4717, 0.4748, 0.4761, 0.4748, 0.4752, 0.4756, 0.4761, 0.4766, 0.4755, 0.4761, 0.4777, 0.4782, 0.4789, 0.479, 0.479, 0.4791, 0.48, 0.4809, 0.4802, 0.4803, 0.4831, 0.4829, 0.4817, 0.4821, 0.4829, 0.484, 0.4839, 0.4853, 0.4858, 0.485, 0.4853, 0.4872, 0.487, 0.4863, 0.4868, 0.4867, 0.4874, 0.4885, 0.4873, 0.4864, 0.4877, 0.4874, 0.4857, 0.4847, 0.4844, 0.4838, 0.4837, 0.484, 0.4791, 0.4843, 0.4924, 0.4923, 0.4938, 0.4954, 0.4943, 0.4936, 0.4966, 0.4983, 0.4966, 0.4968, 0.4961, 0.4974, 0.4989, 0.4977, 0.4998, 0.5011, 0.5003, 0.4989, 0.4999, 0.5, 0.499, 0.5004, 0.4998, 0.5027, 0.5044, 0.5023, 0.5024, 0.5035, 0.502, 0.5013, 0.5042, 0.5039, 0.5049, 0.505, 0.5035, 0.5061, 0.5067, 0.5053, 0.5072, 0.5076, 0.5061, 0.5066, 0.5063, 0.5047, 0.5087, 0.509, 0.5069, 0.5073, 0.5088, 0.5076, 0.5067, 0.5078, 0.508, 0.507, 0.5086, 0.51, 0.5088, 0.5068, 0.5085, 0.5105, 0.5095, 0.5106, 0.5126, 0.5115, 0.5115, 0.513, 0.514, 0.5141, 0.5138, 0.5162, 0.5172, 0.5162, 0.5164, 0.5167, 0.5151, 0.5164, 0.5199, 0.5259, 0.5234, 0.5179, 0.5221, 0.521, 0.5173, 0.5193, 0.5185, 0.5184, 0.5206, 0.5195, 0.5203, 0.5188, 0.5192, 0.5212, 0.521, 0.5207, 0.5202, 0.5197, 0.5199, 0.5227, 0.5232, 0.5216, 0.5234, 0.5268, 0.5276, 0.5257, 0.5247, 0.5256, 0.5249, 0.5225, 0.5233, 0.5255, 0.5253, 0.5246, 0.5233, 0.5247, 0.5256, 0.5249, 0.5277, 0.5275, 0.5251, 0.5259, 0.526, 0.5261, 0.526, 0.5241, 0.5231, 0.5228, 0.5235, 0.5246, 0.5218, 0.521, 0.5217, 0.5199, 0.5186, 0.5184, 0.516, 0.5139, 0.5118, 0.5095, 0.5066, 0.505, 0.5061, 0.5061, 0.5071, 0.5105, 0.5107, 0.5165, 0.5227, 0.5241, 0.5268, 0.529, 0.5315, 0.5331, 0.5327, 0.5351, 0.5373, 0.5355, 0.5337, 0.5312, 0.5277, 0.5255, 0.5258, 0.5298, 0.5355, 0.5338, 0.5273, 0.5265, 0.5312, 0.5329, 0.5335, 0.5342, 0.5357, 0.5366, 0.5379, 0.5387, 0.5391, 0.5402, 0.5407, 0.5421, 0.5431, 0.5446, 0.5448, 0.5466, 0.5471, 0.546, 0.5475, 0.5477, 0.5458, 0.5451, 0.5478, 0.5481, 0.5482, 0.5483, 0.5469, 0.5466, 0.5485, 0.5487, 0.547, 0.546, 0.5461, 0.5478, 0.5478, 0.5472, 0.5462, 0.5441, 0.543, 0.5432, 0.5437, 0.5441, 0.5431, 0.5428, 0.5441, 0.5452, 0.5458, 0.5445, 0.543, 0.5423, 0.5416, 0.5411, 0.5405, 0.5404, 0.5402, 0.5391, 0.5392, 0.5403, 0.5402, 0.5419, 0.5412, 0.5407, 0.5423, 0.5431, 0.5456, 0.546, 0.546, 0.547, 0.5472, 0.5472, 0.5485, 0.5488, 0.5498, 0.5502, 0.5495, 0.5493, 0.5487, 0.5493, 0.5495, 0.5489, 0.5494, 0.5487, 0.5478, 0.5494, 0.549, 0.5473, 0.5455, 0.5449, 0.545, 0.5451, 0.541, 0.5478, 0.5485, 0.5489, 0.5449, 0.5454, 0.5443, 0.5453, 0.5448, 0.5443, 0.5449, 0.5454, 0.5441, 0.5422, 0.5416, 0.5426, 0.5443, 0.5422, 0.5427, 0.5433, 0.5444, 0.5429, 0.542, 0.5434, 0.5439, 0.547, 0.5459, 0.5451, 0.5462, 0.5473, 0.5489, 0.5481, 0.5453, 0.544, 0.5457, 0.5461, 0.5445, 0.5435, 0.5451, 0.548, 0.5483, 0.5454, 0.5443, 0.5432, 0.542, 0.5448, 0.5447, 0.5435, 0.5432, 0.542, 0.5425, 0.5432, 0.5423, 0.5377, 0.5393, 0.5399, 0.5362, 0.5329, 0.5276, 0.5245, 0.5184, 0.5116, 0.5058, 0.5047, 0.4991, 0.4961, 0.4931, 0.4894, 0.484, 0.4836, 0.4756, 0.4673, 0.457, 0.4494, 0.448, 0.4475, 0.4492, 0.4513, 0.4556, 0.462, 0.4654, 0.4686, 0.4771, 0.486, 0.4921, 0.4961, 0.505, 0.5101, 0.5152, 0.5202, 0.5261, 0.5264, 0.5314, 0.5354, 0.5356, 0.5374, 0.5377, 0.5366, 0.54, 0.5435, 0.5431, 0.5456, 0.5453, 0.5443, 0.5472, 0.547, 0.5464, 0.5466, 0.5481, 0.5475, 0.5461, 0.5455, 0.5452, 0.5479, 0.5477, 0.5478, 0.5448, 0.5434, 0.5439, 0.5469, 0.5449, 0.5435, 0.5436, 0.5439, 0.5446, 0.5429, 0.5398, 0.5403, 0.5406, 0.5412, 0.5385, 0.5394, 0.5348, 0.53, 0.5286, 0.5219, 0.5168, 0.5091, 0.504, 0.5024, 0.5032, 0.5125, 0.5201, 0.5228, 0.5185, 0.5055, 0.4968, 0.5063, 0.5127, 0.5156, 0.5161, 0.5158, 0.5212, 0.5259, 0.5244, 0.5241, 0.5265, 0.5273, 0.5306, 0.5323, 0.5336, 0.5331, 0.5328, 0.5322, 0.5312, 0.5333, 0.5319, 0.5369, 0.5195, 0.5329, 0.5367, 0.5282, 0.5309, 0.5309, 0.5294, 0.529, 0.527, 0.5298, 0.5299, 0.5282, 0.5272, 0.5274, 0.5263, 0.5254, 0.5238, 0.5247, 0.5233, 0.5245, 0.523, 0.5234, 0.5242, 0.522, 0.5208, 0.5182, 0.5194, 0.5206, 0.5201, 0.5195, 0.5232, 0.5234, 0.5233, 0.5249, 0.526, 0.5292, 0.531, 0.5332, 0.5328, 0.5345, 0.5378, 0.539, 0.54, 0.5409, 0.5422, 0.5424, 0.5442, 0.5456, 0.5474, 0.5484, 0.549, 0.5476, 0.5474, 0.549, 0.5494, 0.5487, 0.547, 0.5463, 0.5446, 0.5423, 0.5421, 0.5423, 0.5423, 0.5412, 0.5375, 0.5343, 0.5333, 0.5342, 0.5364, 0.5363, 0.5351, 0.5365, 0.5373, 0.5372, 0.5358, 0.5408, 0.538, 0.5305, 0.5418, 0.5435, 0.544, 0.5389, 0.5409, 0.5421, 0.5407, 0.5388, 0.5392, 0.539, 0.5376, 0.5373, 0.5349, 0.5303, 0.5297, 0.5287, 0.5302, 0.5291, 0.528, 0.5267, 0.5243, 0.5231, 0.5216, 0.5203, 0.5195, 0.5181, 0.5133, 0.513, 0.5089, 0.5033, 0.4987, 0.4938, 0.4891, 0.4862, 0.486, 0.4821, 0.48, 0.4775, 0.4782, 0.4849, 0.4865, 0.4831, 0.48, 0.4799, 0.4865, 0.4919, 0.4937, 0.4946, 0.4955, 0.4967, 0.4996, 0.5023, 0.5042, 0.5055, 0.5067, 0.507, 0.5088, 0.5107, 0.5128, 0.5135, 0.5145, 0.5152, 0.5132, 0.5149, 0.5156, 0.5168, 0.5158, 0.5127, 0.5094, 0.5064, 0.5014, 0.4905, 0.4768, 0.4632, 0.4517, 0.4375, 0.4305, 0.4203, 0.4111, 0.4023, 0.385, 0.3692, 0.3578, 0.3568, 0.3559, 0.3497, 0.3419, 0.3405, 0.3323, 0.318, 0.3016, 0.2901, 0.2805, 0.2611, 0.2404, 0.2269, 0.2245, 0.2347, 0.2275, 0.212, 0.2083, 0.2181, 0.2379, 0.257, 0.2579, 0.2621, 0.2698, 0.2807, 0.2965, 0.3135, 0.327, 0.3352, 0.3421, 0.3445, 0.3538, 0.3685, 0.3807, 0.3891, 0.4048, 0.4127, 0.4194, 0.4366, 0.4433, 0.4434, 0.451, 0.4605, 0.4672, 0.4747, 0.4799, 0.4864, 0.4929, 0.4975, 0.5007, 0.5048, 0.508, 0.5122, 0.5144, 0.5133, 0.5163, 0.5188, 0.5206, 0.5226, 0.5232, 0.5228, 0.5238, 0.5218, 0.5228, 0.5227, 0.5253, 0.5248, 0.5235, 0.5221, 0.5207, 0.5213, 0.5216, 0.5196, 0.5214, 0.5197, 0.5187, 0.5202, 0.5202, 0.5201, 0.5159, 0.5163, 0.5156, 0.5144, 0.5133, 0.5106, 0.5099, 0.5106, 0.5097, 0.5078, 0.5056, 0.5033, 0.4993, 0.4955, 0.4939, 0.4888, 0.4856, 0.4915, 0.4855, 0.4784, 0.4815, 0.4816, 0.4781, 0.4717, 0.4732, 0.4632, 0.4628, 0.4727, 0.4713, 0.4747, 0.4707, 0.4698, 0.466, 0.4651, 0.4701, 0.4753, 0.4707, 0.4664, 0.4691, 0.471, 0.4714, 0.4722, 0.4699, 0.4742, 0.4752, 0.4785, 0.4777, 0.4759, 0.4703, 0.4714, 0.4737, 0.4703, 0.4659, 0.4629, 0.4652, 0.4674, 0.4682, 0.4653, 0.4642, 0.4614, 0.4574, 0.4565, 0.4568, 0.4533, 0.4509, 0.4481, 0.4447, 0.4356, 0.4288, 0.4197, 0.4116, 0.3987, 0.3966, 0.398, 0.3978, 0.3939, 0.3868, 0.3752, 0.3625, 0.3476, 0.3534, 0.3639, 0.3682, 0.3571, 0.3489, 0.3412, 0.3387, 0.348, 0.3498, 0.3545, 0.3562, 0.3554, 0.3557, 0.3561, 0.3588, 0.3647, 0.3568, 0.3483, 0.3402, 0.3375, 0.3444, 0.3469, 0.3507, 0.3598, 0.3614, 0.3574, 0.3537, 0.3573, 0.3562, 0.3564, 0.3588, 0.3565, 0.349, 0.3542, 0.3566, 0.3592, 0.3621, 0.3624, 0.37, 0.3753, 0.3749, 0.3838, 0.3949, 0.3887, 0.3864, 0.3962, 0.3991, 0.4028, 0.4081, 0.4083, 0.4095, 0.4162, 0.4228, 0.4222, 0.4242, 0.4277, 0.4344, 0.4393, 0.4409, 0.4433, 0.4458, 0.4435, 0.4394, 0.4405, 0.444, 0.445, 0.4452, 0.4459, 0.4511, 0.4548, 0.4542, 0.4531, 0.453, 0.4546, 0.4584, 0.4596, 0.4593, 0.4581, 0.4614, 0.4641, 0.4631, 0.4652, 0.4669, 0.4697, 0.4716, 0.4765, 0.48, 0.48, 0.4768, 0.4811, 0.4838, 0.4831, 0.4844, 0.4888, 0.4895, 0.4892, 0.4878, 0.485, 0.4885, 0.489, 0.4897, 0.4925, 0.4949, 0.496, 0.4947, 0.4945, 0.493, 0.4909, 0.4898, 0.494, 0.4964, 0.4972, 0.498, 0.4985, 0.499, 0.4969, 0.4963, 0.4971, 0.4971, 0.4963, 0.4961, 0.4964, 0.4965, 0.492, 0.4865, 0.4816, 0.4729, 0.4646, 0.457, 0.4498, 0.4455, 0.4371, 0.4359, 0.4253, 0.4168, 0.3954, 0.389, 0.3978, 0.4152, 0.4227, 0.4185, 0.4058, 0.3801, 0.3693, 0.3708, 0.3984, 0.4326, 0.4382, 0.4327, 0.4219, 0.4155, 0.4211, 0.4331, 0.4408, 0.4458, 0.4487, 0.4495, 0.452, 0.4513, 0.4543, 0.4546, 0.4578, 0.4582, 0.458, 0.4521, 0.4552, 0.4515, 0.4457, 0.4436, 0.4403, 0.4301, 0.4227, 0.4094, 0.3975, 0.3969, 0.3903, 0.3859, 0.3812, 0.3761, 0.3614, 0.3453, 0.3312, 0.3112, 0.2957, 0.2867, 0.2803, 0.2817, 0.273, 0.2623, 0.2488, 0.2246, 0.211, 0.1999, 0.1883, 0.1763, 0.1734, 0.1843, 0.1874, 0.1793, 0.1769, 0.1874, 0.1891, 0.1861, 0.1913, 0.2104, 0.2077, 0.2117, 0.2153, 0.2287, 0.25, 0.2598, 0.2682, 0.2764, 0.2791, 0.2811, 0.2819, 0.2913, 0.3049, 0.3191, 0.33, 0.3353, 0.3377, 0.3452, 0.3483, 0.3416, 0.3388, 0.3403, 0.3435, 0.3432, 0.3405, 0.3363, 0.3305, 0.325, 0.3155, 0.3057, 0.2999, 0.2891, 0.2635, 0.2276, 0.1825, 0.1486, 0.1346, 0.1122, 0.0767, 0.0543, 0.0455, 0.0406, 0.036, 0.0333, 0.0376, 0.0477, 0.0617, 0.0671, 0.0629, 0.0585, 0.0478, 0.0293, 0.0196, 0.0217, 0.0275, 0.0345, 0.0429, 0.0433, 0.0445, 0.0481, 0.0429, 0.0433, 0.0426, 0.0447, 0.0516, 0.0497, 0.0488, 0.0499, 0.0432, 0.038, 0.0404, 0.0426, 0.0434, 0.0409, 0.0411, 0.0514, 0.0648, 0.0776, 0.0799, 0.0777, 0.0761, 0.0691, 0.068, 0.0763, 0.0841, 0.0943, 0.1042, 0.1106, 0.1202, 0.1326, 0.1413, 0.1475, 0.1642, 0.1843, 0.2028, 0.2054, 0.2095, 0.2319, 0.2527, 0.261, 0.2683, 0.2853, 0.2991, 0.3125, 0.3209, 0.3267, 0.3378, 0.3539, 0.3657, 0.3628, 0.3618, 0.358, 0.3497, 0.3494, 0.3555, 0.3599, 0.3488, 0.3452, 0.3606, 0.3752, 0.3797, 0.378, 0.3776, 0.3806, 0.3728, 0.3683, 0.3799, 0.3876, 0.3962, 0.4113, 0.4149, 0.4107, 0.4134, 0.4099, 0.4098, 0.4165, 0.4179, 0.415, 0.4125, 0.4162, 0.4177, 0.4135, 0.4148, 0.4174, 0.4138, 0.4102, 0.4087, 0.4086, 0.4019, 0.397, 0.4036, 0.4127, 0.4206, 0.4191, 0.414, 0.4116, 0.4098, 0.4113, 0.4119, 0.4154, 0.4192, 0.4172, 0.4181, 0.4183, 0.4114, 0.4118, 0.416, 0.428, 0.4273, 0.4276, 0.4265, 0.4343, 0.4451, 0.4311, 0.4287, 0.4404, 0.4383, 0.4348, 0.4428, 0.4417, 0.443, 0.4531, 0.4577, 0.4586, 0.4573, 0.4562, 0.4549, 0.456, 0.4609, 0.4652, 0.4609, 0.4553, 0.4558, 0.458, 0.4561, 0.4566, 0.461, 0.4595, 0.4566, 0.4572, 0.457, 0.4506, 0.4436, 0.4283, 0.4181, 0.4382, 0.4536, 0.4522, 0.4503, 0.4484, 0.4498, 0.4461, 0.4316, 0.4095, 0.4074, 0.4313, 0.448, 0.4471, 0.4482, 0.4463, 0.4421, 0.4423, 0.441, 0.4393, 0.4361, 0.4334, 0.4234, 0.4131, 0.4143, 0.419, 0.4169, 0.4089, 0.3994, 0.3901, 0.3903, 0.3803, 0.3607, 0.3535, 0.3499, 0.3301, 0.3066, 0.2897, 0.2736, 0.2548, 0.2531, 0.257, 0.2465, 0.2302, 0.2252, 0.2219, 0.2209, 0.22, 0.2255, 0.2324, 0.2354, 0.2278, 0.2122, 0.2003, 0.1988, 0.1904, 0.1569, 0.1279, 0.1351, 0.165, 0.188, 0.2017, 0.2071, 0.2006, 0.1885, 0.1793, 0.1797, 0.1869, 0.1968, 0.2059, 0.2097, 0.2039, 0.1966, 0.1868, 0.1782, 0.1799, 0.1871, 0.187, 0.1818, 0.1794, 0.1775, 0.1664, 0.1543, 0.1495, 0.142, 0.1314, 0.1244, 0.1132, 0.0987, 0.098, 0.1019, 0.1016, 0.1069, 0.1147, 0.1161, 0.1125, 0.1102, 0.1046, 0.0931, 0.077, 0.0606, 0.0518, 0.0513, 0.0598, 0.0694, 0.0732, 0.0736, 0.0765, 0.0785, 0.0813, 0.0863, 0.0937, 0.0933, 0.0894, 0.0883, 0.0884, 0.0887, 0.0871, 0.0816, 0.0744, 0.0685, 0.0733, 0.087, 0.0984, 0.1063, 0.107, 0.0956, 0.0852, 0.0894, 0.1052, 0.113, 0.1136, 0.1144, 0.1122, 0.1077, 0.1053, 0.11, 0.1161, 0.1206, 0.1203, 0.1123, 0.0991, 0.0882, 0.0831, 0.0814, 0.078, 0.0751, 0.0784, 0.0824, 0.0807, 0.0794, 0.079, 0.0775, 0.0753, 0.073, 0.0711, 0.0708, 0.07, 0.0691, 0.0682, 0.0671, 0.071, 0.0786, 0.0838, 0.0838, 0.0813, 0.0817, 0.0843, 0.086, 0.0888, 0.0928, 0.1, 0.1083, 0.1139, 0.1152, 0.1206, 0.1294, 0.1336, 0.1317, 0.1302, 0.133, 0.1351, 0.1286, 0.1238, 0.1298, 0.1394, 0.1428, 0.1378, 0.1345, 0.1369, 0.1395, 0.14, 0.1387, 0.1392, 0.1424, 0.1448, 0.1433, 0.1409, 0.1433, 0.1537, 0.1646, 0.1727, 0.1752, 0.1795, 0.1783, 0.1751, 0.1774, 0.1777, 0.1734, 0.1708, 0.1723, 0.1738, 0.1724, 0.1763, 0.1795, 0.1857, 0.196, 0.2048, 0.2067, 0.2047, 0.2053, 0.214, 0.2235, 0.2306, 0.2356, 0.2445, 0.2543, 0.2616, 0.2603, 0.2596, 0.2679, 0.2774, 0.2764, 0.2812, 0.2868, 0.2895, 0.296, 0.301, 0.3088, 0.3188, 0.3271, 0.3281, 0.3274, 0.3355, 0.3407, 0.3492, 0.3546, 0.3628, 0.3732, 0.3731, 0.3797, 0.3823, 0.3907, 0.3974, 0.4007, 0.4039])
19
+
20
+
21
+ saturn_alb_wl = np.array([300.4, 300.8, 301.20000000000005, 301.59999999999997, 302.0, 302.4, 302.8, 303.20000000000005, 303.59999999999997, 304.0, 304.4, 304.8, 305.20000000000005, 305.59999999999997, 306.0, 306.4, 306.8, 307.2, 307.59999999999997, 308.0, 308.40000000000003, 308.8, 309.2, 309.59999999999997, 310.0, 310.40000000000003, 310.8, 311.2, 311.59999999999997, 312.0, 312.40000000000003, 312.8, 313.2, 313.59999999999997, 314.0, 314.40000000000003, 314.8, 315.2, 315.59999999999997, 316.0, 316.40000000000003, 316.8, 317.2, 317.59999999999997, 318.0, 318.40000000000003, 318.79999999999995, 319.2, 319.6, 320.0, 320.40000000000003, 320.79999999999995, 321.2, 321.6, 322.0, 322.40000000000003, 322.79999999999995, 323.2, 323.6, 324.0, 324.40000000000003, 324.79999999999995, 325.2, 325.6, 326.0, 326.40000000000003, 326.79999999999995, 327.2, 327.6, 328.0, 328.40000000000003, 328.79999999999995, 329.2, 329.6, 330.0, 330.40000000000003, 330.79999999999995, 331.2, 331.6, 332.0, 332.4, 332.8, 333.2, 333.6, 334.0, 334.4, 334.8, 335.2, 335.6, 336.0, 336.4, 336.8, 337.2, 337.6, 338.0, 338.4, 338.8, 339.2, 339.6, 340.0, 340.4, 340.8, 341.2, 341.6, 342.0, 342.4, 342.8, 343.2, 343.6, 344.0, 344.4, 344.8, 345.2, 345.6, 346.0, 346.4, 346.8, 347.2, 347.6, 348.0, 348.4, 348.8, 349.2, 349.6, 350.0, 350.4, 350.8, 351.2, 351.6, 352.0, 352.4, 352.8, 353.2, 353.6, 354.0, 354.4, 354.8, 355.2, 355.6, 356.0, 356.4, 356.8, 357.20000000000005, 357.59999999999997, 358.0, 358.4, 358.8, 359.20000000000005, 359.59999999999997, 360.0, 360.4, 360.8, 361.20000000000005, 361.59999999999997, 362.0, 362.4, 362.8, 363.20000000000005, 363.59999999999997, 364.0, 364.4, 364.8, 365.20000000000005, 365.59999999999997, 366.0, 366.4, 366.8, 367.20000000000005, 367.59999999999997, 368.0, 368.4, 368.8, 369.2, 369.59999999999997, 370.0, 370.40000000000003, 370.8, 371.2, 371.59999999999997, 372.0, 372.40000000000003, 372.8, 373.2, 373.59999999999997, 374.0, 374.40000000000003, 374.8, 375.2, 375.59999999999997, 376.0, 376.40000000000003, 376.8, 377.2, 377.59999999999997, 378.0, 378.40000000000003, 378.8, 379.2, 379.59999999999997, 380.0, 380.40000000000003, 380.8, 381.2, 381.6, 382.0, 382.40000000000003, 382.79999999999995, 383.2, 383.6, 384.0, 384.40000000000003, 384.79999999999995, 385.2, 385.6, 386.0, 386.40000000000003, 386.79999999999995, 387.2, 387.6, 388.0, 388.40000000000003, 388.79999999999995, 389.2, 389.6, 390.0, 390.40000000000003, 390.79999999999995, 391.2, 391.6, 392.0, 392.40000000000003, 392.79999999999995, 393.2, 393.6, 394.0, 394.4, 394.8, 395.2, 395.6, 396.0, 396.4, 396.8, 397.2, 397.6, 398.0, 398.4, 398.8, 399.2, 399.6, 400.0, 400.4, 400.8, 401.2, 401.6, 402.0, 402.4, 402.8, 403.2, 403.6, 404.0, 404.4, 404.8, 405.2, 405.6, 406.0, 406.4, 406.8, 407.2, 407.6, 408.0, 408.4, 408.8, 409.2, 409.6, 410.0, 410.4, 410.8, 411.2, 411.6, 412.0, 412.4, 412.8, 413.2, 413.6, 414.0, 414.4, 414.8, 415.2, 415.6, 416.0, 416.4, 416.8, 417.2, 417.6, 418.0, 418.4, 418.8, 419.20000000000005, 419.59999999999997, 420.0, 420.4, 420.8, 421.20000000000005, 421.59999999999997, 422.0, 422.4, 422.8, 423.20000000000005, 423.59999999999997, 424.0, 424.4, 424.8, 425.20000000000005, 425.59999999999997, 426.0, 426.4, 426.8, 427.20000000000005, 427.59999999999997, 428.0, 428.4, 428.8, 429.20000000000005, 429.59999999999997, 430.0, 430.4, 430.8, 431.20000000000005, 431.59999999999997, 432.0, 432.40000000000003, 432.8, 433.2, 433.59999999999997, 434.0, 434.40000000000003, 434.8, 435.2, 435.59999999999997, 436.0, 436.40000000000003, 436.8, 437.2, 437.59999999999997, 438.0, 438.40000000000003, 438.8, 439.2, 439.59999999999997, 440.0, 440.40000000000003, 440.8, 441.2, 441.59999999999997, 442.0, 442.40000000000003, 442.8, 443.2, 443.6, 444.0, 444.40000000000003, 444.79999999999995, 445.2, 445.6, 446.0, 446.40000000000003, 446.79999999999995, 447.2, 447.6, 448.0, 448.40000000000003, 448.79999999999995, 449.2, 449.6, 450.0, 450.40000000000003, 450.79999999999995, 451.2, 451.6, 452.0, 452.40000000000003, 452.79999999999995, 453.2, 453.6, 454.0, 454.40000000000003, 454.79999999999995, 455.2, 455.6, 456.0, 456.4, 456.8, 457.2, 457.6, 458.0, 458.4, 458.8, 459.2, 459.6, 460.0, 460.4, 460.8, 461.2, 461.6, 462.0, 462.4, 462.8, 463.2, 463.6, 464.0, 464.4, 464.8, 465.2, 465.6, 466.0, 466.4, 466.8, 467.2, 467.6, 468.0, 468.4, 468.8, 469.2, 469.6, 470.0, 470.4, 470.8, 471.2, 471.6, 472.0, 472.4, 472.8, 473.2, 473.6, 474.0, 474.4, 474.8, 475.2, 475.6, 476.0, 476.4, 476.8, 477.2, 477.6, 478.0, 478.4, 478.8, 479.2, 479.6, 480.0, 480.4, 480.8, 481.2, 481.59999999999997, 482.0, 482.4, 482.8, 483.20000000000005, 483.59999999999997, 484.0, 484.4, 484.8, 485.20000000000005, 485.59999999999997, 486.0, 486.4, 486.8, 487.20000000000005, 487.59999999999997, 488.0, 488.4, 488.8, 489.20000000000005, 489.59999999999997, 490.0, 490.4, 490.8, 491.20000000000005, 491.59999999999997, 492.0, 492.4, 492.8, 493.20000000000005, 493.59999999999997, 494.0, 494.40000000000003, 494.8, 495.2, 495.59999999999997, 496.0, 496.40000000000003, 496.8, 497.2, 497.59999999999997, 498.0, 498.40000000000003, 498.8, 499.2, 499.59999999999997, 500.0, 500.4, 500.8, 501.2, 501.6, 502.0, 502.4, 502.8, 503.2, 503.6, 504.0, 504.4, 504.8, 505.2, 505.6, 506.0, 506.4, 506.8, 507.2, 507.6, 508.0, 508.4, 508.8, 509.2, 509.6000000000001, 510.0, 510.4, 510.8, 511.2, 511.6000000000001, 512.0, 512.4, 512.8000000000001, 513.2, 513.5999999999999, 514.0, 514.4, 514.8000000000001, 515.2, 515.5999999999999, 516.0, 516.4, 516.8000000000001, 517.2, 517.5999999999999, 518.0, 518.4, 518.8000000000001, 519.2, 519.5999999999999, 520.0, 520.4, 520.8000000000001, 521.2, 521.5999999999999, 522.0, 522.4, 522.8000000000001, 523.2, 523.5999999999999, 524.0, 524.4, 524.8000000000001, 525.2, 525.5999999999999, 526.0, 526.4, 526.8000000000001, 527.2, 527.5999999999999, 528.0, 528.4, 528.8000000000001, 529.2, 529.5999999999999, 530.0, 530.4, 530.8000000000001, 531.2, 531.5999999999999, 532.0, 532.4, 532.8000000000001, 533.2, 533.5999999999999, 534.0, 534.4, 534.8000000000001, 535.2, 535.5999999999999, 536.0, 536.4, 536.8000000000001, 537.2, 537.6, 538.0, 538.4, 538.8, 539.2, 539.6, 540.0, 540.4, 540.8, 541.2, 541.6, 542.0, 542.4, 542.8, 543.2, 543.6, 544.0, 544.4, 544.8, 545.2, 545.6, 546.0, 546.4, 546.8, 547.2, 547.6, 548.0, 548.4, 548.8, 549.2, 549.6, 550.0, 550.4, 550.8, 551.2, 551.6, 552.0, 552.4, 552.8, 553.2, 553.6, 554.0, 554.4, 554.8, 555.2, 555.6, 556.0, 556.4, 556.8, 557.2, 557.6, 558.0, 558.4, 558.8, 559.2, 559.6, 560.0, 560.4, 560.8, 561.2, 561.6, 562.0, 562.4, 562.8, 563.2, 563.6, 564.0, 564.4, 564.8, 565.2, 565.6, 566.0, 566.4, 566.8, 567.2, 567.6, 568.0, 568.4, 568.8, 569.2, 569.6, 570.0, 570.4, 570.8, 571.2, 571.6, 572.0, 572.4, 572.8, 573.2, 573.6, 574.0, 574.4, 574.8, 575.2, 575.6, 576.0, 576.4, 576.8, 577.2, 577.6, 578.0, 578.4, 578.8, 579.2, 579.6, 580.0, 580.4, 580.8, 581.2, 581.6, 582.0, 582.4, 582.8, 583.2, 583.6, 584.0, 584.4, 584.8, 585.2, 585.6, 586.0, 586.4, 586.8, 587.2, 587.6, 588.0, 588.4000000000001, 588.8, 589.1999999999999, 589.6, 590.0, 590.4000000000001, 590.8, 591.1999999999999, 591.6, 592.0, 592.4000000000001, 592.8, 593.1999999999999, 593.6, 594.0, 594.4000000000001, 594.8, 595.1999999999999, 595.6, 596.0, 596.4000000000001, 596.8, 597.1999999999999, 597.6, 598.0, 598.4000000000001, 598.8, 599.1999999999999, 599.6, 600.0, 600.4000000000001, 600.8, 601.1999999999999, 601.6, 602.0, 602.4000000000001, 602.8, 603.1999999999999, 603.6, 604.0, 604.4000000000001, 604.8, 605.1999999999999, 605.6, 606.0, 606.4000000000001, 606.8, 607.1999999999999, 607.6, 608.0, 608.4000000000001, 608.8, 609.1999999999999, 609.6, 610.0, 610.4000000000001, 610.8, 611.1999999999999, 611.6, 612.0, 612.4000000000001, 612.8, 613.1999999999999, 613.6, 614.0, 614.4, 614.8000000000001, 615.1999999999999, 615.6, 616.0, 616.4, 616.8000000000001, 617.1999999999999, 617.6, 618.0, 618.4, 618.8000000000001, 619.1999999999999, 619.6, 620.0, 620.4, 620.8000000000001, 621.1999999999999, 621.6, 622.0, 622.4, 622.8000000000001, 623.1999999999999, 623.6, 624.0, 624.4, 624.8000000000001, 625.1999999999999, 625.6, 626.0, 626.4, 626.8000000000001, 627.1999999999999, 627.6, 628.0, 628.4, 628.8000000000001, 629.1999999999999, 629.6, 630.0, 630.4, 630.8000000000001, 631.1999999999999, 631.6, 632.0, 632.4, 632.8000000000001, 633.1999999999999, 633.6, 634.0, 634.4, 634.8000000000001, 635.1999999999999, 635.6, 636.0, 636.4, 636.8000000000001, 637.2, 637.5999999999999, 638.0, 638.4, 638.8000000000001, 639.2, 639.5999999999999, 640.0, 640.4, 640.8000000000001, 641.2, 641.5999999999999, 642.0, 642.4, 642.8000000000001, 643.2, 643.5999999999999, 644.0, 644.4, 644.8000000000001, 645.2, 645.5999999999999, 646.0, 646.4, 646.8000000000001, 647.2, 647.5999999999999, 648.0, 648.4, 648.8000000000001, 649.2, 649.5999999999999, 650.0, 650.4, 650.8000000000001, 651.2, 651.5999999999999, 652.0, 652.4, 652.8000000000001, 653.2, 653.5999999999999, 654.0, 654.4, 654.8000000000001, 655.2, 655.5999999999999, 656.0, 656.4, 656.8000000000001, 657.2, 657.5999999999999, 658.0, 658.4, 658.8000000000001, 659.2, 659.5999999999999, 660.0, 660.4, 660.8000000000001, 661.2, 661.5999999999999, 662.0, 662.4, 662.8, 663.2, 663.6, 664.0, 664.4, 664.8, 665.2, 665.6, 666.0, 666.4, 666.8, 667.2, 667.6, 668.0, 668.4, 668.8, 669.2, 669.6, 670.0, 670.4, 670.8, 671.2, 671.6, 672.0, 672.4, 672.8, 673.2, 673.6, 674.0, 674.4, 674.8, 675.2, 675.6, 676.0, 676.4, 676.8, 677.2, 677.6, 678.0, 678.4, 678.8, 679.2, 679.6, 680.0, 680.4, 680.8, 681.2, 681.6, 682.0, 682.4, 682.8, 683.2, 683.6, 684.0, 684.4, 684.8, 685.2, 685.6, 686.0, 686.4, 686.8, 687.2, 687.6, 688.0, 688.4, 688.8, 689.2, 689.6, 690.0, 690.4, 690.8, 691.2, 691.6, 692.0, 692.4, 692.8, 693.2, 693.6, 694.0, 694.4, 694.8, 695.2, 695.6, 696.0, 696.4, 696.8, 697.2, 697.6, 698.0, 698.4, 698.8, 699.2, 699.6, 700.0, 700.4, 700.8, 701.2, 701.6, 702.0, 702.4, 702.8, 703.2, 703.6, 704.0, 704.4, 704.8, 705.2, 705.6, 706.0, 706.4, 706.8, 707.2, 707.6, 708.0, 708.4, 708.8, 709.2, 709.6, 710.0, 710.4, 710.8, 711.2, 711.6, 712.0, 712.4, 712.8, 713.1999999999999, 713.6, 714.0, 714.4000000000001, 714.8, 715.1999999999999, 715.6, 716.0, 716.4000000000001, 716.8, 717.1999999999999, 717.6, 718.0, 718.4000000000001, 718.8, 719.1999999999999, 719.6, 720.0, 720.4000000000001, 720.8, 721.1999999999999, 721.6, 722.0, 722.4000000000001, 722.8, 723.1999999999999, 723.6, 724.0, 724.4000000000001, 724.8, 725.1999999999999, 725.6, 726.0, 726.4000000000001, 726.8, 727.1999999999999, 727.6, 728.0, 728.4000000000001, 728.8, 729.1999999999999, 729.6, 730.0, 730.4000000000001, 730.8, 731.1999999999999, 731.6, 732.0, 732.4000000000001, 732.8, 733.1999999999999, 733.6, 734.0, 734.4000000000001, 734.8, 735.1999999999999, 735.6, 736.0, 736.4000000000001, 736.8, 737.1999999999999, 737.6, 738.0, 738.4, 738.8000000000001, 739.1999999999999, 739.6, 740.0, 740.4, 740.8000000000001, 741.1999999999999, 741.6, 742.0, 742.4, 742.8000000000001, 743.1999999999999, 743.6, 744.0, 744.4, 744.8000000000001, 745.1999999999999, 745.6, 746.0, 746.4, 746.8000000000001, 747.1999999999999, 747.6, 748.0, 748.4, 748.8000000000001, 749.1999999999999, 749.6, 750.0, 750.4, 750.8000000000001, 751.1999999999999, 751.6, 752.0, 752.4, 752.8000000000001, 753.1999999999999, 753.6, 754.0, 754.4, 754.8000000000001, 755.1999999999999, 755.6, 756.0, 756.4, 756.8000000000001, 757.1999999999999, 757.6, 758.0, 758.4, 758.8000000000001, 759.1999999999999, 759.6, 760.0, 760.4, 760.8000000000001, 761.1999999999999, 761.6, 762.0, 762.4, 762.8000000000001, 763.2, 763.5999999999999, 764.0, 764.4, 764.8000000000001, 765.2, 765.5999999999999, 766.0, 766.4, 766.8000000000001, 767.2, 767.5999999999999, 768.0, 768.4, 768.8000000000001, 769.2, 769.5999999999999, 770.0, 770.4, 770.8000000000001, 771.2, 771.5999999999999, 772.0, 772.4, 772.8000000000001, 773.2, 773.5999999999999, 774.0, 774.4, 774.8000000000001, 775.2, 775.5999999999999, 776.0, 776.4, 776.8000000000001, 777.2, 777.5999999999999, 778.0, 778.4, 778.8000000000001, 779.2, 779.5999999999999, 780.0, 780.4, 780.8000000000001, 781.2, 781.5999999999999, 782.0, 782.4, 782.8000000000001, 783.2, 783.5999999999999, 784.0, 784.4, 784.8000000000001, 785.2, 785.5999999999999, 786.0, 786.4, 786.8000000000001, 787.2, 787.6, 788.0, 788.4, 788.8, 789.2, 789.6, 790.0, 790.4, 790.8, 791.2, 791.6, 792.0, 792.4, 792.8, 793.2, 793.6, 794.0, 794.4, 794.8, 795.2, 795.6, 796.0, 796.4, 796.8, 797.2, 797.6, 798.0, 798.4, 798.8, 799.2, 799.6, 800.0, 800.4, 800.8, 801.2, 801.6, 802.0, 802.4, 802.8, 803.2, 803.6, 804.0, 804.4, 804.8, 805.2, 805.6, 806.0, 806.4, 806.8, 807.2, 807.6, 808.0, 808.4, 808.8, 809.2, 809.6, 810.0, 810.4, 810.8, 811.2, 811.6, 812.0, 812.4, 812.8, 813.2, 813.6, 814.0, 814.4, 814.8, 815.2, 815.6, 816.0, 816.4, 816.8, 817.2, 817.6, 818.0, 818.4, 818.8, 819.2, 819.6, 820.0, 820.4, 820.8, 821.2, 821.6, 822.0, 822.4, 822.8, 823.2, 823.6, 824.0, 824.4, 824.8, 825.2, 825.6, 826.0, 826.4, 826.8, 827.2, 827.6, 828.0, 828.4, 828.8, 829.2, 829.6, 830.0, 830.4, 830.8, 831.2, 831.6, 832.0, 832.4, 832.8, 833.2, 833.6, 834.0, 834.4, 834.8, 835.2, 835.6, 836.0, 836.4, 836.8, 837.2, 837.6, 838.0, 838.4000000000001, 838.8, 839.1999999999999, 839.6, 840.0, 840.4000000000001, 840.8, 841.1999999999999, 841.6, 842.0, 842.4000000000001, 842.8, 843.1999999999999, 843.6, 844.0, 844.4000000000001, 844.8, 845.1999999999999, 845.6, 846.0, 846.4000000000001, 846.8, 847.1999999999999, 847.6, 848.0, 848.4000000000001, 848.8, 849.1999999999999, 849.6, 850.0, 850.4000000000001, 850.8, 851.1999999999999, 851.6, 852.0, 852.4000000000001, 852.8, 853.1999999999999, 853.6, 854.0, 854.4000000000001, 854.8, 855.1999999999999, 855.6, 856.0, 856.4000000000001, 856.8, 857.1999999999999, 857.6, 858.0, 858.4000000000001, 858.8, 859.1999999999999, 859.6, 860.0, 860.4000000000001, 860.8, 861.1999999999999, 861.6, 862.0, 862.4000000000001, 862.8, 863.1999999999999, 863.6, 864.0, 864.4, 864.8000000000001, 865.1999999999999, 865.6, 866.0, 866.4, 866.8000000000001, 867.1999999999999, 867.6, 868.0, 868.4, 868.8000000000001, 869.1999999999999, 869.6, 870.0, 870.4, 870.8000000000001, 871.1999999999999, 871.6, 872.0, 872.4, 872.8000000000001, 873.1999999999999, 873.6, 874.0, 874.4, 874.8000000000001, 875.1999999999999, 875.6, 876.0, 876.4, 876.8000000000001, 877.1999999999999, 877.6, 878.0, 878.4, 878.8000000000001, 879.1999999999999, 879.6, 880.0, 880.4, 880.8000000000001, 881.1999999999999, 881.6, 882.0, 882.4, 882.8000000000001, 883.1999999999999, 883.6, 884.0, 884.4, 884.8000000000001, 885.1999999999999, 885.6, 886.0, 886.4, 886.8000000000001, 887.2, 887.5999999999999, 888.0, 888.4, 888.8000000000001, 889.2, 889.5999999999999, 890.0, 890.4, 890.8000000000001, 891.2, 891.5999999999999, 892.0, 892.4, 892.8000000000001, 893.2, 893.5999999999999, 894.0, 894.4, 894.8000000000001, 895.2, 895.5999999999999, 896.0, 896.4, 896.8000000000001, 897.2, 897.5999999999999, 898.0, 898.4, 898.8000000000001, 899.2, 899.5999999999999, 900.0, 900.4, 900.8000000000001, 901.2, 901.5999999999999, 902.0, 902.4, 902.8000000000001, 903.2, 903.5999999999999, 904.0, 904.4, 904.8000000000001, 905.2, 905.5999999999999, 906.0, 906.4, 906.8000000000001, 907.2, 907.5999999999999, 908.0, 908.4, 908.8000000000001, 909.2, 909.5999999999999, 910.0, 910.4, 910.8000000000001, 911.2, 911.5999999999999, 912.0, 912.4, 912.8, 913.2, 913.6, 914.0, 914.4, 914.8, 915.2, 915.6, 916.0, 916.4, 916.8, 917.2, 917.6, 918.0, 918.4, 918.8, 919.2, 919.6, 920.0, 920.4, 920.8, 921.2, 921.6, 922.0, 922.4, 922.8, 923.2, 923.6, 924.0, 924.4, 924.8, 925.2, 925.6, 926.0, 926.4, 926.8, 927.2, 927.6, 928.0, 928.4, 928.8, 929.2, 929.6, 930.0, 930.4, 930.8, 931.2, 931.6, 932.0, 932.4, 932.8, 933.2, 933.6, 934.0, 934.4, 934.8, 935.2, 935.6, 936.0, 936.4, 936.8, 937.2, 937.6, 938.0, 938.4, 938.8, 939.2, 939.6, 940.0, 940.4, 940.8, 941.2, 941.6, 942.0, 942.4, 942.8, 943.2, 943.6, 944.0, 944.4, 944.8, 945.2, 945.6, 946.0, 946.4, 946.8, 947.2, 947.6, 948.0, 948.4, 948.8, 949.2, 949.6, 950.0, 950.4, 950.8, 951.2, 951.6, 952.0, 952.4, 952.8, 953.2, 953.6, 954.0, 954.4, 954.8, 955.2, 955.6, 956.0, 956.4, 956.8, 957.2, 957.6, 958.0, 958.4, 958.8, 959.2, 959.6, 960.0, 960.4, 960.8, 961.2, 961.6, 962.0, 962.4, 962.8, 963.1999999999999, 963.6, 964.0, 964.4000000000001, 964.8, 965.1999999999999, 965.6, 966.0, 966.4000000000001, 966.8, 967.1999999999999, 967.6, 968.0, 968.4000000000001, 968.8, 969.1999999999999, 969.6, 970.0, 970.4000000000001, 970.8, 971.1999999999999, 971.6, 972.0, 972.4000000000001, 972.8, 973.1999999999999, 973.6, 974.0, 974.4000000000001, 974.8, 975.1999999999999, 975.6, 976.0, 976.4000000000001, 976.8, 977.1999999999999, 977.6, 978.0, 978.4000000000001, 978.8, 979.1999999999999, 979.6, 980.0, 980.4000000000001, 980.8, 981.1999999999999, 981.6, 982.0, 982.4000000000001, 982.8, 983.1999999999999, 983.6, 984.0, 984.4000000000001, 984.8, 985.1999999999999, 985.6, 986.0, 986.4000000000001, 986.8, 987.1999999999999, 987.6, 988.0, 988.4, 988.8000000000001, 989.1999999999999, 989.6, 990.0, 990.4, 990.8000000000001, 991.1999999999999, 991.6, 992.0, 992.4, 992.8000000000001, 993.1999999999999, 993.6, 994.0, 994.4, 994.8000000000001, 995.1999999999999, 995.6, 996.0, 996.4, 996.8000000000001, 997.1999999999999, 997.6, 998.0, 998.4, 998.8000000000001, 999.1999999999999, 999.6, 1000.0, 1000.4, 1000.8, 1001.2, 1001.6, 1002.0, 1002.4, 1002.8, 1003.2, 1003.6, 1004.0, 1004.4, 1004.8, 1005.2, 1005.6, 1006.0, 1006.4, 1006.8, 1007.2, 1007.6, 1008.0, 1008.4, 1008.8, 1009.2, 1009.6, 1010.0, 1010.4, 1010.8, 1011.2, 1011.6, 1012.0, 1012.4, 1012.8, 1013.2, 1013.6, 1014.0, 1014.4, 1014.8, 1015.2, 1015.6, 1016.0, 1016.4, 1016.8, 1017.2000000000002, 1017.6, 1018.0, 1018.4, 1018.8, 1019.2000000000002, 1019.6, 1020.0, 1020.4, 1020.8, 1021.2000000000002, 1021.6, 1022.0, 1022.4, 1022.8, 1023.2000000000002, 1023.6, 1024.0, 1024.4, 1024.8, 1025.1999999999998, 1025.6000000000001, 1026.0, 1026.4, 1026.8, 1027.1999999999998, 1027.6000000000001, 1028.0, 1028.4, 1028.8, 1029.1999999999998, 1029.6000000000001, 1030.0, 1030.4, 1030.8, 1031.1999999999998, 1031.6000000000001, 1032.0, 1032.4, 1032.8, 1033.1999999999998, 1033.6000000000001, 1034.0, 1034.4, 1034.8, 1035.1999999999998, 1035.6000000000001, 1036.0, 1036.4, 1036.8, 1037.1999999999998, 1037.6000000000001, 1038.0, 1038.4, 1038.8, 1039.1999999999998, 1039.6000000000001, 1040.0, 1040.4, 1040.8, 1041.1999999999998, 1041.6000000000001, 1042.0, 1042.4, 1042.8, 1043.1999999999998, 1043.6000000000001, 1044.0, 1044.4, 1044.8, 1045.1999999999998, 1045.6000000000001, 1046.0, 1046.4, 1046.8, 1047.1999999999998, 1047.6000000000001, 1048.0, 1048.4, 1048.8, 1049.1999999999998, 1049.6000000000001, 1050.0])
22
+
23
+ saturn_alb_rf = np.array([0.2532, 0.2359, 0.2404, 0.2456, 0.2363, 0.2425, 0.2459, 0.2468, 0.2454, 0.2446, 0.2431, 0.2366, 0.2328, 0.2335, 0.2339, 0.232, 0.2293, 0.2303, 0.2322, 0.2315, 0.2316, 0.2313, 0.2306, 0.2341, 0.236, 0.2368, 0.2365, 0.2357, 0.2349, 0.2351, 0.2343, 0.2341, 0.2345, 0.2337, 0.2321, 0.2322, 0.2324, 0.2309, 0.2297, 0.23, 0.2303, 0.2292, 0.2274, 0.2254, 0.2255, 0.2269, 0.2263, 0.2261, 0.2261, 0.2243, 0.2226, 0.2218, 0.2231, 0.2237, 0.2222, 0.2213, 0.222, 0.2213, 0.2206, 0.2209, 0.2202, 0.2201, 0.2205, 0.2195, 0.2177, 0.2176, 0.2175, 0.217, 0.2172, 0.2172, 0.2163, 0.2165, 0.2161, 0.2156, 0.2148, 0.2144, 0.2148, 0.2138, 0.2134, 0.2134, 0.2135, 0.213, 0.2124, 0.2126, 0.2129, 0.2108, 0.2103, 0.2114, 0.2092, 0.21, 0.2105, 0.2084, 0.2105, 0.2089, 0.2082, 0.2085, 0.2065, 0.2071, 0.2081, 0.2059, 0.205, 0.2056, 0.2051, 0.2059, 0.2049, 0.2042, 0.2038, 0.203, 0.2022, 0.2029, 0.2056, 0.2039, 0.2024, 0.2015, 0.2028, 0.2037, 0.2026, 0.2018, 0.2033, 0.2023, 0.201, 0.2003, 0.2007, 0.2022, 0.2009, 0.1985, 0.1983, 0.1992, 0.1997, 0.1974, 0.1998, 0.2004, 0.1982, 0.1983, 0.1977, 0.1973, 0.1963, 0.196, 0.1968, 0.1964, 0.196, 0.1987, 0.2005, 0.1983, 0.1979, 0.1997, 0.2002, 0.1983, 0.1967, 0.1954, 0.1961, 0.1988, 0.1987, 0.1955, 0.1976, 0.1979, 0.1939, 0.1954, 0.1969, 0.1946, 0.1948, 0.1972, 0.1966, 0.1949, 0.1942, 0.1938, 0.1949, 0.1952, 0.1945, 0.1959, 0.1958, 0.197, 0.196, 0.1963, 0.1957, 0.1957, 0.1991, 0.1969, 0.1942, 0.1986, 0.2001, 0.197, 0.1974, 0.2031, 0.2012, 0.1989, 0.2044, 0.2021, 0.1977, 0.2004, 0.2007, 0.2004, 0.1992, 0.1985, 0.1993, 0.1983, 0.2001, 0.201, 0.201, 0.2026, 0.2003, 0.2005, 0.201, 0.2034, 0.2047, 0.2057, 0.2071, 0.207, 0.2061, 0.2074, 0.2057, 0.204, 0.2045, 0.2046, 0.2085, 0.2064, 0.2045, 0.2087, 0.2074, 0.2075, 0.2088, 0.2112, 0.2105, 0.2102, 0.2113, 0.2122, 0.2128, 0.2113, 0.2115, 0.2141, 0.2143, 0.2137, 0.2216, 0.2256, 0.2195, 0.2197, 0.2198, 0.2188, 0.2196, 0.2196, 0.2212, 0.2286, 0.2299, 0.2247, 0.2223, 0.2231, 0.2242, 0.2231, 0.2246, 0.2255, 0.2261, 0.2283, 0.2271, 0.2277, 0.228, 0.2283, 0.2285, 0.2315, 0.2343, 0.232, 0.2351, 0.2381, 0.2346, 0.2369, 0.2372, 0.2395, 0.2404, 0.2405, 0.2433, 0.243, 0.2428, 0.2433, 0.243, 0.2446, 0.2473, 0.249, 0.2473, 0.2494, 0.2494, 0.2509, 0.2521, 0.2525, 0.2556, 0.2568, 0.2546, 0.2587, 0.2595, 0.259, 0.2614, 0.2609, 0.2604, 0.2622, 0.2649, 0.2669, 0.2664, 0.2671, 0.2699, 0.2698, 0.2701, 0.2743, 0.2738, 0.2726, 0.2731, 0.2763, 0.2762, 0.2773, 0.2846, 0.2795, 0.2816, 0.2836, 0.2816, 0.2842, 0.2874, 0.2856, 0.2893, 0.2889, 0.2865, 0.2948, 0.2945, 0.2912, 0.2935, 0.295, 0.2964, 0.2943, 0.3, 0.3024, 0.3055, 0.3093, 0.3012, 0.2977, 0.3075, 0.3088, 0.3017, 0.3028, 0.3103, 0.3107, 0.3055, 0.3104, 0.3127, 0.3114, 0.3117, 0.3126, 0.3149, 0.3168, 0.3145, 0.3237, 0.3243, 0.3197, 0.3209, 0.3207, 0.326, 0.3277, 0.3217, 0.3258, 0.3269, 0.3273, 0.3297, 0.3302, 0.3326, 0.3315, 0.333, 0.3332, 0.3325, 0.3365, 0.3403, 0.3383, 0.339, 0.3404, 0.3378, 0.3423, 0.3442, 0.3418, 0.3452, 0.3455, 0.3466, 0.3459, 0.3452, 0.3475, 0.3495, 0.3489, 0.3494, 0.3543, 0.3554, 0.3546, 0.353, 0.3533, 0.3552, 0.356, 0.3571, 0.3584, 0.3575, 0.3585, 0.3599, 0.3611, 0.3621, 0.3623, 0.3639, 0.3642, 0.3634, 0.3639, 0.3645, 0.367, 0.3707, 0.371, 0.3704, 0.371, 0.3719, 0.3726, 0.3732, 0.3739, 0.3742, 0.375, 0.376, 0.3764, 0.3777, 0.3769, 0.3767, 0.3799, 0.3814, 0.3806, 0.3815, 0.3819, 0.3823, 0.3835, 0.383, 0.3835, 0.3853, 0.3859, 0.3866, 0.3873, 0.3875, 0.3875, 0.3888, 0.3899, 0.3898, 0.3904, 0.3928, 0.3929, 0.3927, 0.3932, 0.3933, 0.3948, 0.3952, 0.3963, 0.397, 0.3966, 0.3968, 0.3988, 0.3993, 0.3986, 0.3995, 0.3998, 0.4006, 0.402, 0.4015, 0.401, 0.4024, 0.4029, 0.4026, 0.4025, 0.4025, 0.4013, 0.4006, 0.4006, 0.3969, 0.4018, 0.4087, 0.4089, 0.4107, 0.4127, 0.4124, 0.4116, 0.4141, 0.4161, 0.4146, 0.4149, 0.4148, 0.4155, 0.4169, 0.4159, 0.4182, 0.4196, 0.4194, 0.4184, 0.4194, 0.4203, 0.4197, 0.4209, 0.4204, 0.4227, 0.4247, 0.4233, 0.4238, 0.425, 0.4239, 0.4232, 0.4264, 0.4262, 0.4272, 0.4278, 0.4269, 0.4291, 0.43, 0.4289, 0.4305, 0.4313, 0.4307, 0.4317, 0.4319, 0.4307, 0.434, 0.4345, 0.4329, 0.4336, 0.4349, 0.4344, 0.434, 0.4351, 0.4353, 0.4347, 0.4362, 0.4378, 0.4373, 0.4357, 0.4366, 0.4387, 0.4383, 0.4395, 0.4414, 0.4408, 0.441, 0.4426, 0.4438, 0.444, 0.4437, 0.446, 0.447, 0.4465, 0.4473, 0.4481, 0.4465, 0.4473, 0.4506, 0.4558, 0.4542, 0.4496, 0.4538, 0.4541, 0.4518, 0.4531, 0.4523, 0.4521, 0.4543, 0.4537, 0.4547, 0.4539, 0.4538, 0.4558, 0.4562, 0.4565, 0.456, 0.4559, 0.4562, 0.4587, 0.4592, 0.4577, 0.4593, 0.4623, 0.4637, 0.4626, 0.4618, 0.4631, 0.463, 0.4616, 0.4623, 0.4643, 0.4644, 0.4645, 0.4639, 0.465, 0.4658, 0.4655, 0.4681, 0.4683, 0.4666, 0.4674, 0.4678, 0.4681, 0.4683, 0.4672, 0.4664, 0.4662, 0.4669, 0.4683, 0.4661, 0.4655, 0.4666, 0.4657, 0.465, 0.4648, 0.4629, 0.4616, 0.4599, 0.4577, 0.4547, 0.4524, 0.4526, 0.4525, 0.4544, 0.4589, 0.4601, 0.4657, 0.4724, 0.4746, 0.4779, 0.4804, 0.4831, 0.4848, 0.4847, 0.4867, 0.4889, 0.4883, 0.488, 0.4877, 0.4871, 0.4869, 0.4868, 0.4883, 0.4917, 0.4917, 0.4896, 0.4893, 0.4917, 0.4924, 0.4925, 0.4929, 0.4941, 0.4949, 0.4958, 0.4964, 0.4965, 0.4979, 0.4985, 0.4995, 0.5008, 0.5024, 0.5028, 0.5045, 0.5048, 0.5042, 0.5059, 0.5064, 0.505, 0.5049, 0.5074, 0.5078, 0.5083, 0.5089, 0.5083, 0.5082, 0.5098, 0.5106, 0.5094, 0.509, 0.5094, 0.5113, 0.5116, 0.5116, 0.5112, 0.5102, 0.5103, 0.5112, 0.5115, 0.512, 0.5118, 0.5113, 0.5113, 0.5128, 0.5139, 0.5134, 0.5125, 0.512, 0.5117, 0.5113, 0.5108, 0.5106, 0.5105, 0.5095, 0.5101, 0.5104, 0.5115, 0.5136, 0.5135, 0.5139, 0.5141, 0.516, 0.519, 0.5202, 0.5211, 0.5221, 0.5229, 0.5243, 0.5254, 0.5263, 0.5284, 0.5284, 0.5294, 0.5297, 0.529, 0.5298, 0.5305, 0.5317, 0.5326, 0.5314, 0.5321, 0.5339, 0.5344, 0.5336, 0.5327, 0.5346, 0.5363, 0.5334, 0.5341, 0.5372, 0.5391, 0.5373, 0.5369, 0.5361, 0.5364, 0.5386, 0.537, 0.5357, 0.5362, 0.5362, 0.5364, 0.5356, 0.5333, 0.5345, 0.5353, 0.5344, 0.5341, 0.5362, 0.5374, 0.5368, 0.5374, 0.537, 0.5392, 0.54, 0.5408, 0.5434, 0.5434, 0.543, 0.5445, 0.5452, 0.5454, 0.545, 0.5444, 0.5457, 0.5458, 0.5455, 0.5473, 0.5484, 0.5476, 0.5475, 0.5473, 0.5442, 0.5443, 0.544, 0.5442, 0.5473, 0.5441, 0.5406, 0.5419, 0.5436, 0.5423, 0.5413, 0.5399, 0.5376, 0.5353, 0.5311, 0.5272, 0.5207, 0.5123, 0.5048, 0.502, 0.4969, 0.4905, 0.4887, 0.4871, 0.4829, 0.4813, 0.476, 0.468, 0.4593, 0.4456, 0.4347, 0.4307, 0.4313, 0.4318, 0.435, 0.4426, 0.4495, 0.4543, 0.46, 0.4709, 0.4826, 0.4902, 0.5001, 0.51, 0.5175, 0.527, 0.5341, 0.5364, 0.5411, 0.5478, 0.5504, 0.5531, 0.5561, 0.5535, 0.5576, 0.5633, 0.5639, 0.5651, 0.5659, 0.5657, 0.5683, 0.5687, 0.569, 0.5696, 0.5699, 0.5698, 0.5717, 0.5721, 0.5696, 0.5679, 0.5692, 0.5743, 0.5745, 0.5718, 0.5681, 0.5677, 0.5723, 0.5716, 0.5707, 0.5706, 0.5707, 0.5702, 0.571, 0.5703, 0.5705, 0.5698, 0.5668, 0.5686, 0.5687, 0.566, 0.5653, 0.5616, 0.5603, 0.5598, 0.5601, 0.5584, 0.5557, 0.5585, 0.5607, 0.5612, 0.5663, 0.5659, 0.5604, 0.557, 0.5581, 0.5617, 0.5627, 0.5599, 0.5626, 0.5649, 0.5631, 0.5639, 0.565, 0.5661, 0.5672, 0.567, 0.5686, 0.5694, 0.5699, 0.568, 0.566, 0.5674, 0.5666, 0.5642, 0.5615, 0.5605, 0.5617, 0.5575, 0.5604, 0.5617, 0.5612, 0.5589, 0.5573, 0.5578, 0.5583, 0.558, 0.557, 0.5556, 0.554, 0.5546, 0.553, 0.5496, 0.5495, 0.5498, 0.5483, 0.5491, 0.5505, 0.5483, 0.5465, 0.5424, 0.5398, 0.5425, 0.5423, 0.5414, 0.5442, 0.5461, 0.5481, 0.5499, 0.5509, 0.5532, 0.5571, 0.5611, 0.5629, 0.5658, 0.5707, 0.5728, 0.5763, 0.5775, 0.5791, 0.5835, 0.5845, 0.5874, 0.5903, 0.5905, 0.5923, 0.5951, 0.5964, 0.5976, 0.598, 0.5966, 0.5959, 0.5951, 0.5944, 0.5939, 0.5886, 0.5855, 0.5853, 0.5852, 0.5851, 0.581, 0.5757, 0.5736, 0.5748, 0.5746, 0.5766, 0.5819, 0.5807, 0.581, 0.5824, 0.5816, 0.584, 0.5808, 0.587, 0.5866, 0.5871, 0.5874, 0.5891, 0.5927, 0.594, 0.5934, 0.5925, 0.5919, 0.5903, 0.5933, 0.5946, 0.592, 0.591, 0.5878, 0.586, 0.5874, 0.5866, 0.5871, 0.5851, 0.5806, 0.5768, 0.5719, 0.5679, 0.5656, 0.5643, 0.5642, 0.5592, 0.5535, 0.5476, 0.5398, 0.5352, 0.5292, 0.5269, 0.5216, 0.517, 0.5117, 0.5071, 0.5102, 0.5168, 0.5242, 0.5212, 0.5136, 0.5156, 0.5272, 0.5383, 0.5427, 0.5432, 0.544, 0.5462, 0.5492, 0.5535, 0.5565, 0.5599, 0.5606, 0.5615, 0.5631, 0.5661, 0.5694, 0.5713, 0.5718, 0.5713, 0.5695, 0.5683, 0.5736, 0.5746, 0.5714, 0.568, 0.5645, 0.5605, 0.5523, 0.5384, 0.5192, 0.5041, 0.484, 0.4607, 0.4464, 0.4351, 0.4262, 0.4169, 0.4001, 0.3807, 0.3676, 0.3666, 0.3686, 0.362, 0.3549, 0.353, 0.3499, 0.3379, 0.3234, 0.312, 0.298, 0.2736, 0.2457, 0.2278, 0.2252, 0.2321, 0.2272, 0.2106, 0.2053, 0.2182, 0.2413, 0.2594, 0.2649, 0.269, 0.2765, 0.2901, 0.3095, 0.3312, 0.3465, 0.3542, 0.3597, 0.3654, 0.3785, 0.3971, 0.4098, 0.421, 0.4329, 0.4419, 0.4525, 0.47, 0.4806, 0.4809, 0.4864, 0.5004, 0.5115, 0.5196, 0.5307, 0.5395, 0.5466, 0.554, 0.5583, 0.5626, 0.5702, 0.5767, 0.5794, 0.5826, 0.5838, 0.5857, 0.5892, 0.5912, 0.5931, 0.5955, 0.5954, 0.594, 0.5935, 0.5952, 0.5944, 0.594, 0.5952, 0.594, 0.5945, 0.5945, 0.5932, 0.5937, 0.5923, 0.588, 0.5916, 0.599, 0.5966, 0.5918, 0.5909, 0.5907, 0.5913, 0.5923, 0.5918, 0.5912, 0.5901, 0.5887, 0.5878, 0.5858, 0.5839, 0.5835, 0.5819, 0.5795, 0.5762, 0.5742, 0.5761, 0.5712, 0.5676, 0.5717, 0.5585, 0.5534, 0.5614, 0.567, 0.5543, 0.5554, 0.5644, 0.5555, 0.5519, 0.5513, 0.5493, 0.5473, 0.5455, 0.5444, 0.5474, 0.551, 0.549, 0.5463, 0.5474, 0.5504, 0.5528, 0.5496, 0.5477, 0.5494, 0.5518, 0.5549, 0.55, 0.5498, 0.5492, 0.5469, 0.5443, 0.5377, 0.5292, 0.5248, 0.5277, 0.532, 0.5326, 0.5303, 0.5256, 0.5236, 0.5213, 0.5179, 0.5159, 0.5118, 0.5072, 0.5063, 0.5039, 0.4987, 0.4892, 0.4784, 0.4615, 0.4466, 0.443, 0.4419, 0.4422, 0.444, 0.4424, 0.4344, 0.4175, 0.4069, 0.4235, 0.4498, 0.4609, 0.4557, 0.448, 0.4455, 0.448, 0.4517, 0.4545, 0.46, 0.4658, 0.468, 0.4707, 0.471, 0.4708, 0.4701, 0.4678, 0.4648, 0.459, 0.4518, 0.4424, 0.4325, 0.4298, 0.4352, 0.4354, 0.4279, 0.42, 0.4128, 0.4124, 0.4137, 0.412, 0.4053, 0.398, 0.3971, 0.4002, 0.4036, 0.4059, 0.4049, 0.4091, 0.4165, 0.4166, 0.4236, 0.4343, 0.4308, 0.4289, 0.4368, 0.4447, 0.4484, 0.4514, 0.4522, 0.4541, 0.4655, 0.4729, 0.4745, 0.4764, 0.4808, 0.4883, 0.4958, 0.4992, 0.5, 0.5064, 0.5059, 0.5003, 0.5025, 0.508, 0.5079, 0.5076, 0.5112, 0.5196, 0.5269, 0.5302, 0.5304, 0.5311, 0.5336, 0.5358, 0.5364, 0.5358, 0.5377, 0.5418, 0.5435, 0.5471, 0.5482, 0.5477, 0.5473, 0.5488, 0.5545, 0.5637, 0.5612, 0.558, 0.568, 0.5703, 0.5659, 0.5642, 0.5679, 0.5703, 0.5727, 0.5699, 0.5694, 0.5716, 0.5713, 0.5717, 0.572, 0.5742, 0.5744, 0.575, 0.5757, 0.5755, 0.5727, 0.5741, 0.5766, 0.5767, 0.5757, 0.5766, 0.5801, 0.5837, 0.5825, 0.5829, 0.5851, 0.585, 0.5814, 0.5771, 0.574, 0.5758, 0.5738, 0.5705, 0.5633, 0.5537, 0.5393, 0.5264, 0.5176, 0.5087, 0.5014, 0.494, 0.484, 0.4677, 0.4449, 0.4328, 0.4399, 0.4558, 0.4663, 0.4599, 0.4402, 0.4173, 0.4053, 0.4121, 0.4436, 0.4808, 0.4901, 0.4799, 0.4648, 0.4575, 0.4713, 0.4897, 0.504, 0.5144, 0.5171, 0.5192, 0.5172, 0.5122, 0.5159, 0.5223, 0.5279, 0.5298, 0.53, 0.5368, 0.5205, 0.5117, 0.5113, 0.5105, 0.5043, 0.4956, 0.4895, 0.4764, 0.4567, 0.4488, 0.4529, 0.4244, 0.4171, 0.4223, 0.4105, 0.3939, 0.3759, 0.3559, 0.3381, 0.3237, 0.3159, 0.314, 0.3115, 0.3031, 0.2856, 0.2647, 0.2499, 0.2387, 0.2229, 0.206, 0.1998, 0.2037, 0.206, 0.2016, 0.2003, 0.206, 0.2109, 0.2126, 0.2145, 0.2256, 0.241, 0.2348, 0.2402, 0.2632, 0.2861, 0.3007, 0.3118, 0.3194, 0.3206, 0.3279, 0.3346, 0.3411, 0.3555, 0.3709, 0.3831, 0.3901, 0.3968, 0.4049, 0.4084, 0.4041, 0.3985, 0.4012, 0.4041, 0.4017, 0.3972, 0.3965, 0.3939, 0.3875, 0.3809, 0.3717, 0.3679, 0.3603, 0.3452, 0.3155, 0.2694, 0.23, 0.2084, 0.1785, 0.1336, 0.0939, 0.0727, 0.0639, 0.0561, 0.0503, 0.0528, 0.0656, 0.0814, 0.0898, 0.0875, 0.0805, 0.0674, 0.0468, 0.0325, 0.0331, 0.0412, 0.052, 0.0614, 0.0638, 0.0636, 0.064, 0.0613, 0.0582, 0.0591, 0.0644, 0.0704, 0.0701, 0.067, 0.0639, 0.0568, 0.051, 0.0522, 0.0562, 0.0564, 0.0529, 0.0535, 0.0647, 0.0837, 0.0985, 0.1014, 0.0986, 0.0949, 0.0893, 0.0884, 0.0964, 0.1092, 0.1244, 0.1348, 0.1427, 0.1566, 0.1727, 0.1804, 0.1872, 0.2059, 0.2339, 0.2545, 0.2585, 0.2638, 0.2861, 0.3098, 0.3193, 0.3311, 0.3494, 0.3651, 0.3769, 0.3842, 0.3919, 0.3994, 0.4166, 0.4295, 0.4328, 0.4317, 0.423, 0.4117, 0.4065, 0.4105, 0.4075, 0.3914, 0.3845, 0.3941, 0.4114, 0.4219, 0.4258, 0.4286, 0.4286, 0.4213, 0.4185, 0.4319, 0.4456, 0.457, 0.4678, 0.4695, 0.47, 0.4732, 0.475, 0.4782, 0.4862, 0.4936, 0.4954, 0.4983, 0.5043, 0.5059, 0.506, 0.5125, 0.5176, 0.519, 0.5227, 0.5275, 0.5278, 0.5231, 0.5213, 0.5247, 0.5299, 0.5341, 0.5319, 0.5288, 0.5319, 0.5368, 0.54, 0.5406, 0.54, 0.5387, 0.5405, 0.5426, 0.5421, 0.5393, 0.537, 0.5381, 0.5408, 0.5395, 0.5383, 0.5403, 0.5449, 0.5456, 0.5466, 0.5472, 0.5449, 0.5459, 0.5454, 0.5429, 0.5411, 0.5433, 0.5521, 0.5571, 0.5588, 0.5591, 0.55, 0.5425, 0.5399, 0.5466, 0.5556, 0.5471, 0.5372, 0.5396, 0.5299, 0.5292, 0.5411, 0.5362, 0.528, 0.5281, 0.5287, 0.5247, 0.5114, 0.4973, 0.4775, 0.4665, 0.4918, 0.5124, 0.5127, 0.5213, 0.5211, 0.5167, 0.5061, 0.4859, 0.4656, 0.4622, 0.4839, 0.5031, 0.5077, 0.5102, 0.5095, 0.5021, 0.5045, 0.5035, 0.4997, 0.4971, 0.4869, 0.4784, 0.4707, 0.4705, 0.4798, 0.4818, 0.4697, 0.458, 0.4494, 0.4427, 0.4373, 0.4274, 0.4158, 0.4044, 0.3875, 0.3659, 0.3442, 0.3216, 0.3025, 0.2985, 0.3019, 0.2913, 0.2722, 0.2602, 0.2534, 0.2493, 0.2478, 0.2542, 0.2667, 0.2751, 0.2702, 0.2558, 0.2468, 0.2461, 0.2353, 0.204, 0.1743, 0.1766, 0.2068, 0.2344, 0.2497, 0.2568, 0.2543, 0.2396, 0.2245, 0.223, 0.2351, 0.2515, 0.2639, 0.2678, 0.27, 0.274, 0.2738, 0.2675, 0.2657, 0.2701, 0.2711, 0.2643, 0.2599, 0.2557, 0.2426, 0.2279, 0.2178, 0.2061, 0.192, 0.1798, 0.1634, 0.1442, 0.1341, 0.1332, 0.1353, 0.1419, 0.1509, 0.1537, 0.1489, 0.1453, 0.1393, 0.1255, 0.1056, 0.0885, 0.0795, 0.0807, 0.0905, 0.1023, 0.109, 0.1111, 0.1115, 0.1148, 0.1176, 0.123, 0.1308, 0.1333, 0.1304, 0.1299, 0.1307, 0.1334, 0.1372, 0.1337, 0.1238, 0.1152, 0.122, 0.142, 0.1636, 0.1788, 0.1767, 0.1621, 0.1516, 0.1572, 0.175, 0.1871, 0.1871, 0.1821, 0.1775, 0.1736, 0.1711, 0.1731, 0.1804, 0.1886, 0.1921, 0.1814, 0.1596, 0.1395, 0.1275, 0.1224, 0.1166, 0.1133, 0.116, 0.117, 0.114, 0.1098, 0.1071, 0.1067, 0.108, 0.1068, 0.1022, 0.0992, 0.0969, 0.0952, 0.0944, 0.0968, 0.0988, 0.1056, 0.1136, 0.1163, 0.1184, 0.1189, 0.1202, 0.1232, 0.1283, 0.1377, 0.1519, 0.1646, 0.1761, 0.1856, 0.1968, 0.2106, 0.2232, 0.2273, 0.2229, 0.2278, 0.2339, 0.2274, 0.221, 0.2297, 0.2469, 0.2543, 0.2479, 0.2434, 0.2513, 0.2531, 0.253, 0.2585, 0.2595, 0.2614, 0.2721, 0.2803, 0.2794, 0.2751, 0.283, 0.2946, 0.3037, 0.3012, 0.3095, 0.3115, 0.3101, 0.3144, 0.3275, 0.3397, 0.3473, 0.3525, 0.3522, 0.3537, 0.3653, 0.3681, 0.3724, 0.3743, 0.3804, 0.3777, 0.3832, 0.3988, 0.4144, 0.4089, 0.4154, 0.4166, 0.4233, 0.4361, 0.4387, 0.4316, 0.4308, 0.4384, 0.455, 0.4609, 0.4572, 0.4528, 0.4525, 0.4692, 0.4792, 0.4898, 0.4827, 0.4878, 0.5029, 0.5137, 0.505, 0.5001, 0.4982, 0.5068, 0.5108, 0.5128, 0.5073, 0.5172, 0.5294, 0.5386, 0.5348, 0.5402, 0.5355])
SkyMapMod/planets.py ADDED
@@ -0,0 +1,131 @@
1
+ from astropy.coordinates import get_body, Galactic
2
+ from astropy.time import Time
3
+ import ephem
4
+ import numpy as np
5
+
6
+ from .zodiac import convolution, integral
7
+ from .albedo_of_planets import venus_alb_wl, venus_alb_rf, mars_alb_wl, mars_alb_rf, jupiter_alb_wl, jupiter_alb_rf, saturn_alb_wl, saturn_alb_rf
8
+ from .solar_spectrum import wavelenght_newguey2003, flux_newguey2003
9
+ from .band_V_data import wavelenght_band_V, trancparency_band_V
10
+
11
+
12
+ def coordinates_of_planet(body, date, time): #'mars', 'venus', 'saturn', 'jupyter'; 2023-11-03; 12:00:00; return galactic coordinates of planet
13
+ datetime = date + ' ' + time
14
+ time = Time(datetime)
15
+ planet = get_body(body, time, ephemeris='builtin').transform_to(Galactic)
16
+ l = planet.l.value
17
+ b = planet.b.value
18
+ return(l, b)
19
+
20
+
21
+ def magnitude_of_planet(body, date, time): #'mars', 'venus', 'saturn', 'jupyter'; 2023-11-03; 12:00:00';
22
+ if body == 'venus':
23
+ planet = ephem.Venus()
24
+ elif body == 'mars':
25
+ planet = ephem.Mars()
26
+ elif body == 'jupiter':
27
+ planet = ephem.Jupiter()
28
+ elif body == 'saturn':
29
+ planet = ephem.Saturn()
30
+
31
+ # Устанавливаем дату наблюдения
32
+ date = date.replace('-', '/')
33
+
34
+ planet.compute(date)
35
+ return planet.mag
36
+
37
+ def mag_to_phot(mag): #в полосе V
38
+ F0 = 10**6
39
+ F1 = F0 * 10 ** (- 0.4 * mag) # фот / (см^2 сек)
40
+ F1 = F1 * 10**4 / (4 * np.pi) # фот / (м^2 сек ср) -- НАДО ЛИ ДЕЛИТЬ НА 4 ПИ?
41
+ return F1
42
+
43
+
44
+ def venus_spectrum(date, time, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, albedo_wl = venus_alb_wl, albedo_rf = venus_alb_rf, V_wl = wavelenght_band_V, V_tr = trancparency_band_V):
45
+ body = 'venus'
46
+ planet_wl, planet_sp = convolution(albedo_wl, albedo_rf, Sun_sp_wl, Sun_sp_fx) #вот тут получили спектр планеты
47
+ planet_V_wl, planet_V_sp = convolution(planet_wl, planet_sp, V_wl, V_tr) #свернули с полосой V, фот / (сек ср м^2 нм)
48
+
49
+ norm = integral(planet_V_wl, planet_V_sp) # фот / (сек ср м^2)
50
+
51
+ #делаем так, чтобы эта норма была равна F1
52
+ mag = magnitude_of_planet(body, date, time)
53
+ F1 = mag_to_phot(mag)
54
+ A = F1 / norm
55
+ planet_sp = A * planet_sp
56
+
57
+ #########################################################################
58
+ # planet_V_wl, planet_V_sp = convolution(planet_wl, planet_sp, V_wl, V_tr)
59
+ # integr = integral(planet_V_wl, planet_V_sp)
60
+ # print('Техническая проверка на совпадение двух чисел:', integr, F1)
61
+ #########################################################################
62
+
63
+ return planet_wl, planet_sp
64
+
65
+ def mars_spectrum(date, time, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, albedo_wl = mars_alb_wl, albedo_rf = mars_alb_rf, V_wl = wavelenght_band_V, V_tr = trancparency_band_V):
66
+ body = 'mars'
67
+ planet_wl, planet_sp = convolution(albedo_wl, albedo_rf, Sun_sp_wl, Sun_sp_fx) #вот тут получили спектр планеты
68
+ planet_V_wl, planet_V_sp = convolution(planet_wl, planet_sp, V_wl, V_tr) #свернули с полосой V, фот / (сек ср м^2 нм)
69
+
70
+ norm = integral(planet_V_wl, planet_V_sp) # фот / (сек ср м^2)
71
+
72
+ #делаем так, чтобы эта норма была равна F1
73
+ mag = magnitude_of_planet(body, date, time)
74
+ F1 = mag_to_phot(mag)
75
+ A = F1 / norm
76
+ planet_sp = A * planet_sp
77
+
78
+ #########################################################################
79
+ # planet_V_wl, planet_V_sp = convolution(planet_wl, planet_sp, V_wl, V_tr)
80
+ # integr = integral(planet_V_wl, planet_V_sp)
81
+ # print('Техническая проверка на совпадение двух чисел:', integr, F1)
82
+ #########################################################################
83
+
84
+ return planet_wl, planet_sp
85
+
86
+ def jupiter_spectrum(date, time, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, albedo_wl = jupiter_alb_wl, albedo_rf = jupiter_alb_rf, V_wl = wavelenght_band_V, V_tr = trancparency_band_V):
87
+ body = 'jupiter'
88
+ planet_wl, planet_sp = convolution(albedo_wl, albedo_rf, Sun_sp_wl, Sun_sp_fx) #вот тут получили спектр планеты
89
+ planet_V_wl, planet_V_sp = convolution(planet_wl, planet_sp, V_wl, V_tr) #свернули с полосой V, фот / (сек ср м^2 нм)
90
+
91
+ norm = integral(planet_V_wl, planet_V_sp) # фот / (сек ср м^2)
92
+
93
+ #делаем так, чтобы эта норма была равна F1
94
+ mag = magnitude_of_planet(body, date, time)
95
+ F1 = mag_to_phot(mag)
96
+ A = F1 / norm
97
+ planet_sp = A * planet_sp
98
+
99
+ #########################################################################
100
+ # planet_V_wl, planet_V_sp = convolution(planet_wl, planet_sp, V_wl, V_tr)
101
+ # integr = integral(planet_V_wl, planet_V_sp)
102
+ # print('Техническая проверка на совпадение двух чисел:', integr, F1)
103
+ #########################################################################
104
+
105
+ return planet_wl, planet_sp
106
+
107
+ def saturn_spectrum(date, time, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, albedo_wl = saturn_alb_wl, albedo_rf = saturn_alb_rf, V_wl = wavelenght_band_V, V_tr = trancparency_band_V):
108
+ body = 'saturn'
109
+ planet_wl, planet_sp = convolution(albedo_wl, albedo_rf, Sun_sp_wl, Sun_sp_fx) #вот тут получили спектр планеты
110
+ planet_V_wl, planet_V_sp = convolution(planet_wl, planet_sp, V_wl, V_tr) #свернули с полосой V, фот / (сек ср м^2 нм)
111
+
112
+ norm = integral(planet_V_wl, planet_V_sp) # фот / (сек ср м^2)
113
+
114
+ #делаем так, чтобы эта норма была равна F1
115
+ mag = magnitude_of_planet(body, date, time)
116
+ F1 = mag_to_phot(mag)
117
+ A = F1 / norm
118
+ planet_sp = A * planet_sp
119
+
120
+ #########################################################################
121
+ # planet_V_wl, planet_V_sp = convolution(planet_wl, planet_sp, V_wl, V_tr)
122
+ # integr = integral(planet_V_wl, planet_V_sp)
123
+ # print('Техническая проверка на совпадение двух чисел:', integr, F1)
124
+ #########################################################################
125
+
126
+ return planet_wl, planet_sp
127
+
128
+ # Внутри основной функции: если планета попадает в поле зрения (через coordinates_of_planet(.....)),
129
+ # ищем ее спектр venus_spectrum(.....)
130
+
131
+
SkyMapMod/star_catalogues.py CHANGED
@@ -5,6 +5,7 @@ from .data.load_data import load_star_temperatures, load_star_brightness
5
5
  from.zodiac import black_body_phot
6
6
 
7
7
  #на вход подаем галактические координаты (l -- долгота, b -- широта) и массив длин волн, по которым строим спектр
8
+ #РАЗОБРАТЬСЯ С КОЭФФИЦИЕНТАМИ А -- НАДО ЛИ НОРМИРОВАТЬ НА C (КАЖЕТСЯ НАДО) А ЕЩЕ ВРОДЕ Я ЕГО ОТНОРМИРОВАЛА И ПОДГРУЗИЛА -- ПРОВЕРИТЬ
8
9
  def star_spectrum(l, b, lmbd):
9
10
  j = int((1800 - round(l*10)) + 3600) #пересчет из галактических к позиции в массиве
10
11
  i = int(900 + round(b*10))
@@ -13,5 +14,6 @@ def star_spectrum(l, b, lmbd):
13
14
  A = star_brightness[i][j]
14
15
  T = star_temperatures[i][j]
15
16
  print('Информация по звезде в этом пикселе: T = ', T, ', A = ', A, sep='')
17
+ print('Техническая строчка: i = ', i, ', j = ', j, sep='')
16
18
  spectrum = A * black_body_phot(lmbd, T)
17
19
  return spectrum
SkyMapMod/zodiac.py CHANGED
@@ -1,37 +1,36 @@
1
1
  import numpy as np
2
- import math
3
2
  import matplotlib.pyplot as plt
4
3
  import random
5
4
 
6
5
  #коэффициент пересчета из единиц S10 в количество фотонов от Солнца
7
6
  #вычислен на основе спектра АЧТ от звезды класса A0, солнечного спектра и факта, что 1 ед. S10 соответствует 100 фотонов от звезды класса A0
8
- k = 100.63910099562007
7
+ N_S10 = 100.63910099562007
9
8
 
10
9
  #переопределяю тригонометрические функции, чтобы работать с градусами
11
10
  def cos(x):
12
- return math.cos(x * math.pi / 180)
11
+ return np.cos(x * np.pi / 180)
13
12
 
14
13
  def sin(x):
15
- return math.sin(x * math.pi / 180)
14
+ return np.sin(x * np.pi / 180)
16
15
 
17
16
  def tan(x):
18
- return math.tan(x * math.pi / 180)
17
+ return np.tan(x * np.pi / 180)
19
18
 
20
19
  def acos(x):
21
- return math.acos(x) * 180 / math.pi
20
+ return np.arccos(x) * 180 / np.pi
22
21
 
23
22
  def asin(x):
24
- return math.asin(x) * 180 / math.pi
23
+ return np.arcsin(x) * 180 / np.pi
25
24
 
26
25
  def atan(x):
27
- return math.atan(x) * 180 / math.pi
26
+ return np.arctan(x) * 180 / np.pi
28
27
 
29
28
  #переопределяю прочие функции для красоты записи
30
29
  def exp(x):
31
- return math.exp(x)
30
+ return np.exp(x)
32
31
 
33
32
  def sqrt(x):
34
- return math.sqrt(x)
33
+ return np.sqrt(x)
35
34
 
36
35
  #функция Хевисайда (в контексте статьи определяется именно так)
37
36
  def u(x):
@@ -95,7 +94,7 @@ def D(lmbd, beta, lmbd_sun):
95
94
  # f = 60 * exp(- eps1 / 16) * (0.15 + 0.85 * exp(-6 * cos(0.6 * h)**4))
96
95
  f = 60 * exp(- eps1 / 16) * (0.85 + 0.15 * exp(-6 * cos(0.6 * h)**4))
97
96
  g = 25 * exp(- acos(cos(2 * lmbd1) * cos(0.7 * beta1)) / 10)
98
- D = A * ((f * math.exp(-d / 10) + g * math.exp(-d / 8)) * u(d) + (f + g) * (1 - u(d)))
97
+ D = A * ((f * np.exp(-d / 10) + g * np.exp(-d / 8)) * u(d) + (f + g) * (1 - u(d)))
99
98
  return D
100
99
 
101
100
  #Компонента E -- остаточная, не учитывалась
@@ -183,27 +182,12 @@ def integral(wl, spec):
183
182
  result += (spec[i] + spec[i+1]) / 2 * (wl[i+1] - wl[i])
184
183
  return result
185
184
 
186
- def N_S10(Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, T_A0 = 9800):
187
- h = 6.63e-34
188
- c = 3e8
189
- A0_sp_wl = Sun_sp_wl
190
- A0_sp_fx = black_body(A0_sp_wl, T_A0)
191
-
192
- A0_V_wl, A0_V_fx = convolution(A0_sp_wl, A0_sp_fx, V_wl, V_tr)
193
- A0_V_fx = A0_V_fx / (np.sum(A0_V_fx * A0_V_wl)) #надо домножать на A0_V_wl при нормировке.
194
185
 
195
- Sun_V_wl, Sun_V_fx = convolution(Sun_sp_wl, Sun_sp_fx, V_wl, V_tr)
196
- Sun_V_fx = Sun_V_fx / (np.sum(Sun_V_fx * Sun_V_wl))
197
-
198
- mean_A0 = np.sum(A0_V_fx * h * c / A0_V_wl * 1e9)
199
- mean_Sun = np.sum(Sun_V_fx * h * c / Sun_V_wl * 1e9)
200
- N_S10_to_Sun = 100 * mean_A0 / mean_Sun
201
- return (N_S10_to_Sun)
202
186
 
203
187
 
204
- def zodiacal_spectrum(lmbd, beta, lmbd_sun, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, T_A0 = 9800):
188
+ def zodiacal_spectrum(lmbd, beta, lmbd_sun, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, T_A0 = 9800, N_S10 = N_S10):
205
189
  S10 = zodiacal_light(lmbd, beta, lmbd_sun)
206
- N_S10_to_Sun = N_S10(Sun_sp_wl, Sun_sp_fx, V_wl, V_tr, T_A0)
190
+ N_S10_to_Sun = N_S10
207
191
  N_phot = round(S10 * N_S10_to_Sun)
208
192
  probabilities = Sun_sp_fx / np.sum(Sun_sp_fx) #или все-таки нормируем на интеграл?
209
193
  counts = np.zeros(Sun_sp_wl.shape[0])
@@ -213,30 +197,3 @@ def zodiacal_spectrum(lmbd, beta, lmbd_sun, Sun_sp_wl = wavelenght_newguey2003,
213
197
  return(Sun_sp_wl, counts)
214
198
 
215
199
 
216
- def N_S10_2(Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, T_A0 = 9800):
217
- h = 6.63e-34
218
- c = 3e8
219
- Sun_V_wl, Sun_V_fx = convolution(Sun_sp_wl, Sun_sp_fx, V_wl, V_tr) #сворачиваю с полосой V
220
-
221
- Sun_norm = integral(Sun_V_wl, Sun_V_fx)#интегрирую спектр, photons/(m^2 sec sr)
222
- Sun_V_fx_norm = Sun_V_fx / Sun_norm #спектр отнормирован на 1 фотон
223
-
224
- Sun_V_fx_norm_Vt = Sun_V_fx_norm * h * c / (Sun_V_wl * 10**(-9)) #перевела в мощности
225
- Sun = np.sum(Sun_V_fx_norm_Vt)
226
-
227
- A0_sp_wl = Sun_sp_wl #нм
228
- A0_sp_fx = black_body_phot(A0_sp_wl, T_A0) #фот/ (м^2 ср сек м)
229
-
230
- A0_V_wl, A0_V_fx = convolution(A0_sp_wl, A0_sp_fx, V_wl, V_tr)
231
- norm_S10 = integral(A0_V_wl * 10**(-9), A0_V_fx)
232
-
233
- A = 100/norm_S10
234
-
235
- A0_V_fx = A * A0_V_fx
236
- A0_norm = integral(A0_V_wl, A0_V_fx)
237
- A0_V_fx_norm = A0_V_fx / A0_norm #photons/(m^2 sec sr)
238
-
239
- A0_V_fx_norm_Vt = A0_V_fx_norm * h * c / (A0_V_wl * 10**(-9)) #Vt/(m^2 sec sr)
240
- A0 = np.sum(A0_V_fx_norm_Vt)
241
- N_S10_to_Sun = 100 * A0 / Sun
242
- return (N_S10_to_Sun, c)
{skymapmod-0.1.16.dist-info → skymapmod-0.1.18.dist-info}/METADATA RENAMED
@@ -1,6 +1,6 @@
1
1
  Metadata-Version: 2.4
2
2
  Name: SkyMapMod
3
- Version: 0.1.16
3
+ Version: 0.1.18
4
4
  Summary: for modeling
5
5
  Project-URL: Homepage, https://istina.msu.ru/profile/aminevann/
6
6
  Author-email: Anna Amineva <amineva.2219@physics.msu.ru>
{skymapmod-0.1.16.dist-info → skymapmod-0.1.18.dist-info}/RECORD RENAMED
@@ -1,18 +1,20 @@
1
- SkyMapMod/__init__.py,sha256=Q0vaTOPF5a5ZVCmJpg6VYOnHAjGpoheLKM_6SewIuQM,271
2
- SkyMapMod/airglow.py,sha256=4Uo4e2sk7Ux2Hc4MwgaZ7UrZwl2qe_ms0crCLOgXHKM,2293
1
+ SkyMapMod/__init__.py,sha256=I_JAagpYc5Z4R0OexeCJd754kRrccDk2gAtcezN53wI,327
2
+ SkyMapMod/airglow.py,sha256=g4NY79WfHya4-5jccCx1VTGJNvb0_DiotH9oYQsWd5Q,2281
3
3
  SkyMapMod/airglow_spectrum.py,sha256=4pY4QSGjIDvjs_jvxrowsZ9WLtIUwR0DesoCK6moWjs,127314
4
+ SkyMapMod/albedo_of_planets.py,sha256=vsGkMwoa2aByRfXML5JShtinIxq3V_rFPWMgvU4mfro,88531
4
5
  SkyMapMod/band_V_data.py,sha256=fOj9H2IdZuqErAxnoLq4A02UXfYMgLNBiLjwrzjfWTs,15223
5
6
  SkyMapMod/modtran_default_kp_transparency.py,sha256=QeLYQA9FNUsKp-eRtOMQFxfq0pdwCnqVIl0X6Vvf7GU,78306
7
+ SkyMapMod/planets.py,sha256=7xzna0cybGgDPOaNzx97hQ08jWYN9Ogbtk89nkId_RE,6702
6
8
  SkyMapMod/solar_radio_flux.py,sha256=I02O6Mltv3lfqxif3BHUNAzeElQ5AZax_2L09A3cXb8,842
7
9
  SkyMapMod/solar_spectrum.py,sha256=LjBPSTI6Yq86pwF7haUja7RjRYH4C8HZL8SjyI7zp6g,24433
8
- SkyMapMod/star_catalogues.py,sha256=IGzzyk7ITioNWbQJVQq8FBskiBwjNV3O-uYCCXD_ICg,849
10
+ SkyMapMod/star_catalogues.py,sha256=hl8-EMMCfJGz1MGBqdwqIY_d1-9G8j8cZGyWmE3OSBA,1166
9
11
  SkyMapMod/transparency.py,sha256=HiSbbVX80neh3n284f64sWnrbxpRzPGKvFyid9sVo7I,376
10
- SkyMapMod/zodiac.py,sha256=EepP00vdMbXu-we0aTOh5yAPrinNhm6IFYUiyeteWl4,10285
12
+ SkyMapMod/zodiac.py,sha256=pgOvKoMoz6P0JFyFDihhxTzIgbLFlM-Ey6D08RnzxYU,8280
11
13
  SkyMapMod/data/__init__.py,sha256=_292qgM72hPa2vxDNSaI9QvPunvi-aSXUsZki99q63o,83
12
14
  SkyMapMod/data/load_data.py,sha256=p-Nq5St713HTTDunSWu9ANrF9V0Slaak17IhD0n_bSY,931
13
15
  SkyMapMod/data/star_brightness.npy,sha256=OVrYGTl3I2wbWAMmFnDkunkRyPLL2-g-nD_JX_0FBdQ,51840128
14
16
  SkyMapMod/data/star_temperatures.npy,sha256=Yjp252bzqmKQz00K0ADTkstjzGFulRMgyW1b8kaDOPE,51840128
15
- skymapmod-0.1.16.dist-info/METADATA,sha256=TWmgczXOq_bUMcAz-1TM24MmZtxgxXmcBHntGqV9kyk,469
16
- skymapmod-0.1.16.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
17
- skymapmod-0.1.16.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
18
- skymapmod-0.1.16.dist-info/RECORD,,
17
+ skymapmod-0.1.18.dist-info/METADATA,sha256=gsBcBZwbqwerABc92sEPE09CbeBIn8NUQRJTYdLo6Z4,469
18
+ skymapmod-0.1.18.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
19
+ skymapmod-0.1.18.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
20
+ skymapmod-0.1.18.dist-info/RECORD,,