# 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  
#
# NEW MESH Dec 2013
#  xic     4.201000000E+05  4.299000000E+05   9.800000000E+03   117450    
#  yic     4.307100000E+06  4.316900000E+06   9.800000000E+03   117450    
#  zic    -2.208199951E+03  1.081999969E+02   2.316399948E+03   117450   
#
#  xval    1.000000000E+00  9.000000000E+00   8.000000000E+00   117450
#  nlayer  1                47                46                117450
#
#  Note this mesh is numbered from top to bottom
#
#
# 50 km X 50 km domain. Four corners for this domain:
# 
# SW x240000 y4043000 mins
# SE x290000 y4043000
# NE x290000 y4093000 maxs
# NW x240000 y4093000
# 

define MO mo

# Define for current data file
define XMIN 420100.
define YMIN 4307100.
define ZMIN -2210.
define ZMIN -5000.

define XMAX 430100.
define YMAX 4317100.
define ZMAX 110.
define ZMAX 5000.

define Z_BOT -4000.
define Z_BOT -2400.

define Z_TOP  1700.
define Z_TOP   300.

# flat and top layers have next numbers after max xval
define Z_BOT_ID  10
define Z_BOT_VAL 10.

define Z_TOP_ID  11
define Z_TOP_VAL 11.

# tag top and bottom flat points
define TOP_LAYER 50
define BOT_LAYER 0


read avs points.inp mo
cmo / printatt / mo / -xyz- / minmax
cmo addatt mo imat VINT/scalar/nnodes
cmo copyatt/ mo mo / imat xval
cmo printatt mo -all- 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
# this helps connect and makes it easier to 
# remove "outside" elements
#########
cmo copy motmp mo
cmo setatt motmp zic Z_BOT 
cmo select motmp
filter 1,0,0
rmpoint compress
  cmo setatt motmp imat Z_BOT_ID
  cmo setatt motmp xval Z_BOT_VAL 
  cmo setatt motmp nlayer BOT_LAYER 

cmo copy motmp2 motmp
  cmo setatt motmp2 zic    Z_TOP 
  cmo setatt motmp2 imat   Z_TOP_ID
  cmo setatt motmp2 xval   Z_TOP_VAL 
  cmo setatt motmp2 nlayer TOP_LAYER 

#### 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

# STOP and Evaluate subset points
finish


#### 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