OMIT.f

OMIT Manual Page

OMIT is an auxillary program to modify TYCHO7 models, primarily to omit from the grid those regions in which little happens, but have short sound travel times and hence are computationally expensive for the explicit hydrodynamic mode (mode = 0). The control file is "omit.in", listed below, should reside locally in the directory in which OMIT is executed.

omit.in

'.....................................................................'
' OMIT: for TYCHO7                                    '
'.....................................................................'
' Display device: /ps, /xwin ................' 'device'   '/xwin'
' Line style.0=solid line, else dots.........' 'linesty'   1
' nopac (-1=thomson,0=opal2,1=opal1).........' 'nopac'     0
' scale opacity (0=no, else yes).............' 'nkscale'   0
' scale thoul diffusion (0=no,else = factor).' 'fthoul'    0.0d0
' (-1=complete ionization,else=compute)......' 'noion'     0
' (0=radiation pressure,else=no).............' 'norad'     0
' (0=H2 molecule,else=no)....................' 'nomole'    0
' (0=H2 coulomb plasma,else=no)..............' 'nocoul'    0
' non-nuclear neutrino flag(-1=no,0=on)......' 'neutro'    0
' flag to force OPAL EOS (0=no)..............' 'nopaleos'  1
' Mixing (0=PAY mix+diff;1=mix;2=diff;3=none)' 'mixmode'   3
' Filename of model..........................' 'cmod'     'imodel'
' Variable flag for x axis..(see genvec.f)...' 'cvarx'     'xm'
' Variable flag for y axis..(see genvec.f....' 'cvary'     'cour'
' Override value for xmin....................' 'xminf'     0.d0
' Override value for xmax (no if = xmin).....' 'xmaxf'     0.0d0
' Override value for ymin....................' 'yminf'     0.0d0
' Override value for ymax (no if = ymin).....' 'ymaxf'     0.0d0
' Keep inside (>=0) or outside (<0)..........' 'inside'    -1
' Choose slice point in solar masses.........' 'slice'     0.42d0
' Explosion energy to be added(ergs).........' 'explode'   5.0d47
' Mass (at bottom) to put it in(grams).......' 'massex'    1.0d32
'.....................................................................'

Meaning of new parameters:

  • "inside<0": a negative integer will cause OMIT to slice of the inside (keeping mantle),
  • "inside = 0 or >0": Zero or a positive integer will cause OMIT to slice off the outside (keeping core).
  • "slice": the mass coordinate (in solar mass units) for the slice boundary.
  • "explode": if this real*8 number is positive, it is the kinetic energy (in ergs) to be added just above "slice". This should imply "inside" < 0. A negative or zero value of "explode" is a signal to do nothing.
  • "massex": a real*8 number denoting the mass (in grams) over which "explode" is spread.

    • OMIT will graph the model named in "cmod", with a vertical gray line showing the slice coordinate, and the usual parameters control the graph.

    Variables which are also in "params.d" need to be set in order to define the EOS, opacities, and mixing mode to use.

    Upon success, OMIT will produce a model file (which can be copied into "imodel") named "newomit".

    Return to main TYCHO page.

    E-mail:darnett@as.arizona.edu

    TYCHO is not public domain software, but is freely available for non-commercial purposes. Absolutely no guarantee is given; use at your own risk...

    Last updated 2-20-02