# create inner cylinder point distribution
# these points are explicitly defined and added
# to rim points for distibution that uses smoothing


# -------------------------------------
# cylinder (well) object parameters
# RLEN is 1st distance from center, 2nd distance
# NRAYS number of rays around the center or sections
# NRADIAl points to distribute in the R direction from R0 to R1
# Ratio factor so radial point will have dr_n = R_RATIO*dr_n-1

define ID_CYL 2
define RLEN 0.0
define NRAYS 8  
define NRADIAL 3 
define R0 0.
define R1 2.
define R_RATIO 1.

cmo / create / mo / / / triplane
createpts/rtz/NRADIAL,NRAYS,1/R0,0.0, 0.0 /R1,360.,RLEN/ &
            1,1,1/1 0 0/ R_RATIO 1.0 1.0

# remove duplicate points
filter / 1 0 0
rmpoint / compress

# set imt value to 2 to distinquish from rest of well points
pset/pring/seq/0,0,0
cmo/setatt/mo/imt1/pset,get,pring/ID_CYL

dump gmv  cylinder_well.gmv mo
dump avs2 cylinder_well.inp mo 

finish