Method of characteristics for pressure transients in cylindrical geometry with axial symmetry: the formulation