initialization function and grid boundaries are encapsulated in a single parameter.
Length of grid specified by initialization array if this is what is passed as input to the code.
Write to log file in out directory
Ability to calculate a diffusion-like derivative.
Trying entirely different artificial viscosity scheme.
Can output at basically equal increments of time. Fixed recently introduced bug where artificial viscosity was not multiplied by the Laplacian of the velocity.
script used in trying to adjust parameters to solve for NGC4551 profile.
power law extrapolations for the boundary conditions
Another adjustment script for ngc4551.
Script to compute profile for ngc4551 from scratch.