We establish boundary gradient estimates for a solution for selfsimilar transonic potential flow on the subsonic region. The gradient estimates enable us to verify the ellipticity of the solution in the region. We consider both convex and general domains. We provide the sufficient conditions to show that the flow becomes transonic.