# subset to smaller bounding box shape area 
#
# FULL MESH
#  xic    1.060500000E+05  4.510500000E+05 3.450000000E+05   1102042  
#  yic    3.869700000E+06  4.442200000E+06 5.725000000E+05   1102042  
#  zic   -4.021226013E+03  1.698531097E+03 5.719757111E+03   1102042  


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

define YMIN 3910500.
define YMAX 4445000.


read avs points.inp mo
cmo / printatt / mo / -xyz- / minmax

#########
# SUBSET
#########
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 printatt mo -all- minmax

#########
# ADD FLAT BOTTOM and FLAT TOP
#########
cmo copy motmp mo
cmo setatt motmp zic -4000.
cmo select motmp
filter 1,0,0
rmpoint compress

cmo setatt motmp imat 10
cmo setatt motmp xval 10.0
cmo setatt motmp nlayer 0 

cmo copy motmp2 motmp
cmo setatt motmp2 zic 1700.
cmo setatt motmp2 imat 11
cmo setatt motmp2 xval 11. 
cmo setatt motmp2 nlayer 50 

#### MERGE #######
cmo printatt mo -all- minmax
cmo printatt motmp -all- minmax
cmo printatt motmp2 -all- minmax

cmo copy motet motmp
copypts/ motet / mo / 0,0/ 1,0,0
copypts/ motet / motmp2 / 0,0/ 1,0,0
cmo printatt motet -all- minmax

dump gmv points_to_connect.gmv motet


#### CONNECT #######

cmo / select motet
cmo / setatt / motet / imt / 1 0 0 / 1
cmo / setatt / motet / itp / 1 0 0 / 0
filter / 1 0 0
rmpoint / compress

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

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

## color elements
cmo setatt motet itetclr 15
cmo select motet
infile color_elements.mlgi

cmo / printatt / motet / -all- / minmax
dump / gmv / connect_tets.gmv / motet
dump / avs2 / connect_tets.inp / motet

#### REMOVE TOP AND BOTTOM #######
rmmat/10/
rmmat/11/
rmpoint compress
resetpts itp

dump gmv tet_trunc.gmv motet
dump avs tet_trunc.inp motet


finish