In [1]:

    

%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
from array_io import *


    

In [2]:

    

fdata = "sedov_taylor.txt"
r, rho, p, u = read_four_arrays(fdata)


    

In [6]:

    

plt.plot(r/r.max(),rho/rho.max())
plt.plot(r/r.max(),p/p.max())
plt.plot(r/r.max(),u/u.max())
plt.xlim([0,1])
plt.ylim([0,1])


    

    
Out[6]:





(0, 1)






    

In [ ]: