s , in a homogeneous magnetic field are analyzed by using a potential model with constituent quarks. To obtain anisotropic wave functions and the corresponding eigenvalues, the cylindrical Gaussian expansion method is applied, where the wave functions for transverse and longitudinal directions in the cylindrical coordinate are expanded by the Gaussian bases separately. Energy level structures in the wide range of magnetic fields are obtained and the deformation of the wave functions is shown, which reflects effects of the spin mixing, the Zeeman splitting and quark Landau levels. The contribution from the magnetic catalysis in heavy-light mesons is discussed as a change of the light constituent quark mass.