* create a 3-D hexahedral mesh object and name it 3dmesh

    cmo/create/3dmesh///hex/

    * create a cylinder inside a cylinder

    * inner cylinder will be replaced by tet grid

    surface/cyl1/reflect/cylinder/0.0,0.0,0.0/0.0,0.0,10.0/5.0/

    surface/top/reflect/plane/0.,0.,10./1.,0.,10./1.,1.,10. /

    surface/bot/reflect/plane/0.,0.,0./1.,0.,0./1.,1.,0. /

    surface/cyl2/intrface/cylinder/0.0,0.0,1.0/0.0,0.0,9.0/4.0/

    surface/cutpl1/intrface/plane/0.,0.,1./1.,0.,1./1.,1.,1./

    surface/cutpl2/intrface/plane/0.,0.,9./1.,0.,9./1.,1.,9./

    *define geometric regions

    region/inner/ lt cyl2 and lt cutpl2 and gt cutpl1 /

    region/outer/ le cyl1 and le top and ge bot and &

    ( ge cyl2 or ge cutpl2 or le cutpl1 )

    * define material regions

    mregion/inner/ lt cyl2 and lt cutpl2 and gt cutpl1 /

    mregion/outer/ le cyl1 and le top and ge bot and &

    ( gt cyl2 or gt cutpl2 or lt cutpl1 )

    * distribute points in cylindrical arrangement

    rz/rtz/11,61,21/0.,0.,0./5.,360.,10./1,1,1/

    * assign material colors to the points

    * identify points that are on material interfaces

    * identify constrained points

    * connect hex mesh

    rzbrick/xyz/61,11,21/1,0,0/connect/

    zq/itp/1,0,0/0

    zq/imt/1,0,0/0

    zq/icr/1,0,0/0

    setpts

    dump/x3d/x3d1

    settets

    zq/itp/100,130,1

    dump/gmv/gmv1

    * remove inner region

    rmmat/inner

    rmpoint/compress

    geniee

    *zq/xic/1,0,0

    *zq/itp/1,0,0

    dump/gmv/gmv2

    * grid tet part

    cmo/create/tetgrid///tet

    rz/rtz/9,61,17/0.,0.,1./1.,360.,9./1,1,1/

    filter/1,0,0/

    search

    zq/imt/1,0,0/1

    zq/xic/1,0,0

    zq/itp/1,0,0

    * add meshes

    dump/x3d/x3d1

    addmesh/pyramid/hybrid/tetgrid/3dmesh/

    surface/cyl1/reflect/cylinder/0.0,0.0,0.0/0.0,0.0,6.0/2.0/

    surface/top/reflect/plane/0.,0.,6./1.,0.,6./1.,1.,6. /

    surface/bot/reflect/plane/0.,0.,0./1.,0.,0./1.,1.,0. /

    surface/cyl2/intrface/cylinder/0.0,0.0,2.0/0.0,0.0,4.0/1.0/

    surface/cutpl1/intrface/plane/0.,0.,2./1.,0.,2./1.,1.,2./

    surface/cutpl2/intrface/plane/0.,0.,4./1.,0.,4./1.,1.,4./

    *define geometric regions

    region/inner/ le cyl2 and le cutpl2 and ge cutpl1 /

    region/outer/ le cyl1 and le top and ge bot and &

    ( gt cyl2 or gt cutpl2 or lt cutpl1 )

    * define material regions

    mregion/inner/ lt cyl2 and lt cutpl2 and gt cutpl1 /

    mregion/outer/ le cyl1 and le top and ge bot and &

    ( gt cyl2 or gt cutpl2 or lt cutpl1 )

    filter/1,0,0/

    dump/x3d/x3d2

    geniee

    zq/imt/1,0,0/0/

    zq/itp/1,0,0/0/

    zq/icr/1,0,0/0/

    setpts

    settets

    * dump mesh to some output form

    dump/gmv/gmv.hybrid/hybrid

    * terminate processing

    finish