reset

import stl

/parameters

$r = 65

$h = 325

$number = 1

$cylinder_number = 1452

vertex create "translation_start" coordinates 0 0 0

vertex create "rotation_strat" coordinates $r 0 $h

vertex create "alignment_strat" coordinates $r 0 -$h

$volume_f = "f_volume."

$volume_r = "volume."

$vertex_f = "f_vertex."

$vertex = "vertex."

do para "$number" init 1 cond ($number .le. $cylinder_number)

$volume_end   =  $volume_f + NTOS($number)

volume delete $volume_end lowertopology onlymesh

enddo

do para "$number" init 1 cond ($number .le. $cylinder_number)

 $volume_start =  $volume_r + NTOS($number+2000)

$volume_end   =  $volume_f + NTOS($number)

volume create $volume_start height (10*$r) radius1 $r radius2 $r radius3 $r zaxis frustum

$rotation_end  = $vertex_f + NTOS($number*2-1)

$alignment_end = $vertex_f + NTOS($number*2)

$translation_end = $vertex + NTOS($number+10000)

vertex create $translation_end centroid volume $volume_end

volume align $volume_start translation "translation_start" $translation_end rotation1 \

 "rotation_strat" $rotation_end rotation2 "alignment_strat" $alignment_end

volume delete $volume_end lowertopology

enddo