Changed the flag for constant gamma from 4 to -1, to reflect the actual select case statement in the code.
Changed the indirect potential to use the position as calculated by wiggle.F (which evolves the position based on gravitational force), rather than usiing indirect routine which assumes position of star is such that center of mass of the system remains fixed.
Changed x and y from reals to double precision.
Commented out line that set OMMAX=zero,  added call to RITE subroutine to print initial data, added if statement so that massflow.xxxxxx writes every 10 steps instead of every iteration, commented out write statements for temperature file, uncomment if you're doing a run with radiation cooling, added flush statement for tcoef file, added if statement so that subroutine CENTMASS is called every IDIAG steps instead of every step.
removed if statements that were included for debugging purposes.
Changed the velocity limiter from two*sound to six*SOUND, because the angular momentum was being clipped on some runs.
Changed readin format of fort.2,added calculation of tmass and total angular momentum for diagnostics, changed calculation of sound speed to use maximum density in disk instead of denny(2,2), changed the normalization of the system, included a place to hardwire mass_star, this will need to be read in eventually, added conditional such that coolhead_full.xxxxxx is not written for non cooling runs, commented out gamma1.xxxxxx file, we don't use it, changed format statements so that mfrac.dat is written with more decimal places.
Changed ppred to reflect the changes made to state.f regarding form of the energy equation.