Input desired special table:
   1=charm only
   2=cc coherent only
   3=nc coherent only
calculating for -n CHARM -v 1
 ERROR IN NEUGEN3:
 INITIALIZE_CONFIGURATION: UNKNOWN CONFIGURATION/VERSION
  NAME = CHARM      VER=  1
  RESET TO DEFAULT
  
  *****  CERN Computer Program Library     -     Reference:  W5051  *****
  *****  PDFLIB Version:  8.04  Released on  2000-04-17  at  12.24  *****
  
   ****  NEUGEN Configuration and Parameters **** 
  NEUGEN CVS $Name: v3_5_5 $  Built: Nov 14 2007 13:13:20

   Configuration :  MODBYRS   
         version :     5
         OK flag :  PASS
  
   Model Choices --> 
      dismodel    :  Bodek-Yang Modified LO
  res/dis overlap :  Exclusive channel reduction   
   Intranuke ver. :  RENUKE                        
    Hadronization :  AGKY-tuned KNO + JETSET       
     piabsorption :  Ransome model                 
   Binding Energy :  Constant amount subtracted    
   formation zone :  SKAT parametrization-OLD      
   Bodek-Yang ver :  Version 2                     
   Pauli blocking :  Fermi gas model               
  CC form factors :  BBBA-2006                     
   Carrot MC errs :  No                            
   Daikon MC errs :  No                            
  
   QEL Model Parameters --> 
            ma   :     0.9900
            mv   :     0.8400
         fa(0)   :    -1.2670
          mu_p   :     2.7930
          mu_n   :    -1.9130
           eta   :     0.1200
  
   Resonance Model Parameters --> 
         ma(res) :     1.1200
         omega   :     1.0500
         Z       :     0.7620
  
   DIS Model Parameters --> 
      PDF Group  :     5
      PDF Set    :    12
      Q2min      :     0.8000
    Scale factor :     1.0320
  
   Bodek-Yang Model Parameters --> 
             a   :     0.5380
             b   :     0.3050
         cv_1u   :     0.2910
         cv_2u   :     0.1890
           csu   :     0.3630
         cv_1d   :     0.2020
         cv_2d   :     0.2550
           csd   :     0.6210
      PDF Group  :     5
        PDF Set  :    12
  
   Overlap parameters --> 
       Wcut_dis  :       1.700
       Wcut_res  :     1.7000
   How Many res  :    17
  
   Coherent Model Parameters --> 
          f_pi   :     0.1307
            r0   :     1.0000
      re/im pi   :     0.3000
       ma(coh)   :     1.0000
  
   Hadronization Model Parameters -->
      jetset_low :        2.3
       jetset_hi :        3.0
         a(im=1) :     0.4000
         a(im=2) :    -0.2000
         a(im=3) :     0.0200
         a(im=4) :     0.8000
         b(im=1) :     1.4200
         b(im=2) :     1.4200
         b(im=3) :     1.2800
         b(im=4) :     0.9500
         c(im=1) :     7.9300
         c(im=2) :     5.2200
         c(im=3) :     5.2200
         c(im=4) :     7.9300
     s-sbar prob :     0.0500
  
   RES / DIS Overlap Parameters -->
     reduction for multiplicity 2,3,4 
      CC (im=1)  :     0.1000   1.0000   1.0000
      CC (im=2)  :     0.3000   1.0000   1.0000
      CC (im=3)  :     0.3000   1.0000   1.0000
      CC (im=4)  :     0.1000   1.0000   1.0000
      NC (im=1)  :     0.1000   1.0000   1.0000
      NC (im=2)  :     0.3000   1.0000   1.0000
      NC (im=3)  :     0.3000   1.0000   1.0000
      NC (im=4)  :     0.1000   1.0000   1.0000
  
    Intranuke Parameters --> 
   Nuclear density :     3.4000
    Nuclear radius :     1.4000
    deBroglie ring :     0.5000
      Minimum step :     0.0000
   Pi abs scale fac:     1.0000
   Ransome fraction:     1.0000
   Binding energy  :     0.0080
  
   Formation Zone Parameters --> 
    formation time :     0.3420
      b (NOMAD par):     0.0000
 i,e =  1  0.100000001
 i,e =  2  0.126896098
 i,e =  3  0.161026224
 i,e =  4  0.204336002
 i,e =  5  0.25929442
 i,e =  6  0.329034537
 i,e =  7  0.417531997
 i,e =  8  0.529831827
 i,e =  9  0.672335923
 i,e =  10  0.85316807
 i,e =  11  1.08263707
 i,e =  12  1.37382424
 i,e =  13  1.74332929
 i,e =  14  2.21221685
 i,e =  15  2.80721712
 i,e =  16  3.56224942
 i,e =  17  4.5203557
 i,e =  18  5.73615551
 i,e =  19  7.27895832
 i,e =  20  9.23671436
 i,e =  21  11.7210312
 i,e =  22  14.8735323
 i,e =  23  18.8739338
 i,e =  24  23.9502869
 i,e =  25  30.391983
 i,e =  26  38.5662422
 i,e =  27  48.9390602
 i,e =  28  62.1017494
 i,e =  29  78.8046875
 i,e =  30  100.
 Calculating for iflav, int_type, im = 1 1 1
 Calculating for iflav, int_type, im = 1 1 2
 Calculating for iflav, int_type, im = 1 1 3
 Calculating for iflav, int_type, im = 1 1 4
 Calculating for iflav, int_type, im = 1 2 1
 Calculating for iflav, int_type, im = 1 2 2
 Calculating for iflav, int_type, im = 1 2 3
 Calculating for iflav, int_type, im = 1 2 4
 Calculating for iflav, int_type, im = 1 3 1
 Calculating for iflav, int_type, im = 1 3 2
 Calculating for iflav, int_type, im = 1 3 3
 Calculating for iflav, int_type, im = 1 3 4
 Calculating for iflav, int_type, im = 1 4 1
 Calculating for iflav, int_type, im = 1 4 2
 Calculating for iflav, int_type, im = 1 4 3
 Calculating for iflav, int_type, im = 1 4 4
 Calculating for iflav, int_type, im = 1 5 1
 Calculating for iflav, int_type, im = 1 5 2
 Calculating for iflav, int_type, im = 1 5 3
 Calculating for iflav, int_type, im = 1 5 4
 Calculating for iflav, int_type, im = 2 1 1
 Calculating for iflav, int_type, im = 2 1 2
 Calculating for iflav, int_type, im = 2 1 3
 Calculating for iflav, int_type, im = 2 1 4
 Calculating for iflav, int_type, im = 2 2 1
 Calculating for iflav, int_type, im = 2 2 2
 Calculating for iflav, int_type, im = 2 2 3
 Calculating for iflav, int_type, im = 2 2 4
 Calculating for iflav, int_type, im = 2 3 1
  
  
  Nucleon PDFs :  GRV98 Set LO Structure Functions
                  Ngroup =  5,   Nset =  12
  --------------------------------------------------------------------------------------
 Calculating for iflav, int_type, im = 2 3 2
 Calculating for iflav, int_type, im = 2 3 3
 Calculating for iflav, int_type, im = 2 3 4
 Calculating for iflav, int_type, im = 2 4 1
 Calculating for iflav, int_type, im = 2 4 2
 Calculating for iflav, int_type, im = 2 4 3
 Calculating for iflav, int_type, im = 2 4 4
 Calculating for iflav, int_type, im = 2 5 1
 Calculating for iflav, int_type, im = 2 5 2
 Calculating for iflav, int_type, im = 2 5 3
 Calculating for iflav, int_type, im = 2 5 4
 Calculating for iflav, int_type, im = 3 1 1
 Calculating for iflav, int_type, im = 3 1 2
 Calculating for iflav, int_type, im = 3 1 3
 Calculating for iflav, int_type, im = 3 1 4
 Calculating for iflav, int_type, im = 3 2 1
 Calculating for iflav, int_type, im = 3 2 2
 Calculating for iflav, int_type, im = 3 2 3
 Calculating for iflav, int_type, im = 3 2 4
 Calculating for iflav, int_type, im = 3 3 1
 Calculating for iflav, int_type, im = 3 3 2
 Calculating for iflav, int_type, im = 3 3 3
 Calculating for iflav, int_type, im = 3 3 4
 Calculating for iflav, int_type, im = 3 4 1
 Calculating for iflav, int_type, im = 3 4 2
 Calculating for iflav, int_type, im = 3 4 3
 Calculating for iflav, int_type, im = 3 4 4
 Calculating for iflav, int_type, im = 3 5 1
 Calculating for iflav, int_type, im = 3 5 2
 Calculating for iflav, int_type, im = 3 5 3
 Calculating for iflav, int_type, im = 3 5 4
 Calculating for iflav, int_type, im = 4 1 1
 Calculating for iflav, int_type, im = 4 1 2
 Calculating for iflav, int_type, im = 4 1 3
 Calculating for iflav, int_type, im = 4 1 4
 Calculating for iflav, int_type, im = 4 2 1
 Calculating for iflav, int_type, im = 4 2 2
 Calculating for iflav, int_type, im = 4 2 3
 Calculating for iflav, int_type, im = 4 2 4
 Calculating for iflav, int_type, im = 4 3 1
 Calculating for iflav, int_type, im = 4 3 2
 Calculating for iflav, int_type, im = 4 3 3
 Calculating for iflav, int_type, im = 4 3 4
 Calculating for iflav, int_type, im = 4 4 1
 Calculating for iflav, int_type, im = 4 4 2
 Calculating for iflav, int_type, im = 4 4 3
 Calculating for iflav, int_type, im = 4 4 4
 Calculating for iflav, int_type, im = 4 5 1
 Calculating for iflav, int_type, im = 4 5 2
 Calculating for iflav, int_type, im = 4 5 3
 Calculating for iflav, int_type, im = 4 5 4
 attempt to write: CHARM1_xsec.dat
 write successful!