This article deals with the degenerate parabolic system with nonlinear bound- ary flux. By constructing the self-similar supersolution and subsolution, we obtain the critical global existence curve and the critical Fujita curve for the problem. Especially for the blow-up case, it is rather technical. It comes from the construction of the so-called Zel'dovich-Kompaneetz-Barenblatt profile.