parker wind working?!
Reynolds number formalism for artificial viscosity and optional parameters for function defining initial conditions
Fixed upadating of auxillary quantities like pressures for ghost zones
Exporting separate movie for each of the quantities of interest (i.e. density, velocity, temperature, mass flux)
Separate sub-step method grid is updated after each substep
Adding mass source term
Calculate and store all necessary time derivs at the beginning of each substep
Changing variables to log_rho
Boundaries are updated by power law extapolations. Artificial viscosity is now set by the min(abs(v), cs) instead of just min(v, cs).
Uncommented code for updating outer boundary. Saved time slices stored in numpy array without time stamps.
Getting rid of large useless blocks of commented code
Improved names for out files
Beginnings of log(r) grid, and modified the cfl function.
Array of deltas for each grid pt
log(r) grid attempt 1
Write solution to file. Saved array included time-stamps and radii
Trying to get outer wind to go supersonic
Tracking driver with analytic solutions for Bondi and Parker.
Decreased Reynold number and doubled outer grid radius in parker.py
Minor changes in grid setup (outer radius, number of grid points)