!go go_kl_sd_feed -------------------------------------- program kl_sd running version: $id: kl_sd.f,v 1.5 1998/05/13 15:24:08 tl exp $ *** selectible level density version *** -------------------------------------- calculation id: z,a of nucleus............: >>> 192-hg(80) decay out of sd only (y/n)[n]: modulus # for normal casc writeout: # of cascades.............: cascade cut off u energies: common entry distribution for normal and sd> entry dist hs2 file name: __ spectrum name........: # times to smooth the input data: specific entry and exit spins for normal and sd> (for decay width precalculation, kl_e?_gen entry spin: exit spin: decay parameters for normal and sd>> paring energy...........: level density parameters: level dens type.........: min level dens for stop.: quadrupole moments @ entry: quadrupole moments @ exit: rotational damping width: rotational damping cut off: n gdr strenght spectrum: __file contains 301 data points spectrum..: gdr_n.top _contains.: gdr strength fun n _id.......: sd gdr strenght spectrum: __file contains 301 data points spectrum..: gdr_s.top _contains.: gdr strength fun sd _id.......: normal yrast line data file: data from: hg192_nd_yrast.dat # data points in file...: 54 modulus for data logging: 6 0. 0. 8.00000 2.21600 18.0000 4.38772 24.0000 6.38096 31.0000 9.91800 43.0000 15.8783 55.0000 23.7458 67.0000 33.5206 79.0000 45.2028 complete array constructed __ data smoothed 5 times spectrum..: kl_n_yrast.top _contains.: normal yrast _id.......: sd yrast line data file: data from: hg192_sd_yrast.dat # data points in file...: 46 modulus for data logging: 5 0. 5.37216 10.00000 6.00000 20.0000 7.69818 30.0000 10.35598 40.0000 13.8519 50.0000 18.0788 60.0000 22.9491 70.0000 28.3927 80.0000 34.3530 90.0000 40.7839 complete array constructed __ data smoothed 0 times spectrum..: kl_s_yrast.top _contains.: sd yrast _id.......: spectrum..: kl_yrast.top _contains.: yrast line _id.......: nd-j1 moment of inertia file: # times to smooth data: data from: hg192_j1_nd.dat # data points in file...: 4 modulus for data logging: 1 0. 60.9800 33.0000 60.9800 36.0000 75.5000 90.0000 75.5000 complete array constructed __ data smoothed 3 times spectrum..: j1_nd.top _contains.: j1_nd _id.......: sd-j1 moment of inertia file: # times to smooth data: data from: hg192_j1_sd.dat # data points in file...: 2 modulus for data logging: 1 0. 120.830 90.0000 120.830 complete array constructed __ data smoothed 3 times spectrum..: j1_sd.top _contains.: j1_sd _id.......: spectrum..: eg_nd.top _contains.: eg_nd _id.......: spectrum..: eg_sd.top _contains.: eg_sd _id.......: rigid body mois for rho (nd,sd): e2 mean energy deviation (%): spectrum..: kl_vn.top _contains.: vn _id.......: spectrum..: kl_vs.top _contains.: vs _id.......: sd barrier hight data file: data from: hg192_bsd.dat # data points in file...: 2 modulus for data logging: 1 0. 1.00000 80.0000 7.50000 complete array constructed __ data smoothed 1 times minimum barrier hight: spectrum..: kl_bsd.top _contains.: (mod)sd barrier _id.......: sd barrier hbar*omega file: data from: homg.dat # data points in file...: 2 modulus for data logging: 1 1.00000 0.600000 80.0000 0.600000 complete array constructed __ data smoothed 1 times spectrum..: kl_homega.top _contains.: hbar*omega _id.......: spectrum..: kl_yrast_comb.top _contains.: combined yrast line _id.......: -- precalcultion of e1&e2 decay widths>> u range for e1 calculation..........: egamma mean range for e2 calculation: mean spins for precalc (n,sd): 9 9 direct calculation precalculating all decay widths (8-total) __done precalculating ** kl_cm2_nout/initialization start version 92/9/24 cm2 data file: /home/tl/KEEP/cm2_1998.dat data file calc id: dummy ** loading data base ** 0.00100 100 0.00126 100 0.00158 100 0.00200 100 0.00251 100 0.00316 100 0.00398 100 0.00501 100 0.00631 100 0.00794 100 0.01000 100 0.01259 100 0.01585 100 0.01995 100 0.02512 100 0.03162 100 0.03981 100 0.05012 100 0.06310 100 0.07943 100 0.10000 100 0.12589 100 0.15849 100 0.19953 100 0.25119 100 0.31623 100 0.39811 100 0.50119 100 0.63096 100 0.79433 100 1.00000 100 1.25892 100 1.58489 100 1.99526 100 2.51188 100 3.16227 100 3.98107 100 5.01187 100 6.30957 100 7.94327 100 9.99999 100 12.58924 100 15.84891 100 19.95260 100 25.11884 100 31.62274 100 39.81068 100 50.11868 100 63.09568 100 79.43275 100 99.99992 100 125.89244 100 # data points read: 51 constructed c3 array if gamma/d less than .10000E-02 using two level if gamma/d greater than 20.000 using wigner otherwise the goe calculation is used ** kl_cm2_nout/initialization done output cascade file name: cascade file: do.tmp __ file marked with id: -------------------------------- setup done, start of monte carlo -------------------------------- kl_entry_exp: setting up total_entry.mascii open... input data smoothed 1 times (smoothed) matrix average spin..: 32.562 (smoothed) matrix average energy: 16.384 kl_entry_exp: setup done processed 10.0 % of the cascades processed 20.0 % of the cascades processed 30.0 % of the cascades processed 40.0 % of the cascades processed 50.0 % of the cascades processed 60.0 % of the cascades processed 70.0 % of the cascades processed 80.0 % of the cascades processed 90.0 % of the cascades processed 100.0 % of the cascades n_raw: 1000 n_fin: 1000 mean raw spin of entry point: 32.695 mean raw energy of entry point: 16.472 mean final spin of entry point: 32.695 mean final energy of entry point: 16.472 **** kl_sd done **** Note: IEEE floating-point exception flags raised: Inexact; Underflow; See the Numerical Computation Guide, ieee_flags(3M) >