# subset to smaller bounding box shape area 

define XMIN 232800.
define YMIN 3863975.
define ZMIN -5000.
define XMAX 239500.
define YMAX 4447925.
define ZMAX 5000.

read avs points.inp mo

**

read avs INFILE mo
cmo / printatt / mo / -xyz- / minmax
pset/pbox/geom/xyz/1,0,0/ XMIN YMIN ZMIN / XMAX YMAX ZMAX 
pset/pdel/not pbox
rmpoint/ pset,get,pdel 
rmpoint/compress
cmo / printatt / mo / -xyz- / minmax
dump gmv subset.gmv mo

cmo / delete / cmotmp
cmo / select mo
cmo / setatt / mo / imt / 1 0 0 / 1
cmo / setatt / mo / itp / 1 0 0 / 0
filter / 1 0 0
rmpoint / compress

# sort / mo / index / ascending / ikey / zic yic xic
# reorder / mo / ikey
# cmo / DELATT / mo / ikey

connect noadd
resetpts / itp
cmo/setatt//ipolydat/no



dump / gmv / connect_tets.gmv / mo
dump / avs2 / connect_tets.inp / mo


finish