corrected poly.f90 to begin loop in num_ghost. solver.py now calculates the scale time step based on self.dt