OBNINSK_MICRO.readme
Code for calculation of microscopic total level densities
(Prepared in Obninsk, submitted by R. Capote 24 September 1997)
***************************************************************
Contents
--------
The file OBNINSK_MICRO.FOR contains a microscopical statistical code for
level density calculation. Collective enhancement of the level density
is considered phenomenologically.
We are including here three test input files as well as the corresponding
output obtained by the OBNINSK_MICRO.FOR. We should remark that files
Z50A120.SPE, Z92A238.SPE and Z40A90.SPE containing single-particle levels
schemes for all the three test cases were retrieved from the Moller-Nix
database (moller_levels.dat) using retrieval code moller_levels.for.
Test files
----------
**************************************************************************
TEST CASE #1: FILE SN120M.INP
_______________________________________________________________________________
Z50A120.SPE
GZ GN TMAX TMIN TSTEP
18.7 12. 4.2 .4 .2
_______________________________________________________________________________
TEST CASE #1: FILE SN120M.OUT
_______________________________________________________________________________
ELEMENT Z=50.000 MASS NUMBER A= 120.000 DEFORMATION BETA=-.0168
FORCES GZ = 18.70000 GN = 12.00000
Z/N G LAM0 DEL0 TCR ECOND E D2N
50.0000 .1558 -8.6878 -.0004 .0000 .0000 -947.7775 .0000
70.0000 .1000 -8.4441 1.2508 .6630 4.0461 -1403.7170 11.1973
T QVIB QROT U LnRO LnST D2A SI2 F|| A A1 A2 A3 M2
.40 1.15 1.00 1.20 .00 1.76 .86 5.38 13.4 3.52 -1.27 -17.78 4.75 6.28
.60 1.49 1.00 5.50 6.47 8.92 3.60 21.27 35.5 9.87 25.80 4.04 10.21 5.91
.80 1.98 1.00 9.86 12.11 14.76 5.75 31.69 39.6 11.82 15.01 9.08 11.68 5.51
1.00 2.50 1.00 13.91 16.02 18.74 7.07 36.81 36.8 11.63 13.62 9.86 11.59 5.21
1.20 2.92 1.00 18.91 20.10 22.89 8.43 42.29 35.2 11.55 12.93 10.32 11.55 5.02
1.40 3.12 1.00 24.89 24.26 27.11 9.80 48.02 34.3 11.52 12.53 10.63 11.54 4.90
1.60 3.04 1.00 31.89 28.47 31.38 11.20 53.99 33.7 11.51 12.29 10.88 11.56 4.82
1.80 2.75 1.00 39.95 32.74 35.71 12.61 60.20 33.4 11.52 12.13 11.08 11.59 4.77
2.00 2.36 1.00 49.08 37.06 40.08 14.04 66.63 33.3 11.55 12.02 11.26 11.63 4.74
2.20 1.98 1.00 59.27 41.44 44.50 15.49 73.24 33.3 11.58 11.95 11.41 11.68 4.73
2.40 1.67 1.00 70.51 45.88 48.99 16.94 80.01 33.3 11.61 11.91 11.54 11.72 4.72
2.60 1.43 1.00 82.79 50.40 53.55 18.40 86.90 33.4 11.64 11.89 11.65 11.77 4.72
2.80 1.27 1.00 96.14 54.99 58.18 19.86 93.88 33.5 11.67 11.87 11.75 11.81 4.73
3.00 1.16 1.00 110.54 59.65 62.88 21.33 100.96 33.7 11.69 11.87 11.83 11.85 4.73
3.20 1.09 1.00 126.01 64.38 67.64 22.79 108.12 33.8 11.71 11.87 11.91 11.89 4.74
3.40 1.05 1.00 142.54 69.16 72.45 24.25 115.37 33.9 11.73 11.87 11.98 11.93 4.76
3.60 1.03 1.00 160.13 73.99 77.31 25.71 122.69 34.1 11.75 11.88 12.04 11.96 4.77
3.80 1.01 1.00 178.75 78.84 82.19 27.18 130.10 34.2 11.76 11.89 12.10 11.99 4.79
4.00 1.01 1.00 198.38 83.71 87.09 28.64 137.58 34.4 11.78 11.90 12.15 12.02 4.80
*********************************************************************************************
_____________________________________________________________________________________________
TEST CASE #2: FILE U238M.INP
_____________________________________________________________________________________________
Z92A238.SPE
GZ GN TMAX TMIN TSTEP
20.5 14. 2. .2 .1
_____________________________________________________________________________________________
TEST CASE #2: FILE U238M.OUT
_____________________________________________________________________________________________
ELEMENT Z=92.000 MASS NUMBER A= 238.000 DEFORMATION BETA= .2050
FORCES GZ = 20.50000 GN = 14.00000
Z/N G LAM0 DEL0 TCR ECOND E D2N
92.0000 .0861 -6.9599 .8751 .5073 .7600 -1563.7030 6.9661
146.0000 .0588 -5.7530 .7076 .4052 .7944 -2722.9580 6.9437
T QVIB QROT U LnRO LnST D2A SI2 F|| A A1 A2 A3 M2
.20 1.06 .66 .22 -3.79 -2.78 .27 1.20 6.0 2.19 -.36 -33.24 3.45 4.51
.30 1.26 4.91 1.23 .83 2.64 1.30 5.96 19.9 7.15 -21.29 -3.65 8.82 4.57
.40 1.57 17.94 3.44 8.20 10.49 3.31 15.47 38.7 13.62 17.77 11.81 14.49 4.67
.50 1.97 28.32 6.06 14.35 16.84 5.13 23.01 46.0 16.87 16.77 18.01 17.38 4.49
.60 2.44 37.32 8.43 18.57 21.19 6.38 30.16 50.3 17.48 17.13 19.11 18.09 4.73
.70 2.90 47.48 11.31 22.85 25.60 7.64 38.37 54.8 17.96 17.50 19.91 18.67 5.02
.80 3.31 58.94 14.77 27.33 30.18 8.98 47.63 59.5 18.46 17.88 20.65 19.22 5.31
.90 3.59 71.40 18.79 31.90 34.85 10.37 57.71 64.1 18.95 18.25 21.28 19.71 5.57
1.00 3.69 84.13 23.35 36.51 39.55 11.79 68.37 68.4 19.40 18.61 21.79 20.14 5.80
1.10 3.62 85.12 28.40 40.99 44.09 13.24 79.38 72.2 19.80 18.93 22.19 20.49 6.00
1.20 3.39 17.40 33.94 43.85 47.02 14.70 90.56 75.5 20.15 19.22 22.49 20.79 6.16
1.30 3.08 1.43 39.93 45.78 49.01 16.17 101.77 78.3 20.46 19.49 22.71 21.04 6.29
1.40 2.72 1.01 46.36 49.82 53.11 17.63 112.92 80.7 20.72 19.73 22.86 21.24 6.40
1.50 2.37 1.00 53.24 54.19 57.52 19.09 123.94 82.6 20.93 19.94 22.97 21.40 6.49
1.60 2.05 1.00 60.53 58.55 61.92 20.54 134.80 84.3 21.12 20.13 23.04 21.53 6.56
1.70 1.79 1.00 68.26 62.91 66.31 21.98 145.50 85.6 21.27 20.30 23.08 21.65 6.62
1.80 1.58 1.00 76.43 67.26 70.71 23.42 156.04 86.7 21.40 20.45 23.11 21.74 6.66
1.90 1.42 1.00 85.02 71.63 75.10 24.84 166.41 87.6 21.51 20.59 23.12 21.82 6.70
*********************************************************************************************
_____________________________________________________________________________________________
TEST CASE #3: FILE ZR40M.INP
_____________________________________________________________________________________________
Z40A90.SPE
GZ GN TMAX TMIN TSTEP
10. 13. 4.2 .4 .2
_____________________________________________________________________________________________
TEST CASE #3: FILE ZR40M.OUT
_____________________________________________________________________________________________
ELEMENT Z=40.000 MASS NUMBER A= 90.000 DEFORMATION BETA=-.0101
FORCES GZ = 10.00000 GN = 13.00000
Z/N G LAM0 DEL0 TCR ECOND E D2N
40.0000 .1111 -6.5230 .0016 .0000 .0000 -696.8522 .0000
50.0000 .1444 -10.4488 .0035 .0000 .0000 -1082.0020 .0000
T QVIB QROT U LnRO LnST D2A SI2 F|| A A1 A2 A3 M2
.40 1.07 1.00 1.03 1.52 3.17 1.10 4.38 10.9 4.53 4.01 6.46 5.09 3.97
.60 1.29 1.00 2.73 2.70 4.64 2.00 7.74 12.9 5.48 5.06 7.59 6.19 3.87
.80 1.62 1.00 5.61 5.64 7.80 3.19 12.00 15.0 6.57 5.93 8.76 7.21 3.76
1.00 2.03 1.00 9.33 9.08 11.40 4.48 16.64 16.6 7.37 6.58 9.33 7.84 3.72
1.20 2.43 1.00 13.73 12.59 15.04 5.74 21.30 17.8 7.87 7.05 9.54 8.20 3.71
1.40 2.75 1.00 18.80 16.07 18.62 6.95 25.87 18.5 8.17 7.39 9.59 8.42 3.72
1.60 2.91 1.00 24.55 19.52 22.15 8.14 30.33 19.0 8.36 7.65 9.59 8.57 3.73
1.80 2.89 1.00 31.04 22.95 25.64 9.30 34.70 19.3 8.50 7.85 9.58 8.68 3.73
2.00 2.70 1.00 38.30 26.37 29.12 10.45 39.02 19.5 8.60 8.02 9.58 8.76 3.73
2.20 2.42 1.00 46.34 29.79 32.59 11.60 43.32 19.7 8.68 8.15 9.57 8.84 3.73
2.40 2.12 1.00 55.18 33.22 36.07 12.75 47.62 19.8 8.74 8.27 9.58 8.90 3.73
2.60 1.83 1.00 64.84 36.69 39.58 13.90 51.95 20.0 8.80 8.37 9.59 8.96 3.74
2.80 1.59 1.00 75.32 40.20 43.13 15.05 56.30 20.1 8.84 8.45 9.61 9.01 3.74
3.00 1.41 1.00 86.65 43.76 46.73 16.21 60.69 20.2 8.89 8.53 9.63 9.06 3.75
3.20 1.27 1.00 98.83 47.38 50.39 17.35 65.14 20.4 8.92 8.60 9.65 9.11 3.75
3.40 1.17 1.00 111.88 51.07 54.11 18.50 69.63 20.5 8.95 8.66 9.68 9.16 3.76
3.60 1.11 1.00 125.79 54.80 57.88 19.65 74.19 20.6 8.98 8.72 9.71 9.20 3.78
3.80 1.06 1.00 140.56 58.59 61.69 20.80 78.82 20.7 9.00 8.77 9.73 9.24 3.79
4.00 1.04 1.00 156.18 62.41 65.54 21.95 83.51 20.9 9.03 8.82 9.76 9.28 3.80