From 0fb2aa3bdc502e0fec29b57d3e8ef08abfc750ef Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Juan=20Manuel=20Franco=20Pati=C3=B1o?= <45358012+juafrapat@users.noreply.github.com> Date: Thu, 5 Dec 2019 08:56:33 +0100 Subject: [PATCH] Add files via upload --- 232Th.exp | 336 ++++++++++++++++++++++++++++++++++++++++++++++++++ 232Th.inp | 62 +++++----- 233U.exp | 173 ++++++++++++++++++++++++++ 233U.inp | 59 ++++----- 235U.exp | 172 ++++++++++++++++++++++++++ 235U.inp | 50 ++++---- 238U.exp | 253 +++++++++++++++++++++++++++++++++++++ 238U.inp | 65 +++++----- 239Pu.exp | 114 +++++++++++++++++ 239Pu.inp | 60 ++++----- exp.py | 38 ++++++ frescoPRC.f90 | 52 +++----- fun.f95 | 13 +- runall.sh | 10 +- script.py | 62 ++++++++-- 15 files changed, 1318 insertions(+), 201 deletions(-) create mode 100644 232Th.exp create mode 100644 233U.exp create mode 100644 235U.exp create mode 100644 238U.exp create mode 100644 239Pu.exp create mode 100644 exp.py diff --git a/232Th.exp b/232Th.exp new file mode 100644 index 0000000..4219de2 --- /dev/null +++ b/232Th.exp @@ -0,0 +1,336 @@ + 0.0496 0 13.613 0.2602 #dY:1.9% + 0.0548 0 12.76 0.2342 #dY:1.8% + 0.0605 0 13.224 0.2095 #dY:1.6% + 0.0669 0 13.178 0.1894 #dY:1.4% + 0.074 0 12.73 0.1727 #dY:1.4% + 0.0817 0 12.926 0.1621 #dY:1.3% + 0.0903 0 12.375 0.1503 #dY:1.2% + 0.0998 0 12.268 0.135 #dY:1.1% + 0.1104 0 11.921 0.1247 #dY:1.0% + 0.122 0 12.104 0.1164 #dY:1.0% + 0.1348 0 11.529 0.1099 #dY:1.0% + 0.149 0 11.352 0.106 #dY:0.9% + 0.1647 0 11.189 0.1003 #dY:0.9% + 0.182 0 11.106 0.096 #dY:0.9% + 0.2011 0 10.778 0.0943 #dY:0.9% + 0.2223 0 10.528 0.0923 #dY:0.9% + 0.2457 0 10.416 0.0915 #dY:0.9% + 0.298 0 9.89 0.1157 #dY:1.2% + 0.3317 0 9.3511 0.0932 #dY:1.0% + 0.358 0 9.418 0.1379 #dY:1.5% + 0.447 0 8.771 0.1226 #dY:1.4% + 0.5 0 8.462 0.0725 #dY:0.9% + 0.54 0 8.023 0.0975 #dY:1.2% + 0.625 0 7.947 0.1016 #dY:1.3% + 0.65 0 7.8165 0.0555 #dY:0.7% + 0.7 0 7.7133 0.0532 #dY:0.7% + 0.75 00 7.4018 0.0541 #dY:0.7% + 0.8 0 7.1835 0.0583 #dY:0.8% + 0.85 0 7.0725 0.0576 #dY:0.8% + 0.9 0 7.0855 0.0581 #dY:0.8% + 0.95 0 6.9535 0.0665 #dY:1.0% + 1 0 7.0827 0.0728 #dY:1.0% + 1.1 0 6.9405 0.0461 #dY:0.7% + 1.2 00 6.8694 0.0442 #dY:0.6% + 1.3 0 6.8833 0.0472 #dY:0.7% + 1.4 0 6.8391 0.0457 #dY:0.7% + 1.5 0 6.882 0.0448 #dY:0.7% + 1.6 0 6.862 0.0455 #dY:0.7% + 1.7 0 6.88 0.044 #dY:0.6% + 1.8 0 6.915 0.045 #dY:0.7% + 1.9 0 6.9853 0.0474 #dY:0.7% + 2 0 7.1262 0.0467 #dY:0.7% + 2.1 0 7.1405 0.0472 #dY:0.7% + 2.2 0 7.2176 0.0464 #dY:0.6% + 2.3 0 7.1759 0.0461 #dY:0.6% + 2.4 0 7.3234 0.0476 #dY:0.6% + 2.5 0 7.4239 0.0491 #dY:0.7% + 2.6 0 7.4965 0.0492 #dY:0.7% + 2.7 0 7.5261 0.0493 #dY:0.7% + 2.8 0 7.5515 0.05 #dY:0.7% + 2.9 0 7.5997 0.051 #dY:0.7% + 3 00 7.6597 0.05 #dY:0.7% + 3.1 0 7.7103 0.0516 #dY:0.7% + 3.2 00 7.7249 0.0531 #dY:0.7% + 3.3 0 7.8685 0.0539 #dY:0.7% + 3.4 0 7.921 0.0544 #dY:0.7% + 3.5 0 7.8328 0.0547 #dY:0.7% + 3.6 0 7.8683 0.055 #dY:0.7% + 3.7 0 7.8076 0.055 #dY:0.7% + 3.8 0 7.7584 0.0551 #dY:0.7% + 3.9 00 7.8143 0.0544 #dY:0.7% + 4 0 7.7437 0.0543 #dY:0.7% + 4.2 0 7.7965 0.0389 #dY:0.5% + 4.4 0 7.7004 0.0394 #dY:0.5% + 4.6 0 7.6114 0.0396 #dY:0.5% + 4.8 00 7.5538 0.0402 #dY:0.5% + 5 0 7.5 0.0416 #dY:0.6% + 5.2 0 7.3675 0.0438 #dY:0.6% + 5.4 0 7.2042 0.0459 #dY:0.6% + 5.6 0 7.2223 0.0503 #dY:0.7% + 5.8 0 7.0443 0.0537 #dY:0.8% + 6 0 6.962 0.0602 #dY:0.9% + 6.2 00 6.8648 0.0632 #dY:0.9% + 6.4 0 6.8016 0.0668 #dY:1.0% + 6.6 0 6.6587 0.0697 #dY:1.0% + 6.8 0 6.6351 0.0731 #dY:1.1% + 7 0 6.593 0.0762 #dY:1.2% + 7.5 0 6.29 0.0536 #dY:0.9% + 8 0 6.0349 0.0575 #dY:1.0% + 8.5 00 5.916 0.0601 #dY:1.0% + 9 0 5.7913 0.0608 #dY:1.0% + 9.5 0 5.6814 0.0604 #dY:1.1% + 10 0 5.5842 0.0583 #dY:1.0% + 10.14 0 5.764 0.01703 #dY:0.3% + 10.24 0 5.752 0.01702 #dY:0.3% + 10.34 00 5.751 0.01707 #dY:0.3% + 10.45 0 5.738 0.01725 #dY:0.3% + 10.55 0 5.731 0.01751 #dY:0.3% + 10.66 0 5.719 0.01781 #dY:0.3% + 10.77 00 5.738 0.01814 #dY:0.3% + 10.87 0 5.719 0.0184 #dY:0.3% + 10.98 0 5.71 0.01859 #dY:0.3% + 11.09 0 5.71 0.01867 #dY:0.3% + 11.2 0 5.704 0.0187 #dY:0.3% + 11.32 00 5.699 0.0187 #dY:0.3% + 11.43 0 5.715 0.01876 #dY:0.3% + 11.55 0 5.727 0.01878 #dY:0.3% + 11.66 0 5.687 0.01873 #dY:0.3% + 11.78 0 5.68 0.01881 #dY:0.3% + 11.9 00 5.704 0.01915 #dY:0.3% + 12.02 0 5.707 0.01951 #dY:0.3% + 12.14 0 5.691 0.01965 #dY:0.3% + 12.26 0 5.696 0.01942 #dY:0.3% + 12.38 0 5.708 0.01929 #dY:0.3% + 12.51 0 5.718 0.01936 #dY:0.3% + 12.63 00 5.708 0.01943 #dY:0.3% + 12.76 0 5.712 0.01959 #dY:0.3% + 12.89 0 5.729 0.01989 #dY:0.3% + 13.02 0 5.73 0.02007 #dY:0.4% + 13.15 0 5.722 0.02006 #dY:0.4% + 13.28 0 5.721 0.02005 #dY:0.4% + 13.41 00 5.763 0.02022 #dY:0.4% + 13.55 0 5.78 0.02031 #dY:0.4% + 13.69 0 5.78 0.02024 #dY:0.4% + 13.82 0 5.773 0.0201 #dY:0.3% + 13.96 0 5.774 0.02004 #dY:0.3% + 14.1 0 5.811 0.02019 #dY:0.3% + 14.24 0 5.826 0.0203 #dY:0.3% + 14.39 00 5.837 0.0204 #dY:0.3% + 14.53 0 5.865 0.02059 #dY:0.4% + 14.68 0 5.882 0.02092 #dY:0.4% + 14.82 0 5.877 0.02113 #dY:0.4% + 14.97 0 5.877 0.02123 #dY:0.4% + 15.12 0 5.887 0.01898 #dY:0.3% + 15.28 0 5.907 0.01377 #dY:0.2% + 15.43 00 5.924 0.01387 #dY:0.2% + 15.58 0 5.936 0.01401 #dY:0.2% + 15.74 0 5.943 0.01415 #dY:0.2% + 15.9 0 5.974 0.01426 #dY:0.2% + 16.06 0 5.986 0.01429 #dY:0.2% + 16.22 0 6 0.01422 #dY:0.2% + 16.38 0 6.005 0.01414 #dY:0.2% + 16.55 0 6.001 0.01409 #dY:0.2% + 16.72 0 6.021 0.01403 #dY:0.2% + 16.88 0 6.05 0.01404 #dY:0.2% + 17.05 00 6.054 0.01406 #dY:0.2% + 17.22 0 6.054 0.01406 #dY:0.2% + 17.4 0 6.088 0.01407 #dY:0.2% + 17.57 0 6.096 0.01407 #dY:0.2% + 17.75 0 6.104 0.01409 #dY:0.2% + 17.93 00 6.108 0.01412 #dY:0.2% + 18.11 0 6.113 0.01416 #dY:0.2% + 18.29 0 6.117 0.01419 #dY:0.2% + 18.47 0 6.144 0.01422 #dY:0.2% + 18.66 0 6.154 0.01425 #dY:0.2% + 18.85 0 6.162 0.01427 #dY:0.2% + 19.04 00 6.192 0.01434 #dY:0.2% + 19.23 0 6.187 0.01436 #dY:0.2% + 19.42 0 6.188 0.01438 #dY:0.2% + 19.62 0 6.187 0.01434 #dY:0.2% + 19.81 0 6.198 0.01426 #dY:0.2% + 20.01 0 6.202 0.01416 #dY:0.2% + 20.21 00 6.201 0.01406 #dY:0.2% + 20.42 0 6.206 0.01395 #dY:0.2% + 20.62 0 6.212 0.01384 #dY:0.2% + 20.83 0 6.2 0.01377 #dY:0.2% + 21.04 0 6.198 0.01372 #dY:0.2% + 21.25 0 6.196 0.01364 #dY:0.2% + 21.46 0 6.196 0.01358 #dY:0.2% + 21.68 0 6.189 0.01351 #dY:0.2% + 21.9 0 6.184 0.01342 #dY:0.2% + 22.12 0 6.207 0.01333 #dY:0.2% + 22.34 0 6.2 0.01325 #dY:0.2% + 22.56 0 6.188 0.01317 #dY:0.2% + 22.79 0 6.171 0.01311 #dY:0.2% + 23.02 0 6.152 0.01305 #dY:0.2% + 23.25 0 6.156 0.01304 #dY:0.2% + 23.48 0 6.142 0.01299 #dY:0.2% + 23.72 0 6.141 0.01293 #dY:0.2% + 23.96 0 6.131 0.01286 #dY:0.2% + 24.2 0 6.108 0.0128 #dY:0.2% + 24.44 0 6.086 0.01273 #dY:0.2% + 24.69 0 6.081 0.01265 #dY:0.2% + 24.94 0 6.055 0.01257 #dY:0.2% + 25.19 0 6.048 0.01251 #dY:0.2% + 25.44 0 6.038 0.01242 #dY:0.2% + 25.7 00 6.007 0.01232 #dY:0.2% + 25.95 0 5.985 0.01222 #dY:0.2% + 26.21 0 5.989 0.01216 #dY:0.2% + 26.48 00 5.974 0.01207 #dY:0.2% + 26.74 0 5.954 0.01199 #dY:0.2% + 27.01 0 5.927 0.01188 #dY:0.2% + 27.28 0 5.909 0.01179 #dY:0.2% + 27.56 0 5.877 0.0117 #dY:0.2% + 27.84 00 5.85 0.01161 #dY:0.2% + 28.12 0 5.829 0.01151 #dY:0.2% + 28.4 0 5.803 0.01141 #dY:0.2% + 28.68 0 5.791 0.01133 #dY:0.2% + 28.97 0 5.757 0.01123 #dY:0.2% + 29.26 0 5.734 0.01114 #dY:0.2% + 29.56 0 5.706 0.01105 #dY:0.2% + 29.85 0 5.678 0.01097 #dY:0.2% + 30.15 00 5.65 0.01089 #dY:0.2% + 1.817 0.009085 7.096 0.06732 #dY:0.9% + 1.854 0.009272 7.131 0.06765 #dY:0.9% + 1.893 0.009465 7.158 0.06791 #dY:0.9% + 1.933 0.009665 7.205 0.06835 #dY:0.9% + 1.974 0.009871 7.204 0.06321 #dY:0.9% + 2.017 0.01008 7.221 0.06336 #dY:0.9% + 2.061 0.0103 7.321 0.06424 #dY:0.9% + 2.107 0.01053 7.36 0.06458 #dY:0.9% + 2.154 0.01077 7.344 0.05966 #dY:0.8% + 2.203 0.01101 7.39 0.06004 #dY:0.8% + 2.253 0.01127 7.482 0.06078 #dY:0.8% + 2.306 0.01153 7.541 0.06126 #dY:0.8% + 2.36 0.0118 7.529 0.06117 #dY:0.8% + 2.417 0.01208 7.528 0.06116 #dY:0.8% + 2.475 0.01238 7.689 0.06247 #dY:0.8% + 2.536 0.01268 7.668 0.0623 #dY:0.8% + 2.582 0.01291 7.778 0.06508 #dY:0.8% + 2.631 0.01315 7.747 0.06294 #dY:0.8% + 2.698 0.01349 7.787 0.06326 #dY:0.8% + 2.767 0.01384 7.873 0.05944 #dY:0.8% + 2.84 0.0142 7.891 0.05958 #dY:0.8% + 2.915 0.01457 7.943 0.05997 #dY:0.8% + 2.993 0.01497 7.948 0.06001 #dY:0.8% + 3.075 0.01537 7.962 0.06011 #dY:0.8% + 3.16 0.0158 7.949 0.06001 #dY:0.8% + 3.249 0.01624 8.021 0.06056 #dY:0.8% + 3.341 0.01671 8.004 0.06043 #dY:0.8% + 3.438 0.01719 8.039 0.06069 #dY:0.8% + 3.539 0.0177 8.079 0.061 #dY:0.8% + 3.645 0.01822 8.025 0.06059 #dY:0.8% + 3.756 0.01878 8.03 0.06063 #dY:0.8% + 3.872 0.01936 8.074 0.06096 #dY:0.8% + 3.993 0.01997 8.022 0.06056 #dY:0.8% + 4.121 0.0206 8.021 0.06516 #dY:0.8% + 4.254 0.02127 8.016 0.07034 #dY:0.9% + 4.395 0.02197 7.998 0.07588 #dY:0.9% + 4.542 0.02271 7.901 0.07496 #dY:0.9% + 4.698 0.02349 7.866 0.07462 #dY:0.9% + 4.862 0.02431 7.84 0.07438 #dY:0.9% + 5.034 0.02517 7.691 0.07881 #dY:1.0% + 5.216 0.02608 7.659 0.07848 #dY:1.0% + 5.408 0.02704 7.552 0.07738 #dY:1.0% + 5.611 0.02806 7.411 0.07594 #dY:1.0% + 5.826 0.02913 7.338 0.07519 #dY:1.0% + 6.054 0.03027 7.149 0.07896 #dY:1.1% + 6.296 0.03148 7.107 0.07882 #dY:1.1% + 6.552 0.03276 6.938 0.08268 #dY:1.2% + 6.825 0.03413 6.827 0.08135 #dY:1.2% + 7.115 0.03558 6.59 0.07853 #dY:1.2% + 7.425 0.03712 6.482 0.08276 #dY:1.3% + 7.755 0.03878 6.254 0.08529 #dY:1.4% + 8.108 0.04054 6.35 0.0866 #dY:1.4% + 8.486 0.04243 6.151 0.08389 #dY:1.4% + 8.891 0.04446 6.155 0.08941 #dY:1.5% + 9.481 0.04741 6.025 0.08478 #dY:1.4% + 10.04 0.1004 5.948 0.09234 #dY:1.6% + 10.56 0.1056 5.916 0.09184 #dY:1.6% + 11.12 0.1112 5.936 0.08684 #dY:1.5% + 11.74 0.1174 5.768 0.07867 #dY:1.4% + 12.4 0.124 5.868 0.07492 #dY:1.3% + 13.12 0.1312 5.815 0.06449 #dY:1.1% + 13.91 0.1391 5.933 0.06108 #dY:1.0% + 14.77 0.1477 6.012 0.05735 #dY:1.0% + 15.71 0.1571 6.038 0.0576 #dY:1.0% + 16.75 0.1675 6.109 0.05828 #dY:1.0% + 17.89 0.1789 6.26 0.05972 #dY:1.0% + 19.15 0.1915 6.264 0.07464 #dY:1.2% + 20.25 0.2025 6.437 0.1426 #dY:2.2% + 0.048 0.008022 13.73 0.2912 #dY:2.1% + 0.063 0.00805 13.25 0.2763 #dY:2.1% + 0.078 0.008022 12.87 0.27 #dY:2.1% + 0.092 0.00804 12.64 0.2659 #dY:2.1% + 0.111 0.008062 12.19 0.2617 #dY:2.1% + 0.137 0.008122 11.78 0.2589 #dY:2.2% + 0.172 0.008246 11.26 0.2555 #dY:2.3% + 0.195 0.008476 10.93 0.1027 #dY:0.9% + 0.223 0.008509 10.64 0.2533 #dY:2.4% + 0.244 0.008443 10.51 0.09984 #dY:0.9% + 0.297 0.008443 9.758 0.08919 #dY:0.9% + 0.35 0.008476 9.454 0.09525 #dY:1.0% + 0.401 0.008443 9.133 0.09058 #dY:1.0% + 0.458 0.008443 8.62 0.06917 #dY:0.8% + 0.51 0.008693 8.34 0.05787 #dY:0.7% + 0.568 0.008654 8.115 0.1035 #dY:1.3% + 0.61 0.008693 7.938 0.05886 #dY:0.7% + 0.67 0.008616 7.693 0.127 #dY:1.7% + 0.715 0.00858 7.584 0.07179 #dY:0.9% + 0.773 0.008616 7.341 0.1083 #dY:1.5% + 0.821 0.008654 7.347 0.06515 #dY:0.9% + 0.924 0.008693 7.163 0.0662 #dY:0.9% + 1.023 0.008544 7.104 0.05399 #dY:0.8% + 1.128 0.008509 6.973 0.06472 #dY:0.9% + 1.23 0.008476 6.869 0.07216 #dY:1.1% + 1.333 0.008443 6.875 0.09678 #dY:1.4% + 1.533 0.008246 6.977 0.05115 #dY:0.7% + 1.631 0.008179 6.867 0.0677 #dY:1.0% + 1.734 0.008223 7.086 0.1439 #dY:2.0% + 1.825 0.008158 7.168 0.07103 #dY:1.0% + 1.931 0.008223 7.233 0.05567 #dY:0.8% + 2.025 0.008139 7.193 0.07193 #dY:1.0% + 2.134 0.0105 7.357 0.06475 #dY:0.9% + 2.279 0.009381 7.435 0.06495 #dY:0.9% + 2.429 0.01124 7.547 0.06527 #dY:0.9% + 2.574 0.009881 7.663 0.06559 #dY:0.9% + 2.791 0.009035 7.719 0.07237 #dY:0.9% + 2.994 0.009381 7.9 0.07201 #dY:0.9% + 3.188 0.009278 8.023 0.07486 #dY:0.9% + 3.387 0.009228 8.063 0.07496 #dY:0.9% + 3.585 0.00913 8.09 0.07757 #dY:1.0% + 3.781 0.008989 8.017 0.07995 #dY:1.0% + 3.974 0.0089 8.064 0.08355 #dY:1.0% + 4.185 0.008773 8.08 0.09067 #dY:1.1% + 4.389 0.008693 8.099 0.1034 #dY:1.3% + 4.598 0.008297 7.828 0.1855 #dY:2.4% + 4.807 0.008297 7.754 0.1037 #dY:1.3% + 0.0065 0.0005 15.17 0.01 + 0.0075 0.0005 16.62 0.01 + 0.0085 0.0005 15.62 0.01 + 0.0095 0.0005 15.5 0.01 + 0.0125 0.0025 14.84 0.01 + 0.0175 0.0025 14.57 0.01 + 0.025 0.005 14.58 0.01 + 0.035 0.005 14.25 0.01 + 0.045 0.005 13.97 0.01 + 0.055 0.005 13.36 0.01 + 0.065 0.005 13.34 0.01 + 0.075 0.005 12.66 0.01 + 0.085 0.005 12.5 0.01 + 0.095 0.005 12.1 0.01 + 0.105 0.005 11.72 0.01 + 0.115 0.005 11.95 0.01 + 0.125 0.005 11.89 0.01 + 0.135 0.005 11.52 0.01 + 0.145 0.005 11.38 0.01 + 0.175 0.025 11.03 0.01 + 0.25 0.05 10.21 0.01 + 0.35 0.05 9.37 0.01 + 0.45 0.05 9.07 0.01 + 0.55 0.05 8.08 0.01 + 0.65 0.05 7.66 0.01 + 0.75 0.05 7.39 0.01 + 0.85 0.05 7.24 0.01 + 0.95 0.05 7.21 0.01 diff --git a/232Th.inp b/232Th.inp index b6fc407..c06cb39 100644 --- a/232Th.inp +++ b/232Th.inp @@ -1,45 +1,47 @@ &Target Z=90.0 A=232.0 eferm=-5.613395 - BETA20=0.211 BETA40=0.063 BETA60=1.8000E-03 + BETA20=0.201 BETA40=0.074 BETA60=3.0000E-05 nstat=18 / -&State Et=0.0000E+00 Jval=0.0 Kval=0.0 NBAND=1 COEFF=0.2110 / -&State Et=4.9400E+01 Jval=2.0 Kval=0.0 NBAND=1 COEFF=0.2110 / -&State Et=1.6210E+02 Jval=4.0 Kval=0.0 NBAND=1 COEFF=0.2110 / -&State Et=3.3330E+02 Jval=6.0 Kval=0.0 NBAND=1 COEFF=0.2110 / -&State Et=5.5690E+02 Jval=8.0 Kval=0.0 NBAND=1 COEFF=0.2110 / -&State Et=8.2680E+02 Jval=10.0 Kval=0.0 NBAND=1 COEFF=0.211 / +&State Et=0.0000E+00 Jval=0.0 Kval=0.0 NBAND=1 COEFF=0.2010 / +&State Et=4.9400E+01 Jval=2.0 Kval=0.0 NBAND=1 COEFF=0.2010 / +&State Et=1.6210E+02 Jval=4.0 Kval=0.0 NBAND=1 COEFF=0.2010 / +&State Et=3.3330E+02 Jval=6.0 Kval=0.0 NBAND=1 COEFF=0.2010 / +&State Et=5.5690E+02 Jval=8.0 Kval=0.0 NBAND=1 COEFF=0.2010 / +&State Et=8.2680E+02 Jval=10.0 Kval=0.0 NBAND=1 COEFF=0.201 / -&State Et=7.1440E+02 Jval=1.0 Kval=0.0 NBAND=-2 COEFF=0.047 / -&State Et=7.7440E+02 Jval=3.0 Kval=0.0 NBAND=-2 COEFF=0.047 / -&State Et=8.8380E+02 Jval=5.0 Kval=0.0 NBAND=-2 COEFF=0.047 / -&State Et=1.0429E+03 Jval=7.0 Kval=0.0 NBAND=-2 COEFF=0.047 / -&State Et=1.2496E+03 Jval=9.0 Kval=0.0 NBAND=-2 COEFF=0.047 / +&State Et=7.1440E+02 Jval=1.0 Kval=0.0 NBAND=-2 COEFF=0.054069 / +&State Et=7.7440E+02 Jval=3.0 Kval=0.0 NBAND=-2 COEFF=0.054069 / +&State Et=8.8380E+02 Jval=5.0 Kval=0.0 NBAND=-2 COEFF=0.054069 / +&State Et=1.0429E+03 Jval=7.0 Kval=0.0 NBAND=-2 COEFF=0.054069 / +&State Et=1.2496E+03 Jval=9.0 Kval=0.0 NBAND=-2 COEFF=0.054069 / -&State Et=7.3060E+02 Jval=0.0 Kval=0.0 NBAND=3 COEFF=0.0100 / -&State Et=7.7420E+02 Jval=2.0 Kval=0.0 NBAND=3 COEFF=0.0100 / +&State Et=7.3060E+02 Jval=0.0 Kval=0.0 NBAND=3 COEFF=0.001005 / +&State Et=7.7420E+02 Jval=2.0 Kval=0.0 NBAND=3 COEFF=0.001005 / +&State Et=8.7390E+02 Jval=4.0 Kval=0.0 NBAND=3 COEFF=0.001005 / -&State Et=1.0786E+03 Jval=0.0 Kval=0.0 NBAND=4 COEFF=0.0080 / -&State Et=1.1217E+03 Jval=2.0 Kval=0.0 NBAND=4 COEFF=0.0080 / -&State Et=1.2221E+03 Jval=4.0 Kval=0.0 NBAND=4 COEFF=0.0080 / +&State Et=1.0786E+03 Jval=0.0 Kval=0.0 NBAND=4 COEFF=0.024924 / +&State Et=1.1217E+03 Jval=2.0 Kval=0.0 NBAND=4 COEFF=0.024924 / -&State Et=7.8530E+02 Jval=2.0 Kval=2.0 NBAND=5 COEFF=0.0080 / -&State Et=8.2960E+02 Jval=3.0 Kval=2.0 NBAND=5 COEFF=0.0080 / +!&State Et=1.2221E+03 Jval=4.0 Kval=0.0 NBAND=4 COEFF=0.024924 / -&Fresco hcm=0.10 rmatch= 20.00 Jmax=20 - Ngrid=41 pythonFlag=0 Nenergy=35 absend= 0.001 - relativ= 2 / +&State Et=7.8530E+02 Jval=2.0 Kval=2.0 NBAND=5 COEFF=0.0017085 / +&State Et=8.2960E+02 Jval=3.0 Kval=2.0 NBAND=5 COEFF=0.0017085 / + +&Fresco hcm=0.10 rmatch= 20.00 Jmax=30 + Ngrid=41 pythonFlag=0 Nenergy=37 absend= 0.001 + relativ= 3 / &Energies elab=0.001,0.002,0.004,0.006,0.008,0.01,0.02,0.04, - 0.06,0.08,0.100,0.200,0.400,0.600,0.800,1.0000,1.40,2.0000, + 0.06,0.08,0.100,0.150,0.200,0.30,0.400,0.600,0.800,1.0000,1.40,2.0000, 2.20,2.70,3.20,3.70,4.2,4.8,5.00, 7.00,10.0000,12.00,15.00,17.00,20.00,23.00,26.00,28.00,30.00 / -&Potential v0a=50.41 v0b=0.0292 lambdhf= 0.00977 - cviso=17.50 vspo=5.64 lambdso=0.005 ccoul=1.36 - av=11.91 bv=81.86 w0=17.70 bs=10.85 wspo=-3.1 - bso=160.0 ea=55.0 alphav=0.355 cs=0.01331 cwiso=29.0 - adv=238.0 rhfa=1.2490 rhfb=-0.00171 ahfa=0.638 ahfb=0.00219 - rv=1.2657 ava=0.6960 avb=-0.00021 rsa=1.1701 rsb=0.0041 - as=0.617 rso=1.1214 aso=0.59 rc=1.1289 ac=0.547 / +&Potential v0a=50.59 v0b=0.0608 lambdhf= 0.00982 + cviso=16.79 vspo=6.83 lambdso=0.005 ccoul=1.43 + av=11.66 bv=81.61 w0a=17.30 w0b=-0.031 bs=10.55 wspo=-3.1 + bso=160.0 ea=55.0 alphav=0.355 cs=0.01343 cwiso=28.96 + adv=238.0 rhfa=1.2426 rhfb=-0.00246 ahfa=0.638 ahfb=0.00306 + rv=1.2704 avv=0.6970 rsa=1.1860 rsb=0.0060 + as=0.616 rso=1.0156 aso=0.549 rc=1.2262 ac=0.674 / / diff --git a/233U.exp b/233U.exp new file mode 100644 index 0000000..acfb527 --- /dev/null +++ b/233U.exp @@ -0,0 +1,173 @@ + 0.058 0.004 13.05 0.1536 #dY:1.2% + 0.063 0.004 13.14 0.1528 #dY:1.2% + 0.074 0.004 12.72 0.1535 #dY:1.2% + 0.103 0.002 12 0.136 #dY:1.1% + 0.122 0.003 11.81 0.1415 #dY:1.2% + 0.139 0.004 11.39 0.1403 #dY:1.2% + 0.147 0.005 11.22 0.1307 #dY:1.2% + 0.166 0.005 11.11 0.1303 #dY:1.2% + 0.188 0.005 11.05 0.1324 #dY:1.2% + 0.193 0.005 10.68 0.1311 #dY:1.2% + 0.228 0.005 10.21 0.1408 #dY:1.4% + 0.24 0.005 10.36 0.1234 #dY:1.2% + 0.297 0.014 9.84 0.1327 #dY:1.3% + 0.328 0.014 9.62 0.1402 #dY:1.5% + 0.352 0.005 9.31 0.1288 #dY:1.4% + 0.383 0.013 9.232 0.1141 #dY:1.2% + 0.435 0.011 8.763 0.104 #dY:1.2% + 0.47 0.011 8.466 0.1469 #dY:1.7% + 0.495 0.005 8.193 0.09758 #dY:1.2% + 0.534 0.011 8.06 0.1136 #dY:1.4% + 0.5855 0.015 7.759 0.1097 #dY:1.4% + 0.6343 0.015 7.628 0.1079 #dY:1.4% + 0.6889 0.015 7.569 0.1136 #dY:1.5% + 0.741 0.015 7.228 0.105 #dY:1.5% + 0.7889 0.015 7.184 0.1016 #dY:1.4% + 0.8421 0.015 6.982 0.09874 #dY:1.4% + 0.8928 0.015 6.858 0.09962 #dY:1.5% + 0.9505 0.015 6.77 0.1035 #dY:1.5% + 0.997 0.01 6.902 0.07928 #dY:1.1% + 1.201 0.015 6.741 0.1011 #dY:1.5% + 1.275 0.015 6.622 0.1139 #dY:1.7% + 1.348 0.015 6.767 0.1182 #dY:1.7% + 1.394 0.015 6.649 0.1023 #dY:1.5% + 1.449 0.015 6.743 0.1188 #dY:1.8% + 1.5 0.015 6.842 0.1039 #dY:1.5% + 1.55 0.015 6.785 0.1215 #dY:1.8% + 1.642 0.015 6.857 0.09697 #dY:1.4% + 1.768 0.015 6.902 0.09761 #dY:1.4% + 1.87 0.015 7.097 0.1013 #dY:1.4% + 1.976 0.015 7.052 0.09973 #dY:1.4% + 2.069 0.015 7.255 0.1026 #dY:1.4% + 2.2 0.015 7.39 0.1045 #dY:1.4% + 2.334 0.015 7.512 0.1062 #dY:1.4% + 2.467 0.015 7.498 0.106 #dY:1.4% + 2.568 0.015 7.517 0.1063 #dY:1.4% + 2.676 0.015 7.743 0.1095 #dY:1.4% + 2.754 0.015 7.609 0.1076 #dY:1.4% + 2.848 0.015 7.832 0.1108 #dY:1.4% + 2.964 0.015 7.978 0.1128 #dY:1.4% + 3.147 0.015 7.897 0.1117 #dY:1.4% + 3.299 0.015 7.948 0.1124 #dY:1.4% + 3.444 0.015 8.096 0.1145 #dY:1.4% + 3.55 0.015 8.028 0.1135 #dY:1.4% + 3.703 0.015 8.099 0.1145 #dY:1.4% + 3.854 0.015 8.109 0.1147 #dY:1.4% + 3.972 0.015 8.02 0.1193 #dY:1.5% + 4.054 0.015 8.041 0.1137 #dY:1.4% + 4.202 0.015 8.058 0.114 #dY:1.4% + 4.326 0.015 8.095 0.1409 #dY:1.7% + 4.432 0.015 7.972 0.1127 #dY:1.4% + 1.817 0.009085 7.053 0.1135 #dY:1.6% + 1.854 0.009272 7.052 0.1135 #dY:1.6% + 1.893 0.009465 7.099 0.1142 #dY:1.6% + 1.933 0.009665 7.15 0.1151 #dY:1.6% + 1.974 0.009871 7.108 0.1144 #dY:1.6% + 2.017 0.01008 7.177 0.1155 #dY:1.6% + 2.061 0.0103 7.26 0.1168 #dY:1.6% + 2.107 0.01053 7.284 0.1152 #dY:1.6% + 2.154 0.01077 7.365 0.1165 #dY:1.6% + 2.203 0.01101 7.396 0.1169 #dY:1.6% + 2.253 0.01127 7.465 0.118 #dY:1.6% + 2.306 0.01153 7.489 0.1184 #dY:1.6% + 2.36 0.0118 7.482 0.1183 #dY:1.6% + 2.417 0.01208 7.604 0.1202 #dY:1.6% + 2.475 0.01238 7.666 0.1212 #dY:1.6% + 2.536 0.01268 7.658 0.1211 #dY:1.6% + 2.582 0.01291 7.748 0.1213 #dY:1.6% + 2.631 0.01315 7.737 0.1206 #dY:1.6% + 2.698 0.01349 7.76 0.121 #dY:1.6% + 2.767 0.01384 7.828 0.122 #dY:1.6% + 2.84 0.0142 7.89 0.123 #dY:1.6% + 2.915 0.01457 7.924 0.1235 #dY:1.6% + 2.993 0.01497 7.947 0.1239 #dY:1.6% + 3.075 0.01537 7.967 0.1242 #dY:1.6% + 3.16 0.0158 7.994 0.1246 #dY:1.6% + 3.249 0.01624 8.007 0.1248 #dY:1.6% + 3.341 0.01671 8.054 0.1255 #dY:1.6% + 3.438 0.01719 8.09 0.1261 #dY:1.6% + 3.539 0.0177 8.056 0.1256 #dY:1.6% + 3.645 0.01822 8.066 0.1257 #dY:1.6% + 3.756 0.01878 8.071 0.1258 #dY:1.6% + 3.872 0.01936 8.074 0.1259 #dY:1.6% + 3.993 0.01997 8.109 0.1264 #dY:1.6% + 4.121 0.0206 8.067 0.1258 #dY:1.6% + 4.254 0.02127 8.06 0.1277 #dY:1.6% + 4.395 0.02197 8.078 0.1303 #dY:1.6% + 4.542 0.02271 8.055 0.1299 #dY:1.6% + 4.698 0.02349 7.901 0.1274 #dY:1.6% + 4.862 0.02431 7.804 0.1285 #dY:1.6% + 5.034 0.02517 7.769 0.1279 #dY:1.6% + 5.216 0.02608 7.613 0.1253 #dY:1.6% + 5.408 0.02704 7.565 0.1245 #dY:1.6% + 5.611 0.02806 7.412 0.122 #dY:1.6% + 5.826 0.02913 7.31 0.1203 #dY:1.6% + 6.054 0.03027 7.218 0.1216 #dY:1.7% + 6.296 0.03148 7.172 0.1209 #dY:1.7% + 6.552 0.03276 6.965 0.1174 #dY:1.7% + 6.825 0.03413 6.808 0.1147 #dY:1.7% + 7.115 0.03558 6.73 0.1134 #dY:1.7% + 7.425 0.03712 6.54 0.1131 #dY:1.7% + 7.755 0.03878 6.484 0.1121 #dY:1.7% + 8.108 0.04054 6.374 0.1133 #dY:1.8% + 8.486 0.04243 6.16 0.1095 #dY:1.8% + 8.891 0.04446 6.092 0.1115 #dY:1.8% + 9.481 0.04741 6.024 0.1083 #dY:1.8% + 10.04 0.1004 6.09 0.1154 #dY:1.9% + 10.56 0.1056 5.94 0.1121 #dY:1.9% + 11.12 0.1112 5.984 0.1095 #dY:1.8% + 11.74 0.1174 5.817 0.1045 #dY:1.8% + 12.4 0.124 5.825 0.1019 #dY:1.7% + 13.12 0.1312 5.868 0.1001 #dY:1.7% + 13.91 0.1391 5.934 0.09751 #dY:1.6% + 14.77 0.1477 5.95 0.09777 #dY:1.6% + 15.71 0.1571 6.053 0.09741 #dY:1.6% + 16.75 0.1675 6.067 0.09764 #dY:1.6% + 17.89 0.1789 6.19 0.1019 #dY:1.6% + 19.15 0.1915 6.215 0.1075 #dY:1.7% + 20.25 0.2025 6.333 0.1415 #dY:2.2% + 0.048 0.008022 13.19 0.2867 #dY:2.2% + 0.063 0.00805 12.79 0.2766 #dY:2.2% + 0.078 0.008022 12.49 0.2711 #dY:2.2% + 0.092 0.00804 12.22 0.2674 #dY:2.2% + 0.111 0.008062 11.78 0.2622 #dY:2.2% + 0.137 0.008122 11.41 0.2583 #dY:2.3% + 0.172 0.008246 10.94 0.2529 #dY:2.3% + 0.195 0.008476 10.71 0.1842 #dY:1.7% + 0.223 0.008509 10.41 0.2477 #dY:2.4% + 0.244 0.008443 10.04 0.1756 #dY:1.7% + 0.297 0.008443 9.7 0.1691 #dY:1.7% + 0.35 0.008476 9.38 0.1676 #dY:1.8% + 0.401 0.008443 8.899 0.1542 #dY:1.7% + 0.458 0.008443 8.426 0.1465 #dY:1.7% + 0.51 0.008693 8.141 0.1305 #dY:1.6% + 0.568 0.008654 7.791 0.1446 #dY:1.9% + 0.61 0.008693 7.776 0.1336 #dY:1.7% + 0.67 0.008616 7.567 0.1602 #dY:2.1% + 0.715 0.00858 7.365 0.1262 #dY:1.7% + 0.773 0.008616 7.23 0.1391 #dY:1.9% + 0.821 0.008654 7.057 0.1217 #dY:1.7% + 0.924 0.008693 6.868 0.1168 #dY:1.7% + 1.023 0.008544 6.756 0.1097 #dY:1.6% + 1.128 0.008509 6.738 0.1224 #dY:1.8% + 1.23 0.008476 6.558 0.1225 #dY:1.9% + 1.333 0.008443 6.644 0.2058 #dY:3.1% + 1.533 0.008246 6.78 0.1089 #dY:1.6% + 1.631 0.008179 6.788 0.1177 #dY:1.7% + 1.734 0.008223 6.966 0.1119 #dY:1.6% + 1.825 0.008158 7.015 0.1221 #dY:1.7% + 1.931 0.008223 7.073 0.1156 #dY:1.6% + 2.025 0.008139 7.21 0.1247 #dY:1.7% + 2.134 0.0105 7.316 0.1302 #dY:1.8% + 2.279 0.009381 7.372 0.121 #dY:1.6% + 2.429 0.01124 7.496 0.1319 #dY:1.8% + 2.574 0.009881 7.582 0.1234 #dY:1.6% + 2.791 0.009035 7.638 0.1312 #dY:1.7% + 2.994 0.009381 7.881 0.1326 #dY:1.7% + 3.188 0.009278 7.978 0.1339 #dY:1.7% + 3.387 0.009228 8.013 0.1348 #dY:1.7% + 3.585 0.00913 8.03 0.1364 #dY:1.7% + 3.781 0.008989 8.048 0.1381 #dY:1.7% + 3.974 0.0089 8.092 0.1401 #dY:1.7% + 4.185 0.008773 7.9 0.1397 #dY:1.8% + 4.389 0.008693 7.894 0.1424 #dY:1.8% diff --git a/233U.inp b/233U.inp index 09debc7..0c6e336 100644 --- a/233U.inp +++ b/233U.inp @@ -1,38 +1,41 @@ &Target Z=92.0 A=233.0 eferm=-6.3036 - BETA20=0.200 BETA40=0.129 BETA60=-1.5200E-02 + BETA20=0.181 BETA40=0.149 BETA60=-2.5000E-02 nstat=16 / -&State Et=0.0000E+00 Jval=2.5 Kval=2.5 NBAND=1 COEFF=0.200 / -&State Et=4.0350E+01 Jval=3.5 Kval=2.5 NBAND=1 COEFF=0.200 / -&State Et=9.2160E+01 Jval=4.5 Kval=2.5 NBAND=1 COEFF=0.200 / -&State Et=1.5523E+02 Jval=5.5 Kval=2.5 NBAND=1 COEFF=0.200 / -&State Et=2.2947E+02 Jval=6.5 Kval=2.5 NBAND=1 COEFF=0.200 / -&State Et=3.1460E+02 Jval=7.5 Kval=2.5 NBAND=1 COEFF=0.200 / -&State Et=4.1170E+02 Jval=8.5 Kval=2.5 NBAND=1 COEFF=0.200 / -&State Et=5.1755E+02 Jval=9.5 Kval=2.5 NBAND=1 COEFF=0.200 / -&State Et=6.3527E+02 Jval=10.5 Kval=2.5 NBAND=1 COEFF=0.20 / +&State Et=0.0000E+00 Jval=2.5 Kval=2.5 NBAND=1 COEFF=0.181 / +&State Et=4.0350E+01 Jval=3.5 Kval=2.5 NBAND=1 COEFF=0.181 / +&State Et=9.2160E+01 Jval=4.5 Kval=2.5 NBAND=1 COEFF=0.181 / +&State Et=1.5523E+02 Jval=5.5 Kval=2.5 NBAND=1 COEFF=0.181 / +&State Et=2.2947E+02 Jval=6.5 Kval=2.5 NBAND=1 COEFF=0.181 / +&State Et=3.1460E+02 Jval=7.5 Kval=2.5 NBAND=1 COEFF=0.181 / +&State Et=4.1170E+02 Jval=8.5 Kval=2.5 NBAND=1 COEFF=0.181 / +&State Et=5.1755E+02 Jval=9.5 Kval=2.5 NBAND=1 COEFF=0.181 / +&State Et=6.3527E+02 Jval=10.5 Kval=2.5 NBAND=1 COEFF=0.181 / -&State Et=3.9850E+02 Jval=0.5 Kval=0.5 NBAND=2 COEFF=0.034 / -&State Et=4.1580E+02 Jval=1.5 Kval=0.5 NBAND=2 COEFF=0.034 / -&State Et=4.5610E+02 Jval=2.5 Kval=0.5 NBAND=2 COEFF=0.034 / +&State Et=3.9850E+02 Jval=0.5 Kval=0.5 NBAND=2 COEFF=0.001041 / +&State Et=4.1580E+02 Jval=1.5 Kval=0.5 NBAND=2 COEFF=0.001041 / +&State Et=4.5610E+02 Jval=2.5 Kval=0.5 NBAND=2 COEFF=0.001041 / -&State Et=2.9880E+02 Jval=2.5 Kval=2.5 NBAND=-3 COEFF=0.065 / -&State Et=3.2080E+02 Jval=3.5 Kval=2.5 NBAND=-3 COEFF=0.065 / -&State Et=3.5380E+02 Jval=4.5 Kval=2.5 NBAND=-3 COEFF=0.065 / -&State Et=3.9760E+02 Jval=5.5 Kval=2.5 NBAND=-3 COEFF=0.065 / +&State Et=2.9880E+02 Jval=2.5 Kval=2.5 NBAND=-3 COEFF=0.043078 / +&State Et=3.2080E+02 Jval=3.5 Kval=2.5 NBAND=-3 COEFF=0.043078 / +&State Et=3.5380E+02 Jval=4.5 Kval=2.5 NBAND=-3 COEFF=0.043078 / +&State Et=3.9760E+02 Jval=5.5 Kval=2.5 NBAND=-3 COEFF=0.043078 / -&Fresco hcm=0.10 rmatch= 20.00 Jmax=20 - Ngrid=41 pythonFlag=0 Nenergy=2 absend= 0.001 - relativ= 2 / +&Fresco hcm=0.10 rmatch= 20.00 Jmax=30 + Ngrid=41 pythonFlag=0 Nenergy=37 absend= 0.001 + relativ= 3 / &Energies - elab=0.001,2.00/ + elab=0.001,0.002,0.004,0.006,0.008,0.01,0.02,0.04, + 0.06,0.08,0.100,0.150,0.200,0.30,0.400,0.600,0.800,1.0000,1.40,2.0000, + 2.20,2.70,3.20,3.70,4.2,4.8,5.00, + 7.00,10.0000,12.00,15.00,17.00,20.00,23.00,26.00,28.00,30.00 / -&Potential v0a=50.41 v0b=0.0292 lambdhf= 0.00977 - cviso=17.5 vspo=5.64 lambdso=0.005 ccoul=1.36 - av=11.91 bv=81.86 w0=17.70 bs=10.85 wspo=-3.1 - bso=160.0 ea=55.0 alphav=0.355 cs=0.01331 cwiso=29.0 - adv=238.0 rhfa=1.2490 rhfb=-0.00171 ahfa=0.638 ahfb=0.00219 - rv=1.2657 ava=0.6960 avb=-0.00021 rsa=1.1701 rsb=0.0041 - as=0.617 rso=1.1214 aso=0.59 rc=1.1289 ac=0.547 / +&Potential v0a=50.59 v0b=0.0608 lambdhf= 0.00982 + cviso=16.79 vspo=6.83 lambdso=0.005 ccoul=1.43 + av=11.66 bv=81.61 w0a=17.30 w0b=-0.031 bs=10.55 wspo=-3.1 + bso=160.0 ea=55.0 alphav=0.355 cs=0.01343 cwiso=28.96 + adv=238.0 rhfa=1.2426 rhfb=-0.00246 ahfa=0.638 ahfb=0.00306 + rv=1.2704 avv=0.6970 rsa=1.1860 rsb=0.0060 + as=0.616 rso=1.0156 aso=0.549 rc=1.2262 ac=0.674 / / diff --git a/235U.exp b/235U.exp new file mode 100644 index 0000000..7dfab58 --- /dev/null +++ b/235U.exp @@ -0,0 +1,172 @@ + 0.048 0.008022 13.22 0.1692 #dY:1.3% + 0.063 0.00805 12.77 0.2134 #dY:1.7% + 0.078 0.008022 12.31 0.1581 #dY:1.3% + 0.092 0.00804 12.08 0.1564 #dY:1.3% + 0.111 0.008062 11.74 0.1545 #dY:1.3% + 0.137 0.008122 11.38 0.1534 #dY:1.3% + 0.172 0.008246 10.89 0.1512 #dY:1.4% + 0.195 0.008476 10.54 0.1022 #dY:1.0% + 0.223 0.008509 10.32 0.1497 #dY:1.5% + 0.244 0.008443 10.15 0.1172 #dY:1.2% + 0.297 0.008443 9.667 0.09413 #dY:1.0% + 0.35 0.008476 9.256 0.09831 #dY:1.1% + 0.401 0.008443 8.795 0.101 #dY:1.1% + 0.458 0.008443 8.532 0.0682 #dY:0.8% + 0.51 0.008693 8.261 0.06014 #dY:0.7% + 0.568 0.008654 8.071 0.1286 #dY:1.6% + 0.61 0.008693 7.808 0.05733 #dY:0.7% + 0.67 0.008616 7.64 0.1398 #dY:1.8% + 0.715 0.00858 7.488 0.0758 #dY:1.0% + 0.773 0.008616 7.324 0.144 #dY:2.0% + 0.821 0.008654 7.211 0.06465 #dY:0.9% + 0.924 0.008693 7.179 0.07635 #dY:1.1% + 1.023 0.008544 6.927 0.05566 #dY:0.8% + 1.128 0.008509 6.82 0.06234 #dY:0.9% + 1.23 0.008476 6.809 0.07015 #dY:1.0% + 1.333 0.008443 6.781 0.09517 #dY:1.4% + 1.533 0.008246 6.85 0.054 #dY:0.8% + 1.631 0.008179 6.892 0.07044 #dY:1.0% + 1.734 0.008223 7.068 0.05695 #dY:0.8% + 1.825 0.008158 7.074 0.07433 #dY:1.1% + 1.931 0.008223 7.16 0.05736 #dY:0.8% + 2.025 0.008139 7.156 0.0738 #dY:1.0% + 2.134 0.0105 7.404 0.08124 #dY:1.1% + 2.279 0.009381 7.46 0.06564 #dY:0.9% + 2.429 0.01124 7.602 0.08024 #dY:1.1% + 2.574 0.009881 7.638 0.06565 #dY:0.9% + 2.791 0.009035 7.768 0.07758 #dY:1.0% + 2.994 0.009381 7.801 0.07532 #dY:1.0% + 3.188 0.009278 7.818 0.07538 #dY:1.0% + 3.387 0.009228 7.941 0.07584 #dY:1.0% + 3.585 0.00913 8.014 0.07847 #dY:1.0% + 3.781 0.008989 7.944 0.08061 #dY:1.0% + 3.974 0.0089 8.145 0.08375 #dY:1.0% + 4.185 0.008773 7.98 0.08648 #dY:1.1% + 4.389 0.008693 7.883 0.09038 #dY:1.1% + 4.598 0.008297 7.972 0.09581 #dY:1.2% + 4.807 0.008297 7.773 0.1113 #dY:1.4% + 1.817 0.009085 7.046 0.05895 #dY:0.8% + 1.854 0.009272 7.07 0.05915 #dY:0.8% + 1.893 0.009465 7.12 0.05957 #dY:0.8% + 1.933 0.009665 7.139 0.05973 #dY:0.8% + 1.974 0.009871 7.207 0.0603 #dY:0.8% + 2.017 0.01008 7.137 0.05971 #dY:0.8% + 2.061 0.0103 7.31 0.05709 #dY:0.8% + 2.107 0.01053 7.283 0.05688 #dY:0.8% + 2.154 0.01077 7.377 0.05762 #dY:0.8% + 2.203 0.01101 7.325 0.05721 #dY:0.8% + 2.253 0.01127 7.468 0.05833 #dY:0.8% + 2.306 0.01153 7.458 0.05825 #dY:0.8% + 2.36 0.0118 7.533 0.05883 #dY:0.8% + 2.417 0.01208 7.51 0.05865 #dY:0.8% + 2.475 0.01238 7.608 0.05942 #dY:0.8% + 2.536 0.01268 7.639 0.05613 #dY:0.7% + 2.582 0.01291 7.757 0.05805 #dY:0.7% + 2.631 0.01315 7.688 0.0565 #dY:0.7% + 2.698 0.01349 7.759 0.05702 #dY:0.7% + 2.767 0.01384 7.781 0.05718 #dY:0.7% + 2.84 0.0142 7.841 0.05762 #dY:0.7% + 2.915 0.01457 7.9 0.05805 #dY:0.7% + 2.993 0.01497 7.89 0.05798 #dY:0.7% + 3.075 0.01537 7.89 0.05798 #dY:0.7% + 3.16 0.0158 7.92 0.0582 #dY:0.7% + 3.249 0.01624 7.978 0.05863 #dY:0.7% + 3.341 0.01671 7.982 0.05866 #dY:0.7% + 3.438 0.01719 8 0.05879 #dY:0.7% + 3.539 0.0177 8.001 0.0588 #dY:0.7% + 3.645 0.01822 8.001 0.05601 #dY:0.7% + 3.756 0.01878 8.016 0.05611 #dY:0.7% + 3.872 0.01936 8.041 0.05629 #dY:0.7% + 3.993 0.01997 8.019 0.05613 #dY:0.7% + 4.121 0.0206 8.027 0.05619 #dY:0.7% + 4.254 0.02127 8.06 0.06032 #dY:0.7% + 4.395 0.02197 7.98 0.06483 #dY:0.8% + 4.542 0.02271 7.917 0.06671 #dY:0.8% + 4.698 0.02349 7.961 0.06708 #dY:0.8% + 4.862 0.02431 7.819 0.06588 #dY:0.8% + 5.034 0.02517 7.722 0.06993 #dY:0.9% + 5.216 0.02608 7.61 0.06891 #dY:0.9% + 5.408 0.02704 7.502 0.06793 #dY:0.9% + 5.611 0.02806 7.311 0.0662 #dY:0.9% + 5.826 0.02913 7.253 0.06568 #dY:0.9% + 6.054 0.03027 7.154 0.06973 #dY:1.0% + 6.296 0.03148 7.028 0.0685 #dY:1.0% + 6.552 0.03276 6.928 0.06753 #dY:1.0% + 6.825 0.03413 6.81 0.06638 #dY:1.0% + 7.115 0.03558 6.688 0.06519 #dY:1.0% + 7.425 0.03712 6.48 0.06796 #dY:1.0% + 7.755 0.03878 6.307 0.06704 #dY:1.1% + 8.108 0.04054 6.37 0.07263 #dY:1.1% + 8.486 0.04243 6.266 0.07144 #dY:1.1% + 8.891 0.04446 6.219 0.07091 #dY:1.1% + 9.481 0.04741 6.068 0.06592 #dY:1.1% + 10.04 0.1004 5.95 0.07871 #dY:1.3% + 10.56 0.1056 5.859 0.07271 #dY:1.2% + 11.12 0.1112 5.762 0.07311 #dY:1.3% + 11.74 0.1174 5.783 0.06335 #dY:1.1% + 12.4 0.124 5.766 0.05708 #dY:1.0% + 13.12 0.1312 5.888 0.05739 #dY:1.0% + 13.91 0.1391 5.973 0.05409 #dY:0.9% + 14.77 0.1477 5.937 0.05003 #dY:0.8% + 15.71 0.1571 6.037 0.05087 #dY:0.8% + 0.00015 5e-5 45.95 1 #dY:2.2% + 0.00025 5e-5 44.71 1 #dY:2.2% + 0.00035 5e-5 31.11 0.18 #dY:0.6% + 0.00045 5e-5 29.73 0.19 #dY:0.6% + 0.00055 5e-5 31.24 0.21 #dY:0.7% + 0.00065 5e-5 27.15 0.21 #dY:0.8% + 0.00075 5e-5 28.57 0.22 #dY:0.8% + 0.00085 5e-5 23.69 0.24 #dY:1.0% + 0.00095 5e-5 24.64 0.25 #dY:1.0% + 0.00105 5e-5 25.93 0.3 #dY:1.2% + 0.00115 5e-5 25.39 0.3 #dY:1.2% + 0.00125 5e-5 22.86 0.3 #dY:1.3% + 0.00135 5e-5 24.08 0.3 #dY:1.2% + 0.00145 5e-5 21.78 0.3 #dY:1.4% + 0.00155 5e-5 21.2 0.3 #dY:1.4% + 0.00165 5e-5 21.87 0.3 #dY:1.4% + 0.00175 5e-5 21.27 0.3 #dY:1.4% + 0.00185 5e-5 20.37 0.3 #dY:1.5% + 0.00195 5e-5 23.09 0.3 #dY:1.3% + 0.0025 0.0005 19.75 0.1 #dY:0.5% + 0.0035 0.0005 18.8 0.14 #dY:0.7% + 0.0045 0.0005 18.1 0.16 #dY:0.9% + 0.0065 0.0005 17 0.18 #dY:1.1% + 0.0075 0.0005 16.65 0.19 #dY:1.1% + 0.0085 0.0005 17.05 0.2 #dY:1.2% + 0.0095 0.0005 16.45 0.2 #dY:1.2% + 0.015 0.005 15.45 0.1 #dY:0.6% + 0.025 0.005 14.5 0.1 #dY:0.7% + 0.035 0.005 14 0.1 #dY:0.7% + 0.045 0.005 13.55 0.1 #dY:0.7% + 0.055 0.005 13.3 0.1 #dY:0.8% + 0.065 0.005 12.9 0.1 #dY:0.8% + 0.075 0.005 12.6 0.1 #dY:0.8% + 0.085 0.005 12.65 0.1 #dY:0.8% + 0.095 0.005 12.3 0.1 #dY:0.8% + 0.105 0.005 12.1 0.1 #dY:0.8% + 0.115 0.005 11.9 0.1 #dY:0.8% + 0.125 0.005 11.8 0.1 #dY:0.8% + 0.135 0.005 11.85 0.1 #dY:0.8% + 0.145 0.005 11.6 0.1 #dY:0.9% + 0.155 0.005 11.7 0.1 #dY:0.9% + 0.165 0.005 11.4 0.1 #dY:0.9% + 0.175 0.005 11.05 0.1 #dY:0.9% + 0.185 0.005 11 0.1 #dY:0.9% + 0.195 0.005 10.6 0.1 #dY:0.9% + 0.225 0.025 10.5 0.1 #dY:1.0% + 0.275 0.025 10 0.1 #dY:1.0% + 0.325 0.025 9.6 0.1 #dY:1.0% + 0.375 0.025 9.2 0.1 #dY:1.1% + 0.425 0.025 8.95 0.1 #dY:1.1% + 0.475 0.025 8.8 0.1 #dY:1.1% + 0.525 0.025 8.25 0.1 #dY:1.2% + 0.575 0.025 8 0.1 #dY:1.3% + 0.65 0.05 7.8 0.1 #dY:1.3% + 0.75 0.05 7.5 0.1 #dY:1.3% + 0.85 0.05 7.2 0.1 #dY:1.4% + 0.95 0.05 7 0.1 #dY:1.4% + 16.75 0.1675 6.085 0.05127 #dY:0.8% + 17.89 0.1789 6.201 0.05225 #dY:0.8% + 19.15 0.1915 6.179 0.06117 #dY:1.0% + 20.25 0.2025 6.298 0.1065 #dY:1.7% diff --git a/235U.inp b/235U.inp index 6927945..d6dd8df 100644 --- a/235U.inp +++ b/235U.inp @@ -1,5 +1,5 @@ &Target Z=92.0 A=235.0 eferm=-5.92151 - BETA20=0.220 BETA40=0.109 BETA60=-5.700E-02 + BETA20=0.220 BETA40=0.102 BETA60=-6.200E-02 nstat=21 / &State Et=0.0000E+00 Jval=3.5 Kval=3.5 NBAND=-1 COEFF=0.220 / @@ -11,37 +11,37 @@ &State Et=4.3940E+02 Jval=9.5 Kval=3.5 NBAND=-1 COEFF=0.220 / &State Et=5.5120E+02 Jval=10.5 Kval=3.5 NBAND=-1 COEFF=0.22 / -&State Et=4.4560E+02 Jval=3.5 Kval=3.5 NBAND=2 COEFF=0.0455 / -&State Et=5.1060E+02 Jval=4.5 Kval=3.5 NBAND=2 COEFF=0.0455 / -&State Et=5.7880E+02 Jval=5.5 Kval=3.5 NBAND=2 COEFF=0.0455 / -&State Et=6.9020E+02 Jval=6.5 Kval=3.5 NBAND=2 COEFF=0.0455 / +&State Et=4.4560E+02 Jval=3.5 Kval=3.5 NBAND=2 COEFF=0.000209 / +&State Et=5.1060E+02 Jval=4.5 Kval=3.5 NBAND=2 COEFF=0.000209 / +&State Et=5.7880E+02 Jval=5.5 Kval=3.5 NBAND=2 COEFF=0.000209 / +&State Et=6.9020E+02 Jval=6.5 Kval=3.5 NBAND=2 COEFF=0.000209 / -&State Et=3.9320E+02 Jval=1.5 Kval=1.5 NBAND=3 COEFF=0.040 / -&State Et=4.2670E+02 Jval=2.5 Kval=1.5 NBAND=3 COEFF=0.040 / -&State Et=4.7380E+02 Jval=3.5 Kval=1.5 NBAND=3 COEFF=0.040 / -&State Et=5.3320E+02 Jval=4.5 Kval=1.5 NBAND=3 COEFF=0.040 / -&State Et=6.0820E+02 Jval=5.5 Kval=1.5 NBAND=3 COEFF=0.040 / +&State Et=3.9320E+02 Jval=1.5 Kval=1.5 NBAND=3 COEFF=0.00682 / +&State Et=4.2670E+02 Jval=2.5 Kval=1.5 NBAND=3 COEFF=0.00682 / +&State Et=4.7380E+02 Jval=3.5 Kval=1.5 NBAND=3 COEFF=0.00682 / +&State Et=5.3320E+02 Jval=4.5 Kval=1.5 NBAND=3 COEFF=0.00682 / +&State Et=6.0820E+02 Jval=5.5 Kval=1.5 NBAND=3 COEFF=0.00682 / -&State Et=6.3780E+02 Jval=1.5 Kval=1.5 NBAND=-4 COEFF=0.030 / -&State Et=6.6645E+02 Jval=2.5 Kval=1.5 NBAND=-4 COEFF=0.030 / -&State Et=7.0120E+02 Jval=3.5 Kval=1.5 NBAND=-4 COEFF=0.030 / -&State Et=7.5500E+02 Jval=4.5 Kval=1.5 NBAND=-4 COEFF=0.030 / +&State Et=6.3780E+02 Jval=1.5 Kval=1.5 NBAND=-4 COEFF=0.01496 / +&State Et=6.6645E+02 Jval=2.5 Kval=1.5 NBAND=-4 COEFF=0.01496 / +&State Et=7.0120E+02 Jval=3.5 Kval=1.5 NBAND=-4 COEFF=0.01496 / +&State Et=7.5500E+02 Jval=4.5 Kval=1.5 NBAND=-4 COEFF=0.01496 / -&Fresco hcm=0.10 rmatch= 20.00 Jmax=20 - Ngrid=41 pythonFlag=0 Nenergy=35 absend= 0.001 - relativ= 2 / +&Fresco hcm=0.10 rmatch= 20.00 Jmax=30 + Ngrid=41 pythonFlag=0 Nenergy=37 absend= 0.001 + relativ= 3 / &Energies elab=0.001,0.002,0.004,0.006,0.008,0.01,0.02,0.04, - 0.06,0.08,0.100,0.200,0.400,0.600,0.800,1.0000,1.40,2.0000, + 0.06,0.08,0.100,0.150,0.200,0.30,0.400,0.600,0.800,1.0000,1.40,2.0000, 2.20,2.70,3.20,3.70,4.2,4.8,5.00, 7.00,10.0000,12.00,15.00,17.00,20.00,23.00,26.00,28.00,30.00 / -&Potential v0a=50.41 v0b=0.0292 lambdhf= 0.00977 - cviso=17.50 vspo=5.64 lambdso=0.005 ccoul=1.36 - av=11.91 bv=81.86 w0=17.70 bs=10.85 wspo=-3.1 - bso=160.0 ea=55.0 alphav=0.355 cs=0.01331 cwiso=29.0 - adv=238.0 rhfa=1.2490 rhfb=-0.00171 ahfa=0.638 ahfb=0.00219 - rv=1.2657 ava=0.6960 avb=-0.00021 rsa=1.1701 rsb=0.0041 - as=0.617 rso=1.1214 aso=0.59 rc=1.1289 ac=0.547 / +&Potential v0a=50.59 v0b=0.0608 lambdhf= 0.00982 + cviso=16.79 vspo=6.83 lambdso=0.005 ccoul=1.43 + av=11.66 bv=81.61 w0a=17.30 w0b=-0.031 bs=10.55 wspo=-3.1 + bso=160.0 ea=55.0 alphav=0.355 cs=0.01343 cwiso=28.96 + adv=238.0 rhfa=1.2426 rhfb=-0.00246 ahfa=0.638 ahfb=0.00306 + rv=1.2704 avv=0.6970 rsa=1.1860 rsb=0.0060 + as=0.616 rso=1.0156 aso=0.549 rc=1.2262 ac=0.674 / / diff --git a/238U.exp b/238U.exp new file mode 100644 index 0000000..e9bbf36 --- /dev/null +++ b/238U.exp @@ -0,0 +1,253 @@ + 0.048 0.008022 13.3 0.25 #dY:1.9% + 0.063 0.00805 12.85 0.2282 #dY:1.8% + 0.078 0.008022 12.42 0.216 #dY:1.7% + 0.092 0.00804 12.39 0.2098 #dY:1.7% + 0.111 0.008062 11.92 0.2019 #dY:1.7% + 0.137 0.008122 11.64 0.1954 #dY:1.7% + 0.172 0.008246 11.21 0.1901 #dY:1.7% + 0.195 0.008476 10.81 0.09693 #dY:0.9% + 0.223 0.008509 10.52 0.185 #dY:1.8% + 0.244 0.008443 10.31 0.09768 #dY:0.9% + 0.297 0.008443 9.842 0.08628 #dY:0.9% + 0.35 0.008476 9.466 0.08682 #dY:0.9% + 0.401 0.008443 9.09 0.08107 #dY:0.9% + 0.458 0.008443 8.68 0.0635 #dY:0.7% + 0.51 0.008693 8.486 0.0542 #dY:0.6% + 0.568 0.008654 8.202 0.08541 #dY:1.0% + 0.61 0.008693 8.087 0.05306 #dY:0.7% + 0.67 0.008616 7.868 0.1032 #dY:1.3% + 0.715 0.00858 7.745 0.06112 #dY:0.8% + 0.773 0.008616 7.508 0.08805 #dY:1.2% + 0.821 0.008654 7.382 0.05472 #dY:0.7% + 0.924 0.008693 7.249 0.05706 #dY:0.8% + 1.023 0.008544 7.17 0.04588 #dY:0.6% + 1.128 0.008509 7.118 0.05448 #dY:0.8% + 1.23 0.008476 6.985 0.05889 #dY:0.8% + 1.333 0.008443 6.975 0.08208 #dY:1.2% + 1.533 0.008246 7.092 0.04238 #dY:0.6% + 1.631 0.008179 7.097 0.05606 #dY:0.8% + 1.734 0.008223 7.146 0.1105 #dY:1.5% + 1.825 0.008158 7.234 0.05788 #dY:0.8% + 1.931 0.008223 7.299 0.04606 #dY:0.6% + 2.025 0.008139 7.338 0.0588 #dY:0.8% + 2.134 0.0105 7.349 0.05296 #dY:0.7% + 2.279 0.009381 7.408 0.05382 #dY:0.7% + 2.429 0.01124 7.62 0.05427 #dY:0.7% + 2.574 0.009881 7.668 0.06179 #dY:0.8% + 2.791 0.009035 7.73 0.05929 #dY:0.8% + 2.994 0.009381 7.815 0.05857 #dY:0.7% + 3.188 0.009278 7.926 0.06043 #dY:0.8% + 3.387 0.009228 7.951 0.06043 #dY:0.8% + 3.585 0.00913 7.973 0.06295 #dY:0.8% + 3.781 0.008989 8.014 0.06644 #dY:0.8% + 3.974 0.0089 8.039 0.06819 #dY:0.8% + 4.185 0.008773 8.065 0.07441 #dY:0.9% + 4.389 0.008693 8.039 0.08342 #dY:1.0% + 4.598 0.008297 7.97 0.0973 #dY:1.2% + 4.807 0.008297 7.735 0.08379 #dY:1.1% + 5.293 0 7.495 0.0227 #dY:0.3% + 5.346 0 7.489 0.02349 #dY:0.3% + 5.399 0 7.472 0.02435 #dY:0.3% + 5.454 0 7.449 0.02241 #dY:0.3% + 5.509 0 7.399 0.02185 #dY:0.3% + 5.564 00 7.384 0.02172 #dY:0.3% + 5.62 0 7.385 0.02169 #dY:0.3% + 5.676 00 7.319 0.02115 #dY:0.3% + 5.733 0 7.294 0.02074 #dY:0.3% + 5.791 0 7.278 0.0205 #dY:0.3% + 5.849 0 7.241 0.02028 #dY:0.3% + 5.908 0 7.215 0.02023 #dY:0.3% + 5.967 0 7.168 0.02001 #dY:0.3% + 6.027 0 7.124 0.01971 #dY:0.3% + 6.088 0 7.097 0.01964 #dY:0.3% + 6.149 0 7.09 0.01977 #dY:0.3% + 6.211 00 7.065 0.02057 #dY:0.3% + 6.273 0 7.047 0.02259 #dY:0.3% + 6.336 00 7.011 0.02251 #dY:0.3% + 6.4 0 6.932 0.01925 #dY:0.3% + 6.464 0 6.909 0.0182 #dY:0.3% + 6.529 0 6.886 0.01763 #dY:0.3% + 6.595 0 6.852 0.01687 #dY:0.2% + 6.661 0 6.825 0.01674 #dY:0.2% + 6.728 0 6.778 0.01655 #dY:0.2% + 6.796 00 6.763 0.01643 #dY:0.2% + 6.864 0 6.738 0.01619 #dY:0.2% + 6.933 0 6.717 0.01595 #dY:0.2% + 7.003 0 6.695 0.01587 #dY:0.2% + 7.073 0 6.659 0.01582 #dY:0.2% + 7.144 0 6.614 0.01588 #dY:0.2% + 7.216 0 6.594 0.01648 #dY:0.2% + 7.289 0 6.569 0.01782 #dY:0.3% + 7.362 0 6.556 0.01973 #dY:0.3% + 7.436 0 6.503 0.02005 #dY:0.3% + 7.51 0 6.46 0.01958 #dY:0.3% + 7.586 0 6.45 0.01931 #dY:0.3% + 7.662 0 6.442 0.01987 #dY:0.3% + 7.739 0 6.422 0.02128 #dY:0.3% + 7.817 0 6.376 0.02123 #dY:0.3% + 7.896 0 6.352 0.0204 #dY:0.3% + 7.975 0 6.332 0.01977 #dY:0.3% + 8.055 0 6.313 0.01946 #dY:0.3% + 8.136 0 6.291 0.0189 #dY:0.3% + 8.218 0 6.25 0.0175 #dY:0.3% + 8.3 0 6.219 0.01654 #dY:0.3% + 8.384 0 6.204 0.01608 #dY:0.3% + 8.468 0 6.168 0.01574 #dY:0.3% + 8.553 0 6.168 0.01568 #dY:0.3% + 8.639 00 6.132 0.01555 #dY:0.3% + 8.726 0 6.11 0.01556 #dY:0.3% + 8.814 0 6.102 0.0156 #dY:0.3% + 8.902 0 6.089 0.01559 #dY:0.3% + 8.992 0 6.052 0.01569 #dY:0.3% + 9.082 0 6.054 0.01607 #dY:0.3% + 9.173 0 6.04 0.01633 #dY:0.3% + 9.266 0 5.994 0.01634 #dY:0.3% + 9.359 0 5.993 0.01632 #dY:0.3% + 9.453 0 5.966 0.01613 #dY:0.3% + 9.548 0 5.95 0.016 #dY:0.3% + 9.644 0 5.967 0.01611 #dY:0.3% + 9.741 0 5.942 0.01609 #dY:0.3% + 9.838 0 5.911 0.016 #dY:0.3% + 9.937 0 5.904 0.0159 #dY:0.3% + 10.04 0 5.895 0.01588 #dY:0.3% + 10.14 0 5.867 0.01581 #dY:0.3% + 10.24 0 5.858 0.01579 #dY:0.3% + 10.34 0 5.85 0.01581 #dY:0.3% + 10.45 00 5.835 0.01596 #dY:0.3% + 10.55 0 5.832 0.0162 #dY:0.3% + 10.66 0 5.829 0.01649 #dY:0.3% + 10.77 0 5.861 0.01688 #dY:0.3% + 10.87 0 5.84 0.01708 #dY:0.3% + 10.98 0 5.813 0.01718 #dY:0.3% + 11.09 0 5.797 0.0172 #dY:0.3% + 11.2 0 5.773 0.01715 #dY:0.3% + 11.32 0 5.781 0.01719 #dY:0.3% + 11.43 0 5.78 0.0172 #dY:0.3% + 11.55 0 5.769 0.01716 #dY:0.3% + 11.66 0 5.775 0.01722 #dY:0.3% + 11.78 0 5.795 0.01738 #dY:0.3% + 11.9 0 5.797 0.01765 #dY:0.3% + 12.02 0 5.803 0.01801 #dY:0.3% + 12.14 0 5.789 0.01812 #dY:0.3% + 12.26 0 5.778 0.01786 #dY:0.3% + 12.38 0 5.77 0.01768 #dY:0.3% + 12.51 0 5.768 0.01772 #dY:0.3% + 12.63 0 5.79 0.01789 #dY:0.3% + 12.76 0 5.81 0.01811 #dY:0.3% + 12.89 0 5.83 0.01843 #dY:0.3% + 13.02 0 5.797 0.01846 #dY:0.3% + 13.15 0 5.785 0.01842 #dY:0.3% + 13.28 0 5.803 0.01849 #dY:0.3% + 13.41 0 5.81 0.01859 #dY:0.3% + 13.55 0 5.821 0.01868 #dY:0.3% + 13.69 0 5.828 0.01865 #dY:0.3% + 13.82 0 5.845 0.0186 #dY:0.3% + 13.96 0 5.843 0.01853 #dY:0.3% + 14.1 0 5.848 0.01862 #dY:0.3% + 14.24 0 5.851 0.01869 #dY:0.3% + 14.39 0 5.878 0.01888 #dY:0.3% + 14.53 0 5.916 0.01915 #dY:0.3% + 14.68 0 5.911 0.01941 #dY:0.3% + 14.82 0 5.915 0.01963 #dY:0.3% + 14.97 0 5.942 0.01984 #dY:0.3% + 15.12 0 5.956 0.01757 #dY:0.3% + 15.28 0 5.959 0.0122 #dY:0.2% + 15.43 0 5.974 0.01231 #dY:0.2% + 15.58 0 5.994 0.01247 #dY:0.2% + 15.74 0 5.991 0.01258 #dY:0.2% + 15.9 0 5.997 0.01265 #dY:0.2% + 16.06 0 6.018 0.01272 #dY:0.2% + 16.22 0 6.029 0.01266 #dY:0.2% + 16.38 0 6.035 0.01261 #dY:0.2% + 16.55 0 6.035 0.01257 #dY:0.2% + 16.72 0 6.048 0.01253 #dY:0.2% + 16.88 0 6.066 0.01254 #dY:0.2% + 17.05 0 6.08 0.01259 #dY:0.2% + 17.22 0 6.095 0.01263 #dY:0.2% + 17.4 0 6.103 0.01262 #dY:0.2% + 17.57 0 6.124 0.01266 #dY:0.2% + 17.75 00 6.128 0.01268 #dY:0.2% + 17.93 0 6.144 0.01275 #dY:0.2% + 18.11 0 6.159 0.01283 #dY:0.2% + 18.29 0 6.165 0.01286 #dY:0.2% + 18.47 0 6.149 0.01282 #dY:0.2% + 18.66 0 6.164 0.01287 #dY:0.2% + 18.85 0 6.188 0.01295 #dY:0.2% + 19.04 0 6.188 0.01297 #dY:0.2% + 19.23 0 6.189 0.013 #dY:0.2% + 19.42 0 6.183 0.013 #dY:0.2% + 19.62 0 6.201 0.01302 #dY:0.2% + 19.81 0 6.21 0.01297 #dY:0.2% + 20.01 00 6.206 0.01287 #dY:0.2% + 20.21 0 6.217 0.0128 #dY:0.2% + 20.42 00 6.228 0.01273 #dY:0.2% + 20.62 0 6.228 0.01262 #dY:0.2% + 20.83 0 6.221 0.01256 #dY:0.2% + 21.04 0 6.199 0.01247 #dY:0.2% + 21.25 00 6.201 0.0124 #dY:0.2% + 21.46 0 6.199 0.01234 #dY:0.2% + 21.68 0 6.204 0.01229 #dY:0.2% + 21.9 0 6.221 0.01226 #dY:0.2% + 22.12 0 6.221 0.01215 #dY:0.2% + 22.34 00 6.201 0.01203 #dY:0.2% + 22.56 0 6.209 0.012 #dY:0.2% + 22.79 0 6.207 0.01195 #dY:0.2% + 23.02 0 6.192 0.01189 #dY:0.2% + 23.25 0 6.172 0.01182 #dY:0.2% + 23.48 0 6.152 0.01174 #dY:0.2% + 23.72 0 6.149 0.01168 #dY:0.2% + 23.96 0 6.147 0.01163 #dY:0.2% + 24.2 0 6.129 0.01156 #dY:0.2% + 24.44 0 6.122 0.01151 #dY:0.2% + 24.69 0 6.125 0.01145 #dY:0.2% + 24.94 0 6.111 0.01137 #dY:0.2% + 25.19 0 6.085 0.01126 #dY:0.2% + 25.44 0 6.063 0.01115 #dY:0.2% + 25.7 0 6.055 0.01107 #dY:0.2% + 25.95 0 6.039 0.01097 #dY:0.2% + 26.21 00 6.028 0.01089 #dY:0.2% + 26.48 0 6.012 0.01079 #dY:0.2% + 26.74 0 5.984 0.01067 #dY:0.2% + 27.01 0 5.96 0.01056 #dY:0.2% + 27.28 0 5.94 0.01045 #dY:0.2% + 27.56 0 5.917 0.01036 #dY:0.2% + 27.84 0 5.899 0.01027 #dY:0.2% + 28.12 0 5.874 0.01015 #dY:0.2% + 28.4 0 5.866 0.01007 #dY:0.2% + 28.68 0 5.842 0.009961 #dY:0.2% + 28.97 0 5.811 0.009849 #dY:0.2% + 29.26 0 5.784 0.009734 #dY:0.2% + 29.56 0 5.766 0.009649 #dY:0.2% + 29.85 0 5.736 0.009554 #dY:0.2% + 5.5E-03 0.5E-03 19.13 0.4 #dY:2.1% + 6.5E-03 0.5E-03 18.44 0.71 #dY:3.9% + 7.5E-03 0.5E-03 16.16 0.55 #dY:3.4% + 8.5E-03 0.5E-03 15.7 0.5 #dY:3.2% + 9.5E-03 0.5E-03 14.7 0.26 #dY:1.8% + 11E-03 1E-03 13.94 0.3 #dY:2.2% + 13E-03 1E-03 14.68 0.41 #dY:2.8% + 15E-03 1E-03 13.84 0.27 #dY:2.0% + 17E-03 1E-03 14.44 0.49 #dY:3.4% + 19E-03 1E-03 13.84 0.17 #dY:1.2% + 21E-03 1E-03 14.18 0.31 #dY:2.2% + 23E-03 1E-03 13.92 0.26 #dY:1.9% + 25E-03 1E-03 13.4 0.16 #dY:1.2% + 27E-03 1E-03 14.08 0.3 #dY:2.1% + 29E-03 1E-03 14.1 0.29 #dY:2.1% + 31E-03 1E-03 13.55 0.14 #dY:1.0% + 33E-03 1E-03 12.94 0.13 #dY:1.0% + 35E-03 1E-03 13.24 0.21 #dY:1.6% + 37E-03 1E-03 13.32 0.31 #dY:2.3% + 39E-03 1E-03 13.19 0.18 #dY:1.4% + 41.5E-03 1.5E-03 13.31 0.27 #dY:2.0% + 44.5E-03 1.5E-03 12.99 0.22 #dY:1.7% + 47.5E-03 1.5E-03 12.68 0.13 #dY:1.0% + 50.5E-03 1.5E-03 12.9 0.2 #dY:1.6% + 54E-03 2E-03 13.04 0.14 #dY:1.1% + 58E-03 2E-03 12.98 0.25 #dY:1.9% + 62E-03 2E-03 12.82 0.33 #dY:2.6% + 66E-03 2E-03 12.86 0.21 #dY:1.6% + 70E-03 2E-03 12.37 0.2 #dY:1.6% + 74E-03 2E-03 12.3 0.23 #dY:1.9% + 78E-03 2E-03 12.33 0.2 #dY:1.6% + 82E-03 2E-03 12.1 0.16 #dY:1.3% diff --git a/238U.inp b/238U.inp index cfddb7c..b54e4c6 100644 --- a/238U.inp +++ b/238U.inp @@ -1,42 +1,45 @@ &Target Z=92.0 A=238.0 eferm=-5.4803 - BETA20=0.231 BETA40=0.062 BETA60=-9.6000E-03 + BETA20=0.223 BETA40=0.075 BETA60=-6.2000E-03 nstat=18 / -&State Et=0.0000E+00 Jval=0.0 Kval=0.0 NBAND=1 COEFF=0.2310 / -&State Et=4.4900E+01 Jval=2.0 Kval=0.0 NBAND=1 COEFF=0.2310 / -&State Et=1.4840E+02 Jval=4.0 Kval=0.0 NBAND=1 COEFF=0.2310 / -&State Et=3.0740E+02 Jval=6.0 Kval=0.0 NBAND=1 COEFF=0.2310 / -&State Et=5.1740E+02 Jval=8.0 Kval=0.0 NBAND=1 COEFF=0.2310 / -&State Et=7.7590E+02 Jval=10.0 Kval=0.0 NBAND=1 COEFF=0.231 / +&State Et=0.0000E+00 Jval=0.0 Kval=0.0 NBAND=1 COEFF=0.2230 / +&State Et=4.4900E+01 Jval=2.0 Kval=0.0 NBAND=1 COEFF=0.2230 / +&State Et=1.4840E+02 Jval=4.0 Kval=0.0 NBAND=1 COEFF=0.2230 / +&State Et=3.0740E+02 Jval=6.0 Kval=0.0 NBAND=1 COEFF=0.2230 / +&State Et=5.1740E+02 Jval=8.0 Kval=0.0 NBAND=1 COEFF=0.2230 / +&State Et=7.7590E+02 Jval=10.0 Kval=0.0 NBAND=1 COEFF=0.223 / -&State Et=6.7980E+02 Jval=1.0 Kval=0.0 NBAND=-2 COEFF=0.062 / -&State Et=7.3130E+02 Jval=3.0 Kval=0.0 NBAND=-2 COEFF=0.062 / -&State Et=8.2670E+02 Jval=5.0 Kval=0.0 NBAND=-2 COEFF=0.062 / -&State Et=9.6640E+02 Jval=7.0 Kval=0.0 NBAND=-2 COEFF=0.062 / -&State Et=1.1507E+03 Jval=9.0 Kval=0.0 NBAND=-2 COEFF=0.062 / +&State Et=6.7980E+02 Jval=1.0 Kval=0.0 NBAND=-2 COEFF=0.053966 / +&State Et=7.3130E+02 Jval=3.0 Kval=0.0 NBAND=-2 COEFF=0.053966 / +&State Et=8.2670E+02 Jval=5.0 Kval=0.0 NBAND=-2 COEFF=0.053966 / +&State Et=9.6640E+02 Jval=7.0 Kval=0.0 NBAND=-2 COEFF=0.053966 / +&State Et=1.1507E+03 Jval=9.0 Kval=0.0 NBAND=-2 COEFF=0.053966 / -&State Et=9.9300E+02 Jval=0.0 Kval=0.0 NBAND=3 COEFF=0.0240 / -&State Et=1.0373E+03 Jval=2.0 Kval=0.0 NBAND=3 COEFF=0.0240 / +&State Et=9.9300E+02 Jval=0.0 Kval=0.0 NBAND=3 COEFF=0.014049 / +&State Et=1.0373E+03 Jval=2.0 Kval=0.0 NBAND=3 COEFF=0.014049 / -&State Et=9.2700E+02 Jval=0.0 Kval=0.0 NBAND=4 COEFF=0.0110 / -&State Et=9.6630E+02 Jval=2.0 Kval=0.0 NBAND=4 COEFF=0.0110 / -&State Et=1.0564E+03 Jval=4.0 Kval=0.0 NBAND=4 COEFF=0.0110 / +&State Et=9.2700E+02 Jval=0.0 Kval=0.0 NBAND=4 COEFF=0.0016056 / +&State Et=9.6630E+02 Jval=2.0 Kval=0.0 NBAND=4 COEFF=0.0016056 / +&State Et=1.0564E+03 Jval=4.0 Kval=0.0 NBAND=4 COEFF=0.0016056 / -&State Et=1.0603E+03 Jval=2.0 Kval=2.0 NBAND=5 COEFF=0.0700 / -&State Et=1.1057E+03 Jval=3.0 Kval=2.0 NBAND=5 COEFF=0.0700 / +&State Et=1.0603E+03 Jval=2.0 Kval=2.0 NBAND=5 COEFF=0.004906 / +&State Et=1.1057E+03 Jval=3.0 Kval=2.0 NBAND=5 COEFF=0.004906 / -&Fresco hcm=0.10 rmatch= 20.00 Jmax=20 - Ngrid=41 pythonFlag=0 Nenergy=2 absend= 0.001 - relativ= 2 / +&Fresco hcm=0.10 rmatch= 20.00 Jmax=30 + Ngrid=41 pythonFlag=0 Nenergy=37 absend= 0.001 + relativ= 3 / &Energies - elab=1.000,2.0000 / - -&Potential v0a=50.41 v0b=0.0292 lambdhf= 0.00977 - cviso=17.50 vspo=5.64 lambdso=0.005 ccoul=1.36 - av=11.91 bv=81.86 w0=17.70 bs=10.85 wspo=-3.1 - bso=160.0 ea=55.0 alphav=0.355 cs=0.01331 cwiso=29.0 - adv=238.0 rhfa=1.2490 rhfb=-0.00171 ahfa=0.638 ahfb=0.00219 - rv=1.2657 ava=0.6960 avb=-0.00021 rsa=1.1701 rsb=0.0041 - as=0.617 rso=1.1214 aso=0.59 rc=1.1289 ac=0.547 / + elab=0.001,0.002,0.004,0.006,0.008,0.01,0.02,0.04, + 0.06,0.08,0.100,0.150,0.200,0.30,0.400,0.600,0.800,1.0000,1.40,2.0000, + 2.20,2.70,3.20,3.70,4.2,4.8,5.00, + 7.00,10.0000,12.00,15.00,17.00,20.00,23.00,26.00,28.00,30.00 / + +&Potential v0a=50.59 v0b=0.0608 lambdhf= 0.00982 + cviso=16.79 vspo=6.83 lambdso=0.005 ccoul=1.43 + av=11.66 bv=81.61 w0a=17.30 w0b=-0.031 bs=10.55 wspo=-3.1 + bso=160.0 ea=55.0 alphav=0.355 cs=0.01343 cwiso=28.96 + adv=238.0 rhfa=1.2426 rhfb=-0.00246 ahfa=0.638 ahfb=0.00306 + rv=1.2704 avv=0.6970 rsa=1.1860 rsb=0.0060 + as=0.616 rso=1.0156 aso=0.549 rc=1.2262 ac=0.674 / / diff --git a/239Pu.exp b/239Pu.exp new file mode 100644 index 0000000..948f933 --- /dev/null +++ b/239Pu.exp @@ -0,0 +1,114 @@ + 0.048 0.008022 13.11 0.2857 #dY:2.2% + 0.063 0.00805 12.56 0.2708 #dY:2.2% + 0.078 0.008022 12.37 0.266 #dY:2.2% + 0.092 0.00804 12.44 0.2653 #dY:2.1% + 0.111 0.008062 11.8 0.2598 #dY:2.2% + 0.137 0.008122 11.46 0.2569 #dY:2.2% + 0.172 0.008246 11.04 0.2546 #dY:2.3% + 0.195 0.008476 10.75 0.2075 #dY:1.9% + 0.223 0.008509 10.46 0.2515 #dY:2.4% + 0.244 0.008443 10.17 0.2296 #dY:2.3% + 0.297 0.008443 10.11 0.1894 #dY:1.9% + 0.35 0.008476 9.767 0.1896 #dY:1.9% + 0.401 0.008443 9.392 0.1873 #dY:2.0% + 0.458 0.008443 8.957 0.1144 #dY:1.3% + 0.51 0.008693 8.526 0.1024 #dY:1.2% + 0.568 0.008654 8.31 0.1607 #dY:1.9% + 0.61 0.008693 8.123 0.09793 #dY:1.2% + 0.67 0.008616 7.912 0.171 #dY:2.2% + 0.715 0.00858 7.674 0.1084 #dY:1.4% + 0.773 0.008616 7.653 0.1756 #dY:2.3% + 0.821 0.008654 7.365 0.0974 #dY:1.3% + 0.924 0.008693 7.084 0.1048 #dY:1.5% + 1.023 0.008544 7.054 0.08922 #dY:1.3% + 1.128 0.008509 7.033 0.0938 #dY:1.3% + 1.23 0.008476 6.945 0.09969 #dY:1.4% + 1.333 0.008443 7.013 0.1126 #dY:1.6% + 1.533 0.008246 7.099 0.08775 #dY:1.2% + 1.631 0.008179 7.083 0.1021 #dY:1.4% + 1.734 0.008223 7.097 0.08914 #dY:1.3% + 1.825 0.008158 7.098 0.1049 #dY:1.5% + 1.931 0.008223 7.297 0.09155 #dY:1.3% + 2.025 0.008139 7.279 0.1084 #dY:1.5% + 2.134 0.0105 7.284 0.09446 #dY:1.3% + 2.279 0.009381 7.377 0.1078 #dY:1.5% + 2.429 0.01124 7.519 0.09771 #dY:1.3% + 2.574 0.009881 7.627 0.1112 #dY:1.5% + 2.791 0.009035 7.74 0.1102 #dY:1.4% + 2.994 0.009381 7.971 0.1122 #dY:1.4% + 3.188 0.009278 7.823 0.1155 #dY:1.5% + 3.387 0.009228 8.056 0.1174 #dY:1.5% + 3.585 0.00913 7.975 0.1194 #dY:1.5% + 3.781 0.008989 7.814 0.1209 #dY:1.5% + 3.974 0.0089 7.984 0.1265 #dY:1.6% + 4.185 0.008773 7.792 0.1451 #dY:1.9% + 4.389 0.008693 8.187 0.1933 #dY:2.4% + 4.598 0.008297 7.874 0.1264 #dY:1.6% + 4.807 0.008297 7.944 0.1382 #dY:1.7% + 1.818 0.00909 7.193 0.08691 #dY:1.2% + 1.856 0.009278 7.26 0.08772 #dY:1.2% + 1.894 0.009472 7.256 0.08767 #dY:1.2% + 1.935 0.009673 7.298 0.08818 #dY:1.2% + 1.976 0.009881 7.366 0.089 #dY:1.2% + 2.019 0.0101 7.381 0.08918 #dY:1.2% + 2.063 0.01032 7.393 0.08933 #dY:1.2% + 2.109 0.01055 7.449 0.08782 #dY:1.2% + 2.157 0.01078 7.518 0.08864 #dY:1.2% + 2.206 0.01103 7.495 0.08836 #dY:1.2% + 2.257 0.01129 7.578 0.08934 #dY:1.2% + 2.31 0.01155 7.606 0.08967 #dY:1.2% + 2.365 0.01182 7.604 0.08965 #dY:1.2% + 2.422 0.01211 7.654 0.09024 #dY:1.2% + 2.48 0.0124 7.711 0.09091 #dY:1.2% + 2.542 0.01271 7.749 0.09136 #dY:1.2% + 2.605 0.01303 7.786 0.0918 #dY:1.2% + 2.671 0.01336 7.757 0.09145 #dY:1.2% + 2.74 0.0137 7.87 0.0911 #dY:1.2% + 2.811 0.01406 7.87 0.0911 #dY:1.2% + 2.886 0.01443 7.947 0.09199 #dY:1.2% + 2.963 0.01482 7.971 0.09227 #dY:1.2% + 3.044 0.01522 7.936 0.09187 #dY:1.2% + 3.128 0.01564 7.99 0.09249 #dY:1.2% + 3.216 0.01608 8.039 0.09306 #dY:1.2% + 3.308 0.01654 8.05 0.09143 #dY:1.1% + 3.403 0.01702 8.012 0.09275 #dY:1.2% + 3.503 0.01752 8.064 0.09335 #dY:1.2% + 3.608 0.01804 8.054 0.09148 #dY:1.1% + 3.717 0.01859 8.075 0.09171 #dY:1.1% + 3.832 0.01916 8.099 0.09199 #dY:1.1% + 3.952 0.01976 8.084 0.09182 #dY:1.1% + 4.078 0.02039 8.083 0.09357 #dY:1.2% + 4.21 0.02105 8.007 0.0944 #dY:1.2% + 4.349 0.02174 7.983 0.09479 #dY:1.2% + 4.495 0.02247 7.943 0.09598 #dY:1.2% + 4.648 0.02324 7.89 0.09823 #dY:1.2% + 4.81 0.02405 7.851 0.09774 #dY:1.2% + 4.981 0.0249 7.748 0.09646 #dY:1.2% + 5.161 0.0258 7.693 0.09578 #dY:1.2% + 5.351 0.02675 7.636 0.09507 #dY:1.2% + 5.551 0.02776 7.615 0.09481 #dY:1.2% + 5.764 0.02882 7.522 0.09365 #dY:1.2% + 5.989 0.02995 7.276 0.09059 #dY:1.2% + 6.228 0.03114 7.153 0.09216 #dY:1.3% + 6.482 0.03241 7.087 0.09131 #dY:1.3% + 6.751 0.03376 6.994 0.09011 #dY:1.3% + 7.038 0.03519 6.882 0.08867 #dY:1.3% + 7.344 0.03672 6.77 0.08723 #dY:1.3% + 7.671 0.03836 6.652 0.089 #dY:1.3% + 8.02 0.0401 6.526 0.08731 #dY:1.3% + 8.394 0.04197 6.38 0.0856 #dY:1.3% + 8.795 0.04398 6.367 0.08542 #dY:1.3% + 9.226 0.04613 6.265 0.08749 #dY:1.4% + 9.69 0.04845 6.216 0.0868 #dY:1.4% + 10.19 0.1019 6.133 0.0893 #dY:1.5% + 10.73 0.1073 6.051 0.0881 #dY:1.5% + 11.31 0.1131 5.994 0.0837 #dY:1.4% + 11.95 0.1195 5.965 0.08003 #dY:1.3% + 12.64 0.1264 5.987 0.07737 #dY:1.3% + 13.39 0.1339 5.968 0.0743 #dY:1.2% + 14.21 0.1421 6.052 0.07535 #dY:1.2% + 15.11 0.1511 6.114 0.07388 #dY:1.2% + 17.19 0.1719 6.142 0.07421 #dY:1.2% + 18.39 0.1839 6.22 0.07744 #dY:1.2% + 19.73 0.1973 6.205 0.08325 #dY:1.3% + 20.91 0.2091 6.295 0.1193 #dY:1.9% diff --git a/239Pu.inp b/239Pu.inp index 1313345..b62069a 100644 --- a/239Pu.inp +++ b/239Pu.inp @@ -1,44 +1,44 @@ &Target Z=94.0 A=239.0 eferm=-6.0902 - BETA20=0.236 BETA40=0.086 BETA60=-3.1000E-02 + BETA20=0.226 BETA40=0.092 BETA60=-5.0000E-03 nstat=19 / -&State Et=0.0000E+00 Jval=0.5 Kval=0.5 NBAND=1 COEFF=0.236 / -&State Et=7.8600E+00 Jval=1.5 Kval=0.5 NBAND=1 COEFF=0.236 / -&State Et=5.7270E+01 Jval=2.5 Kval=0.5 NBAND=1 COEFF=0.236 / -&State Et=7.5700E+01 Jval=3.5 Kval=0.5 NBAND=1 COEFF=0.236 / -&State Et=1.6376E+02 Jval=4.5 Kval=0.5 NBAND=1 COEFF=0.236 / -&State Et=1.9280E+02 Jval=5.5 Kval=0.5 NBAND=1 COEFF=0.236 / -&State Et=3.1850E+02 Jval=6.5 Kval=0.5 NBAND=1 COEFF=0.236 / -&State Et=3.5810E+02 Jval=7.5 Kval=0.5 NBAND=1 COEFF=0.236 / -&State Et=5.1930E+02 Jval=8.5 Kval=0.5 NBAND=1 COEFF=0.236 / +&State Et=0.0000E+00 Jval=0.5 Kval=0.5 NBAND=1 COEFF=0.226 / +&State Et=7.8600E+00 Jval=1.5 Kval=0.5 NBAND=1 COEFF=0.226 / +&State Et=5.7270E+01 Jval=2.5 Kval=0.5 NBAND=1 COEFF=0.226 / +&State Et=7.5700E+01 Jval=3.5 Kval=0.5 NBAND=1 COEFF=0.226 / +&State Et=1.6376E+02 Jval=4.5 Kval=0.5 NBAND=1 COEFF=0.226 / +&State Et=1.9280E+02 Jval=5.5 Kval=0.5 NBAND=1 COEFF=0.226 / +&State Et=3.1850E+02 Jval=6.5 Kval=0.5 NBAND=1 COEFF=0.226 / +&State Et=3.5810E+02 Jval=7.5 Kval=0.5 NBAND=1 COEFF=0.226 / +&State Et=5.1930E+02 Jval=8.5 Kval=0.5 NBAND=1 COEFF=0.226 / -&State Et=2.8551E+02 Jval=2.5 Kval=2.5 NBAND=2 COEFF=0.025 / -&State Et=3.3010E+02 Jval=3.5 Kval=2.5 NBAND=2 COEFF=0.025 / -&State Et=3.8740E+02 Jval=4.5 Kval=2.5 NBAND=2 COEFF=0.025 / -&State Et=4.6200E+02 Jval=5.5 Kval=2.5 NBAND=2 COEFF=0.025 / +&State Et=2.8551E+02 Jval=2.5 Kval=2.5 NBAND=2 COEFF=0.01695 / +&State Et=3.3010E+02 Jval=3.5 Kval=2.5 NBAND=2 COEFF=0.01695 / +&State Et=3.8740E+02 Jval=4.5 Kval=2.5 NBAND=2 COEFF=0.01695 / +&State Et=4.6200E+02 Jval=5.5 Kval=2.5 NBAND=2 COEFF=0.01695 / -&State Et=4.6980E+02 Jval=0.5 Kval=0.5 NBAND=-3 COEFF=0.062 / -&State Et=4.9210E+02 Jval=1.5 Kval=0.5 NBAND=-3 COEFF=0.062 / -&State Et=5.0560E+02 Jval=2.5 Kval=0.5 NBAND=-3 COEFF=0.062 / -&State Et=5.5620E+02 Jval=3.5 Kval=0.5 NBAND=-3 COEFF=0.062 / -&State Et=5.8300E+02 Jval=4.5 Kval=0.5 NBAND=-3 COEFF=0.062 / -&State Et=6.6110E+02 Jval=5.5 Kval=0.5 NBAND=-3 COEFF=0.062 / +&State Et=4.6980E+02 Jval=0.5 Kval=0.5 NBAND=-3 COEFF=0.051076 / +&State Et=4.9210E+02 Jval=1.5 Kval=0.5 NBAND=-3 COEFF=0.051076 / +&State Et=5.0560E+02 Jval=2.5 Kval=0.5 NBAND=-3 COEFF=0.051076 / +&State Et=5.5620E+02 Jval=3.5 Kval=0.5 NBAND=-3 COEFF=0.051076 / +&State Et=5.8300E+02 Jval=4.5 Kval=0.5 NBAND=-3 COEFF=0.051076 / +&State Et=6.6110E+02 Jval=5.5 Kval=0.5 NBAND=-3 COEFF=0.051076 / -&Fresco hcm=0.10 rmatch= 20.00 Jmax=20 +&Fresco hcm=0.10 rmatch= 20.00 Jmax=30 Ngrid=41 pythonFlag=0 Nenergy=37 absend= 0.001 - relativ= 2 / + relativ= 3 / &Energies elab=0.001,0.002,0.004,0.006,0.008,0.01,0.02,0.04, - 0.06,0.08,0.100,0.200,0.400,0.600,0.800,1.0000,1.3,1.40,1.5,2.0000, + 0.06,0.08,0.100,0.150,0.200,0.30,0.400,0.600,0.800,1.0000,1.40,2.0000, 2.20,2.70,3.20,3.70,4.2,4.8,5.00, 7.00,10.0000,12.00,15.00,17.00,20.00,23.00,26.00,28.00,30.00 / -&Potential v0a=50.41 v0b=0.0292 lambdhf= 0.00977 - cviso=17.5 vspo=5.64 lambdso=0.005 ccoul=1.36 - av=11.91 bv=81.86 w0=17.70 bs=10.85 wspo=-3.1 - bso=160.0 ea=55.0 alphav=0.355 cs=0.01331 cwiso=29.0 - adv=238.0 rhfa=1.2490 rhfb=-0.00171 ahfa=0.638 ahfb=0.00219 - rv=1.2657 ava=0.6960 avb=-0.00021 rsa=1.1701 rsb=0.0041 - as=0.617 rso=1.1214 aso=0.59 rc=1.1289 ac=0.547 / +&Potential v0a=50.59 v0b=0.0608 lambdhf= 0.00982 + cviso=16.79 vspo=6.83 lambdso=0.005 ccoul=1.43 + av=11.66 bv=81.61 w0a=17.30 w0b=-0.031 bs=10.55 wspo=-3.1 + bso=160.0 ea=55.0 alphav=0.355 cs=0.01343 cwiso=28.96 + adv=238.0 rhfa=1.2426 rhfb=-0.00246 ahfa=0.638 ahfb=0.00306 + rv=1.2704 avv=0.6970 rsa=1.1860 rsb=0.0060 + as=0.616 rso=1.0156 aso=0.549 rc=1.2262 ac=0.674 / / diff --git a/exp.py b/exp.py new file mode 100644 index 0000000..6be0a59 --- /dev/null +++ b/exp.py @@ -0,0 +1,38 @@ +import numpy as np +import matplotlib.pyplot as plt + + + +filename= "235U.exp" +filename2 = "239Pu.new.xsec_3d_no_j_var" +filename3 = "232Th.new.xsec_3d_no_j_var" + +x,dx,y,dy = np.loadtxt(filename,usecols=(0,1,2,3),unpack=True) +ener, xsec = np.loadtxt(filename2,usecols=(0,4),unpack=True) +a, b = np.loadtxt(filename3,usecols=(0,4),unpack=True) +for i in range(len(x)): + y[i]=y[i]*1000 # a mb + dy[i]=dy[i]*1000 +fig, ax = plt.subplots() + +ratio=[] + +for i in range(len(ener)): + ratio.append(2*(xsec[i]-b[i])/(xsec[i]+b[i])) +#ax.plot(ener,xsec,'--b',zorder=4) +#ax.plot(a,b,"-r",zorder=3) +#ax.errorbar(x,y,xerr=dx,yerr=dy,fmt='.g',zorder=1) +ax.plot(ener,ratio,'--r') +#---------------------------------- +ax.set_xlim([3*10**-2,25]) +ax.set_ylim([-0.07,0.07]) +ax.set_xscale('log') +#ax.set_yscale('log') +ax.set_title(r"R($^{239}$Pu, $^{232}$Th)") +ax.set_xlabel("Energy (MeV)") +ax.set_ylabel(r"R($^{239}$Pu, $^{232}$Th)") +ax.legend(["FRESCO V3.3"]) +ax.grid() +#------------------------------- +#fig.savefig("R-239Pu-232Th.pdf",format="pdf") +plt.show() \ No newline at end of file diff --git a/frescoPRC.f90 b/frescoPRC.f90 index 3201982..c8f2139 100644 --- a/frescoPRC.f90 +++ b/frescoPRC.f90 @@ -9,7 +9,7 @@ CHARACTER*20 input_file CHARACTER*2 SYMBOL(mxsym) !&Target and &fresco - INTEGER Nenergy,nstat,NBAND,Ngrid,Jmax, relativ!,absend + INTEGER Nenergy,nstat,NBAND,Ngrid,Jmax, relativ REAL absend REAL Z,A,eferm REAL BETA20,BETA40,BETA60 @@ -17,8 +17,8 @@ REAL Et,Jval,Kval,COEFF ! Potential REAL v0a,v0b,lambdhf,cviso,vspo,lambdso,ccoul - REAL av,bv,w0,bs,wspo,bso,ea,alphav,cs,cwiso,adv - REAL rhfa,rhfb,ahfa,ahfb,rv,ava,avb,rsa,rsb,as,rso,aso,rc,ac + REAL av,bv,w0a,w0b,bs,wspo,bso,ea,alphav,cs,cwiso,adv + REAL rhfa,rhfb,ahfa,ahfb,rv,avv,rsa,rsb,as,rso,aso,rc,ac !!!!! INTEGER i,sum_neg,sum_pos,gv,ii, j_cal DOUBLE PRECISION k_number, nuc_rad @@ -41,8 +41,8 @@ NAMELIST /energies/ elab NAMELIST /state/ Et,Jval,Kval,NBAND,COEFF NAMELIST /potential/ v0a,v0b,lambdhf,cviso,vspo,lambdso,ccoul, & - av,bv,w0,bs,wspo,bso,ea,alphav,cs,cwiso,adv, & - rhfa,rhfb,ahfa,ahfb,rv,ava,avb,rsa,rsb,as, & + av,bv,w0a,w0b,bs,wspo,bso,ea,alphav,cs,cwiso,adv, & + rhfa,rhfb,ahfa,ahfb,rv,avv,rsa,rsb,as, & rso,aso,rc,ac !------------------------------------------------------ data pottype / 'REAL_VOLUME', 'REAL_VOLUME', 'IMAG_VOLUME','REAL_SURFACE', & @@ -57,20 +57,14 @@ 1021 FORMAT (300(A2,1X)) 450 CONTINUE call get_command_argument(1,input_file) - !WRITE(6,*) 'Write the name of the input file:' - !READ(*,*) input_file OPEN(40,STATUS='old',FILE=input_file,IOSTAT=err) IF(err .NE. 0) THEN WRITE(6,*) 'Error reading the name of the input file.' + WRITE(6,*) "if your input is 'input.inp', you generate FRESCO's inputs typing in the terminal './frescoPRC input.inp'" WRITE(6,*) 'Please, check the name of the file [it must include the extension].' - !GO TO 450 STOP ENDIF - !Reading input - !Old input style commented - !READ(40,'(F5.1,F10.5,F7.4)') Z,A,eferm - !READ(40,'(I2)') nstat READ (40, NML=target, END=761, IOSTAT=ios, ERR=761 ) 761 IF (ios .ne. 0) then WRITE(*,*) 'Input read error while reading Target: ', ios @@ -80,7 +74,6 @@ KBAND(nstat),or_val(nstat),indexx(nstat), STAT=err) CALL error(err,1) DO i=1,nstat - !READ(40,'(E12.5,F5.2,F4.2,I3,F7.5)') Ener_levels(i),J_val(i),KBAND(i),BAND(i),BETA_EFF(i) READ (40, NML=state, END=762, IOSTAT=iosss, ERR=762 ) 762 IF (iosss .ne. 0) THEN WRITE(*,*) 'Input read error while reading States: ', iosss @@ -91,7 +84,6 @@ indexx(i) = i ENDDO Ener_levels = Ener_levels/1000 !Reading in KeV but FRESCO reads it in MeV. - !READ(40,'(I3,F6.2,F7.2,I4,I3)') Jmax,hcm,rmatch,Ngrid,Nenergy READ (40, NML=fresco, END=763, IOSTAT=ioss, ERR=763 ) 763 IF (ioss .ne. 0) then WRITE(*,*) 'Input read error while reading Fresco: ', ioss @@ -100,10 +92,9 @@ ALLOCATE(elab(Nenergy)) IF (relativ.NE.1 .AND. relativ.NE.2 .AND. relativ.NE.3 ) THEN - WRITE(6,*) 'ERROR: relativ flag invalid, use: 1,2 o 3' + WRITE(6,*) 'ERROR: relativ flag invalid, use: 1 (No relativistic correction),2 (fresco3.3) o 3 (frx6j)' STOP ENDIF - !READ(40,'(5E12.5)') (elab(i), i=1, Nenergy) READ (40, NML=energies, END=764, IOSTAT=iosss, ERR=764 ) 764 IF (iosss .ne. 0) then WRITE(*,*) 'Input read error while reading Energies: ', iosss @@ -118,14 +109,6 @@ WRITE(*,*) 'Input read error while reading potential: ', iosss STOP ENDIF - !READ(40,'(7E12.4)') v0a,v0b,lambdhf,cviso,vspo,lambdso,ccoul - !READ(40,'(6E12.4)') av,bv,w0,bs,Wso0,BBso - !READ(40,'(5E12.4)') Ea,alpha,CCs,Cwiso,Ades - !READ(40,'(7E12.4)') rHFl,rHFdep,aHFl,aHFdep,rv,avl,avdep - !READ(40,'(3E12.4)') rsl,rsdep,as - !!READ(40,'(4E12.4)') rso,aso,rc,ac - !READ(40,'(3E12.4)') BETA20,BETA40,BETA60 - !READ(40,'(I1)') pythonFlag BETA_EFF = BETA_EFF/BETA20 ! input like OPTMAN, with x\beta_{20} factor. !Reading done @@ -181,7 +164,7 @@ ALLOCATE(jexc(n_exc),index_exc(n_exc),exc_Bparity(n_exc), & exc_KBAND(n_exc),BETA_PAR(n_exc),STAT=err) CALL error(1,err) - gv = 1; ii = 1 !indices + gv = 1; ii = 1 DO i=1,nstat IF(ABS(BAND(i))==1) THEN ! G.S band MUST be |NBAND|=1 in the input. po=>jgsval(gv) @@ -220,9 +203,9 @@ CALL dispers2(A,Z,NTYPE,E,VR,RR,AR, dv,drv,dav, VD,RVD,AVD, & W,RW,AW, WD,RD,AD, VSO,RSO,ASO, dvso, WSO,WRSO,WASO, & v0a,v0b,lambdhf,cviso,vspo,lambdso,ccoul, & - av,bv,w0,bs,cs,cwiso,wspo,bso, & + av,bv,w0a,w0b,bs,cs,cwiso,wspo,bso, & ea,alphav,eferm,adv, & - rhfa,rhfb,ahfa,ahfb,rv,ava,avb, & + rhfa,rhfb,ahfa,ahfb,rv,avv, & rsa,rsb,as, & rso,aso,rc,ac) RVOL = ACroot * RR @@ -230,15 +213,15 @@ RSURF = ACroot * RD k_number = k_val(E,A) nuc_rad = 1.3d0*((NINT(A) + 1)**(1.0/3.0)) - j_cal = NINT(k_number*nuc_rad*2.2) !!! Solving problem with absend in FRESCO V3.3 + j_cal = NINT(k_number*nuc_rad*2.4) !!! Solving the problem with absend in FRESCO V3.3?? - IF (j_cal .LT. 10) THEN - Jmax = 10 - ELSE IF ( j_cal .GE. 10 ) THEN - Jmax = j_cal - ENDIF + !IF (j_cal .LE. 8) THEN + ! Jmax = 14 + !ELSE IF ( j_cal .GT. 8 ) THEN + ! Jmax = j_cal + !ENDIF - WRITE(66,*) E, k_number*nuc_rad,NINT(k_number*nuc_rad*2.2) + WRITE(66,*) E, k_number*nuc_rad,NINT(k_number*nuc_rad*2.4) WRITE(10,10) E,VR,RR,AR, dv,drv,dav, W,RW,AW, VD,RVD,AVD, WD,RD,AD,VSO,RSO,ASO, dvso,WSO,WRSO,WASO, RC,AC 10 FORMAT(f7.3, 6(f8.3,2f6.3),2f8.3,2f6.3,2f6.3) fname = 'fresco-00-'//POTL//'-s'//CHAR(ICHAR('0')+nexe)//',o'//CHAR(ICHAR('0')+sum_neg)//'-E0000000.in' @@ -279,7 +262,6 @@ ENDIF ENDIF !!!!!!!!! - !WRITE(1,'(a,i4,a,i2)') ' jtmin= 0.0 jtmax=',Jmax,' absend= ',absend 888 CONTINUE WRITE(1,755) Jmax, absend 755 FORMAT (' jtmin= 0.0 jtmax=',i4,' absend= ',f10.6) diff --git a/fun.f95 b/fun.f95 index 66266d9..fbae701 100644 --- a/fun.f95 +++ b/fun.f95 @@ -7,7 +7,7 @@ MODULE modulo !* Subroutine wigner: Compute matrix element of Wigner's functions accoding to B.14 and B.19 * !* from PRC 94 6 (2016), 064605 (even and odd nuclei). * !* * - !* Dispersive functions pack: Calculation of Analytical dispersive integrals. * + !* Dispersive functions pack: Calculation of Analytical dispersive integrals by Roberto Capote * !* * !* Subroutine FORMFACT: Numerical calculation of form factors using * !* Gauss-Legendre quadrature. Those form factors are part of * @@ -15,7 +15,7 @@ MODULE modulo !* * !* Subroutine steps: Build the couplings between G.S band and excited bands for * !* even and odd nuclei according to equations B.12 and B.16 from * - !* PRC 94 6 (2016), 064605. * * + !* PRC 94 6 (2016), 064605. * !*********************************************************************************************** CONTAINS @@ -414,9 +414,9 @@ SUBROUTINE dispers2(A,Z,k,eopt, & v,rvv,avv, dv,drv,dav, dvs,drs,das, w,rw,aw, wd,rwd,awd, & vso,rvso,avso,dvso, wso,rwso,awso, & Vlin,Vdep,lambdaHF,Cviso,Vso0,lambdaso,Ccoul, & - AAv,BBv,W0,BBs,CCs,Cwiso,Wso0,BBso, & + AAv,BBv,W0l,W0dep,BBs,CCs,Cwiso,Wso0,BBso, & Ea,alpha,eferm,Ades, & - rHFl,rHFdep,aHFl,aHFdep,rv,avl,avdep, & + rHFl,rHFdep,aHFl,aHFdep,rv,avpot, & rsl,rsdep,as, & rso,aso,rc,ac) @@ -441,7 +441,8 @@ SUBROUTINE dispers2(A,Z,k,eopt, & V0 = Vlin + Vdep*Au rHF = rHFL + rHFdep * Au aHF = aHFl + aHFdep * Au - av = avl + avdep * Au + W0 = W0l + W0dep*Au + av = avpot rs = rsl + rsdep * Au AAHF = V0 * (1 + (-1)**k * Cviso*asym/V0) AAs = W0 * (1 + (-1)**k * Cwiso*asym/W0) @@ -809,7 +810,7 @@ SUBROUTINE Graphs CLOSE(94) END SUBROUTINE Graphs - DOUBLE PRECISION FUNCTION k_val(E,M) + DOUBLE PRECISION FUNCTION k_val(E,M) !wave number: trying to fix absend=-1.00 problem in FRESCO v3.3 DOUBLE PRECISION hc REAL E, M, mu diff --git a/runall.sh b/runall.sh index 9929645..c1f135a 100644 --- a/runall.sh +++ b/runall.sh @@ -5,7 +5,7 @@ # .in files are FRESCO's inputs and .form files are external form factors used by FRESCO. # #-----------------------------------------------------------------------------------------# # -FRESCO=fresco #Calling FRESCO +FRESCO=frescox #Calling FRESCO GRACE=python destino=fort.4 # .form file for each energy must be named fort.4 when running the correspondent .in file for that energy. XSECS2=xsec_states.out # Data from FRESCO's output fort.13 [Xsec to all excited states one by one]. @@ -16,8 +16,8 @@ lista=lista.txt # names of .in and .form files are here -> generated by frescoPR grace_file=graphs.py outfile=outfile.out # all output files for all energies are saved here. -archivo=prueba.txt #.txt with the name of all NAMELIST inputss -generador=frescoPRC #Modified generator to meet new input reading format (frescoPRC input.inp in console) +archivo=prueba.txt #.txt with the name of all NAMELIST inputs for the tool +generador=frescoPRC_no_j_var #calling frescoPRC tool to generate the inputs for each target in $archivo rm -f $destino; rm -f $XSECS; rm -f $XSECS2; rm -f $outfile; @@ -26,7 +26,7 @@ rm -f $destino; rm -f $XSECS; rm -f $XSECS2; rm -f $outfile; while read este; do ini=${este%%.*} - final=$ini.xsec4 + final=$ini.new.xsec_3d_no_j_var #echo $este >> $final #echo "Energy Elastic Absorption Reaction Total (MeV/mb)" > $final $generador "$este" @@ -57,7 +57,7 @@ while read este; do done < $lista #looping for each pair .in/.form in $lista naming the .form fort.4 so FRESCO identify it. rm -f $lista rm -f fort* - rm -f fresco-* #Borrar todos los inputs entre target y target + rm -f fresco-* #Borrar todos los inputs entre un target y el siguiente rm -f DOMEIC16* done < $archivo diff --git a/script.py b/script.py index 033e4ef..1fcbd6b 100644 --- a/script.py +++ b/script.py @@ -14,10 +14,10 @@ def graph(name, x, y): ax.set_yscale('log') ax.legend(["FRESCO 3.3"])#,loc='center right') ax.grid() - fig.savefig(title + ".pdf", format='pdf') + #fig.savefig(title + ".pdf", format='pdf') #ax.set(xlabel='Time (s)', ylabel='Voltage (mV)', # title='Plot example') - #plt.show() + plt.show() def graph2(name1, x, y, xx, yy): title = name1.replace('.xsec','') fig,ax = plt.subplots() @@ -27,12 +27,48 @@ def graph2(name1, x, y, xx, yy): ax.set_ylabel("Cross section (mb)") ax.set_xscale('log') ax.set_yscale('log') - ax.legend(["FRESCO frxy6j", "FRESCO 3.3"])#,loc='center right') + #ax.legend(["FRESCO 3.3", "FRESCO frx6j"])#,loc='center right') ax.grid() #ax.set(xlabel='Time (s)', ylabel='Voltage (mV)', # title='Plot example') - #plt.show() - fig.savefig(title + ".pdf", format='pdf') + plt.show() + #fig.savefig(title + ".pdf", format='pdf') +def graph3(name1, x, y, xx, yy, xxx, yyy): + title = name1.replace('.xsec','') + fig,ax = plt.subplots() + ax.plot(x, y, '--r',xx, yy, '-.b', xxx, yyy, '+k') + ax.set_title('n + ' + title + '\n Total cross section') + ax.set_xlabel("Energy (MeV)") + ax.set_ylabel("Cross section (mb)") + ax.set_xscale('log') + ax.set_yscale('log') + ax.legend(["FRESCO 3.3 abs=0.001 j=jmax", "FRESCO frxy6j abs=0.001", "FRESCO 3.3 abs=-1.000 j=j_cal"])#,loc='center right') + ax.grid() + #ax.set(xlabel='Time (s)', ylabel='Voltage (mV)', + # title='Plot example') + plt.show() + #fig.savefig(title + ".pdf", format='pdf') +def graph4(name1, x, y, xx, yy, xxx, yyy): + title = name1.replace('.new.xsec','') + for i in range(len(xxx)): + xxx[i]=xxx[i]*1000 # a mb + yyy[i]=yyy[i]*1000 + + fig,ax = plt.subplots() + ax.plot(x, y, '--k',zorder=2) + ax.errorbar(xx,xxx,xerr=yy,yerr=yyy,fmt='.c',zorder=1) + ax.set_title('n + ' + title + '\n Total cross section') + ax.set_xlabel("Energy (MeV)") + ax.set_ylabel("Cross section (mb)") + ax.set_xscale('log') + ax.set_yscale('log') + ax.set_xlim([np.amin(xx),30]) + ax.legend(["FRESCO V3.3","Experimental data"])#,loc='center right') + ax.grid() + #ax.set(xlabel='Time (s)', ylabel='Voltage (mV)', + # title='Plot example') + plt.show() + #fig.savefig(title + ".pdf", format='pdf') #fig, ax = plt.subplots() filename = "prueba.txt" #reading input's names @@ -42,14 +78,18 @@ def graph2(name1, x, y, xx, yy): #print(inputs) for input in inputs: - outfile = input.replace('.inp','.xsec') - outfile2 = input.replace('.inp','.xsec2') + outfile = input.replace('.inp','.new.xsec_3d_no_j_var') + outfile2 = input.replace('.inp','.exp') + #outfile2 = input.replace('.inp','.xsec4') + #outfile3 = input.replace('.inp','.xsec3') energies, total_xsec = np.loadtxt(outfile,usecols=(0,4),unpack=True) - energies2, total_xsec2 = np.loadtxt(outfile2,usecols=(0,4),unpack=True) + x,dx,y,dy = np.loadtxt(outfile2,usecols=(0,1,2,3),unpack=True) + #energies2, total_xsec2 = np.loadtxt(outfile2,usecols=(0,4),unpack=True) + #energies3, total_xsec3 = np.loadtxt(outfile3,usecols=(0,4),unpack=True) #graph(outfile, energies, total_xsec) - graph2(outfile, energies, total_xsec, energies2, total_xsec2) - - + #graph2(outfile, energies, total_xsec, energies2, total_xsec2) + #graph3(outfile, energies, total_xsec, energies2, total_xsec2, energies3, total_xsec3) + graph4(outfile,energies,total_xsec,x,dx,y,dy) #print(energies, total_xsec) #print(outfile) #np.rea(filename,)