*get, xk, kp, Nummer, loc, x
*get, yk, kp, Nummer, loc, y
*get, zk, kp, Nummer, loc, z
nsel, s,,, node(xk, yk, zk)
esln, s
nsle, s

*get, element, elem, 0, num, min
drauf = 0
*dowhile, element

	*get, x0, node, nelem(element, 1), loc, x
	*get, y0, node, nelem(element, 1), loc, y
	*get, xa, node, nelem(element, 2), loc, x
	*get, ya, node, nelem(element, 2), loc, y
	*get, xb, node, nelem(element, 4), loc, x
	*get, yb, node, nelem(element, 4), loc, y
	a1 = xa - x0
	a2 = ya - y0
	b1 = xb - x0
	b2 = yb - y0
	c1 = xk - x0
	c2 = yk - y0
	alpha = (c1*b2 - b1*c2)/(a1*b2 - a2*b1)
	beta  = (a1*c2 - c1*a2)/(a1*b2 - a2*b1)

	*if, alpha, ge, 0, and, beta, ge, 0, then
		*get, x0, node, nelem(element, 3), loc, x
		*get, y0, node, nelem(element, 3), loc, y
		*get, xa, node, nelem(element, 2), loc, x
		*get, ya, node, nelem(element, 2), loc, y
		*get, xb, node, nelem(element, 4), loc, x
		*get, yb, node, nelem(element, 4), loc, y
		a1 = xa - x0
		a2 = ya - y0
		b1 = xb - x0
		b2 = yb - y0
		c1 = xk - x0
		c2 = yk - y0
		alpha = (c1*b2 - b1*c2)/(a1*b2 - a2*b1)
		beta  = (a1*c2 - c1*a2)/(a1*b2 - a2*b1)
		*if, alpha, ge, 0, and, beta, ge, 0, then
			drauf = 1
		*endif	
	*endif
	
	*get, element, elem, element, nxth
	*if, drauf, eq, 1, then
		element = 0
	*endif
*enddo