##############################
import rhinoscriptsyntax as rs
import random, time, sys   ###
sys.path.append("P:/")     ###
import DM_lib as dm        ### reload(dm)
##############################      

rs.UnitSystem(4)                                        # km = 5, meters = 4, cm = 3 etc
rs.ShowGrid(None, 0)                                    # grid > 0 = off
rs.ShowGridAxes(None, 1)                                # y/y/z axen display > 0/1 = off/on
rs.ViewDisplayMode(rs.CurrentView(), "wireframe")
rs.Command("cplane w t enter", 0)						# cPlane World Top
dm.PointRadius(displayModeX=0, rad=3, styl=3)
dm.printDisplay(0)                                      # nomen est omen
rs.EnableRedraw(0) ### MAC >> uncomment !




floors = H  = dm.H = 40                     # default=40 / incl roof
slabs  = L  = dm.L = 11                     # default=11
depth  = D  = dm.D =  4                     # default= 4 / division in building_depth
floorHeight = fH = dm.fH = 4.0              # default= 4.0 / 4.0*(H-1) = 156 meters
                                            #
#############                               # get gridCoords L*D*H = 1760
UnoGridCoords = dm.getUnoGridCoords()       # get from DM_lib as dm
#################___________________________#


dm.newEmptyLayer("UNO::setUp::points", [150,150,200])
dm.newEmptyLayer("UNO::setUp::geo", [150, 20,150])
dm.newEmptyLayer("UNO::setUp::crv", [200, 100, 100])
dm.newEmptyLayer("UNO::setUp", [120,120,140])


##################
depthVec  = dVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 2, 0) ))
lengthVec = lVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) ))
##################

if 1: 
    if 1:
        rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) )
        rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 1, 0) )
        rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 0, 1) )
        rs.ObjectColor(rs.AllObjects()[2], [0, 222, 0] )
        rs.ObjectColor(rs.AllObjects()[1], [222, 0, 0] )
        rs.ObjectColor(rs.AllObjects()[0], [0, 0, 222] )
        rs.ObjectPrintWidth( rs.AllObjects()[0:3], 1.0 )
        
        p0 = dm.getUnoCoord(  0, 0, 0)
        p1 = dm.getUnoCoord( 10, 0, 0)
        p2 = dm.getUnoCoord( 10, 3, 0)
        p3 = dm.getUnoCoord(  0, 3, 0)
        unoBaseCoords = [p0, p1, p2, p3, p0]
        unoBaseCrv = rs.AddCurve( unoBaseCoords, 1)
        unoTopCoords = [dm.getUnoCoord(  0, 0, 39),
                        dm.getUnoCoord( 10, 0, 39),
                        dm.getUnoCoord( 10, 3, 39),
                        dm.getUnoCoord(  0, 3, 39),
                        dm.getUnoCoord(  0, 0, 39),
                       ]
        unoTopCrv = rs.AddCurve( unoTopCoords, 1)
        rs.ZoomExtents()
    if 0:
        ##############################
        rs.CurrentLayer("UNO::setUp::points")
        rs.AddPoints( UnoGridCoords )
        for i in range(L): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(i, 0, 0) ), [0, 222, 0] )
        for i in range(D): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(0, i, 0) ), [222, 0, 0] )
        for i in range(H): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(0, 0, i) ), [0, 0, 222] )
        
        ##############################
    if 0:
        rs.CurrentLayer("UNO::setUp::crv")
        for i in range(8):
            p0 = dm.getUnoCoord( random.randint(0,L-1), random.randint(0,D-1), random.randint(0,H-1) )
            p1 = dm.getUnoCoord( random.randint(0,L-1), random.randint(0,D-1), random.randint(0,H-1) )
            rs.AddCurve( [p0, p1] )
            dm.textDots(  [p0, p1] )
        rs.ZoomExtents()
    if 0:
        coords = UnoGridCoords[:] ### make COPY - not equal !
        #random.shuffle( coords )
        rs.AddCurve( coords[0:], 2 )


#################
### abbreveation:
def gC (l=1, d=2, h=3):
    cor = dm.getUnoCoord(l, d, h)
    return cor
## rs.AddPoint( gC(1,0,1) )
#################


#_____________________________here you go:
dm.newEmptyLayer("myPROJ", [100,100,100])

#cor = dm.getUnoCoord(5,0,0)
#rs.AddPoint(cor)
#rs.AddPoint( dm.getUnoCoord(5,0,2))
#rs.AddPoint(dm.getUnoCoord(5,3,2))


#for z in range (2,40):
#    rs.AddPoint(dm.getUnoCoord(5,3,z))
#    if z%4==0:
#        rs.AddPoint(dm.getUnoCoord(7,3,z))
#        #rs.AddLine(gC(7,3,z),gC(7,3,z))

p0= gC(0,2,0)
p1= gC(10,3,10)
p2=rs.VectorAdd(p1, rs.VectorScale(dVec,80))
p3= gC(2,1,15)
p4=rs.VectorAdd(p3, rs.VectorScale(dVec,-40))
p5= gC(9,3,20)
p6= rs.VectorAdd(p5, rs.VectorScale(dVec,50))
p7= gC(5,0,25)
p8=rs.VectorAdd(p7, rs.VectorScale(dVec,-40))
p9= gC(8,3,28)
p10=rs.VectorAdd(p9, rs.VectorScale(dVec,20))
p11= gC(4,0,30)
p12=gC(5,3,32)
p13=gC(7,3,35)
p14=gC(6,0,38)
p15=gC(5,2,39)


coords= [p0,p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15]
#rs.AddCurve(coords, 1)
crv= rs.AddCurve(coords, 2)
anzCircs=64
divCircs=18
coords= rs.DivideCurve(crv, anzCircs, create_points=0)

allCircCoordsLists= []

for cor in coords:
    rs.AddCircle(cor,5.5)
    paraX= rs.CurveClosestPoint(crv,cor)
    planeX=rs.CurvePerpFrame(crv,paraX)
    circ= rs.AddCircle(planeX,5.5)
    coords= rs.DivideCurve(circ,divCircs, 1)
    rs.DeleteObject(circ)
    allCircCoordsLists.append(coords)




dom=rs.CurveDomain(crv)[1]
print "dom", dom
paraX= dom*0.25
cor= rs.EvaluateCurve(crv,paraX)
rs.AddPoint(cor)
planeX= rs.CurvePerpFrame(crv,paraX)
#rs.AddCircle(planeX,5.5)
rs.AddRectangle(planeX,15,15)

i=0
for i in range(divCircs):
    coords=[]
    for list in allCircCoordsLists:
        p0=list[i]
        #rs.AddPoint(p0)
        coords.append(p0)
    rs.AddCurve(coords,3)
    rs.AddInterpCurve(coords,degree=3,knotstyle=3)
    


#____________________________:here you end

###################################
rs.ZoomExtents()
rs.EnableRedraw(1) ### 4 the MACs
dm.printDisplay(1)
dm.newEmptyLayer("Default")