# 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/addatt/cmotet/cc_neg/vint/scalar/nnodes/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

# check to see if neg ccoefs are on outside of mesh
quality
pset ncoef attribute ccoef 1,0,0 lt -0.001
pset pout attribute itp 1,0,0 eq 10.
pset pin not pout
pset pBAD inter pin ncoef
cmo setatt cmotet cc_neg pset,get,ncoef 1

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


# Extract surface mesh for viewing voronoi polygons
cmo select cmotet
extract/surfmesh/1 0 0/cmos/cmotet/external

cmo copyatt / cmos cmos / itetclr itetclr1
cmo setatt cmos ipolydat yes
dump gmv tet_surfmesh.gmv cmos



finish