from netica import Netica
import numpy as np

# initialize class
ntc = Netica()

# create new environment
env = ntc.newenv()
# initialize environment
ntc.initenv(env)
# create new net
net_p = ntc.newnet('BNZeebrugge', env)

# define nodes
nodeBC1 = ntc.newnode('PeakWaterLevel', 5, net_p)
ntc.setnodelevels(nodeBC1, 5,
                  np.asarray([6.35, 7.1, 7.4, 7.8, 7.9], dtype='float64'))
ntc.setnodetitle(nodeBC1, 'Peak Water Level (m)')

nodeBC2 = ntc.newnode('MaxSignificantWaveHeight', 5, net_p)
ntc.setnodelevels(nodeBC2, 5,
                  np.asarray([5.2, 5.7, 5.9, 6.08, 6.2], dtype='float64'))
ntc.setnodetitle(nodeBC2, 'Max. significant wave height (m)')

nodeR1 = ntc.newnode('Location_Houses', 4, net_p)
ntc.setnodetitle(nodeR1, 'Houses - Location')
ntc.setnodestatenames(nodeR1, "Area1, Area2, Area3, Area4")
ntc.setnodestatetitle(nodeR1, 0, 'Area 1 (has 283 houses)')
ntc.setnodestatetitle(nodeR1, 1, 'Area 2 (has 759 houses)')
ntc.setnodestatetitle(nodeR1, 2, 'Area 3 (has 383 houses)')
ntc.setnodestatetitle(nodeR1, 3, 'Area 4 (has 273 houses)')
from netica import Netica
import numpy as np
 
# initialize class
ntc = Netica()
 
# create new environment
env = ntc.newenv()
# initialize environment
ntc.initenv(env)
# create new net
net_p = ntc.newnet('BNZeebrugge', env)
 
# define nodes
nodeBC1 = ntc.newnode('PeakWaterLevel', 5, net_p)
ntc.setnodelevels(nodeBC1, 5, np.asarray([6.35, 7.1, 7.4, 7.8, 7.9], dtype='float64'))
ntc.setnodetitle(nodeBC1,'Peak Water Level (m)')

nodeBC2 = ntc.newnode('MaxSignificantWaveHeight', 5, net_p)
ntc.setnodelevels(nodeBC2, 5, np.asarray([5.2, 5.7, 5.9, 6.08, 6.2], dtype='float64'))
ntc.setnodetitle(nodeBC2,'Max. significant wave height (m)')

nodeR1 = ntc.newnode('Location_Houses', 4, net_p)
ntc.setnodetitle(nodeR1,'Houses - Location')
ntc.setnodestatenames(nodeR1, "Area1, Area2, Area3, Area4")
ntc.setnodestatetitle(nodeR1, 0, 'Area 1 (has 283 houses)')
ntc.setnodestatetitle(nodeR1, 1, 'Area 2 (has 759 houses)')
ntc.setnodestatetitle(nodeR1, 2, 'Area 3 (has 383 houses)')
ntc.setnodestatetitle(nodeR1, 3, 'Area 4 (has 273 houses)')

nodeH1_R1 = ntc.newnode('MaxInundation_Houses', 0, net_p)