예제 #1
0
def cylinder(center, R, H, N=100, ntype='STRUCT'):
    """Create a cylinder of center C, radius R and hieght H.
    Usage: cylinder((xc,yc,zc), R, H, N)"""
    A = Geom.cylinder(center, R, H, N, ntype)
    if ntype == 'STRUCT':
        return [
            C.convertArrays2ZoneNode('cyl-part1', [A[0]]),
            C.convertArrays2ZoneNode('cyl-part2', [A[1]]),
            C.convertArrays2ZoneNode('cyl-part3', [A[2]]),
            C.convertArrays2ZoneNode('cyl-part4', [A[3]]),
            C.convertArrays2ZoneNode('cyl-part5', [A[4]]),
            C.convertArrays2ZoneNode('cyl-part6', [A[5]]),
            C.convertArrays2ZoneNode('cyl-part7', [A[6]]),
            C.convertArrays2ZoneNode('cyl-part8', [A[7]]),
            C.convertArrays2ZoneNode('cyl-part9', [A[8]]),
            C.convertArrays2ZoneNode('cyl-part10', [A[9]]),
            C.convertArrays2ZoneNode('cyl-part11', [A[10]])
        ]
    else:
        return C.convertArrays2ZoneNode('cyl', [A])
예제 #2
0
# - cylinder (array) -
import Geom as D
import Converter as C

a = D.cylinder((0,0,0), 1., 10.)
C.convertArrays2File(a, 'out.plt')