 199  40  10  15  15   3   1   1   1   1   1   1
   1   5   1   5   0   0   0   0   0   0   0   0
 100   1   1   1   01000   0
 1.0E-06 1.0E-02 1.0E-07 0.0E+00 0.5E+00
  6  7  0  3  0 80  0  0  0 99
  1  2  4  5  3  9 10 11 12 15  8  7  6  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
  7  8  9 10 11 12 14  4  2  1  3  5  6  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
  4  5  6  7  8  9 10  2  3  1  2  3  1  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
  0  0  0  0  0154  0  0  0 99
  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
 17  2  4  5 18  8  9 10 11 12 13 14 24 25 26 30 16 34
  2  3  4  5  6  8  9 10 17 18 19 20 21 28  7
  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
 0.80 1.05 9.99 0.00 0.05 0.50 0.15 0.02 0.45 1.0E-04 1.0E+15 3.0E+19
 2.00E-02 2.000 0.700 0.173 0.053 0.482 0.099 0.038 0.080 0.072
 1.00E+08 0.00E+00 0.00E-01 0.00E-07 0.00E+00 3.00E-01 1.00E-03
 1.00E+03 0.00E+00 0.00E+00 1.00E+00 1.00E-04 0.00E+00 0.00E+00
 1 1 0.00E+02 1 1.00E+01 0.50 0.00
 0 700.0 0   0.0 0.00 0 0 1.00E-02 0.00E+00

 Above is for normal evolution with mesh spaced according to EQU(6)
 and central BC EQU(1)

Mesh spacing coefficients for normal evolution
 0.80 1.20 9.99 0.00 0.05 0.50 0.15 0.02 0.45 1.0E-04 1.0E+15 3.0E+19
Mesh spacing coefficients for AGB evolution
 0.80 1.20 3.00 3.00 0.05 0.01 0.60 0.35 0.45 1.0E-04 1.0E-01 3.0E-01

Data block for binaries with all the trimmings
 14 14  0  9  1102  0  0  0 99
  1  2  4 16 17 19 13 14 29  5  3  9 10 11 12 15 20 18 24 25 26 27 30  8  7  6 23 22 21  0
  7  8  9 10 11 12 14 22 23 24 25 26 27 29  4  2  1  3  5  6 19 17 16 13 28 18 20 21  0  0
  4  5  6  7  8  9 10 19 20 21 22 23 24 25  2  3  1 17 18 16  2  3  1  4  5 17 18 16 20  0

Data block for single stars -- with no orbital equations being solved
  6  7  0  3  0 76  0  0  0 99
  1  2  4  5  3  9 10 11 12 15  8  7  6  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
  7  8  9 10 11 12 14  4  2  1  3  5  6  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
  4  5  6  7  8  9 10  2  3  1  2  3  1  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0

--------------------------------------------------------------------------------
The following template identifies the numbers used in the above input file:
 NM2 IT1 IT2 JIN OUT NCH  JP ITH  IX  IY  IZ IMODE
 ICL ION IAM IOPINUC IBC ICN IML1 IML2 ISGTH IMO IDIFF
 NWRT1 NWRT2 NWRT3 NWRT4 NWRT5 NSAV  NMONT
  EPS     DEL     DH0     DT3     DD 
 NE1 NE2 NE3 NB NEV NF J1 J2 IH JH
  ID(11-110) - 7 lines
  ISX(1-45)  - 3 lines 
  DT1  DT2  CT(1) ............................................ CT(10) 
  ZS      ALPHA  CH    CC    CN    CO    CNE   CMG   CSI   CFE 
  RCD       OS      RML      RMG     ECA       XF       DR
  RMT      RHL       AC      AK1     AK2      ECT      TRB
  IRAM IRS1 VROT1 IRS2 VROT2 FMAC FAM
  IVMC TRC1 IVMS TRC2 MWTS IAGB ISGFAC FACSGMIN SGTHFAC
--------------------------------------------------------------------------------
NWRT1=100    -  print internal details at every 100th timestep only
NWRT2=2      -  print these details at every 2nd meshpoint only
NWRT3=1,2,3  -  print 1, 2 or 3 `pages' of these details
NWRT4=4      -  print a short summary of every 4th model only
NWRT5=5      -  print convergence details only after 5 iterations at each timestep
Each `page' has a selection of 15 output columns (first page) or 11 (second 
or third page; more decimal places), chosen from 31 variables computed in PRINTB.F.
These are:
             1  psi   2  P     3  rho  4  T      5  kappa  6  grada  7  grad 
 8  gr-ga    9  m    10  H1   11  He4 12  C12   13  N14   14  O16   15  Ne20
16  Mg24    17  r    18  L    19  Eth 20  Enuc  21  Eneu  22  dm    23  k**2
24  n/n+1   25  Uhom 26  Vhom 27  U   28  S     29  L/Edd 30  w.l   31  unused at present 
            33  Ne0  34  WT
k**2 is the square of the radius of gyration; 24, 25 and 26 are homology 
invariants, dlog rho/dlog P, dlog r/dlog P and dlog m/dlog P ; L/Edd is the 
ratio of local luminosity to local Eddington luminosity; w.l is convective 
velocity times mixing-length.  The 3 lines of ISX above give the selections
to be printed.
--------------------------------------------------------------------------------
