Proposal and analysis of a new numerical method for the treatment of hydrodynamical shock problems