Parametrization of Riemann-measurable selections for multifunctions of two variables with application to differential inclusions