def get_fortrantextspectra(ndur,nloc): # 生のテキストデータを読み込む import numpy as np import spectra filename,N=spectra.get_fn_and_N(ndur,nloc) data = np.loadtxt(filename+'_spectra.dat',delimiter=',') omega=data[:,0] power=data[:,1] return [omega,power]
def get_textdata(ndur,nloc): # 生のテキストデータを読み込む import numpy as np import spectra filename,N=spectra.get_fn_and_N(ndur,nloc) data = np.loadtxt(filename+'.dat',delimiter=' ') t=data[:,0] f=data[:,1] f[np.where(f==-0.999)]=np.average(f) return [t,f]
def get_fortranspectra(ndur,nloc): import numpy as np import spectra filename,N=spectra.get_fn_and_N(ndur,nloc) format_of_data="<"+str(N/2+1)+"f" dt=np.dtype([("data",format_of_data)]) fd=open(filename+"_spectra.out","r") chunk=np.fromfile(fd,dtype=dt,count=2) omega=chunk[0]['data'] power=chunk[1]['data'] return [omega,power]
def get_data(ndur,nloc): # fortranで変換したバイナリデータを読み込む import numpy as np import spectra filename,N=spectra.get_fn_and_N(ndur,nloc) format_of_data="<"+str(N)+"f" dt=np.dtype([("data",format_of_data)]) fd=open(filename+".out","r") chunk=np.fromfile(fd,dtype=dt,count=2) tt=chunk[0]['data'] ff=chunk[1]['data'] return [tt,ff]
def get_test(ndur): import numpy as np import spectra nloc=1 # 何でもいい hoge,N=spectra.get_fn_and_N(ndur,nloc) format_of_data="<"+str(N)+"f" dt=np.dtype([("data",format_of_data)]) fd=open("test.out","r") chunk=np.fromfile(fd,dtype=dt,count=2) test1=chunk[0]['data'] test2=chunk[1]['data'] return [test1,test2]