'---------------------------------------------------------------'
'Input/Output and Initialization                                '
'---------------------------------------------------------------'
'Name of initialization and dump file..'   'basenm'     'VINEdmp'
'IO library: Ftn:0, XDR=1, HDF=2.......'   'io_fmt'      0
'Dump number for restart (neg=tempdump)'   'idump'       0
'Make checkpoint dumps if makechkpnt=1.'   'makechkpnt'  0
'Override dumpfile value of tnextdump..'   'override'    F
'Initial time step for all particles...'   'dtinit'      1d-3
'Maximum allowable time step ..........'   'dtmax'       .5d0 
'Minimum allowable time step ..........'   'dtmin'       0d0 
'Maximum time in code units to end ....'   'tstop'       200d0
'Max time steps allowed................'   'nstepmax'    100000000
'Code time between output dumps .......'   'dt_out'      5d0
'Code time between on the fly analysis.'   'dt_ana'      0.1d0
'Max time steps between restart dumps..'   'nout_max'    100000
'Max time steps between analysis output'   'nana_max'    100
'Debugging turn on if idbg>0...........'   'idbg'        1
'Print timing information if itme=1....'   'itme'        1
'---------------------------------------------------------------'
'I/O unit conversion settings (see also subr. setunits in io.F).'
'---------------------------------------------------------------'
'Units 0=scalefree 1=phys 2=cosmo 3=?..'   'iunits'      0
'Mass unit M=1vineunit=xmas*physmass...'   'xmas'        1d0
'Length unit D=1vineunit=ylen*physlen..'   'ylen'        10d0
'---------------------------------------------------------------'
'Integrator, boundary and timestep criteria settings            '
'---------------------------------------------------------------'
'Individual time steps.................'   'indts'       1
'Use timestepbins for global timestep..'   'iusebin'     1
'Extrap all parts to now, (indts only).'   'ifullextrap' 1
'Boundary: 0=free, 1=periodic, 2=refl..'   'ibctype'     0
'Random box shift for periodic bound...'   'ishift'      0
'Minimum x-coordinate..................'   'xmin'       -1d0
'Maximum x-coordinate..................'   'xmax'        1d0
'Minimum y-coordinate..................'   'ymin'       -1d0
'Maximum y-coordinate..................'   'ymax'        1d0
'Minimum z-coordinate..................'   'zmin'       -1d0
'Maximum z-coordinate..................'   'zmax'        1d0
'Error tolerance coeff. for RK integ...'   'tol'         1.d-4
'PM error tolerance coeff. for RK integ'   'tolpm'       1.d-6
'Frac of force tstep criterion   ......'   'ftol'        0.3d0
'Frac of veloc tstep criterion   ......'   'vtol'        0.3d0
'PM frac of force tstep criterion   ...'   'ftolpm'      0.003d0
'PM frac of veloc tstep criterion   ...'   'vtolpm'      0.003d0
'Fraction of CFL timestep.......... ...'   'cfl'         0.3d0
'Use smoothing length ts criterion.....'   'i_hts'       1
'Max fractional h change...............'   'dhmax'      .15d0
'---------------------------------------------------------------'
'Thermodynamic and SPH settings                                 '
'---------------------------------------------------------------'
'Minimum smoothing length .............'   'hmin'        1d-5
'Maximum smoothing length .............'   'hmax'        5d-1
'Minimum specific internal energy......'   'umin'        1d-8
'EOS: 0=ideal, 1=LagISO 2=EulISO.......'   'ieos'        0
'Gas pressure if ipres=1...............'   'ipres'       1
'Heating from PdV work  if ipdv =1.....'   'ipdv'        1
'Heating from viscosity if ishok=1.....'   'ishok'       1
'Artificial viscosity if ivisc=1.......'   'ivisc'       1
'Balsara viscosity if ibals=1..........'   'ibals'       1
'Time varying viscosity if ivtime=1....'   'ivtime'      0
'Artificial viscosity decay length (h).'   'vlength'     2d0
'Max AV constant value for alpha.......'   'alphamax'    1d0
'Artificial viscosity constant alpha...'   'alphastar'   1d0
'Artificial viscosity constant beta....'   'betastar'    2d0
'External cooling .....................'   'icool'       0
'Cooling time constant ................'   'taucool'     15d0
'External heating .....................'   'iheat'       0
'---------------------------------------------------------------'
'Gravity and tree settings                                      '
'---------------------------------------------------------------'
'Self gravity if igrav=1...............'   'igrav'       1
'Use GRAPE: igrape=1, else 0...........'   'igrape'      0
'Direct summation if 1, else tree......'   'idirect'     0
'Timesteps between tree builds (globts)'   'maxbuild'    30
'Full(1) or partial(0) tree revisions..'   'irevise'     1
'Tree rebuild: if h_clmp>clfac*h_clmp0.'   'clhfrac'     2.0d0
'MAC:.angle(1).abserrWS(2).gadget(3)...'   'imac'        3
'MAC: opening criterion (angle or err).'   'treeacc'     5e-3
'Softening:0=varkern,1=fixkern,2=Plumr.'   'isoft'       0
'Plummer eps or fixed h for SPH grav...'   'eps'         1d-1
'Max sep between clumps ...............'   'sepmax'      10d0
'Max sep between grape bunches.........'   'bsepmax'     10d0
'Minimum particles to send to GRAPE....'   'igrapelim'   50
'Max node mass for interaction list....'   'gmasslim'    1d100
'---------------------------------------------------------------'
'Physical model settings (non-cosmo)                            '
'---------------------------------------------------------------'
'Particle removal:0=dont,1=eat,2=bleed.'   'irem'        0
'Dust evap on for full vert. column....'   'idustev'     0
'Fraction of column evap before op-thin'   'colfrac'     .005 
'Turn on disk heating from star........'   'istrheat'    0
'Assume outer radius of disk is .......'   'rdisk'       50d0
'---------------------------------------------------------------'
'Cosmology parameter settings                                   '
'---------------------------------------------------------------'
'cosmological expansion if iexpand=1...'   'iexpand'     0
'Number of cosmo dumps.................'   'ncdumps'     11
'Redshifts for cosmo dumps.............'    20 15 12 10 8 4 3 2 1 0.5 0 
'Hubble parameter -- 0<h_0<1 -- .......'   'h_0'         0.65d0
'Omega_matter..........................'   'omegam'      .35d0

