added separate field for c, (Do not put c in second space in k[2], that will lead to tons and tons of problems!!!) and modified the shift-dv and psi_to_v functions accordingly.