To make this figure, retrieve the data:
To obtain the cell central radius for the zone of interest, compute the average of the cell outer radius (column 4 of mass_below_cell_mass_cell_outer_radius_cell_density.txt) and cell inner radius, which is also the cell outer radius for the zone directly below the zone of interest and thus easily acquired from column 4 of mass_below_cell_mass_cell_outer_radius_cell_density.txt as well. Such an average can be computed via an IDL routine, Microsoft Excel, or other computation software. Because the inner radius for the first zone is unavailable, the central radius for the first zone can be computed by combining the following two equations to eliminate Rinner and solve for Rcentral in terms of the cell mass, Router, and ρ:
![]() | (1) |
![]() | (2) |
![]() | (3) |
The cell mass, cell outer radius, and cell density for the first zone are found in columns 3, 4, and 5, respectively, of the first row of mass_below_cell_mass_cell_outer_radius_cell_density.txt
Save the values of the cell central radius for each zone in an appropriately-named text file and then graph the proper combination of the cell central radius and column 5 of mass_below_cell_mass_cell_outer_radius_cell_density.txt versus column 2 of mass_below_cell_mass_cell_outer_radius_cell_density.txt (divided by 1.9891 × 1033 to convert grams to solar masses).