# after tri top is created, stack into 3D

define TRI_IN top_tri_smooth.inp
define TRI_IN top_tri_recon.inp

read avs TRI_IN mo
cmo printatt mo -xyz- minmax

# transpose y and z axis and stack
# switch z to y and zic to (zero)
cmo copyatt mo mo yic zic                                             
cmo setatt mo zic 0.                                                       
resetpts itp 
dump avs stack_bot.inp mo

cmo prinatt mo -xyz- minmax
cmo copy mo_top mo
cmo setatt mo_top zic 1000.
dump avs stack_top.inp mo_top

# stack layers from 0 to 1000 every 25m
cmo create mo_stack
stack/layers/avs/ stack_bot.inp stack_top.inp 1 39 / flip
dump gmv tmp_tri_layers.gmv mo_stack 

# fill layers with prism elements
stack/ fill / mo_pri / mo_stack
cmo setatt mo_pri itetclr 1
cmo setatt mo_pri imt 1
resetpts itp
dump gmv tmpxy_stack.gmv mo_pri

# transpose y and z axis and write
cmo/addatt/mo_pri/zic_sav/VDOUBLE/scalar/nnodes/linear/
cmo printatt mo_pri -xyz- minmax
cmo copyatt mo_pri mo_pri zic_sav zic      

# switch y to z and z to y
cmo copyatt mo_pri mo_pri zic yic      
cmo copyatt mo_pri mo_pri yic zic_sav      
resetpts itp
cmo printatt mo_pri -xyz- minmax

# Warning, these prisms are inside out and useless (almost)
# write points for use with base points to connect 
dump gmv tri_top_prism.gmv mo_pri
dump avs2 tri_top_prism_pts.inp mo_pri 1 0 0 0
cmo printatt mo_pri -xyz- minmax

# write volume fill version for finding tets outside
# This is a prism mesh with no interior points
cmo delete mo
read avs TRI_IN mo                                                  
extrude/mo2/mo/const/1000./volume/0. 1. 0.                                     

dump gmv tri_top_extrude.gmv mo2
dump avs tri_top_extrude.inp mo2

finish