########### Knoten laden ###################
INKNOTEN = open('knoten.txt', 'r')
knoten = []
for line in INKNOTEN:
	label,x,y = line.strip('\n').split('\t')
	label = label.strip()
	x = x.strip()
	y = y.strip()
	tmp = [int(label),float(x),float(y)]
	knoten.append(tmp)
	
#print knoten

INKNOTEN.close()

############ Elemente laden ##################
INELE = open('elemente.txt','r')
elem = []
for line in INELE:
	label,k1,k2,k3,k4 = line.strip(' \n').split('\t')
	tmp = [int(label), int(k1), int(k2), int(k3), int(k4)]
	elem.append(tmp)

#print elem

INELE.close()

############ Knotenverschiebungen #############
INU = open('feld.txt','r')
NL = []
U = []
for line in INU:
	label, ux, uy = line.strip(' \n').split('\t')
	NL.append(int(label))
	tmp=[float(ux), float(uy)]
	U.append(tmp)

############ In ODB schreiben #################
from odbAccess import *

odbneu = Odb(name='BAMSNEU', analysisTitle='derived data', description='test problem',
    path='BAMSNEU.odb')


part1 = odbneu.Part(name='BLOCK', embeddedSpace=TWO_D_PLANAR, type=DEFORMABLE_BODY)
part1.addNodes(nodeData=knoten)
part1.addElements(elementData=elem, type='S4')

a = odbneu.rootAssembly
instance1 = a.Instance(name='BLOCK-1', object=part1)

s = odbneu.Step(name='BAMSNEU', description='', domain=TIME, timePeriod=1.0)
f = s.Frame(incrementNumber=1,frameValue=0.1, description='')
uField = f.FieldOutput(name='U', description='Displacements', type=VECTOR)
uField.addData(position=NODAL, instance=instance1, labels=NL, data=U)
s.setDefaultDeformedField(uField)

#odbneu = openOdb(path='BAMSNEU.odb')
#odbneu.parts['BLOCK'].addNodes(nodeData=knoten)
#odbneu.parts['BLOCK'].addElements(elementData=elem, type='S4')
#a = odbneu.rootAssembly
#instance1 = a.instances['BLOCK-1']
#s = odbneu.steps['BAMSNEU']
#f = s.frames[-1]
#uField = f.fieldOutputs['U']

odbneu.save()