###################################
### DM2_w24 # AGruber@tugraz.at ###
### ue_05 UN_headquaters  NYC   ###
###################################


##############################
import rhinoscriptsyntax as rs
import random, time, sys   ###
sys.path.append("P:/WWW/lsb24802/dm2/")     ###
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(state=1) ### nomen est omen / printDisplay(state=1, scale=1000, thickness=1, color="Display")
rs.EnableRedraw(0)       ### MAC >> problem !



###_________________________________________#
### basic settings for grid to fit UN_slab  #
###                                         # ! no need 2 change !
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]  )


################
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) )
################

if 1: ### setup gebaeude
    if 1:
        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 1:
        rs.CurrentLayer("UNO::setUp::crv")
        dm.UN_slab(showAll=1)
        rs.ZoomExtents()

#################
### 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])

p000 = rs.VectorAdd(gC(10,0,1), rs.VectorScale(-dVec,10))
p00 = rs.VectorAdd(gC(9,0,1), rs.VectorScale(-dVec,10))
p0 = rs.VectorAdd(gC(5,0,20), rs.VectorScale(-dVec, 10))
p1 = rs.VectorAdd(p0,rs.VectorScale(lVec, 75))
p2 = rs.VectorAdd(p1,rs.VectorScale(dVec, 50))
p3 = rs.VectorAdd(p2,rs.VectorScale(-lVec, 150))
p4 = rs.VectorAdd(p3,rs.VectorScale(-dVec, 20))
p5 = rs.VectorAdd(gC(10,0,25), rs.VectorScale(-lVec, 15))
p6 = rs.VectorAdd(p5, rs.VectorScale(-dVec, 100))
coords =[p000,p00,p0,p1,p2,p3,p4,p5,p6]

#rs.AddCurve(coords,1)
crv = rs.AddCurve(coords,2)

anzCircs = 60
divCircs = 30
coords= rs.DivideCurve(crv,anzCircs, create_points=0)

if 1: ###toggle tube
    allCircCoordsLists = []
    for cor in coords:
        pass
        paraX=rs.CurveClosestPoint(crv, cor)
        planeX = rs.CurvePerpFrame(crv,paraX)
        circ = rs.AddCircle(planeX,5)
        coords=rs.DivideCurve(circ, divCircs,0)
        rs.DeleteObject(circ)
        allCircCoordsLists.append(coords)
    
    for i in range (divCircs):
        coords=[]
        for list in allCircCoordsLists:
            p0=list[i]
           #rs.AddPoint(p0)
            coords.append(p0)
            
        #rs.AddCurve(coords,1)
        rs.AddInterpCurve(coords, degree=3, knotstyle=3)
#supports
anzSprt= 10
ankTube= rs.DivideCurve(crv, anzSprt-1, create_points=1)
ankBoden= rs.DivideCurve(unoBaseCrv, anzSprt,create_points=1)
for i in range(anzSprt):
        rs.AddLine(ankTube[i], ankBoden[i])
        


#supports ueber curveclosestpoint, definitions

#____________________________:here you end

###################################
dm.zA( proz=0.95 )
rs.EnableRedraw(1) ### 4 the MACs
#dm.printDisplay(1)
dm.newEmptyLayer("Default")