*
* Create a hexahedral icosahedral sphere mesh
* Number of nodes per shell is 642
*
define / n_icos / 642

cmo/create/cmo1/ / / hex
cmo / geometry / cmo1 / geom_sph_shells
cmo / constraint / cmo1 / dummy
createpts/sphere/8/2/n_icos/r1/r2/0.0,0.0,0.0/1,0,0.0/
createpts/sphere/8/3/n_icos/r2/r3/0.0,0.0,0.0/1,0,0.0/
filter / 1 0 0
rmpoint / compress
setpts
settets
dump/gmv/icosahedron_642.gmv/cmo1/
cmo / status / cmo1
cmo / delete / cmo1

finish