The magnetic fields of the giant planets, Jupiter and Saturn, deviate significantly from a pure magnetic dipole and the cold plasma is mostly centrifugally confined near the equator. The additional contribution of the azimuthal currents leads to the stretching of the magnetic field and the formation of a characteristic, disc-type structure known as a magnetodisc. We present here an updated version of a numerical implementation of Caudal's iterative scheme, used to create models of the magnetosphere. In particular, we include newer equatorial density, temperature and hot plasma profiles obtained from Galileo data. Finally, we describe and use an algorithm to update the angular velocity profile after the end of the iterative process, using information from the magnetodisc. We also present comparisons between the azimuthal current and plasma flow predicted by our model and those derived from spacecraft observations.