nuker solver attempt 1
Fixing conervation checks and some sign issues.
Packaged at lot of the nuker functions into galaxy class.
Tying to convert nuker to cgs
Circularized surface density profile in Nuker consistent units for radius/density
Commenting out converged stop condition
Parameter which specifies new units from the q source function.
Non-isothermal evolution does not have to be preceded by isothermal evolution.
New method to allow resetting wind velocity in the middle of a calculation.
Velocity dispersion term now defined in terms of Mtot.
Incorporating progress bar.
Solve and gradually adjust a parameter in the process.
Changing the target time after each parameter adjustment.
Fixed typos in solve_adjust method.
Got rid of vel dispersion in vw. Revert method.
Dummy arrays in cell class which track the overall vw/phi arrays of the solver. Specific heating rate is now calculated from the grid level vw array.
going back to using specific heating rate within each grid cell, rather than grid level method to calculate heating rate.
Artificial viscosity in entropy evolution equation
Adding ipython notebook from which solver run to git repository.
More detailed output for saved array. Also write to npz file. Can specifiy index of saved array to go back to in the revert method.