**********************************************************
* Example driver for interpolate_mesh_to_points.lgi LaGriT
* control file.
* 1) Build a simple 3D, tet mesh in a unit cube.
* 2) Create some node attributes using math functions and
*    the x,y,z coordinates.
* 3) Call the interpolate_mesh_to_points.lgi control file
*    to demonstrate interpolation from the mesh just created
*    to a point distribution.
*
* Carl Gable
* gable@lanl.gov
*
**********************************************************
*
cmo create cmo
* Create a uniform point distribution on the unit cube.
*
createpts/xyz/11 11 11/0. 0. 0./1. 1. 1./1 1 1
cmo/setatt/cmo/imt/1 0 0/ 1
cmo/setatt/cmo/itp/1 0 0/ 0
define/ATTRIBUTE1/xdisplace
define/ATTRIBUTE2/ydisplace
define/ATTRIBUTE3/zdisplace
*
* Connect points into Delaunay tetrahedra.
*
connect
*
* Turn off some attributes that are usually output by default.
*
cmo/modatt/cmo/imt/ioflag/l
cmo/modatt/cmo/itp/ioflag/l
cmo modatt cmo isn ioflag l
cmo modatt cmo icr ioflag l
*
* Add some attributes and fill them with various floating
* point values using math functions and x,y,z coordinates.
*
cmo / addatt / cmo / ATTRIBUTE1 / vdouble / scalar / nnodes
cmo / addatt / cmo / ATTRIBUTE2 / vdouble / scalar / nnodes
cmo / addatt / cmo / ATTRIBUTE3 / vdouble / scalar / nnodes

cmo / copyatt / cmo cmo ATTRIBUTE1  xic
cmo / copyatt / cmo cmo ATTRIBUTE2  yic
cmo / copyatt / cmo cmo ATTRIBUTE3  zic
*
* Modify the attributes with the math functions.
*
* f = e**x
*
math/exp/cmo/ATTRIBUTE1/1,0,0/cmo/ATTRIBUTE1
*
* f = y**2
*
math/power/cmo/ATTRIBUTE2/1 0 0/cmo/ATTRIBUTE2/2.0
*
* f = sin(z)
*
math/sin/cmo/ATTRIBUTE3/1 0 0/cmo/ATTRIBUTE3
*
dump / avs / unit_cube_tet_11_11_11.inp / cmo
dump / gmv / unit_cube_tet_11_11_11.gmv / cmo / ascii
dump / avs / input.inp / cmo

cmo / delete / cmo
*
* Run interpolation steps.
*
infile interpolate_mesh_to_points.lgi

finish