An efficient technique for building plasma equilibria with preset or constrained q profiles by solving the nonlinear Grad–Shafranov equation iteratively is proposed. The technique uses a set of dynamic virtual poloidal field probes in the plasma core region. Numerical examples are shown for monotonic q profiles and negative central magnetic shear equilibria which are obtained with this technique. As a proof of principle, the technique is applied on the EAST tokamak, to constrain the plasma equilibrium reconstruction with the experimental q = 1 surface information from the ECE diagnostic in addition to the external magnetic measurements.