Numerical modeling of the constriction of the dc positive column in rare gases