# read subset truncated at top and bottom

read avs tet_trunc.inp cmotet
resetpts itp
cmo printatt cmotet -all- minmax

********* Write FEHM setup files ********

cmo/addatt/cmotet/matsav/vint/scalar/nelements/linear/permanent/afgx/0.
cmo copyatt/ cmotet cmotet/ matsav itetclr
cmo setatt cmotet itetclr 1
cmo copyatt/ cmotet cmotet/ imt imat

filter 1,0,0
rmpoint compress
resetpts itp
quality

* Must be single defined mesh
cmo select cmotet
resetpts/parent
rmpoint compress
cmo status cmotet
cmo setatt cmotet ipolydat no

dump fehm tet cmotet ascii
cmo copyatt cmotet cmotet itetclr matsav
dump gmv tet_fehm.gmv cmotet

* translate to 0,0 for better voronoi volume precision
trans 1,0,0 zero
dump stor tet cmotet ascii

cmo select cmotet
extract/surfmesh/1 0 0/cmos/cmotet/external
cmo setatt cmos ipolydat yes
dump gmv tet_surfmesh.gmv cmos



finish