Added test of convective stability. Changed datafolder to 'data2' to separate from oblig1.
Renamed program to better name.
New inital parameters set to reflect entire sun.
Changed to dynamic dm and thus new loop type.
Convection implemente. Model is in theory complete, but something is of course wrong.
Simplified code by removing length-2 arrays for R,P,L,T. By black magic/voodoo this removed divide by zero error. It does also cause Forward Euler to be used instead of Euler Chromer.
Added functions for finding fluxes.
Fixed bug where two variables were both named R.
Moved printing to function for more elegance. Fixed bug with flux functions having wrong indentation.
Calculate and print fluxes.
Avoid saving last point twice.Extra linebreak between printouts.
Dynamic mass step update should happen AFTER parameter updates, or else a modified and wrong dm is used.
Dynamic resolution.
Made dm smarter and stricter when handling too large differences.
Set F_con=0 if not convective unstable. This shows that there is almost no convection.
Used safer method to find the real root of xi. Langtangen taught us not to compare floats with '=='.
Added fluxplotting. Plots in 3 figures.
Added update of fluxes when simulation stops.
Fixed logical error which picked wrong root for xi. Now it is just like before.
Fixed algebraic typo.