Procedure from library primdec.lib (see primdec_lib).
primdec.lib
equidimMax(i); i ideal
ideal of equidimensional locus (of maximal dimension) of i.
Example:
LIB "primdec.lib"; ring r = 32003,(x,y,z),dp; ideal i = intersect(ideal(z),ideal(x,y),ideal(x2,z2),ideal(x5,y5,z5)); equidimMax(i); → _[1]=z