suOlsonP1MUAB = []
suOlsonDiffMUB = []
suOlsonDiffAsym = []
suOlsonWF = []
flux = []
a1 = []
x1, t1, suOlsonWF, a1, flux = dh.WaveFront(
"../data/Temp/Back_1500_WaveFront.txt")
#p1Analytic = GetFromFile("../data/P1AnalyticData.txt")
#EFM,x1 = SuOlsonMyNumericSolution("../data/Temp/SuOlsonEddingtonFactorMinerbo.txt");
#EFM2,x1 = SuOlsonMyNumericSolution("../data/SuOlsonEddingtonFactorMinerbo2.txt");
#suOlsonDiffNumerit,x2 = SuOlsonMyNumericSolution("../data/Temp/diffsuindt1.txt")
#suOlsonDiffMUB,x1 =SuOlsonMyNumericSolution1("../data/Temp/Diff_TH5_DT02.txt")
#suOlsonP1Numerit,x1 = SuOlsonMyNumericSolution("../data/Temp/Diff_TH5_DT01.txt")
#p1Analytic,x1 = SuOlsonMyNumericSolution3("../data/Temp/Diff_TH5_DT001.txt")
x2, suOlsonDiffAsym, a1, t = dh.extract_data("../data/Temp/SuOlsonData.txt",
t0)
x2, y1, a1, t = dh.extract_data("../data/SuOlsonData.txt", t0)
#t0 = 20;
#f = open("dataset.csv","r")
#f1 = open("dataset1.csv","w+")
#data1 = f.readlines()
#f.close()
#f = open("dataset.csv","r")
#X = []
#for i in range(len(data1)):
# data = f.readline()
# x = float(data.split(',')[0])
# y = float(data.split(',')[1])
# X.append(str(x) + "\n" + str(y) + "\n")
#f1.writelines(X)
#f1.close()
tm_2 = []
tm_third = []
EV = 11605
#x1,t1,suOlsonWF, a1, flux = dh.WaveFront("../data/Temp/Back_1500_WaveFront.txt")
#p1Analytic = GetFromFile("../data/P1AnalyticData.txt")
#EFM,x1 = SuOlsonMyNumericSolution("../data/Temp/SuOlsonEddingtonFactorMinerbo.txt");
#EFM2,x1 = SuOlsonMyNumericSolution("../data/SuOlsonEddingtonFactorMinerbo2.txt");
#suOlsonDiffNumerit,x2 = SuOlsonMyNumericSolution("../data/Temp/diffsuindt1.txt")
#suOlsonDiffMUB,x1 =SuOlsonMyNumericSolution1("../data/Temp/Diff_TH5_DT02.txt")
#suOlsonP1Numerit,x1 = SuOlsonMyNumericSolution("../data/Temp/Diff_TH5_DT01.txt")
#p1Analytic,x1 = SuOlsonMyNumericSolution3("../data/Temp/Diff_TH5_DT001.txt")
#x2 ,y_half,a1,t = dh.extract_data("../data/SuOlsonData.txt", 0.5)
#x2 ,y_one,a1,t = dh.extract_data("../data/SuOlsonData.txt", 1)
#x2 ,y_two,a1,t = dh.extract_data("../data/SuOlsonData.txt", 2)
#x2 ,y_three,a1,t = dh.extract_data("../data/SuOlsonData.txt", t0)
x2, y_third, a1, t = dh.extract_data("../data/SuOlsonData.txt", 0.25)
x2, y_2, a1, t = dh.extract_data("../data/SuOlsonData.txt", 0.01)
#x1 ,tm_half,a1,t = dh.extract_data("../data/Temp/SuOlsonData.txt", 0.5)
#x1 ,tm_one,a1,t = dh.extract_data("../data/Temp/SuOlsonData.txt", 1)
#x2 ,tm_two,a1,t = dh.extract_data("../data/Temp/SuOlsonData.txt", 2)
#x2 ,tm_three,a1,t = dh.extract_data("../data/Temp/SuOlsonData.txt", t0)
x1, tm_third, a1, t = dh.extract_data("../data/Temp/SuOlsonData.txt", 0.25)
x1, tm_2, a1, t = dh.extract_data("../data/Temp/SuOlsonData.txt", 0.01)
#plt.xlim(0,2)
for i in range(0, 2000):
# x2[i] = x2[i] * 10
# y_half[i] = y_half[i]/EV
#y_one[i] = y_one[i]/EV
# y_two[i] = y_two[i]/EV
suOlsonP1MUAB = []
suOlsonDiffMUB = []
suOlsonDiffAsym = []
suOlsonWF = []
flux = []
a1 = []
#x1,t1,suOlsonWF, a1, flux = WaveFront("../data/Temp/Back_1500_WaveFront.txt")
#p1Analytic = GetFromFile("../data/P1AnalyticData.txt")
#EFM,x1 = SuOlsonMyNumericSolution("../data/Temp/SuOlsonEddingtonFactorMinerbo.txt");
#EFM2,x1 = SuOlsonMyNumericSolution("../data/SuOlsonEddingtonFactorMinerbo2.txt");
#suOlsonDiffNumerit,x2 = SuOlsonMyNumericSolution("../data/Temp/diffsuindt1.txt")
#suOlsonDiffMUB,x1 =SuOlsonMyNumericSolution1("../data/Temp/Diff_TH5_DT02.txt")
#suOlsonP1Numerit,x1 = SuOlsonMyNumericSolution("../data/Temp/Diff_TH5_DT01.txt")
#p1Analytic,x1 = SuOlsonMyNumericSolution3("../data/Temp/Diff_TH5_DT001.txt")
plt.xlim(0, 2)
x2, y1, y2, t1, t2 = dh.extract_data("../data/SuOlsonData.txt", t0)
line16, = plt.plot(x2, y1[0:2000], 'g', label="weff")
#legend
#x = [0.5,1,2,3,4]
#y = [0,0.001,0.01,0.1,1,3]
#plt.yscale('log')
#plt.xscale('log')
#plt.ylim(0,1)
#plt.axis([0.3,8,0.001,3]);
#plt.xscale('log')
#ticks = [0.3,0.5,1,3,4];
#ticks2 = [0.001,0.01,0.1,1]
#labels = ['0.3','0.5','1','3','4']
#labels2 = ['0.001','0.01','0.1','1']
#plt.xticks(ticks,labels);
#plt.yticks(ticks2,labels2)
suOlsonP1MUAB = []
suOlsonDiffMUB = []
suOlsonDiffAsym = []
suOlsonWF = []
flux = []
a1 = []
#x1,t1,suOlsonWF, a1, flux = dh.WaveFront("../data/Temp/Back_1500_WaveFront.txt")
#p1Analytic = GetFromFile("../data/P1AnalyticData.txt")
#EFM,x1 = SuOlsonMyNumericSolution("../data/Temp/SuOlsonEddingtonFactorMinerbo.txt");
#EFM2,x1 = SuOlsonMyNumericSolution("../data/SuOlsonEddingtonFactorMinerbo2.txt");
#suOlsonDiffNumerit,x2 = SuOlsonMyNumericSolution("../data/Temp/diffsuindt1.txt")
#suOlsonDiffMUB,x1 =SuOlsonMyNumericSolution1("../data/Temp/Diff_TH5_DT02.txt")
#suOlsonP1Numerit,x1 = SuOlsonMyNumericSolution("../data/Temp/Diff_TH5_DT01.txt")
#p1Analytic,x1 = SuOlsonMyNumericSolution3("../data/Temp/Diff_TH5_DT001.txt")
#x2, suOlsonDiffAsym, a1, t = dh.extract_data("../data/Temp/SuOlsonData.txt", t0)
x2, y1, y2, t1, t2 = dh.extract_data("../data/OlsonEnergy.txt", t0)
#t0 = 20;
#f = open("dataset.csv","r")
#f1 = open("dataset1.csv","w+")
#data1 = f.readlines()
#f.close()
#f = open("dataset.csv","r")
#X = []
#for i in range(len(data1)):
# data = f.readline()
# x = float(data.split(',')[0])
# y = float(data.split(',')[1])
# X.append(str(x) + "\n" + str(y) + "\n")
#f1.writelines(X)
#f1.close()
#f.close()
a1 = []
y_half = []
y_one = []
y_two = []
y_three = []
y_third = []
y_2 = []
y_1 = []
tm_half = []
tm_one = []
tm_two = []
tm_three = []
tm_2 = []
tm_third = []
EV = 11605
x1, t1, suOlsonWF = dh.WaveFront("../data/Temp/Back_1500_WaveFront.txt")
x1, y1, y2, t1, t2 = dh.extract_data("../data/BackTemp.txt", t0)
x3, rad_1, rad_2, t1, t2 = dh.extract_data("../data/BackEnergy.txt", t0)
#plt.xlim(0,2)
for i in range(0, 1500):
rad_1[i] = rad_1[i] / EV
x1[i] = x1[i] * 10
line16, = plt.plot(x1[0:1500], y1[0:1500])
line1, = plt.plot(x1[0:1500], rad_1[0:1500])
plt.title("For t = " + str(t1))
#plt.ylabel('Radation tempreture Density - T(x,t)');
#plt.xlabel('x');
plt.show()