changes added to workspace, flux1, flux1_2, reconstruct.f90 and solver to accomodate for the inclusion of polynomial interpolation in lim_type 4