TYCHO Users Manual

TYCHO-6.9: Auxillary programs for data analysis

TYCHO has a variety of auxillary programs to analyze the data it produces. There are three types of data:
  • model "dumps", with file names in the form of a prefix of 2 characters and a suffix of 5 numbers, and which are restartable by simply copying the dump file to imodel) and represent a snapshot of the conditions in the star at that time. TYCHO does not overwrite dumps. Most analysis programs use these data.
  • HR time lines, which give the surface conditions and abundances at every time step since the beginning of the computation of the sequence. These are analyzed with hrplot which gives a conventional Hertzsprung-Russell diagram, and hrt which gives surface variables versus time. TYCHO will overwrite these files unless you restart from the highest model number (last model) in a sequence. Change your prefix sequence identifiers as needed to avoid this overwriting.
  • CV time lines, which give the convective structure and the most active burning zones as a function of time. TYCHO will overwrite these files unless you restart from the highest model number (last model) in a sequence. Change your prefix sequence identifiers as needed to avoid this overwriting. These data are analyzed with cvplot.

    Each of the analysis programs has a input data file (such as genplot.in) which is read by the program initially, and allows the user control over some aspects of the display.

    Genplot

    genplot is one of the most versatile programs for examining model dumps. It can plot a variety of variables, pairwise; see /tycho-6.9/src/getvec.f for a listing (you can hack this subroutine to define additional variables). The figure shown is a graph of convective speed (as estimated by mixing length theory) versus mass coordinate for model m600702 (6 solar masses). The source code is in /tycho-6.9/src/genplot.f; see also /tycho-6.9/genplot.in.

    Genrate

    genrate shows the most active reaction links between nuclei, by breadth of line and by color, in a proton number Z and neutron number N space. For the graph shown (the inner zone K=2 of m600702), the CNO cycle dominates, but the pp chain contributes, and weak NeNaMg cycling occurs. Dotted lines represent weak interactions (beta decay, positron emission, positron decay, beta emission). Some detailed numerical data are also produced; the souce code is in /tycho-6.9/src/genrate.f. See also /tycho-6.9/genrate.in.

    Gennuc

    With 177 nuclei in the standard network, TYCHO produces a lot of information concerning abundances. gennuc plots abundance by mass versus mass coordinate (or radius, log radius, or zone index k), suppressing the representation of those abundances that do not change through the star (for clarity).Clear labelling is a problem for all but the simplest networks. The source code is in /tycho-6.9/src/gennuc.f, see also /tycho-6.9/gennuc.in .

    Hrplot

    hrplot shows the history of the star in the HR (log luminosity-log effective temperature) diagram, and can plot several sequences simultaneously. Some data for well-observed stars is ploted also. Sequences m6, x2, m2, n1, m1, p7 are shown. The source code is in /tycho-6.9/src/hrplot.f; see also /tycho-6.9/hrplot.in

    Cvplot

    cvplot shows the history of convection and nuclear burning for a given star (one sequence). The convective region is shown in red; the line of most vigorous burning in yellow. The horizontal axis shown is step number; log time is also available. The m6 sequence is shown. Originally the entire star is convective (on the Hayashi track). Deuterium burning proceeds all the way to the surface. CNO burning creates a convective core, which is still active at the last plotted model. The source code is in /tycho-6.9/src/cvplot.f; see also /tycho-6.9/cvplot.in.

    E-mail:darnett@as.arizona.edu hi

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