!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) 
>