*
* Carl W. Gable
* gable@lanl.gov
* LaGriT Control File
* http://lagrit.lanl.gov
* http://meshing.lanl.gov
*
* Demonstrate 2D mesh smoothing and mesh to mesh interpolation.
*
read / gmv / smooth.gmv / cmo
*
* Reset material ID numbers so that there are only
* two material types, slab and not_slab
*
* Slab is material number 6
*
resetpts / parents
rmpoint / compress
eltset / e_slab / itetclr / eq / 6
cmo / setatt / cmo / itetclr / 1 0 0 / 1
cmo / setatt / cmo / itetclr / eltset get e_slab / 2
resetpts / itp

loop / do / NLOOP 1 10 1 / loop_end &
smooth / position / esug / 1 0 0 ; recon

dump / gmv / smooth_merge_material.gmv / cmo
*
* Now interpolate element and node properties from the 
* original multi-material mesh onto the smoothed two 
* material mesh.
*
read / gmv / smooth.gmv / cmo_src
*
cmo / addatt / cmo / id_elem / vint / scalar / nelements
cmo / addatt / cmo / id_node / vint / scalar / nnodes

interpolate / map / cmo / id_node / 1 0 0 / cmo_src itetclr
interpolate / map / cmo / id_elem / 1 0 0 / cmo_src itetclr

dump / gmv / smooth_merge_interpolate.gmv / cmo

cmo / select / cmo
refine2d
resetpts / itp

loop / do / NLOOP 1 10 1 / loop_end &
smooth / position / esug / 1 0 0 ; recon

interpolate / map / cmo / id_node / 1 0 0 / cmo_src itetclr
interpolate / map / cmo / id_elem / 1 0 0 / cmo_src itetclr

dump / gmv / smooth_merge_ref_interpolate.gmv / cmo

finish