Buoyancy

Now that the density of the cell can be calculated at every point (assuming constant temperature), I am now able to calculate  the following:
1) The acceleration due to gravity at every point in the Sun.
2) The force of buoyancy acting upon the cell at every point along the Radius.
3) The net acceleration of the cell at every point.

The code used for this is "buoyancy.py". It is similar to "bubble.py" in that it asks for the mass, initial position and temperature of the cell. For a cell of 50E6 kg and 3E6 K, graphs like the following are given:

The second graph is independent of the initial criteria given for the cell, as it is the gravity due to the Sun itself.

The next step is to calculate the motion of the cell. Since I now have acceleration for the cell at every point (up to 0.95 Solar Radii), I can calculate the velocity and position of the cell via integration, provided I fit a function to the acceleration. There exists a command in python which can achieve this - "curve_fit".