

In [2]:

    
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
#%matplotlib auto

Read and plot physical quantities



In [3]:

    
file_path = "../build/output/"
rho_data = np.loadtxt(file_path+"rho_file.txt");
sigma_data = np.loadtxt(file_path+"sigma_file.txt");
kappa_data = np.loadtxt(file_path+"kappa_file.txt");
data = [rho_data, sigma_data, kappa_data]

emission_data = np.loadtxt(file_path+"emission_file.txt");
nottingham_data = np.loadtxt(file_path+"nottingham_file.txt");



In [4]:

    
for d in data:
    f, (a1, a2) = plt.subplots(1, 2, figsize=(10, 2))
    a1.plot(d[:,0], d[:,1], "-")
    a2.plot(d[:,0], d[:,2], "-")
    plt.show()



In [5]:

    
plt.plot(emission_data[:, 0], emission_data[:, 2])
plt.yscale('log')
plt.xscale('log')
plt.show()



In [6]:

    
plt.plot(nottingham_data[:, 0], nottingham_data[:, 2])
plt.xscale('log')
plt.show()



In [7]:

    
nottingham_data









    Out[7]:





array([[  1.00000000e-02,   5.00000000e+02,   8.69824938e+01],
       [  2.00000000e-02,   5.00000000e+02,   8.26941806e+01],
       [  3.00000000e-02,   5.00000000e+02,   7.95082532e+01],
       ..., 
       [  9.57000000e+00,   5.00000000e+02,   6.96139358e-01],
       [  9.58000000e+00,   5.00000000e+02,   6.95856368e-01],
       [  9.59000000e+00,   5.00000000e+02,   6.95570000e+00]])



In [ ]: