In Hele-Shaw flows at vanishing surface tension, the boundary of a viscous fluid develops cusp-like singularities. In recent papers [1,2] we have showed that singularities trigger viscous shocks propagating through the viscous fluid. Here we show that the weak solution of the Hele-Shaw problem describing viscous shocks is equivalent to a semiclassical approximation of a special real solution of the Painlevé I equation. We argue that the Painlevé I equation provides an integrable deformation of the Hele-Shaw problem which describes flow passing through singularities. In this interpretation shocks appear as Stokes level-lines of the Painlevé linear problem.Harmonic measure is the probability of a Brownian excursion starting at a marked point (in this case it is a drain) to exit the domain through that boundary element. The probability density of Brownian excursion is the Poisson excursion kernel, equal to the gradient |∇p| of a harmonic function p with a source at a marked point and Dirichlet boundary condition.A well-known stochastic realization of the Darcy law is DLA -diffusion-limited aggregation [4]. It is realized through Brownian excursions of particles with a small size emanating at a constant rate from a distant point, with absorbing boundary, such that the stopping-set cluster grows. The Darcy law emerges at vanishing particle size, → 0.Computer experiments with DLA [4] show that a boundary, initially featureless, very quickly develops into a branching graph with a width controlled by the size of one particle . This signals that the limit → 0 is impossible.Darcy's law (1) also indicates that Hele-Shaw flow tends to a microscale. Non-linear differential equation (1) is ill-defined. A smooth initial boundary first evolves into a fingering pattern, then fingers at a finite, critical time develop cusp singularities [5,6,7,8]. At that point the differential form of the Darcy law stops making sense. This phenomenon constitutes the major problem in the field.Unlike in fluid dynamics, DLA processes possess a dimensional parameter -the particle size . It plays the role of a minimal area, which regularizes singularities emerging in fluid dynamics.Comparing Hele-Shaw flows to DLA suggests that, after a singularity is reached, the flow should feature shocks -a graph of curved lines where pressure (and velocity) suffer discontinuities. Such solutions of ill-defined non-linear differential equations are known as weak solutions [9].In Refrs. [1,2] we developed the weak solution of Hele-Shaw problem, which we believe to be a solution of the problem of Hele-Shaw singularities. Within this solution, a singularity triggers a branching graph of shocks or cracks propagating through the viscous fluid. Shocks at small Reynolds number is a new phenomena, not yet being observed experimentally. We refer them as viscous shocks. An emerging pattern of viscous shocks is reminiscent DLA patterns.