In this work, some numerical solutions of magnetohydrodynamics equations are investigated in the presence of differential rotation with the use of previously developed algorithm. This algorithm includes the thin shell approximation and a special separation of variables which were used to obtain the radial and latitudinal variations of physical parameters in spherical coordinates. The magnetic field profile is chosen to produce comparable magnetic fluxes found in previous works. The sphericity and density shape parameters relevant to model is determined by using two different known differential rotation profiles. It is found that the shape of variations in physical parameters is strongly dependent to magnetic field profile and there is a considerable change in density with respect to reference model. It is as well shown that the spherical symmetric distributions of physical parameters are broken for the region of study. (C) 2008 Elsevier B.V. All rights reserved.