This article presents a magnetohydrodynamic two-dimensional numerical model of diamagnetic bubble equilibria in an axisymmetric open trap. The theoretical model consists of the Grad-Shafranov equilibrium equation and the transport equation obtained within the resistive single-fluid magnetohydrodynamics with isotropic pressure. Found are the numerical solutions corresponding to the diamagnetic confinement mode. In particular, the equilibria of the diamagnetic bubble in the GDMT are calculated. We investigate the effect of magnetic field corrugation on the equilibrium; the corrugation of the vacuum field is shown to lead to a rather moderate corrugation of the bubble boundary if the period of corrugation is sufficiently small. A valuable numerical result is the distribution of the diamagnetic field, which would be useful for optimizing the position of the wall-stabilization plates.