Commit 1c9ce165 authored by Paola Bacigaluppi's avatar Paola Bacigaluppi
Browse files

deleted some obsolete files

parent 05eaa873
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.axes_grid1.inset_locator import zoomed_inset_axes
from mpl_toolkits.axes_grid1.inset_locator import mark_inset
testnr=3
conv=3
orders=["3"]
styles = ['--^', '-s', ':o']
msize=[7,6,7]
meve=[5,10,15]
lwd=[1,1,2.5]
vars=["rho","v","p" ]#, "Momentum", "Energy"]
#limits=[[[0.45,0.65], [0.37, 0.47]]] Sod first area
#limits=[[[0.63,0.83], [0.1, 0.3]]] #sod second area
limits=[[[0.,1], [-0.01, 1.05]]]#, [[0.,1.],[-0.05,1.4]], [[0.,1.],[0.07,1.1]]]
mesh=["100"]
#or N in [int(2**k) for k in range()]:
sols=[]
det_sol=[]
flux_sol=[]
#dashes = [10, 5, 100, 5]
fig, ax = plt.subplots()
for ord, ordname in enumerate(orders):
sols.append(np.loadtxt("Test"+str(testnr)+"/Test_1D_Sod/B"+ordname+"/sol_prim_1.dat"))
flux_sol.append(np.loadtxt("Test"+str(testnr)+"/Test_1D_Sod/B"+ordname+"/FluxIndicator_100_1.dat"))
det_sol.append(np.loadtxt("Test"+str(testnr)+"/Test_1D_Sod/B"+ordname+"/DiagIndicator_100_1.dat"))
ax.plot(flux_sol[ord][:,0],(flux_sol[ord][:,1]*0.1+0.1), 'x',markersize=4,label="Flux Indicator")
ax.plot(sols[ord][:,0],sols[ord][:,1],'r',linewidth=2, label="Density")
#axins.plot(sols[ord][:,0],sols[ord][:,var+1])
ax.set_xlim(limits[0][0])
ax.set_ylim(limits[0][1])
ax.legend(loc=1)
ax.legend(frameon=False)
#mark_inset(ax, axins, loc1=2, loc2=4, fc="none", ec="0.5")
plt.savefig("Test_1D_Sod_detectionflux_B3_rho_100_1st.pdf")
# plt.show()
plt.clf()
plt.close
fig, ax = plt.subplots()
ax.plot(det_sol[ord][:,0],det_sol[ord][:,1]*0.1, 'x',markersize=4,label="Criteria Indicator")
ax.plot(sols[ord][:,0],sols[ord][:,1],'r',linewidth=2, label="Density")
#axins.plot(sols[ord][:,0],sols[ord][:,var+1])
ax.set_xlim(limits[0][0])
ax.set_ylim(limits[0][1])
ax.legend(loc=1)
ax.legend(frameon=False)
#mark_inset(ax, axins, loc1=2, loc2=4, fc="none", ec="0.5")
plt.savefig("Test_1D_Sod_detectionCriteria_B3_rho_100_1st.pdf")
# plt.show()
plt.clf()
plt.close
limits=[[[0.,1], [-0.01, 1.05]]]#, [[0.,1.],[-0.05,1.4]], [[0.,1.],[0.07,1.1]]]
mesh=["100"]
#or N in [int(2**k) for k in range()]:
sols=[]
det_sol=[]
flux_sol=[]
#dashes = [10, 5, 100, 5]
fig, ax = plt.subplots()
for ord, ordname in enumerate(orders):
sols.append(np.loadtxt("Test"+str(testnr)+"/Test_1D_Sod/B"+ordname+"/sol_prim_11.dat"))
flux_sol.append(np.loadtxt("Test"+str(testnr)+"/Test_1D_Sod/B"+ordname+"/FluxIndicator_100_11.dat"))
det_sol.append(np.loadtxt("Test"+str(testnr)+"/Test_1D_Sod/B"+ordname+"/DiagIndicator_100_11.dat"))
ax.plot(flux_sol[ord][:,0],(flux_sol[ord][:,1]*0.1+0.1), 'x',markersize=4,label="Flux Indicator")
ax.plot(sols[ord][:,0],sols[ord][:,1],'r',linewidth=2, label="Density")
#axins.plot(sols[ord][:,0],sols[ord][:,var+1])
ax.set_xlim(limits[0][0])
ax.set_ylim(limits[0][1])
ax.legend(loc=1)
ax.legend(frameon=False)
#mark_inset(ax, axins, loc1=2, loc2=4, fc="none", ec="0.5")
plt.savefig("Test_1D_Sod_detectionflux_B3_rho_100_10th.pdf")
# plt.show()
plt.clf()
plt.close
fig, ax = plt.subplots()
ax.plot(det_sol[ord][:,0],det_sol[ord][:,1]*0.1, 'x',markersize=4,label="Criteria Indicator")
ax.plot(sols[ord][:,0],sols[ord][:,1],'r',linewidth=2, label="Density")
#axins.plot(sols[ord][:,0],sols[ord][:,var+1])
ax.set_xlim(limits[0][0])
ax.set_ylim(limits[0][1])
ax.legend(loc=1)
ax.legend(frameon=False)
#mark_inset(ax, axins, loc1=2, loc2=4, fc="none", ec="0.5")
plt.savefig("Test_1D_Sod_detectionCriteria_B3_rho_100_10th.pdf")
# plt.show()
plt.clf()
plt.close
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment