#######################
#-*- coding: utf-8 -*-

#########################
### AGruber@tugraz.at ###
### DEM demo  2024/12 ###
#########################


##############################
import rhinoscriptsyntax as rs
import sys
sys.path.append("P:/WWW/tkonrad/DM2/")  
import DEM_lib as dem               ### "dem" is alias 4 calling defs @ OSM_lib
##############################      reload (dem)


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

#######################
#HGTpath = "D:/DEM/HGT/DEM1/"
HGTpath = "P:/hgtData/"
HGTpath = "L:/hgtData/"
#######################

if 1:
    #dem.eAA()
    # lampedusa 35° 30′ 56″ N, 12° 34′ 23″ E
    dem.makeDEM(35, 28, 4, 12, 32, 4, HGTpath=HGTpath, west=0, DEMsec=1, deg=1, eleFact=1.0, move2Guinea=1)
    # wien urania 48 12 42 N 16 23 1 E
    #dem.makeDEM(48, 11, 3, 16, 21, 4, HGTpath=HGTpath, west=0, DEMsec=1, deg=1, eleFact=1.0, move2Guinea=1)

    
    # NYC
    #dem.makeDEM(40, 44, 2,  74,  2, 2, HGTpath=HGTpath, west=1)
    #dem.makeDEM(40, 42, 6,  74,  0, 3, HGTpath=HGTpath, west=1)
    
    #innsbruck
    #dem.makeDEM(47, 16, 4,  11, 23, 4, HGTpath=HGTpath, west=0, DEMsec=1, deg=1, eleFact=1.0, move2Guinea=0)
    
    pass
    #dem.makeDEM(37, 3, 3, 41, 12, 5, HGTpath=HGTpath, west=0, DEMsec=3, deg=1, eleFact=1.0, move2Guinea=0)