We present a numerical model used to compute steady flow solutions of the venous circulation of the leg. The network topology is based on clinical data and the flow is assumed to be steady, incompressible, and one-dimensional. We develop a non Newtonian approach to a one-dimensional flow because the blood viscosity depends on the velocity profile, and we demonstrate theoretically the pertinence of a phenomenological law of equivalent viscosity. Clinical experiments observe hemodynamical variables (i.e. venous pressure, venous area, blood velocity) only at the accessible places. In contrast the numerical model results are not limited to particular locations but can be evaluated on every point of the network. It provides important help to the definition of a clinical protocol. The model was designed to quantify a compression level of elastic compression stockings and to plan clinical studies. We validate the numerical approach using a published clinical trial, where the diameter of superficial and deep veins were measured at different compression pressures. We show also that the viscosity variations in a bed-rest position as a consequence of the application of a European Class II compression stockings. These variations could prevent the hyper-coagulability and the stasis of the blood.