Allow specifying q=const at BField init.
Made Circular Tokamak Model. Fixed cap_r errors. Close to desired poincare plot