From 20528e1c0df4396203e6e26431fe20ce3070be24 Mon Sep 17 00:00:00 2001 From: Christian Tai Udovicic Date: Thu, 6 Jul 2023 10:07:02 -0700 Subject: [PATCH] Add diviner data and clean up docs --- .gitignore | 5 +- data/diviner/div_t2r.txt | 441 ++++++++++++++++++++++ data/diviner/diviner_filter_functions.hdf | Bin 0 -> 18264 bytes roughness/config.py | 2 +- roughness/diviner.py | 93 ++++- 5 files changed, 525 insertions(+), 16 deletions(-) create mode 100755 data/diviner/div_t2r.txt create mode 100644 data/diviner/diviner_filter_functions.hdf diff --git a/.gitignore b/.gitignore index 0d721e7..577b930 100644 --- a/.gitignore +++ b/.gitignore @@ -143,6 +143,9 @@ tmp/ scratch/ # Data folder except !items -/data/* +/data/** +!data/**/ !data/__init__.py !data/m3/m3_wls.json +!data/diviner/div_t2r.txt +!data/diviner/diviner_filter_functions.hdf \ No newline at end of file diff --git a/data/diviner/div_t2r.txt b/data/diviner/div_t2r.txt new file mode 100755 index 0000000..a688b80 --- /dev/null +++ b/data/diviner/div_t2r.txt @@ -0,0 +1,441 @@ +T c3 c4 c5 c6 c7 c8 c9 +0.00000 0 0 0 0 0 0 0 +1.00000 0 0 0 8.3784E-237 1.2191E-115 2.79231E-42 2.38004E-17 +2.00000 0 0 0 2.3241E-119 1.57451E-59 9.22189E-24 3.90897E-11 +3.00000 4.7847E-213 4.4446E-213 5.5714E-213 4.21515E-80 1.1124E-40 1.66728E-17 7.58855E-09 +4.00000 2.2946E-160 2.2051E-160 2.6208E-160 1.97657E-60 3.37424E-31 2.56898E-14 1.32797E-07 +5.00000 9.9113E-129 9.6066E-129 1.1085E-128 1.32683E-48 1.75272E-25 2.33707E-12 8.43872E-07 +6.00000 1.269E-107 1.2235E-107 1.3987E-107 1.0625E-40 1.17243E-21 5.22945E-11 3.16349E-06 +7.00000 1.55348E-92 1.48061E-92 1.70227E-92 4.88771E-35 6.48449E-19 5.41408E-10 8.68285E-06 +8.00000 3.27283E-81 3.07719E-81 3.5967E-81 8.97348E-31 7.54432E-17 3.57093E-09 1.94981E-05 +9.00000 2.12699E-72 1.97356E-72 2.3623E-72 1.92854E-27 3.11043E-15 1.76081E-08 3.81466E-05 +10.00000 2.41138E-65 2.21177E-65 2.72334E-65 9.24032E-25 6.22729E-14 6.99951E-08 6.75525E-05 +11.00000 1.43698E-59 1.3061E-59 1.65799E-59 1.48536E-22 7.40899E-13 2.33278E-07 0.000110967 +12.00000 9.39327E-55 8.48431E-55 1.11093E-54 1.05103E-20 5.99709E-12 6.6967E-07 0.000171894 +13.00000 1.12107E-50 1.00926E-50 1.36228E-50 3.94921E-19 3.62569E-11 1.69362E-06 0.000253995 +14.00000 3.51566E-47 3.16416E-47 4.39742E-47 9.00835E-18 1.74833E-10 3.84771E-06 0.000360986 +15.00000 3.78388E-44 3.41457E-44 4.88097E-44 1.37613E-16 7.04579E-10 7.98413E-06 0.000496524 +16.00000 1.70786E-41 1.54963E-41 2.27796E-41 1.51558E-15 2.45423E-09 1.53431E-05 0.000664106 +17.00000 3.76816E-39 3.44743E-39 5.2185E-39 1.27336E-14 7.57283E-09 2.76204E-05 0.000866978 +18.00000 4.57912E-37 4.23641E-37 6.62603E-37 8.52974E-14 2.10839E-08 4.70192E-05 0.001108065 +19.00000 3.36774E-35 3.16086E-35 5.1374E-35 4.71737E-13 5.37222E-08 7.62813E-05 0.001389922 +20.00000 1.61647E-33 1.54511E-33 2.63009E-33 2.21528E-12 1.26696E-07 0.000118697 0.001714709 +21.00000 5.38313E-32 5.26529E-32 9.47303E-32 9.03689E-12 2.79112E-07 0.000178097 0.00208418 +22.00000 1.30771E-30 1.31672E-30 2.52722E-30 3.26289E-11 5.78828E-07 0.00025882 0.002499691 +23.00000 2.41568E-29 2.52234E-29 5.20529E-29 1.05905E-10 1.13743E-06 0.000365671 0.002962217 +24.00000 3.51331E-28 3.83728E-28 8.55692E-28 3.13042E-10 2.12985E-06 0.000503859 0.003472379 +25.00000 4.1429E-27 4.77926E-27 1.1529E-26 8.51935E-10 3.81907E-06 0.000678933 0.004030476 +26.00000 4.06155E-26 4.99917E-26 1.30062E-25 2.15459E-09 6.58576E-06 0.000896709 0.004636521 +27.00000 3.38259E-25 4.48597E-25 1.25028E-24 5.10411E-09 1.09629E-05 0.001163192 0.005290275 +28.00000 2.43843E-24 3.5142E-24 1.03933E-23 1.14035E-08 1.76745E-05 0.001484498 0.00599128 +29.00000 1.54659E-23 2.4378E-23 7.56651E-23 2.41704E-08 2.76784E-05 0.001866788 0.006738897 +30.00000 8.75519E-23 1.51498E-22 4.87832E-22 4.88519E-08 4.22115E-05 0.00231619 0.007532336 +31.00000 4.47958E-22 8.5149E-22 2.8132E-21 9.45749E-08 6.28364E-05 0.002838743 0.008370679 +32.00000 2.09441E-21 4.36267E-21 1.46413E-20 1.76065E-07 9.14891E-05 0.003440334 0.009252912 +33.00000 9.03399E-21 2.05144E-20 6.93326E-20 3.16286E-07 0.000130526 0.004126654 0.010177944 +34.00000 3.6245E-20 8.90588E-20 3.00937E-19 5.4995E-07 0.000182768 0.004903152 0.011144627 +35.00000 1.36203E-19 3.5885E-19 1.20534E-18 9.2808E-07 0.000251544 0.005775003 0.012151774 +36.00000 4.82195E-19 1.34858E-18 4.48222E-18 1.52377E-06 0.000340729 0.006747076 0.013198173 +37.00000 1.61607E-18 4.74808E-18 1.55609E-17 2.4393E-06 0.00045478 0.007823918 0.014282597 +38.00000 5.14799E-18 1.5727E-17 5.06914E-17 3.81476E-06 0.00059876 0.009009737 0.015403821 +39.00000 1.56385E-17 4.91984E-17 1.55665E-16 5.83828E-06 0.000778367 0.010308392 0.016560621 +40.00000 4.54309E-17 1.45883E-16 4.52509E-16 8.75796E-06 0.000999948 0.011723388 0.017751791 +41.00000 1.26516E-16 4.11413E-16 1.24999E-15 1.28955E-05 0.00127051 0.013257881 0.018976141 +42.00000 3.38456E-16 1.10696E-15 3.29269E-15 1.8662E-05 0.001597728 0.014914675 0.020232505 +43.00000 8.71466E-16 2.84999E-15 8.29747E-15 2.65746E-05 0.001989939 0.016696234 0.021519745 +44.00000 2.16352E-15 7.04027E-15 2.00615E-14 3.72761E-05 0.002456137 0.01860469 0.022836751 +45.00000 5.18751E-15 1.67289E-14 4.66623E-14 5.15556E-05 0.00300596 0.020641854 0.024182447 +46.00000 1.20319E-14 3.83262E-14 1.04671E-13 7.03703E-05 0.003649674 0.022809229 0.025555789 +47.00000 2.70366E-14 8.48436E-14 2.26945E-13 9.48697E-05 0.004398142 0.025108024 0.026955766 +48.00000 5.89446E-14 1.8185E-13 4.76603E-13 0.00012642 0.005262804 0.027539177 0.028381404 +49.00000 1.2486E-13 3.78089E-13 9.71329E-13 0.000166632 0.006255643 0.030103355 0.029831763 +50.00000 2.57323E-13 7.63861E-13 1.9245E-12 0.000217384 0.007389145 0.032800988 0.031305934 +51.00000 5.16618E-13 1.50203E-12 3.71296E-12 0.000280856 0.008676268 0.035632269 0.032803047 +52.00000 1.01166E-12 2.87896E-12 6.98604E-12 0.000359551 0.010130397 0.038597189 0.034322264 +53.00000 1.9346E-12 5.38638E-12 1.28368E-11 0.000456328 0.011765303 0.041695529 0.035862778 +54.00000 3.6168E-12 9.84974E-12 2.30657E-11 0.00057443 0.013595097 0.044926902 0.037423816 +55.00000 6.6176E-12 1.76257E-11 4.05769E-11 0.000717513 0.015634187 0.048290741 0.039004637 +56.00000 1.18621E-11 3.08995E-11 6.99656E-11 0.000889671 0.017897236 0.05178634 0.040604527 +57.00000 2.08508E-11 5.31251E-11 1.18368E-10 0.001095468 0.020399108 0.055412843 0.042222806 +58.00000 3.59736E-11 8.96641E-11 1.96678E-10 0.001339964 0.023154833 0.059169279 0.043858817 +59.00000 6.09704E-11 1.48699E-10 3.21248E-10 0.001628739 0.026179553 0.063054559 0.045511934 +60.00000 1.01598E-10 2.42517E-10 5.16246E-10 0.001967922 0.029488489 0.067067491 0.047181557 +61.00000 1.66576E-10 3.89289E-10 8.16867E-10 0.00236421 0.033096886 0.071206799 0.048867109 +62.00000 2.6892E-10 6.15499E-10 1.27364E-09 0.002824896 0.037019987 0.075471125 0.050568039 +63.00000 4.27776E-10 9.59214E-10 1.95815E-09 0.003357889 0.041272976 0.07985904 0.052283819 +64.00000 6.70927E-10 1.47444E-09 2.97051E-09 0.003971733 0.045870962 0.084369052 0.054013942 +65.00000 1.03817E-09 2.23683E-09 4.44907E-09 0.004675623 0.050828921 0.088999627 0.055757924 +66.00000 1.5858E-09 3.35112E-09 6.58284E-09 0.005479428 0.056161678 0.09374917 0.057515302 +67.00000 2.39252E-09 4.96066E-09 9.62714E-09 0.006393701 0.061883871 0.098616061 0.059285626 +68.00000 3.56713E-09 7.25955E-09 1.39233E-08 0.007429693 0.068009919 0.103598647 0.061068476 +69.00000 5.25838E-09 1.05078E-08 1.99233E-08 0.008599367 0.074553996 0.108695246 0.06286344 +70.00000 7.66759E-09 1.50506E-08 2.82194E-08 0.009915405 0.081530009 0.113904157 0.06467013 +71.00000 1.10645E-08 2.13414E-08 3.95813E-08 0.011391217 0.088951568 0.119223666 0.066488161 +72.00000 1.58072E-08 2.99709E-08 5.50001E-08 0.013040946 0.096831974 0.124652052 0.06831719 +73.00000 2.23666E-08 4.17019E-08 7.57416E-08 0.014879474 0.105184194 0.130187581 0.070156856 +74.00000 3.13568E-08 5.75113E-08 1.0341E-07 0.016922425 0.114020844 0.135828527 0.072006839 +75.00000 4.35722E-08 7.86402E-08 1.40021E-07 0.019186162 0.123354176 0.141573154 0.073866818 +76.00000 6.00314E-08 1.06654E-07 1.88092E-07 0.02168779 0.133196066 0.147419739 0.07573649 +77.00000 8.20316E-08 1.4351E-07 2.50742E-07 0.024445151 0.143557999 0.153366566 0.077615557 +78.00000 1.11212E-07 1.91646E-07 3.31812E-07 0.02747682 0.154451061 0.159411925 0.079503749 +79.00000 1.49629E-07 2.54067E-07 4.35998E-07 0.030802101 0.165885956 0.165554118 0.081400793 +80.00000 1.99848E-07 3.34462E-07 5.6901E-07 0.034441015 0.177872925 0.171791467 0.083306424 +81.00000 2.65041E-07 4.3733E-07 7.37747E-07 0.038414297 0.190421835 0.178122305 0.085220409 +82.00000 3.49116E-07 5.68121E-07 9.50498E-07 0.042743379 0.203542124 0.184544981 0.087142494 +83.00000 4.56851E-07 7.33404E-07 1.21716E-06 0.047450386 0.217242816 0.191057866 0.089072458 +84.00000 5.94055E-07 9.41055E-07 1.54952E-06 0.052558116 0.231532464 0.197659351 0.091010087 +85.00000 7.67755E-07 1.20046E-06 1.96148E-06 0.058090024 0.246419284 0.204347848 0.092955153 +86.00000 9.86399E-07 1.52275E-06 2.46941E-06 0.064070225 0.261910974 0.211121789 0.094907473 +87.00000 1.2601E-06 1.92107E-06 3.0925E-06 0.070523443 0.278014893 0.217979626 0.096866835 +88.00000 1.60089E-06 2.41085E-06 3.8531E-06 0.077475033 0.294737962 0.224919841 0.098833061 +89.00000 2.02304E-06 3.01016E-06 4.77716E-06 0.084950933 0.312086653 0.231940934 0.100805963 +90.00000 2.54336E-06 3.73999E-06 5.89466E-06 0.09297765 0.330067077 0.239041434 0.102785374 +91.00000 3.18159E-06 4.62473E-06 7.24011E-06 0.101582265 0.348684945 0.246219883 0.104771118 +92.00000 3.96081E-06 5.6925E-06 8.85309E-06 0.110792365 0.367945526 0.253474869 0.106763044 +93.00000 4.9079E-06 6.97564E-06 1.07788E-05 0.120636068 0.387853752 0.260804972 0.108760988 +94.00000 6.05398E-06 8.51123E-06 1.30686E-05 0.131141977 0.408414159 0.268208849 0.1107648 +95.00000 7.43502E-06 1.03416E-05 1.57809E-05 0.142339152 0.429630886 0.275685108 0.112774343 +96.00000 9.0924E-06 1.25148E-05 1.89815E-05 0.15425711 0.451507751 0.283232453 0.114789468 +97.00000 1.10735E-05 1.50853E-05 2.27447E-05 0.166925777 0.474048175 0.290849567 0.116810048 +98.00000 1.34326E-05 1.81149E-05 2.71537E-05 0.180375481 0.497255257 0.298535202 0.118835954 +99.00000 1.62311E-05 2.16727E-05 3.23017E-05 0.194636918 0.521131748 0.30628804 0.120867057 +100.00000 1.95392E-05 2.58368E-05 3.82928E-05 0.209741132 0.545680059 0.314106932 0.122903234 +101.00000 2.34359E-05 3.06942E-05 4.52426E-05 0.225719488 0.570902298 0.321990609 0.124944373 +102.00000 2.80102E-05 3.6342E-05 5.32796E-05 0.242603654 0.596800236 0.329937922 0.126990362 +103.00000 3.33624E-05 4.28886E-05 6.25458E-05 0.260425571 0.623375372 0.337947722 0.129041086 +104.00000 3.96048E-05 5.04537E-05 7.31978E-05 0.279217418 0.650628879 0.346018839 0.131096453 +105.00000 4.68629E-05 5.91704E-05 8.54081E-05 0.299011604 0.67856169 0.354150196 0.133156343 +106.00000 5.52763E-05 6.9185E-05 9.93661E-05 0.319840751 0.707174401 0.362340669 0.135220675 +107.00000 6.50005E-05 8.06591E-05 0.000115279 0.341737632 0.736467393 0.370589234 0.137289342 +108.00000 7.62077E-05 9.37698E-05 0.000133374 0.364735186 0.766440789 0.378894815 0.139362262 +109.00000 8.90884E-05 0.000108711 0.000153897 0.388866473 0.79709442 0.38725641 0.141439333 +110.00000 0.000103853 0.000125696 0.000177117 0.414164656 0.828427935 0.395672969 0.143520485 +111.00000 0.000120731 0.000144956 0.000203327 0.440662984 0.860440723 0.404143564 0.145605623 +112.00000 0.000139979 0.000166742 0.000232841 0.468394751 0.893131945 0.412667218 0.147694665 +113.00000 0.000161872 0.000191329 0.000266001 0.497393282 0.926500577 0.421242978 0.149787543 +114.00000 0.000186716 0.000219013 0.000303176 0.52769196 0.960545383 0.429869915 0.151884171 +115.00000 0.000214842 0.000250115 0.000344762 0.559324072 0.99526491 0.43854715 0.153984486 +116.00000 0.00024661 0.000284982 0.000391186 0.592322892 1.030657556 0.44727378 0.156088408 +117.00000 0.000282413 0.000323988 0.000442905 0.626721747 1.066721519 0.456048949 0.158195868 +118.00000 0.000322675 0.000367533 0.000500409 0.662553766 1.103454814 0.464871778 0.160306803 +119.00000 0.000367856 0.000416051 0.000564222 0.699852028 1.140855317 0.473741485 0.162421144 +120.00000 0.000418453 0.000470002 0.000634903 0.738649468 1.178920753 0.48265724 0.164538835 +121.00000 0.000475002 0.000529883 0.000713047 0.778978972 1.217648664 0.491618237 0.166659806 +122.00000 0.000538078 0.000596221 0.00079929 0.820873179 1.257036498 0.500623693 0.168784001 +123.00000 0.000608302 0.000669581 0.000894304 0.864364535 1.297081538 0.509672852 0.170911361 +124.00000 0.000686338 0.000750564 0.000998804 0.909485485 1.337780959 0.51876498 0.173041831 +125.00000 0.000772898 0.000839808 0.001113546 0.956268035 1.379131782 0.527899321 0.175175357 +126.00000 0.000868743 0.000937992 0.001239332 1.004744093 1.421130968 0.537075166 0.177311884 +127.00000 0.000974685 0.001045835 0.001377008 1.054945322 1.46377531 0.546291806 0.179451359 +128.00000 0.001091591 0.001164099 0.001527464 1.106903141 1.507061556 0.555548583 0.181593736 +129.00000 0.001220383 0.00129359 0.001691641 1.160648677 1.550986299 0.564844787 0.183738965 +130.00000 0.001362043 0.001435157 0.001870529 1.216212761 1.595546134 0.574179783 0.185886993 +131.00000 0.00151761 0.001589699 0.002065165 1.273625935 1.640737457 0.583552912 0.188037775 +132.00000 0.001688189 0.00175816 0.002276641 1.332918492 1.686556756 0.592963563 0.190191271 +133.00000 0.001874948 0.001941534 0.002506101 1.394120338 1.733000059 0.602411077 0.192347429 +134.00000 0.002079123 0.002140866 0.002754741 1.457261034 1.780063855 0.611894892 0.194506215 +135.00000 0.002302019 0.002357253 0.003023815 1.522369851 1.827744325 0.621414372 0.196667577 +136.00000 0.002545011 0.002591843 0.003314631 1.589475714 1.876037346 0.630968955 0.198831481 +137.00000 0.002809551 0.002845841 0.003628554 1.658607112 1.924939253 0.640558081 0.200997887 +138.00000 0.003097164 0.003120506 0.003967009 1.72979219 1.974445769 0.650181187 0.203166755 +139.00000 0.003409455 0.003417153 0.004331479 1.803058751 2.024552924 0.659837686 0.205338045 +140.00000 0.003748108 0.003737156 0.004723506 1.878434207 2.075256748 0.66952709 0.207511722 +141.00000 0.004114888 0.004081948 0.005144695 1.955945482 2.126553116 0.679248863 0.209687746 +142.00000 0.004511648 0.004453018 0.005596713 2.035619207 2.178437601 0.68900249 0.211866089 +143.00000 0.004940323 0.004851922 0.006081286 2.117481475 2.230906383 0.698787434 0.21404671 +144.00000 0.005402939 0.005280272 0.006600207 2.201558037 2.283954882 0.708603231 0.216229576 +145.00000 0.005901611 0.005739747 0.007155332 2.287874252 2.337579128 0.718449393 0.218414658 +146.00000 0.006438547 0.006232086 0.00774858 2.376454944 2.391774537 0.728325431 0.22060192 +147.00000 0.007016048 0.006759095 0.008381935 2.467324592 2.446537141 0.738230905 0.22279133 +148.00000 0.007636511 0.007322642 0.00905745 2.56050719 2.501862509 0.748165327 0.224982858 +149.00000 0.008302429 0.007924665 0.009777239 2.656026289 2.557746214 0.758128257 0.227176481 +150.00000 0.009016396 0.008567163 0.010543486 2.753905001 2.614183978 0.768119256 0.229372154 +151.00000 0.009781105 0.009252208 0.01135844 2.854165951 2.671171374 0.778137907 0.231569865 +152.00000 0.010599353 0.009981934 0.012224418 2.956831322 2.728704277 0.788183796 0.233769575 +153.00000 0.011474038 0.010758546 0.013143804 3.061922858 2.786778107 0.798256459 0.235971259 +154.00000 0.012408164 0.011584317 0.014119047 3.169461864 2.84538874 0.808355553 0.23817489 +155.00000 0.013404842 0.012461589 0.015152666 3.279469107 2.904531442 0.818480615 0.240380439 +156.00000 0.014467289 0.013392772 0.016247248 3.391964963 2.964202241 0.828631302 0.242587888 +157.00000 0.015598832 0.014380346 0.017405444 3.506969319 3.024396557 0.838807224 0.244797199 +158.00000 0.016802908 0.01542686 0.018629975 3.624501525 3.085110113 0.849007989 0.247008358 +159.00000 0.018083061 0.016534936 0.01992363 3.744580638 3.146338481 0.859233232 0.249221338 +160.00000 0.019442952 0.017707261 0.021289262 3.86722503 3.208077385 0.869482586 0.251436116 +161.00000 0.02088635 0.018946595 0.022729793 3.992452831 3.270322547 0.879755733 0.253652663 +162.00000 0.022417142 0.020255767 0.024248211 4.120281484 3.333069694 0.890052282 0.255870961 +163.00000 0.024039323 0.021637677 0.02584757 4.250728139 3.396314241 0.900371892 0.258090983 +164.00000 0.025757008 0.023095292 0.02753099 4.383809412 3.460052373 0.910714245 0.260312711 +165.00000 0.027574425 0.024631652 0.029301656 4.519541426 3.524279354 0.921078999 0.262536123 +166.00000 0.029495916 0.026249863 0.031162819 4.657939964 3.588991214 0.931465861 0.264761194 +167.00000 0.031525941 0.027953103 0.03311779 4.799020175 3.654183524 0.941874464 0.266987903 +168.00000 0.033669077 0.029744616 0.03516995 4.942797059 3.719852314 0.952304539 0.269216238 +169.00000 0.035930016 0.031627716 0.037322741 5.08928464 3.785993154 0.96275577 0.27144617 +170.00000 0.038313566 0.033605787 0.039579665 5.238496843 3.852601921 0.973227838 0.273677683 +171.00000 0.040824652 0.035682276 0.041944282 5.390447352 3.919674645 0.98372045 0.275910752 +172.00000 0.043468319 0.037860696 0.044420229 5.54514887 3.987206897 0.994233339 0.278145368 +173.00000 0.04624972 0.040144637 0.04701118 5.702614102 4.055194858 1.004766209 0.280381506 +174.00000 0.049174133 0.04253774 0.04972089 5.862855264 4.123634101 1.015318793 0.282619144 +175.00000 0.052246952 0.045043723 0.052553155 6.025883593 4.192520959 1.025890798 0.284858271 +176.00000 0.055473675 0.04766636 0.055511841 6.191710819 4.261851004 1.036481955 0.287098863 +177.00000 0.058859934 0.050409495 0.058600863 6.360347689 4.331620571 1.047091995 0.28934091 +178.00000 0.062411455 0.053277031 0.061824195 6.531803978 4.401825536 1.057720674 0.291584388 +179.00000 0.066134099 0.056272933 0.065185862 6.706090434 4.472461776 1.068367748 0.293829286 +180.00000 0.070033822 0.059401231 0.068689948 6.883216343 4.543525626 1.079032947 0.29607558 +181.00000 0.074116706 0.062666012 0.072340582 7.063190499 4.615012962 1.089716004 0.298323255 +182.00000 0.078388933 0.066071424 0.076141951 7.246022187 4.686919966 1.100416723 0.300572303 +183.00000 0.08285681 0.069621669 0.08009829 7.431719224 4.759242974 1.111134834 0.302822697 +184.00000 0.087526741 0.07332101 0.084213878 7.620289429 4.831977862 1.121870119 0.305074431 +185.00000 0.092405247 0.07717377 0.088493047 7.811741109 4.905121117 1.132622357 0.307327487 +186.00000 0.09749896 0.08118432 0.092940177 8.006081104 4.978668768 1.14339128 0.309581843 +187.00000 0.102814607 0.08535709 0.097559687 8.203315767 5.052617304 1.154176692 0.311837498 +188.00000 0.108359031 0.08969656 0.102356045 8.40345145 5.126962753 1.164978398 0.314094424 +189.00000 0.11413918 0.094207265 0.107333758 8.606494505 5.201701603 1.175796153 0.31635261 +190.00000 0.120162097 0.09889379 0.112497379 8.812450796 5.276830188 1.186629763 0.318612044 +191.00000 0.126434931 0.103760764 0.117851494 9.02132472 5.352344997 1.197479006 0.320872714 +192.00000 0.132964934 0.108812879 0.123400735 9.233122141 5.428242212 1.208343689 0.323134604 +193.00000 0.139759459 0.114054855 0.129149771 9.447846969 5.504518471 1.21922359 0.325397702 +194.00000 0.146825945 0.119491469 0.135103301 9.6655036 5.581170111 1.230118564 0.327661991 +195.00000 0.154171937 0.125127543 0.141266067 9.886095944 5.658193924 1.241028365 0.329927459 +196.00000 0.161805072 0.130967941 0.147642831 10.10962742 5.735586093 1.2519528 0.332194095 +197.00000 0.169733082 0.137017564 0.154238405 10.33610145 5.813343256 1.262891745 0.334461888 +198.00000 0.177963785 0.143281361 0.161057618 10.56552048 5.891462209 1.273844956 0.33673082 +199.00000 0.186505086 0.149764319 0.168105336 10.79788695 5.969939437 1.284812262 0.339000885 +200.00000 0.19536499 0.156471454 0.175386445 11.03320379 6.048771581 1.295793493 0.341272061 +201.00000 0.204551575 0.163407831 0.182905867 11.27147198 6.127955434 1.306788477 0.343544348 +202.00000 0.214073007 0.170578543 0.19066854 11.51269396 6.207487636 1.317797043 0.345817729 +203.00000 0.223937532 0.177988721 0.198679436 11.75687023 6.287364981 1.32881902 0.348092191 +204.00000 0.234153482 0.185643523 0.206943537 12.00400224 6.36758426 1.339854238 0.350367724 +205.00000 0.244729263 0.193548147 0.215465855 12.25409097 6.448142115 1.350902524 0.352644316 +206.00000 0.255673351 0.201707809 0.224251418 12.50713594 6.529035644 1.361963709 0.354921956 +207.00000 0.266994309 0.210127764 0.233305276 12.76313813 6.610261639 1.37303767 0.357200632 +208.00000 0.278700757 0.218813288 0.24263249 13.02209654 6.691816893 1.384124237 0.359480333 +209.00000 0.290801397 0.227769684 0.25223814 13.2840112 6.773698353 1.395223237 0.361761053 +210.00000 0.303304993 0.237002284 0.262127321 13.5488816 6.855903116 1.406334525 0.364042775 +211.00000 0.316220383 0.246516429 0.272305138 13.81670677 6.938427975 1.417457978 0.366325499 +212.00000 0.32955645 0.2563175 0.282776709 14.08748524 7.021270182 1.428593401 0.368609201 +213.00000 0.34332216 0.266410887 0.293547155 14.36121556 7.104426683 1.439740696 0.370893877 +214.00000 0.357526521 0.276801998 0.304621614 14.63789575 7.187894422 1.450899693 0.373179521 +215.00000 0.372178613 0.28749625 0.316005245 14.91752437 7.271670651 1.462070268 0.375466121 +216.00000 0.387287563 0.298499124 0.327703205 15.20009895 7.355752468 1.473252276 0.377753665 +217.00000 0.402862546 0.309816023 0.339720597 15.48561707 7.440137124 1.484445571 0.380042143 +218.00000 0.418912798 0.321452487 0.352062636 15.77407579 7.524821718 1.495650053 0.382331548 +219.00000 0.435447593 0.333413912 0.364734452 16.06547315 7.609803348 1.506865553 0.384621869 +220.00000 0.452476301 0.345705868 0.377741204 16.35980476 7.695079418 1.518091947 0.386913094 +221.00000 0.470008219 0.358333836 0.391088081 16.65706866 7.780647332 1.529329115 0.389205219 +222.00000 0.488052812 0.371303327 0.40478021 16.95725997 7.866504035 1.540576933 0.391498236 +223.00000 0.5066195 0.384619825 0.418822781 17.26037576 7.952647084 1.55183528 0.393792135 +224.00000 0.52571779 0.398288898 0.433220894 17.5664121 8.039073882 1.563104034 0.396086905 +225.00000 0.545357232 0.412316028 0.447979765 17.87536511 8.125781681 1.574383072 0.398382539 +226.00000 0.565547332 0.426706813 0.463104467 18.18722939 8.212767885 1.585672272 0.400679026 +227.00000 0.586297724 0.441466707 0.478600185 18.50200154 8.300030049 1.596971513 0.40297636 +228.00000 0.607618044 0.456601278 0.494472051 18.81967617 8.387565578 1.608280696 0.405274531 +229.00000 0.629517882 0.472116065 0.510725135 19.14024889 8.475371876 1.619599699 0.407573537 +230.00000 0.652006917 0.488016582 0.527364566 19.46371481 8.563446652 1.630928424 0.409873361 +231.00000 0.675094867 0.504308338 0.544395415 19.79006807 8.651787309 1.642266726 0.412173999 +232.00000 0.698791453 0.520996873 0.561822812 20.11930379 8.740391405 1.65361453 0.414475444 +233.00000 0.723106394 0.538087699 0.579651798 20.45141609 8.829256495 1.664971714 0.416777691 +234.00000 0.74804941 0.555586298 0.597887414 20.78639911 8.918380289 1.676338158 0.419080721 +235.00000 0.773630349 0.573498182 0.616534702 21.12424798 9.007760343 1.68771381 0.421384542 +236.00000 0.799858888 0.591828862 0.635598674 21.46495584 9.09739452 1.699098526 0.423689136 +237.00000 0.826744832 0.610583791 0.655084372 21.80851683 9.187280223 1.710492207 0.425994497 +238.00000 0.854297985 0.629768423 0.674996719 22.15492461 9.277415314 1.721894755 0.42830062 +239.00000 0.882528196 0.649388241 0.6953407 22.50417281 9.367797502 1.733306098 0.430607498 +240.00000 0.911445184 0.669448671 0.716121266 22.85625509 9.458424497 1.744726112 0.43291512 +241.00000 0.941058797 0.689955165 0.737343314 23.21116412 9.549294313 1.756154701 0.435223487 +242.00000 0.97137884 0.710913062 0.759011767 23.56889355 9.640404354 1.767591792 0.437532586 +243.00000 1.002415118 0.732327818 0.781131491 23.92943702 9.73175294 1.779037285 0.439842407 +244.00000 1.034177435 0.754204798 0.803707351 24.29278622 9.823337476 1.790491084 0.44215295 +245.00000 1.066675513 0.776549343 0.826744126 24.65893481 9.915155974 1.801953092 0.444464207 +246.00000 1.09991924 0.799366792 0.850246652 25.02787546 10.00720645 1.813423234 0.446776164 +247.00000 1.133918295 0.822662482 0.874219706 25.39960034 10.09948677 1.824901437 0.44908883 +248.00000 1.168682524 0.846441668 0.898668007 25.77410261 10.19199494 1.83638758 0.451402183 +249.00000 1.204221647 0.870709631 0.923596275 26.15137397 10.28472868 1.847881588 0.453716224 +250.00000 1.240545343 0.895471652 0.949009198 26.5314071 10.37768615 1.859383389 0.456030945 +251.00000 1.277663415 0.920732927 0.974911468 26.91419465 10.47086536 1.87089291 0.458346337 +252.00000 1.31558554 0.94649868 1.001307685 27.29972784 10.56426418 1.882410052 0.460662398 +253.00000 1.354321396 0.972774049 1.028202452 27.68799884 10.65788093 1.893934743 0.462979123 +254.00000 1.393880618 0.999564201 1.055600341 28.07900033 10.75171347 1.905466884 0.465296507 +255.00000 1.43427284 1.026874273 1.083505866 28.47272351 10.84575982 1.917006403 0.467614538 +256.00000 1.475507696 1.054709319 1.111923569 28.86916008 10.94001829 1.92855325 0.46993321 +257.00000 1.517594778 1.083074419 1.140857878 29.26830221 11.03448705 1.940107328 0.472252529 +258.00000 1.560543635 1.111974623 1.170313248 29.67014112 11.12916397 1.951668589 0.474572472 +259.00000 1.604363775 1.1414149 1.200294074 30.07466848 11.22404736 1.963236908 0.476893044 +260.00000 1.64906466 1.171400242 1.230804724 30.48187552 11.31913539 1.974812239 0.479214241 +261.00000 1.694655841 1.201935616 1.261849537 30.89175392 11.41442624 1.986394532 0.48153605 +262.00000 1.741146696 1.233025872 1.293432792 31.30429537 11.50991821 1.997983665 0.483858472 +263.00000 1.788546561 1.26467592 1.325558767 31.71949011 11.60560948 2.009579638 0.486181495 +264.00000 1.836864899 1.29689061 1.358231684 32.13733031 11.70149836 2.021182305 0.48850512 +265.00000 1.88611092 1.329674707 1.391455706 32.55780669 11.79758303 2.032791664 0.490829339 +266.00000 1.936294002 1.363033034 1.425234993 32.98090997 11.89386196 2.044407595 0.493154149 +267.00000 1.98742331 1.396970298 1.459573651 33.40663183 11.99033331 2.056030073 0.495479542 +268.00000 2.03950801 1.431491235 1.494475752 33.834963 12.0869954 2.067659 0.497805513 +269.00000 2.092557353 1.466600466 1.529945342 34.26589469 12.18384656 2.079294327 0.500132057 +270.00000 2.146580333 1.502302643 1.565986378 34.69941711 12.28088526 2.090936006 0.502459174 +271.00000 2.201586075 1.538602356 1.602602875 35.13552196 12.37810997 2.102583938 0.504786852 +272.00000 2.257583572 1.575504143 1.639798673 35.57419897 12.47551885 2.1142381 0.50711509 +273.00000 2.31458182 1.613012537 1.67757773 36.01544034 12.57311039 2.125898393 0.509443878 +274.00000 2.372589643 1.651132017 1.715943796 36.45923532 12.67088305 2.137564793 0.511773216 +275.00000 2.431615952 1.689867032 1.754900741 36.90557558 12.7688353 2.149237203 0.514103104 +276.00000 2.491669613 1.729221945 1.794452259 37.35445137 12.86696563 2.160915597 0.51643353 +277.00000 2.55275928 1.769201147 1.8346021 37.8058534 12.96527249 2.172599879 0.518764488 +278.00000 2.614893692 1.809808917 1.8753539 38.25977189 13.06375437 2.184290048 0.521095978 +279.00000 2.678081547 1.851049587 1.916711322 38.71619805 13.16240974 2.195985983 0.523427995 +280.00000 2.74233137 1.892927378 1.958677942 39.1751221 13.26123707 2.207687682 0.525760534 +281.00000 2.807651729 1.935446452 2.001257278 39.63653428 13.36023513 2.219395049 0.528093593 +282.00000 2.874051109 1.978611 2.044452905 40.10042579 13.45940226 2.231108059 0.530427162 +283.00000 2.941537908 2.022425125 2.088268195 40.56678637 13.55873707 2.242826638 0.53276124 +284.00000 3.010120481 2.066892905 2.132706664 41.03560624 13.65823834 2.254550738 0.535095822 +285.00000 3.079807184 2.112018358 2.177771596 41.5068771 13.75790439 2.26628031 0.537430902 +286.00000 3.150606159 2.157805474 2.22346639 41.9805882 13.85773401 2.278015306 0.539766486 +287.00000 3.222525676 2.204258213 2.2697943 42.45672976 13.95772581 2.289755675 0.542102551 +288.00000 3.295573793 2.251380423 2.316758609 42.93529349 14.05787827 2.30150137 0.544439114 +289.00000 3.36975861 2.299176036 2.364362512 43.41626813 14.15819032 2.313252318 0.546776153 +290.00000 3.445088055 2.34764884 2.412609144 43.89964538 14.25866043 2.325008468 0.549113677 +291.00000 3.521570101 2.396802568 2.461501672 44.38541499 14.35928738 2.336769797 0.551451672 +292.00000 3.599212592 2.446641008 2.511043144 44.87356669 14.46006996 2.348536256 0.553790141 +293.00000 3.678023286 2.497167778 2.561236608 45.36409217 14.56100662 2.360307772 0.55612908 +294.00000 3.758009942 2.548386579 2.612085055 45.8569802 14.66209616 2.372084345 0.558468477 +295.00000 3.839180234 2.600300971 2.663591443 46.35222247 14.7633375 2.383865852 0.560808344 +296.00000 3.921541749 2.652914543 2.715758646 46.84980871 14.86472926 2.395652318 0.563148657 +297.00000 4.005101992 2.706230767 2.768589566 47.34972867 14.96627023 2.40744367 0.565489428 +298.00000 4.089868422 2.760253088 2.822086987 47.85197355 15.06795903 2.419239859 0.567830651 +299.00000 4.175848413 2.81498498 2.876253723 48.35653309 15.16979475 2.43104086 0.570172315 +300.00000 4.2630493 2.8704298 2.9310925 48.863398 15.271776 2.4428466 0.57251442 +301.00000 4.351478158 2.926590907 2.986605926 49.37255949 15.37390142 2.45465703 0.574856977 +302.00000 4.441142874 2.983471344 3.042796728 49.88400291 15.4761704 2.466472102 0.577199935 +303.00000 4.532049416 3.041074842 3.099667837 50.39772824 15.5785814 2.478291816 0.579543351 +304.00000 4.624205884 3.099404272 3.15722101 50.91371595 15.6811329 2.49011617 0.581887168 +305.00000 4.717618247 3.158462791 3.215459767 51.43196115 15.7838249 2.501945166 0.584231443 +306.00000 4.812294177 3.21825327 3.274385571 51.95245407 15.88665435 2.51377856 0.586576118 +307.00000 4.908240496 3.27877944 3.334001647 52.47518004 15.98962278 2.52561635 0.588921252 +308.00000 5.005463599 3.340043884 3.394310048 53.00013907 16.09272559 2.537458537 0.591266786 +309.00000 5.103969879 3.402049473 3.455313118 53.5273165 16.19596585 2.549305366 0.593612664 +310.00000 5.203766584 3.464799365 3.517013494 54.05669767 16.29933898 2.561156348 0.595959 +311.00000 5.304860109 3.528296143 3.579412936 54.58828257 16.40284496 2.573011483 0.598305737 +312.00000 5.407256421 3.592543251 3.642513789 55.12205167 16.50648382 2.584871014 0.600652874 +313.00000 5.510962769 3.657542698 3.706318689 55.65800497 16.61025401 2.596734943 0.603000412 +314.00000 5.615984694 3.72329793 3.770829397 56.1961278 16.71415401 2.608603025 0.605348351 +315.00000 5.72232859 3.789811242 3.836047964 56.73641039 16.81818229 2.620475015 0.607696691 +316.00000 5.830000853 3.857085506 3.901976735 57.27884297 16.92233886 2.632351402 0.610045431 +317.00000 5.93900745 3.925123303 3.968617468 57.82342066 17.02662218 2.644231698 0.612394515 +318.00000 6.049354349 3.993927218 4.035972215 58.3701239 17.13103226 2.656116147 0.614744 +319.00000 6.16104752 4.063499548 4.104043027 58.9189527 17.23556604 2.668004504 0.617093828 +320.00000 6.27409293 4.133843162 4.172831665 59.46989728 17.34022505 2.67989677 0.619444114 +321.00000 6.388496121 4.204960069 4.242339886 60.02294299 17.44500623 2.691793188 0.621794686 +322.00000 6.504263061 4.276853141 4.312569741 60.57808494 17.54990806 2.703693516 0.624145717 +323.00000 6.621399293 4.349524386 4.383523283 61.13530848 17.65493359 2.715597507 0.626497034 +324.00000 6.739910358 4.422976389 4.455201685 61.6946087 17.76007824 2.727505407 0.628848751 +325.00000 6.859801799 4.497211444 4.527606704 62.25597096 17.86534048 2.739417216 0.631200869 +326.00000 6.981078731 4.572231275 4.600740393 62.81939526 17.97072337 2.751332933 0.633553274 +327.00000 7.103747122 4.648038752 4.674604217 63.38486693 18.07622233 2.76325207 0.635906079 +328.00000 7.227812088 4.72463531 4.749199349 63.95237133 18.18183735 2.775175116 0.638259228 +329.00000 7.353278318 4.80202382 4.824528133 64.52190843 18.28756844 2.787101826 0.64061272 +330.00000 7.480151354 4.880205716 4.900591449 65.0934636 18.39341559 2.7990322 0.642966613 +331.00000 7.608436313 4.959183583 4.977390762 65.66703194 18.49937423 2.810966238 0.645320793 +332.00000 7.738137883 5.038958855 5.054927832 66.24259879 18.6054474 2.822903697 0.647675315 +333.00000 7.869261179 5.119533829 5.133204122 66.82015927 18.71163206 2.83484482 0.650030239 +334.00000 8.001811319 5.200910226 5.212220514 67.39969871 18.81792668 2.846789607 0.652385449 +335.00000 8.135792137 5.28308977 5.291978765 67.98121224 18.92433277 2.858737569 0.654741002 +336.00000 8.271209177 5.366074757 5.372479755 68.56469496 19.03084883 2.870689196 0.657096899 +337.00000 8.408066701 5.449866622 5.453725243 69.15012733 19.13747332 2.882644488 0.659453139 +338.00000 8.546368973 5.534467086 5.535716107 69.73750935 19.2442047 2.894602955 0.661809666 +339.00000 8.686121108 5.619877873 5.61845352 70.32682637 19.351043 2.906564842 0.664166536 +340.00000 8.827326942 5.706100704 5.701938656 70.91807348 19.45798819 2.918529904 0.66652375 +341.00000 8.969991592 5.793137016 5.786172978 71.51124093 19.56503876 2.930498631 0.668881249 +342.00000 9.114118042 5.880988816 5.871156781 72.10631894 19.6721947 2.942470534 0.671239093 +343.00000 9.259711407 5.969657254 5.956891823 72.70329773 19.77945296 2.954445612 0.67359728 +344.00000 9.406775524 6.059144051 6.04337869 73.30216753 19.88681508 2.96642411 0.675955753 +345.00000 9.555314656 6.149450356 6.130618848 73.90292347 19.99427951 2.978405784 0.678314512 +346.00000 9.705332213 6.240577891 6.21861259 74.50555087 20.10184626 2.990390634 0.680673615 +347.00000 9.856832886 6.332527804 6.307361088 75.11004974 20.20951228 3.002378659 0.683033004 +348.00000 10.00981966 6.425301817 6.396865221 75.71640053 20.3172791 3.014369616 0.685392737 +349.00000 10.16429679 6.518900792 6.48712587 76.32460324 20.42514518 3.026363993 0.687752756 +350.00000 10.3202677 6.613326451 6.57814362 76.93464811 20.533109 3.038361301 0.690113061 +351.00000 10.47773621 6.708579654 6.669919643 77.54652047 20.64117209 3.050361785 0.692473653 +352.00000 10.63670575 6.804661551 6.762454233 78.16021543 20.74933139 3.062365445 0.694834588 +353.00000 10.79718014 6.901573576 6.85574827 78.77572811 20.8575869 3.074372036 0.697195752 +354.00000 10.95916194 6.999316877 6.949802631 79.39304384 20.96593862 3.086381558 0.699557259 +355.00000 11.12265543 7.097892316 7.044617612 80.01215776 21.07438503 3.098394012 0.701919053 +356.00000 11.28766314 7.197301041 7.140193797 80.63306007 21.1829246 3.110409397 0.704281133 +357.00000 11.45418851 7.2975442 7.236532359 81.25574589 21.29155885 3.122427958 0.706643499 +358.00000 11.62223493 7.398622654 7.333633005 81.88020058 21.40028626 3.134449207 0.709006094 +359.00000 11.79180538 7.500537264 7.431496907 82.50641922 21.5091053 3.146473386 0.711369033 +360.00000 11.96290287 7.603288891 7.53012436 83.13439207 21.61801597 3.158500497 0.713732258 +361.00000 12.13552993 7.706878684 7.629515654 83.76410934 21.72701827 3.170530295 0.716095712 +362.00000 12.30968957 7.811307504 7.729671378 84.39557103 21.83611067 3.182563025 0.718459452 +363.00000 12.48538517 7.916575924 7.830591824 85.0287576 21.94529165 3.194598686 0.720823479 +364.00000 12.66261847 8.022684806 7.932277285 85.66366905 22.05456274 3.206636789 0.723187792 +365.00000 12.84139327 8.129634724 8.034728348 86.30029072 22.16392087 3.218677825 0.725552391 +366.00000 13.02171131 8.237426826 8.137945304 86.9386226 22.27336758 3.230721791 0.727917219 +367.00000 13.20357512 8.346061399 8.241928156 87.57865005 22.38290134 3.242768201 0.730282333 +368.00000 13.3869877 8.455539018 8.346677195 88.22036328 22.49252215 3.254817297 0.732647677 +369.00000 13.57195161 8.56586083 8.452193007 88.8637623 22.60222695 3.266869325 0.735013306 +370.00000 13.75846854 8.677026835 8.558475594 89.50883734 22.7120188 3.278923796 0.737379222 +371.00000 13.94654106 8.789037895 8.665524954 90.15557373 22.82189465 3.290980954 0.739745367 +372.00000 14.1361713 8.901894583 8.773341674 90.80396659 22.93185449 3.303040555 0.742111798 +373.00000 14.32736138 9.015597187 8.881925168 91.45401103 23.0418968 3.315103087 0.744478458 +374.00000 14.52011388 9.130146281 8.991276316 92.1056924 23.15202311 3.327167818 0.746845347 +375.00000 14.71443049 9.245542152 9.10139453 92.75901069 23.26223035 3.33923548 0.749212523 +376.00000 14.91031291 9.361785374 9.212280104 93.41395613 23.37252006 3.351305341 0.751579984 +377.00000 15.10776372 9.478876233 9.323933331 94.07051895 23.48289224 3.363377889 0.753947618 +378.00000 15.30678417 9.596815017 9.436353625 94.72869426 23.59334231 3.375452879 0.756315537 +379.00000 15.50737642 9.715602301 9.549541573 95.38846741 23.70387484 3.387530313 0.758683686 +380.00000 15.70954215 9.835237796 9.663496587 96.04983839 23.81448526 3.399610189 0.761052121 +381.00000 15.91328308 9.955722365 9.778218961 96.71279742 23.92517509 3.411692753 0.763420785 +382.00000 16.11860091 10.07705601 9.893708402 97.37733475 24.03594281 3.423777515 0.76578962 +383.00000 16.32549693 10.19923872 10.00996491 98.04344548 24.14678842 3.43586472 0.768158799 +384.00000 16.5339724 10.3222708 10.12698819 98.71111983 24.25771038 3.447954368 0.77052815 +385.00000 16.74402947 10.44615252 10.24477795 99.38035293 24.36871023 3.460046214 0.77289773 +386.00000 16.9556694 10.57088389 10.36333449 100.051135 24.47978644 3.472140503 0.775267539 +387.00000 17.16889306 10.69646519 10.48265692 100.7234612 24.59093748 3.484237236 0.777637634 +388.00000 17.38370257 10.82289614 10.60274554 101.3973167 24.70216335 3.496336166 0.780007901 +389.00000 17.60009836 10.95017674 10.72359976 102.0727066 24.81346405 3.508437295 0.782378455 +390.00000 17.81808172 11.07830727 10.84521929 102.7496113 24.92483959 3.520540867 0.78474918 +391.00000 18.03765433 11.20728774 10.96760412 103.4280307 25.03628843 3.532646638 0.787120134 +392.00000 18.25881664 11.33711785 11.09075368 104.1079551 25.14781057 3.544754852 0.789491374 +393.00000 18.48156992 11.4677979 11.21466796 104.7893796 25.25940449 3.556865019 0.791862786 +394.00000 18.70591545 11.59932731 11.33934579 105.4722945 25.37107019 3.56897763 0.794234427 +395.00000 18.93185365 11.73170608 11.46478776 106.1566899 25.4828092 3.581092439 0.796606297 +396.00000 19.15938538 11.8649345 11.59099269 106.8425658 25.59461845 3.593209202 0.798978396 +397.00000 19.38851192 11.99901198 11.71796058 107.5299126 25.70649795 3.605328409 0.801350666 +398.00000 19.61923369 12.13393826 11.84569056 108.2187204 25.81844771 3.617449569 0.803723166 +399.00000 19.85155154 12.26971332 11.97418262 108.9089843 25.93046771 3.629572928 0.806095895 +400.00000 20.08546591 12.40633717 12.10343618 109.6006946 26.04255644 3.641698486 0.808468853 +401.00000 20.32097764 12.54380894 12.23345065 110.2938512 26.15471389 3.653825998 0.810841982 +402.00000 20.55808759 12.68212892 12.36422544 110.9884395 26.26694006 3.665955708 0.813215341 +403.00000 20.79679533 12.82129655 12.49575998 111.6844546 26.37923343 3.678087618 0.815588928 +404.00000 21.03710257 12.96131123 12.62805366 112.3818917 26.49159552 3.690221481 0.817962688 +405.00000 21.27900888 13.10217327 12.76110621 113.0807458 26.60402328 3.702357298 0.820336676 +406.00000 21.52251468 13.24388151 12.89491674 113.7810071 26.71651824 3.714495315 0.822710893 +407.00000 21.7676204 13.38643624 13.02948466 114.482666 26.82907886 3.726635041 0.825085282 +408.00000 22.01432648 13.5298366 13.16480939 115.1857223 26.94170516 3.73877721 0.827459843 +409.00000 22.2626329 13.67408229 13.30089034 115.8901664 27.05439712 3.750921089 0.829834633 +410.00000 22.5125401 13.81917305 13.43772693 116.5959884 27.16715322 3.763066922 0.832209652 +411.00000 22.76404849 13.965108 13.57531828 117.3031835 27.27997347 3.775214954 0.834584842 +412.00000 23.01715724 14.11188685 13.7136635 118.0117516 27.39285786 3.787364696 0.836960205 +413.00000 23.27186719 14.25950932 13.8527623 118.7216732 27.50580638 3.799516636 0.839335796 +414.00000 23.52817834 14.40797455 13.99261352 119.4329533 27.61881753 3.811670287 0.841711616 +415.00000 23.78609026 14.55728197 14.13321685 120.145582 27.73189128 3.823825891 0.844087551 +416.00000 24.04560339 14.70743099 14.27457113 120.8595544 27.84502765 3.83598345 0.846463715 +417.00000 24.30671687 14.85842134 14.41667577 121.5748559 27.95822511 3.848142964 0.848840051 +418.00000 24.56943112 15.01025215 14.55953019 122.2914914 28.07148518 3.860304187 0.851216615 +419.00000 24.83374571 15.16292256 14.7031332 123.0094462 28.18480634 3.872467364 0.853593352 +420.00000 25.09966066 15.31643258 14.84748423 123.7287203 28.29818706 3.884632496 0.85597026 +421.00000 25.36717553 15.47078075 14.9925821 124.4492991 28.41162887 3.896799338 0.858347397 +422.00000 25.6362899 15.6259671 15.13842623 125.1711873 28.52513024 3.90896789 0.860724648 +423.00000 25.90700333 15.78199017 15.28501603 125.8943705 28.63869116 3.921138396 0.863102129 +424.00000 26.17931584 15.93884997 15.43235003 126.6188437 28.75231012 3.933310856 0.865479781 +425.00000 26.45322656 16.09654536 15.58042765 127.344602 28.86599017 3.945484782 0.867857606 +426.00000 26.7287355 16.25507575 15.72924801 128.0716405 28.97972672 3.957660662 0.870235601 +427.00000 27.00584138 16.41444029 15.87880994 128.7999494 29.09352283 3.969838253 0.872613826 +428.00000 27.28454419 16.57463812 16.02911284 129.5295239 29.20737545 3.982017797 0.874992166 +429.00000 27.56484352 16.73566837 16.18015556 130.260364 29.32128457 3.994198808 0.877370734 +430.00000 27.84673851 16.89753047 16.3319372 130.992455 29.43525173 4.006381528 0.879749474 +431.00000 28.1302283 17.06022327 16.48445689 131.7257969 29.54927539 4.018566202 0.882128386 +432.00000 28.41531248 17.2237462 16.63771346 132.4603799 29.66335555 4.030752343 0.884507413 +433.00000 28.70199018 17.3880984 16.79170603 133.1961991 29.7774907 4.042940437 0.886886668 +434.00000 28.99026056 17.55327901 16.94643342 133.9332497 29.89168235 4.055129997 0.889266095 +435.00000 29.28012276 17.71928658 17.10189475 134.6715268 30.00592898 4.067321268 0.891645694 +436.00000 29.57157636 17.88612084 17.25808915 135.4110206 30.12023059 4.079514248 0.894025408 +437.00000 29.86462007 18.05378064 17.41501545 136.1517311 30.23458565 4.091708938 0.89640535 +438.00000 30.15925305 18.22226511 17.57267247 136.8936485 30.34899568 4.103905338 0.898785464 +439.00000 30.45547444 18.3915731 17.73105933 137.636763 30.46345917 4.116103204 0.901165693 diff --git a/data/diviner/diviner_filter_functions.hdf b/data/diviner/diviner_filter_functions.hdf new file mode 100644 index 0000000000000000000000000000000000000000..0fa13e27111b39bba8fd1af81bb1ae46f8332f8e GIT binary patch literal 18264 zcmeHu1yEeix8~r%AvgpH1PBrwf(?XVA$V|i26uNt@ZiBU5ZomU?hL^p*bsbh4>FLT zOYFV+y{iAdt=iq%dT(>;-v0WWKKIP+ujcgBnd%OIAt8m0L5>0VA@}Y9kO6oAq@T*K z=JvN}O8-;;OzCI$x1HvWO~0aGIvqaS@6BV#k0JM=K} zuynq|KXI4u`u`n8^_8ULk6HfV`49e?>oW4M9_b&$U;O`y0R#Zz&m5rnJ!|~C{QQv< z0lyQ&@VEUBAOPc&ubz}h6-I%=-=%vY57x#z!@c*a$iu1>&qM#^& zaJQoUv5x(I|3BmZ5%GQ-Ebjluk$>uc^6{SFKkffrwIAyg!X4Xx)~7%CiSmDK#?J`< z>G!X`SHG9Le>(rw{xa~Ffxis=W#BIZe;N48z+VRbGVqsy|56Mf0MLHbYQGD%f2&Io z|J0?w>Q;igS#jOn2zPbv5Bc}{@~3)N$Nj2be?JDe_rv*4euVwe{nvFQ>d!jT#?sF0 z$Nx(IlRENuX`*@O^;iBg20SwMGeAnJgT#koD5{UbXj3=?Y3|W-G&i5i(P4z;7^4#6 zAxp|I3%umWF;wR&$V)+yOaYSOMMoH`V};Q^a6T*Pd$8Cz*W_$fYhJRLe7>SKYw&nS zqgBPdQonWDKPqD2@T~dT?v_2>qNB@Di-XIIuZbLNKL0I;j%(p-vRvQS!~n{7$j&U~ zR{P{G(|~$QwYIu`w=~vrF=-uG+)=eaUQ=n-)<7 zj$JYAw`F?d8sEO#SuhtyS|JwkD16&d`4)fMn*PNV8$%BNsOhL zFij%0ukKy2_WXh2<$k7wVo5z$b!7`eWQ%scOSVU~57->p5QR*0vc(5yGvRwaulSC5 z1PeLCltW`HB|`1+&Pb|)&yeJFJK_+R2k$23nVeuX3#t1g6~V{=ZG5#z3!d8-~;;C^oaEybOa4 zggb&D0a zLeEMNIL_lv7*6e%TZ{&rPYhD^Jd#Xx45V8=f56A_ruu~zJM;M`)S+)7K^FZ^Qgy;z zR8ttcJdNdH4ih* zb$uMscic!?sHe}%O27J~46VH+=)6|993r59vafL_o!(kb<_6aoU3-OSQ~I_w zePFv)^I>bMI#xx0iQ@Xkep^q@hSgVAq+O&21oMg1+6&qaNaY^nbb_=`LNpUQ1x`H^ z+nMRKtv_`W;3_Ccjxz7V=IggCVz`5LkkcsJlHwY3i#+rdLf&)c5$qyzjyzO99_M*! zf@I>Nhg8LMO19dzo?aF@Wz<1j1TSxRPr;1o{t|(;)fs!nc@tkNw9f=#@DwAK%JSPI za@TR9D!A=X3TH)asxI;+QKJ5HX8h?P@;dyg`*wI%p`u&@F#`Z0p>Ie8+qy$MrrnLF zRTJ3FnV#>yS7Vi-<-%5(1m|)_<|mr~nwzS0<`L408zb|s0jXrr*^NRJI&cGO2ZYqx z2^+;=T;0EFQ3e``QnM0n_4u3=sv@JlY^XK5do+Y1?a>a?=-0K;Z@O99t03(dBKDgwS>3uiV>+&$MwZJw`-TI02I>*cAM!{a31 zGp~B~yj+N@Ik6x`r-LJ85o=y>0*GG9Vq1|^4e8V-*Au8cA3|H%EZ=^ z5+uP4h}3cAai* z=+xOh&1yfTbKaT5{nS?3hfdCFbxxFlzm;3P-e9<7WWM_TFo1q|onz)bX6+gO*Y6fB zYxcVN`VqqSe3G_sNV5gD!|Jz?L?hY8BKp`N>9mLMvjR~1>hQ{wRLx%itF(AHnzFDo zGy9{m%P#gMKz z=*D?^?t$QS>kmDSbs7UnZ$WI!bmNjh&w{#ERXd-5~NAwrr<}{N=S4 z0~L8`6plX@vfOAvIe)GKNh{Snd9e~vH^b3Q7hl6D-*U4IxmR+eC~jt`{3p8(oM7vh zA>}R|#FE+pDc>{yydS}mzM94V#^RG z(QA3=W{d&gK7n4v$d?j>S^*;XbFN4HX%3*Gh6=2`VjG3X(`yy;;{Rk$0Dp|!r$1X~i& zLX_RNBT08KX+8!cC*FM(5XtLHKTFm6Mqv)IjqJ!V`%!ADNcIRe@+N?5XIn}pS^Cg& z?GdnJkFK?6DxA^M2ldZt$Eh|$wE@3v+Gm^$ zqj0f8DjTg&@=cEs%g{bwx25O_D<(4xht5ns-Ox-{JhyB?(KPAdFnKc9YC|36hUWKZ z8@`I}$yprCwugHLFg9!s%4krQ76*9)3MkwY2;bz7&w{-W3RD^ck-`3wAlAU^CVnCO zOAYCwU6N)0``wNqJyv~wR(VK16rTrc!cWhxq~Lmiw=x%+wv5b!J@I{+Ih9K+{Dc}Y zZi$44qNv&Q`)nh9cJW%^Gz9BPZvZ2N^hHx@_esRg!DE{*=k*N`Bq(047Z{(H8cM5h z9EnzTw}f2eUOyN9u>Wlut6g}CSg4iT(4ZWJ4;f4-(u0|f4Q=hW1zdh!qMmDM>0^jS zZ9}z8-wBIHt1y)eS{As9+IJj*J?ENpBlbY^*-;kYvf7qmdw^Ok>~joxdb(pwH0zD= zmfFt=_Sh~P99dTP`ntkhNiy^}l;D=P8)5aes&XWf-ED1bfZ=MxO5w!A2pQyk%Y< zjRJpp=KZ-VZ42~NWZ=?=qDo?vxkcNIa;N2fnf%735Bj1JaaTDIVFUAvX9*eg*>_f1 zTogY_n6UIRa9dBn*103)b%vh@$nl}c8B))(84eNz7H=nG@rDn&$}D0J42)!F?Q&K- z47_|IU>QAQx&&FA;s|NpwZT(xYYUHnsQ_xTj|oCs{DP*30Dg}W1mCmpiFl#=dLUoX zkMD!MtD?TE?s*;>Swk#2c=;8rS@msEeMDxV0j4b#T%_rqbHi1%a_ww z%^RU6@;q2XL{kRr`!mX;zMG$t?1 zLtCUl4Qong?0*Yh>}g^^)8XDOIE)@@!lZVX+3m<))`#kvZSaG|LT=b~YG0nYlrPXF zhz*@xE?+SYxgQNmyv7~>kLTb^Y2WvYgM-OAB(ZHryBHWzd+shK-h!-EOXlD{ta$H zCR4W}jv= zXkzn{Y8yutSX63l2DMtK#z+RVjCsd|K550r9SdYb$6tKE_9A358+SfqGQ(`@Wj5<{eR(b~ zD6Q5QbsB>xs4j5QuV%+sK=L>{h)RKmniiCx0f_pLA-|VoEKlU>(Oa z@2;o=mlQz?FgO~N&XWtnH=D(<_POZMP|5DIcV=Kc=pjdCdv^$0ABGyRO8K%z zJIbbWsZT!Ep0_BxrRQRm?3hlKle>(fI7KZQxTjxaV)JrU9(WF;4z%t>P^Dwoh)^1& ze2oN>1=%s)^YPxP3pPkEGH+slHSG~{T*5|#FC7(d%Yn2h8ib0gzNquj1$2tefw;1^ z8^_xBTq>cI#UOO22QHdyhqZY-^A8vfHQf=ovnPjY-G}$_TUlBPW_t*4{DP9%rKZ#o zTZwsr*VYf_pNNke`3LGpYG;u3a~C$%eljr!A=E~!_*%FOM~MfOJjK+atQ?inI}wj7 zZZ>SsTcbG4%Za~!HMyVv3{zHmjp5{+LID^%Vjal-a^cIdHp1IxA{r9gPOVslt-(uz zqUr^}I|}-*eMKNn|7D3;EiM<+f({(zTmBpWGP+?glBei(CN7JpBDf)zo!Bl=xt!EO z<%RF?z{RsN@rem9Oih=?*S9xjH=iZth`gV7z`!Tg`LMbrSq3i2D#g!lv_@|)HOTiV zdWD5)zs6YXS?ce8l4Q3K+ZI4v%90{EazriIH1VH2K1m9Yv^etSFLCq<7E%^#eSrP_ zqp-uMwaIv>oY+Xr$kDQM#?~COVx1Q47_rLXg@Lu1i8#+OFG)sT5o4dfB4LUsGSLCS zgVXIt27YDZocpdtMvc!3pY<)Zc8d-=Y&xCt8Y~Zvjc>HL5uz2P8N;JrbF)xRMODUh zMDKsevkyY?W_J}Hyv&An7L$~eofFhhw6sTdU(QRnCBD86h>kz*vfEqm7pP%_GW#T% z%y^o!dWLRf;#-_eNxwzM=g7vP7nT))h9N;PmW-MN+BH_L zCi|c#$@@v%=MaM#RxjHSX;ciDl1tr_N^!_8+qz*Mg+FFWH$hK{Mt1Kw{e%Mrbs$4( z4FTG!=&e^M+K#QnQc^iF$F3kr3)D)eiOLnr-Ah!iLxQ?*IofWw`XDSy^SrHN^nC%6 zUNkHEc-JDG>$oDtdr_lK;i4hoE9wU}H~5(@v|WA-eKAm}z3%NN_usjS^P{1O%zBtU zOTZ9uUT4xiR#K{!iK0qc{cOh_Zly%^vLF9jZ=J>>{p>_n_qSbqRLV>7W$TF%djI$m zBcg_@JfQF(@HlB<@$Ib)-*Ip&Cam9pnQS!X;_7{pdl1`eA zh;y8&*s8=@;pIJ^Nmv-(RCcrC|6<*w106kilx~&yc3S?LV7UH%OY@rVu;kgY=IiO9 zxfq6qduI;)&6o19HU$=?{trIycr^`2AK^NAMGz?BN-@r(KkH z6UjV9zjE)l1k?7%;%sk))?x_q%W^P2B9Y2|_eS)nTTdD}$Cbh4&Gj_u_NqXH| zX_u5|8c$|TYbB0#H`(0IpwW}2h@xsv4#96rAn#W4*)A-u(8Qj^n3~PHf4ST+)^*oaj?Q=UqZIkNAPMG;16qv&wfy3=FT#_Z)SL6z~;D|`??q`XFTLe>vnN5V} zNvu^KE2(PT8;PCFYvyiTH4ajL6!?t}UQtFrb6=g)x-aMfDt4)X*hnUI56S^TwLZE3 zj(`BS_Hczr%M(L%+vK+cNcP>xOFn{0A#1A@QP4Fr`HVrR-;Ud@+yG zmqsNXQU)ePGGyC4OlOW%o+ru?&5xu)iO)?A{4|HW)?!hk0A!lq3;9xk%88Sv%biZh zSHfw)@Y>5iP3~#Jl^Zzra9JdFKR!0f^7dOQ`LeKV0`*#4Mtq=!S~5idJu+Y1>v+*v z?;+}d-aS|P@B6+fg&<-;y~CAf&HI{+{xig*XM&X<=(KA2@GX6i(gX#Q1te>@rt>mv zaw_wk2c_eCFE-=$4v_*-qN$$KM)Bo*55Gh>@!69@^!EX=auiwLt#HJ#TG_%4UF6?} z3tXrbyNsRloA*|HQk5-U!lK`++}yz0C9m-FuB(swPfF( z4?cofvz{DdgjIePP&tS*qnP0+1$Sx-w&RJ`ob4!^pS$xXx0nN+w z7tiE_IXF+|r9G~yA(FAGi%0p}9^1B{k#Aj&)OKcyv}qT{QxAmz1R!c!8GA+}N)R3> zZ-@i4rsv?Hoet(q-NZhdspBl2syAr`S~}a*BeQOkcd`UZ zgS;?5@i%gE0sC?DlcZ7*IBC%Rk=MXyPLJ%LlDWHZXbX}}D!t^h)MMaIvX%71?Kn(+ z@!UW?{_7WQ(`M`*XL>#a|2T)Km>1N-6ubKmndYltUzB7o@`N0*dh{MM*n8KHBVuVijZ>>FHgize=u zmR-Hbjktud9At{M3uoScMN%1UVE$U$G~Vzgj<@vNW~>PcWW{2R7Cl{WO=i;D2hrx^ zp;fipwZi8k?4pOoWnlt`QY%-&qV%DAp+id`Y|I>+dQE?Qsa{(k=X; zp7D~ki*SRg*iN0$F{6URK!PkA)AWmJu&2JDXN>4kJt|yXIK4F&GnG%RQDw$ zKM)uViZdDwFHL&xiyGxaCJ!MET9R1Nn#{ci)BM2kRNBN11@DNy_dGJe3a>!fDC;#O zUxUD{p-Y!7D$XJqnA zv0`p^2}dP-RKSf@#Py=^Ip4982%}C32Z`c#wOC?|IrHQ*<)aPu#8FYEA+U<$K-|Fa zeI=7`yJ?|W19s=?Dd0SuzT_2OuHYs*L9V$BRcQkapWzG=8)eY)`PU(~wKV8?p|QF8 zgiP+rJ*zyA=|Uac{?DEfePYcJna@zP4CC~~_+ z5o42AXl@{O};uTP8?++lj&kp8qc)qd-H(IlXA2u-wF=!KvQkD#iXy&WiPk% z5iR+41&79PbVZX)-U87Vr|Cv{`l#;!t9B2*m%;&3g@H{LS+3S8*P~SLl~yISYv2|( z-od+pYGe)R4A}i@J!8A0I+hV3UEjM(G<0)iqA|{;_S{(BLOg?(_^%A|IH_``RYr2# zSZ1;%f*;S+N$jgubVr_{v=<-Qx1%g)(1DCphz`_*)5&l~dYv{7N_B%)pww53Q5jKk zq#xcjqA<;*cr2i7N(o?Im8ovZE~q40$D^h&&98b)AJx3Nj{TtUL6w^RVM2%2hQE9n z1FEb2EqMPH;Hzdt^4!SIl!hDQC-D3N-e17?3!eT0{$C(~0r1Emns$As*Nf111a9zl zbaw}@4Ii!?92~?y#~`e^)7{+_`LyR-dy4^uF9bRe6q>KsVg#0a&JN7a7lho+twuGN zKeqFb!lhJ#W{$P>2dzD^D_3}1LNrdS^oL75i7H`gEfE^WFZKIbJW(o_iCSJ$Gz8vV z$tOHG1MnHM;649J6_>n>Fa{i8Hz)ARG!hvB?^)^&}pJX$ePyj^(oYxjj)wm&t^K z3Qt_JMkYMT6;~~VqY95zvj$W=F%(yrgo7DROtXd(Jnt*QWQ4;Rk0r7`5qKghE~5(v z_MLeA_>2cJtB%6aeaD(qgC?F>F)N(HA<`!nRKt0mgfTE>;RxyD7gYUpp2#uF_`68tDXj30_`%cfo7w6Nxo$U{%lr(mGTN_BuU0NB> z*Rtk=yd7{4Iu{y++&2&AqnaI`9zgUOS=`qv=7XCZNDjI(8@b&#M&`qX92pOKC>!bA z*P`bGg&fchI`$iZ;LWr7=s8CryFS~-C*bw^`H(pW3cK#|Mqcp7^n8SlBa2<{lSX>* zT9WQ4{O+3XYS-Cv&Kh^VU==O(-vG@Yf~FI*D=JHY<%<iA`U!8ya_qef^#NUS+#n zqZ{}?zx_0D5*G|0DlCOQJf@zBLl@^W>(?+K!uno6F4bQDrPzl{uuZ& zW|9tm5BMB+rW?)yl*648hvNYyaHoUdG{6_QlSps?@Y(sywX5iY;`x-XEA@i(`E;u* z`-1HGB+wOiLEL?2!S(5=(bI)Q_bELOl7&?F=}Z^K1!eb1(Jh>feVpBax!tp-OuNFl zo%*IIyF8uU>83=xe4U-7reM4G>AM?EsdfeFI|WVQcDWS0kfu1h4-`8gO@VefySo=n zSrvt@+J!{rs?3%=H=wJTZ>`mw?Ys7n%EYFtuJ*0btAr!E+oD@ezpwB^xX`N3>JsE> zg@2v_g+b#D`*q}iDHH)EQ0`)&Oh2i#xU~4 zQ9b-bCPN1H;65mZ-@b>8ZiL7&DeOc%gAA5|2MXfH?IEK%yN2$X_{RW&XGhRxllEw! z;MqDf(&RP@`1EWRs$(J=1>`v!hO(IWzXfuhbwV$V+arN&XCI-1#_%Z`?U7yO=Av$tr1ozjLmB%_i%czjJ#f zN?$E}@Sp7{`aAyrf`R|j=f-|-^ZT{u?`}f>yXVG8g?_Q}Bkqs;|Mt1De|6r8|Ifj; zU$65oTfyJ^68|{-D?a|7clyK7U;N9!Uk3g%@RxzV4E$x_F9UxW_{+e5T?YOKCIs}7 literal 0 HcmV?d00001 diff --git a/roughness/config.py b/roughness/config.py index 898afb2..6b8fafb 100644 --- a/roughness/config.py +++ b/roughness/config.py @@ -38,7 +38,7 @@ ROUGHNESS_DIR = ROOT_DIR / "roughness" DATA_DIR = ROOT_DIR / "data" DATA_DIR_M3 = DATA_DIR / "m3" -DATA_DIR_DIVINER = DATA_DIR / "Diviner" / "lev4" +DATA_DIR_DIVINER = DATA_DIR / "diviner" EXAMPLES_DIR = ROOT_DIR / "examples" FIGURES_DIR = ROOT_DIR / "figures" FORTRAN_DIR = ROUGHNESS_DIR / "fortran" diff --git a/roughness/diviner.py b/roughness/diviner.py index 09f7572..56a13cb 100755 --- a/roughness/diviner.py +++ b/roughness/diviner.py @@ -4,9 +4,9 @@ import numpy as np import pandas as pd import xarray as xr -import roughness.config as cfg import roughness.helpers as rh import roughness.emission as re +from roughness.config import DATA_DIR_DIVINER DIV_C = ("c3", "c4", "c5", "c6", "c7", "c8", "c9") @@ -19,8 +19,8 @@ FILT_SCALE = np.array( [1.0005, 1.0003, 1.0010, 1.0018, 1.00600, 0.71939, 0.79358] ) -FDIV_T2R = "/home/ctaiudovicic/projects/roughness/tmp/div_t2r.txt" -FDIV_FILT = "/nfs/data/josh/home/projects/diviner/diviner_filter_functions.hdf" +FDIV_T2R = DATA_DIR_DIVINER / "div_t2r.txt" +FDIV_FILT = DATA_DIR_DIVINER / "diviner_filter_functions.hdf" def lev4hourly2xr( @@ -41,14 +41,13 @@ def lev4hourly2xr( If savefile: save result to netCDF (.nc) file for quick I/O into xarray, e.g. with xr.open_dataarray(savefile). - Parameters - ---------- - fgrds () - tres (num): time resolution - interp_method (str): interpolation method (see xarray.interpolate_na) - interp_wrap (bool): interpolate around time axis (wrap around 0/24h) - ext (list): extent of ROI - savefile (str): Path to save result to netCDF (.nc) file + Parameters: + fgrds (list of str): list of paths to .grd files (e.g. from glob) + tres (num): time resolution + interp_method (str): interpolation method (see xarray.interpolate_na) + interp_wrap (bool): interpolate around time axis (wrap around 0/24h) + ext (list): extent of ROI + savefile (str): Path to save result to netCDF (.nc) file """ if not tres: # Assume 8 bands (ch3 - ch9 + tbol), diurnal [24 h], infer tres [hr] @@ -105,6 +104,10 @@ def lev4hourly2xr( def fit_poly_daytime(Tday, deg=2): """ Return polynomial fit of order deg to daytime temperature data. + + Parameters: + Tday (np.array): daytime temperature data + deg (int): order of polynomial fit """ tloc = np.linspace(6, 18, len(Tday) + 1)[:-1] nan = np.isnan(Tday) @@ -117,6 +120,10 @@ def fit_poly_daytime(Tday, deg=2): def smooth_daytime(T_xarr, savefile=None): """ Return smoothed daytime T from 6 AM to 6 PM with 2nd order polyfit. + + Parameters: + T_xarr (xr.DataArray): Diviner lev4 hourly temperature data + savefile (str): Path to save result to netCDF (.nc) file """ tres = T_xarr.tloc[1].values - T_xarr.tloc[0].values Tday = T_xarr.sel(tloc=np.arange(6, 18, tres)) @@ -135,6 +142,10 @@ def smooth_daytime(T_xarr, savefile=None): def add_wls_diviner(xarr, bands=("t3", "t4", "t5", "t6", "t7", "t8", "t9")): """ Return xarr with wavelength coordinate. + + Parameters: + xarr (xr.DataArray): Diviner lev4 hourly temperature data + bands (tuple): bands to include in output """ wl = get_diviner_wls().values out = xarr.sel(band=slice("t3", "t9")).assign_coords( @@ -147,6 +158,9 @@ def add_wls_diviner(xarr, bands=("t3", "t4", "t5", "t6", "t7", "t8", "t9")): def get_diviner_wls(bres=1): """ Return diviner c3-c9 wavelength arr with band_res values in each bandpass. + + Parameters: + bres (int): number of wavelength values in each bandpass """ if bres == 1: wls = np.mean((DIV_WL_MIN, DIV_WL_MAX), axis=0) @@ -161,6 +175,10 @@ def get_diviner_wls(bres=1): def div_integrated_rad(emission, units="W/m^2/sr/um"): """ Return integrated radiance of diviner bands. + + Parameters: + emission (xr.DataArray): Diviner radiated emission + units (str): units of output """ out = [] for dwlmin, dwlmax in zip(DIV_WL_MIN, DIV_WL_MAX): @@ -179,10 +197,19 @@ def load_div_lev4( smoothday=False, invert_y=False, load_cached=True, - divdir=cfg.DATA_DIR_DIVINER, + divdir=DATA_DIR_DIVINER / "lev4", ): """ Return Diviner lev4 data as xarray. + + Parameters: + roi (str): region of interest + ext (list): extent of ROI + savefile (str): Path to save result to netCDF (.nc) file + smoothday (bool): smooth daytime temperatures with 2nd order polyfit + invert_y (bool): invert y axis (sometimes Diviner data inverted) + load_cached (bool): load cached data if available + divdir (str): path to Diviner data directory """ roi_str = roi.replace("'", "_").lower() savepath = Path(divdir) / roi_str / savefile @@ -210,7 +237,15 @@ def load_div_lev4( def div_tbol(divbt, wl1=DIV_WL1, wl2=DIV_WL2, tmin=DIV_TMIN, iters=ITERS): """ - Return tbol from Diviner C3-C9 brightness temperatures. + Return bolometric temperature from Diviner C3-C9 brightness temperatures + (s.o.m, Paige et al., 2010, in Science). + + Parameters: + divbt (np.array): Diviner brightness temperatures + wl1 (np.array): lower wavelength bound of each band + wl2 (np.array): upper wavelength bound of each band + tmin (np.array): minimum temperature of each band + iters (np.array): number of iterations for each band """ def f(wl, bt, n_iter): @@ -230,6 +265,12 @@ def cfst(bt, wl1, wl2, n_iter=100): Return fraction of Planck radiance between wl1 and wl2 at bt. Translated from JB, via DAP via Houghton 'Physics of Atmospheres' + + Parameters: + bt (float): brightness temperature + wl1 (float): lower wavelength bound + wl2 (float): upper wavelength bound + n_iter (int): number of iterations """ return 1.53989733e-1 * (f(wl2, bt, n_iter) - f(wl1, bt, n_iter)) @@ -261,6 +302,12 @@ def load_div_filters( ): """ Return Diviner filter functions + + Parameters: + fdiv_filt (str): path to Diviner filter functions + scale (bool): scale filter functions by FILT_SCALE + wlunits (bool): convert wavenumber to wavelength + bands (tuple): bands to include in output """ hdf = xr.load_dataset(fdiv_filt) wn = hdf.xaxis.values.squeeze() @@ -284,6 +331,9 @@ def load_div_filters( def load_div_t2r(fdiv_t2r=FDIV_T2R): """ Return Diviner temperature to radiance lookup table. + + Parameters: + fdiv_t2r (str): path to Diviner t2r lookup table """ return pd.read_csv(fdiv_t2r, index_col=0, header=0, delim_whitespace=True) @@ -291,6 +341,10 @@ def load_div_t2r(fdiv_t2r=FDIV_T2R): def divfilt_rad(wnrad, div_filt=None): """ Return radiance of Diviner bands (convolve rad with Diviner filter funcs). + + Parameters: + wnrad (xr.DataArray): Diviner radiance in wavenumber space + div_filt (xr.DataArray): Diviner filter functions """ if div_filt is None: div_filt = load_div_filters() @@ -301,6 +355,10 @@ def divfilt_rad(wnrad, div_filt=None): def divrad2bt(divrad, fdiv_t2r=FDIV_T2R): """ Return brightness temperatures from Diviner radiance (e.g. from div_filt). + + Parameters: + divrad (xr.DataArray): Diviner radiance + fdiv_t2r (str): path to Diviner t2r lookup table """ t2r = load_div_t2r(fdiv_t2r) # cached lookup rad = divrad.values @@ -323,7 +381,14 @@ def divrad2bt(divrad, fdiv_t2r=FDIV_T2R): def emission2tbol_xr(emission, wls=None, div_filt=None): - """Return tbol from input emission array at wls.""" + """ + Return tbol from input emission array at wls. + + Parameters: + emission (xr.DataArray): Diviner radiated emission + wls (xr.DataArray): wavelengths to integrate over + div_filt (xr.DataArray): Diviner filter functions + """ if wls is None: wls = emission.wavelength elif not isinstance(wls, xr.DataArray):