# tag and write a zone file for the center points

read avs cyl_tet_delaunay.inp mo

# overwrite previous zones for imt 1 and 2
# we need node values for well, so change background to 5
pset/p2/attribute imt/1,0,0/ eq 2
pset/p1/attribute imt/1,0,0/ eq 1
cmo setatt mo imt pset,get,p2 5
cmo setatt mo imt pset,get,p1 2


# pset select the center node
pset/px/attribute xic/1,0,0/ eq 0.
pset/pcenter/ attribute yic /pset,get,px/ eq 0.

cmo setatt mo imt pset,get,pcenter 1
cmo printatt mo imt minmax

dump zone_imt cyl_center mo 1 
dump zone_imt cyl_well_sides mo 2 
dump zone_imt cyl_background mo 5 

dump gmv cyl_zones.gmv mo

finish