A method based on the utilization of magnetic field and field gradient distributions for calculating the gradient of the magnetic surface function, ∇ψ, in stellarator systems is treated. The distributions of ∇ψ and associated quantities over the magnetic surfaces of torsatrons obtained by this method are discussed. The method can also be applied to other toroidal systems which confine plasma in an asymmetrical magnetic field.