Skip to content
Snippets Groups Projects
Commit 0c483e40 authored by Yngve Levinsen's avatar Yngve Levinsen
Browse files

fixes in ess.TraceWin.field_map

parent 165f9722
No related branches found
No related tags found
No related merge requests found
...@@ -952,16 +952,27 @@ class field_map: ...@@ -952,16 +952,27 @@ class field_map:
raise ValueError("Cannot find file {}".format(filename)) raise ValueError("Cannot find file {}".format(filename))
fin = file(filename,'r') fin = file(filename,'r')
l=fin.readline().split() l=fin.readline().split()
self.size=[] self.start=[]
self.end=[]
numindexes=[] numindexes=[]
while len(l)==2: while len(l)>1:
numindexes.append(int(l[0])+1) numindexes.append(int(l[0])+1)
self.size.append(float(l[1])) if len(l)==2:
self.start.append(0.0)
self.end.append(float(l[1]))
else:
self.start.append(float(l[1]))
self.end.append(float(l[2]))
l=fin.readline().split() l=fin.readline().split()
self.z=numpy.arange(self.start[0],self.end[0],(self.end[0]-self.start[0])/numindexes[0])
if len(self.start)>1:
self.y=numpy.arange(self.start[1],self.end[1],(self.end[1]-self.start[1])/numindexes[1])
if len(self.start)>2:
self.x=numpy.arange(self.start[2],self.end[2],(self.end[2]-self.start[2])/(numindexes[2]))
self.norm=float(l[0]) self.norm=float(l[0])
self.map=numpy.loadtxt(fin).reshape(numindexes) self.map=numpy.loadtxt(fin).reshape(numindexes)
print self.map.shape
def savemap(self, filename): def savemap(self, filename):
fout=file(filename,'w') fout=file(filename,'w')
...@@ -974,3 +985,4 @@ class field_map: ...@@ -974,3 +985,4 @@ class field_map:
data=self.map.reshape(l) data=self.map.reshape(l)
for j in data: for j in data:
fout.write('{}\n'.format(j)) fout.write('{}\n'.format(j))
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment