Exemplo n.º 1
0
def compare(f1, f2):

    s = open(f1)
    table = ft.read(s)
    time = table[:, 0]
    q = table[:, 1:22]
    s2 = open(f2)
    table2 = ft.read(s2)
    time2 = table2[:, 0]
    q2 = table2[:, 1:22]
    return time, time2, q, q2
Exemplo n.º 2
0
def plot_rover(filename):

    s = open(filename)
    table = ft.read(s)
    time = table[:, 0]
    q = table[:, 1:22]
    velo = table[:, 22:25]
    y = table[:, 25:31]
    lamb = table[:, 31]

    plt.subplot(2, 1, 1)
    plt.plot(time, q[:, 0:3], 'o-')
    plt.subplot(2, 1, 2)
    plt.plot(time, y[:, 0:3], 'o-')
    plt.show()
Exemplo n.º 3
0
def analyze_filedata():
    """
    Reads ``f C r`` from the command line, where ``f`` is the name of
    a file, ``C`` and ``r`` are the initial guesses of the parameters
    in the error model E = C*h^r. A 4th command-line argument can be
    the plot title.
    """
    # read filename and initial guess of C and r in error formula E=C*h^r
    f = open(sys.argv[1], 'r')
    C = float(sys.argv[2])
    r = float(sys.argv[3])
    try:
        plot_title = sys.argv[4]
    except:
        plot_title = ''
    from scitools.filetable import read
    data = read(f)
    print data
    OneDiscretizationPrm.analyze(data[:,0], data[:,1],
                                 initial_guess=(C,r),
                                 plot_title=plot_title)
Exemplo n.º 4
0
def analyze_filedata():
    """
    Reads ``f C r`` from the command line, where ``f`` is the name of
    a file, ``C`` and ``r`` are the initial guesses of the parameters
    in the error model E = C*h^r. A 4th command-line argument can be
    the plot title.
    """
    # read filename and initial guess of C and r in error formula E=C*h^r
    f = open(sys.argv[1], 'r')
    C = float(sys.argv[2])
    r = float(sys.argv[3])
    try:
        plot_title = sys.argv[4]
    except:
        plot_title = ''
    from scitools.filetable import read
    data = read(f)
    print data
    OneDiscretizationPrm.analyze(data[:, 0],
                                 data[:, 1],
                                 initial_guess=(C, r),
                                 plot_title=plot_title)
Exemplo n.º 5
0
f=open('mz'+str(bandindex)+'.txt','r')
B=float(f.readline())
C=float(f.readline())
f.close()

f=open('mxy'+str(bandindex)+'.txt','r')
A=float(f.readline())
f.close()

def func(kx,ky,kz,p):
     A, B, C = p
     return A*(kx**2+ky**2)+B*kz**2+C
  

s=open('go7CIISeso.dat','r')
table=ft.read(s)
s.close()

end=int(datsf)

kx = table[0:end, 0]
ky = table[0:end, 1]
kz = table[0:end, 2]

e1=(table[0:end,bandindex]+table[0:end,bandindex-1])/2*ry

r4=sqrt(kx**2+ky**2+(kz_factor*kz)**2)

plot(r4,e1,'r-')
hold('on')	
 
Exemplo n.º 6
0
fileListLambda2=(('d129_4'),('d129_5'),('drag_06'),('d129_7'),('d129_12'))
fileListLambda3=(('d129_5'),('d257_5'))#,('d257_7'))
fileListLambda4=(('drag_06'),('d257_6'))#,('d257_7'))

legendLayer=('layer=0.1', 'layer=0.2','layer=0.3')
legendLambda=('lambda=1e5','lambda=1e6','lambda=1e7','lambda=1e9','lambda=1e11', 'Ref from Folke : 0.6726')
legendLambda2=('lambda=1e4','lambda=1e5','lambda=1e6','lambda=1e7','lambda=1e12','Ref from Folke : 0.6726')
legendLambda3=('lambda=1e5','257 - lambda=1e5','Ref from Folke : 0.6726')
legendLambda3=('lambda=1e6','257 - lambda=1e6','Ref from Folke : 0.6726')
plt.hold('off')
plt.xlabel('time')
plt.hold('on')
plt.ylabel('drag')
plt.axis([0,70,0.3,1])
plt.grid('on')

for filename in fileListLambda3:
	print ("my file is ", filename)
	file=open(filename)
	table=ft.read(file)
	time=table[:,0]
	drag=table[:,1]
	file.close()
	plt.plot(time,drag,'--')
plt.axhline(y=0.6726,xmin=0,xmax=22,color='r')
plt.legend(legendLambda3)
#plt.hold('on')

plt.savefig('DragRe133_CompLambda3.pdf')
plt.show()
Exemplo n.º 7
0
def loaddata(filename):
    
    datafile=open(filename,'r')
    data=ft.read(datafile)
    datafile.close()
    return data