# Select center cylinder set of nodes and cells for well
# MUST DEFINE MAT_ID, WELL_ID, CMO

cmo setatt CMO imt MAT_ID 
cmo setatt CMO itetclr MAT_ID 

# pset select the center node
pset/px/attribute xic/1,0,0/ eq 0.
pset/pcenter/ attribute yic /pset,get,px/ eq 0.
# eltset select cells connected to center node
eltset/ewell/ inclusive pset,get,pcenter
# pset select all node vertices of well cells
pset/pwell/ eltset / ewell

# overwrite material MAT_ID with WELL_ID for well
cmo setatt CMO imt pset,get,pwell WELL_ID 
cmo setatt CMO itetclr eltset,get,ewell WELL_ID 

finish