##############################
import rhinoscriptsyntax as rs
import random, time, sys   ###
sys.path.append("C:\Weronika\erasmus Graz\Digitale Methoden - Python\DM_lib")     ###
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(1)

#delete everything start from scratch
allobj = rs.AllObjects()
rs.DeleteObjects(allobj)

###_________________________________________#
### basic settings for grid to fit UN_slab  #
###                                         # 
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 from DM_lib as dm:
UnoGridCoords  = dm.UnoGridCoords           # get gridCoords L*D*H = 1760
UnoPanelCoords = dm.UnoPanelCoords          # get panelCoords [frontPanels, backPanels, sidePanels, upSidePanels] / default arguments s.u.
################____________________________#
################
lengthVec = lVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) )) ## rs.AddPoint( dm.getUnoCoord(1, 0, 0) )
depthVec  = dVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 1, 0), dm.getUnoCoord(0, 0, 0) )) ## rs.AddPoint( dm.getUnoCoord(0, 1, 0) )


dm.newEmptyLayer("UNO::setUp", [120,120,240])



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()[1], [200,200,200] )
        rs.ObjectColor(rs.AllObjects()[2], [200,200,200] )
        rs.ObjectColor(rs.AllObjects()[0], [200,200,200] )
        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()
        


coords = []
dirL = 1
Lval = 0
for L in range(40):
    d=0
    if L>5 and L%10 == 0:
        dirL *= -1
    
    
    d = L%1 
    cor = dm.getUnoCoord(Lval, d, L*1)
    coords.append( cor )
    Lval += dirL 
crv = rs.AddCurve( coords, 2)
dom = rs.CurveDomain( crv )[1]

a = dom*3
cor = rs.EvaluateCurve( crv, a )
rs.AddPoint( cor )
rad = 6
coords = rs.DivideCurve( crv, 30, create_points=1)
allCoords = []
anzCirc = 1



for cor in coords:
    a = rs.CurveClosestPoint(crv, cor)
    normEbene = rs.CurvePerpFrame(crv, a)
    circ = rs.AddCircle( normEbene, rad)