Create a tet sphere mesh with icosahedral point distribution
* Number of nodes per shell is 642
*

cmo/create/cmo1/ / / tet

infile lagrit.input_geom

cmo / geometry / cmo1 / geom_sph_shells
cmo / constraint / cmo1 / dummy

define / NPOINT / 162

createpts/sphere/2/3/NPOINT/0./r1/0.0,0.0,0.0/1,0,0.0/
createpts/sphere/2/3/NPOINT/r1/r2/0.0,0.0,0.0/1,0,0.0/
createpts/sphere/2/3/NPOINT/r2/r3/0.0,0.0,0.0/1,0,0.0/
filter / 1 0 0
rmpoint / compress
setpts

connect
settets
rmregion / r_r1
rmpoint / compress

extract / surfmesh / 1 0 0  / cmo_tri / cmo1

cmo / delete / cmo1
cmo / select / cmo_tri

eltset / e_1 / facecol / eq / 1
rmpoint / element / eltset get e_1
eltset / e_2 / facecol / eq / 2
rmpoint / element / eltset get e_2

rmpoint / compress

finish