ionnames = ["He2+", "O6+", "O7+"]
parkerwinkel = pi * 0.2
marks = [
] #15.35,29.38,58.2,91.,117.7,127.6,138.5,144.,165.,172.3,180.8,192.,207.5,210.7,219.,222.5,227.,243.8,249.7,257.9,263.5,271.,298.5,324.5,345.,351.3]
mres = 240
pres = 720
sires = 3600
bintime = 4 * 3600
year = 2007
timeframe = [mark - 4, mark + 4]
#Multiplotklasse initialisieren
mp = multiplot(6, 0)
mp.setgv(0)
mp.setxrange(timeframe[0], timeframe[1])
mpi = 0
#vertikale Marks setzen
for i in range(len(marks)):
mp.addplotmark("vline")
mp.plotmarks[i].setvlinepos(marks[i])
# Anfang: Magdaten auswerten
mdata = magdata(mres, year, timeframe[0], timeframe[1], "/data/ivar/ace/mag/")
mdata.load()
binres = bintime / mres
mbin = []
for i in range(len(mdata.magbx)):
# Anfang: Winkel ausrechnen
jahr = []
for i in range(1):
jahr.append([])
#jahr[0] = [2004,[[1,30]]]#120],[150,360]]]
#jahr[1] = [2005,[[1,30]]]#36],[64,103],[105,122],[135,191],[193,230]]]
#jahr[2] = [2006,[[1,30]]]#345],[347,360]]]
jahr[0] = [2007, [[1, 360]]]
gp = Gnuplot()
ionnames = ["He1+"] #,"O6+","O7+","O8+"]
intervalrange = 20
intervalrange = 0.1
maxrange = 5
mp = multiplot(len(ionnames))
mp.setxrange(0, maxrange)
mp.setxlabel("vel/swvel")
datamin = []
datamax = []
for i in range(len(ionnames)):
datamin.append(1e30)
datamax.append(0)
for jh in range(len(jahr)):
print "Jahr " + str(jahr[jh][0])
pathh1 = "/data/wimmer/ace/puis/1hdata/" + str(jahr[jh][0]) + "/"
pathh2 = "/data/wimmer/ace/1hdata/" + str(jahr[jh][0]) + "/"
data5 = zeros([100, 100], float)
for i in range(100):
data1[i][0] = float(i)
data2[i][0] = float(i)
data3[i][0] = float(i)
data4[i][0] = float(i)
data5[i][0] = float(i)
data1[i][1] = float(i)
data2[i][1] = float(i)**2
data3[i][1] = sqrt(float(i))
data4[i][1] = float(i)**3
data5[i][1] = float(i)**3
panels = 4
mp = multiplot(panels)
mp.addplotmark("vline")
mp.addplotmark("vline")
#mp.addplotmark("box")
#mp.addplotmark("box")
mp.plotmarks[0].setvlinepos(10.)
mp.plotmarks[1].setvlinepos(15.)
mp.plotmarks[1].setpanel([0, 2])
#mp.plotmarks[2].setpanel([-1])
#mp.plotmarks[2].setboxrange(0.,10.)
#mp.plotmarks[3].setpanel([-1])
#mp.plotmarks[3].setboxrange(15.,40.)
for panel in range(panels):
mp.panel[panel].setylabel("Test")
if (panel == 1):
#temp1 = [[20,164.6],[25,172.3],[30,178.4],[35,183.1],[40,186.6],[45,189.2],[50,191.1]]
#temp2 = [[20,173.3],[25,181.5],[30,188],[35,195],[40,199],[45,202],[50,205]]
#quot = []
#for i in range(len(temp2)):
# quot.append([temp[i][0],temp[i][1]/temp2[i][1]])
#print quot
pardiff = []
for i in range(len(temp2)):
pardiff.append([temp1[i][0], sqrt(abs(temp1[i][1]**2 - temp2[i][1]**2))])
print pardiff
mp = multiplot(2)
pdata = Data(temp1)
mp.panel[0].adddata(pdata)
pdata = Data(temp2)
mp.panel[0].adddata(pdata)
#pdata = Data(quot)
#mp.panel[0].adddata(pdata)
mp.panel[0].setyrange(140.0, 190.)
mp.panel[1].setlogscale(0)
mp.panel[1].setylabel("Abs")
pdata = Data(pardiff)
mp.panel[1].adddata(pdata)
mp.panel[1].setyrange(20., 70.)
mp.panel[1].setlogscale(0)
phistarr[26].set_option_colonsep("title", "'Si^{9+}'")
phistarr[26].set_option_colonsep("with", "steps lt 5 lw 2")
phistarr2[-1].set_option_colonsep("with", "steps lt 1 lw 2")
phistarr2[36].set_option_colonsep("with", "steps lt 3 lw 2")
phistarr2[9].set_option_colonsep("with", "steps lt 4 lw 2")
phistarr2[26].set_option_colonsep("with", "steps lt 5 lw 2")
pmeandv[-1].set_option_colonsep("with", "l lt 1 lw 2")
pmeandv[36].set_option_colonsep("with", "l lt 3 lw 2")
pmeandv[9].set_option_colonsep("with", "l lt 4 lw 2")
pmeandv[26].set_option_colonsep("with", "l lt 5 lw 2")
pmeandv2[-1].set_option_colonsep("with", "l lt 1 lw 2")
pmeandv2[36].set_option_colonsep("with", "l lt 3 lw 2")
pmeandv2[9].set_option_colonsep("with", "l lt 4 lw 2")
pmeandv2[26].set_option_colonsep("with", "l lt 5 lw 2")
mp = multiplot(2, 1)
mp.setname("paperdvhistfinal")
mp.setxlabel("|v_{ip}|/C_{A}")
mp.setxrange(0, 2)
mp.panel[1].setylabel("# occurences")
mp.panel[1].setyrange(0, 180)
mp.panel[1].key = "set key"
mp.panel[1].adddata(phistarr[-1])
mp.panel[1].adddata(phistarr[36])
mp.panel[1].adddata(phistarr[26])
mp.panel[1].adddata(phistarr[9])
mp.panel[1].adddata(pmeandv[-1])
mp.panel[1].adddata(pmeandv[36])
mp.panel[1].adddata(pmeandv[26])
mp.panel[1].adddata(pmeandv[9])
mp.panel[0].setylabel("# occurences")
break
integral_sw = integral(data.densspec[takt], sw_ende, sw_anfang)
integral_spth = integral(data.densspec[takt], 0, sw_ende)
if (integral_sw > 0):
vergleich2.append(integral_spth / integral_sw)
else:
vergleich2.append(0)
#print sw_anfang
#print sw_ende
#gp("set style data lines")
#gp.plot(vergleich)
mp = multiplot(3)
mp.setxrange(timeframe[0][0], timeframe[-1][1])
#marks = [data.densspec[takt][sw_anfang][0],data.densspec[takt][sw_ende][0], data.densspec[takt][maximum][0]]
#for i in range(len(marks)):
# mp.addplotmark("vline")
# mp.plotmarks[i].setvlinepos(marks[i])
mp.panel[0].adddata(Data(zip(data.time, data.vel)))
mp.panel[0].setyrange(0., 1000.)
mp.panel[0].setlogscale(0)
mp.panel[0].setylabel("Geschwindigkeit")
mp.panel[1].adddata(Data(zip(data.time, vergleich)))
mp.panel[1].setyrange(0., 0.04)
mp.panel[1].setlogscale(0)
vyb2 = tan(syncatt.phi) * (hdata.vel + valf * 0.5) - (
hdata.vely + valf * 0.5 * tan(tmpmagphi))
vzb2 = tan(syncatt.theta) * (hdata.vel + valf * 0.5) - (
hdata.velz + valf * 0.5 * tan(magdata.theta))
dvb2 = sqrt(vyb2**2 + vzb2**2)
counts = sum(hedata.countspec[:, :, 1], axis=1)
#vspec=[]
#for i in range(len(hedata.time)):
# vd=vdetector(func=twobimaxmag,para=[hdata.vel[i]+0.2*valf[i],0.,0.,20.,30.,30.,0.03,0.8,0.8*valf[i],tmpmagphi[i],magdata.theta[i]],aspphi=-syncatt.phi[i],asptheta=syncatt.theta[i])
#vd=vdetector(func=twobimaxmag,para=[hdata.velx[i],hdata.vely[i],hdata.velz[i],20.,30.,30.,0.02,0.5,valf[i],tmpmagphi[i],magdata.theta[i]],aspphi=-syncatt.phi[i],asptheta=syncatt.theta[i])
# vspec.append(array(vd.vspec))
mp = multiplot(4)
mp.setgv(0)
mp.setname("sp_%3.i_%3.i" % (int(timeframe[0, 0]), int(timeframe[-1, 1])))
mp.setxrange(timeframe[0, 0], timeframe[-1, 1])
mp.panel[0].adddata(Data(magdata.time, thetaarr, with_="l lt 1"))
mp.panel[0].adddata(
Data(array([[timeframe[0, 0], pi / 2.], [timeframe[-1, 1], pi / 2.]]),
with_="l lt 0"))
mp.panel[0].setyrange(0, pi)
mp.panel[0].setytics(
"('0' 0,'' %f 1,'' %f 1,'' %f 1,'{/Symbol p}/4' %f,'' %f 1,'' %f 1,'' %f 1,'{/Symbol p}/2' %f,'' %f 1,'' %f 1,'' %f 1,'{/Symbol p}3/4' %f,'' %f 1,'' %f 1,'' %f 1,'{/Symbol p}' %f)"
%
(pi / 16., pi / 8., pi * 3 / 16., pi / 4., pi * 5. / 16., pi * 6 / 16.,
pi * 7 / 16., pi / 2., pi * 9. / 16., pi * 10. / 16., pi * 11. / 16.,
pi * 3. / 4., pi * 13. / 16., pi * 14. / 16., pi * 15. / 16., pi))
mp.panel[1].adddata(Data(magdata.time, magdata.phi, with_="l lt 1"))
[224.5, 227.2, 0, 0], [242.3, 244.3, 0, 0],
[296., 297.6, 298.3, 298.5]]
ionnames = ["He2+", "O6+", "O7+"]
parkerwinkel = pi * 0.2
mres = 240
pres = 720
sires = 3600
bintime = 4 * 3600
year = 2007
#Multiplotklasse initialisieren
mp = multiplot(7, 1)
mp.setgv(1)
mp.setxlabel("Tage relativ zum Ereignis")
mp.setxrange(-4, +4)
mpi = 0
# parameteuebergabe zum plotten eines einzelnen ereignisses
if (len(sys.argv) > 1):
del mark
mark = [float(sys.argv[1])]
mp.setgv(0)
mp.settitle("Ereignis " + str(mark[0]))
#dips einzeichnen
for i in range(len(event)):
mp.addplotmark("vline")
for p in range(len(hist)):
for i in range(len(data[0][w])):
bin = int(data[0][w][i][1]/binsize)
sum[p,bin] += (data[p+1][w][i][1]-hist[p,w,bin,1])**2
for i in range(len(sum[p])):
hist[p,w,i,2] = sqrt(sum[p,i] / ((hist[p,w,i,3])**2))
# Ende: Standartabweichung bestimmen
# Anfang: Geschwindigkeiten anfügen
for p in range(len(hist)):
hist[p,w,:,0] = linspace(0,(bincount-1)*binsize,bincount) + binsize/2
# Ende: Geschwindigkeiten anfügen
mp = multiplot(plotcount+1,1)
mpi = 0
mp.setgv(1)
mp.setxlabel("Geschwindigkeitsbins")
mp.setxrange(200,bincount*binsize)
#Anfang: Panels plotten
for p in range(len(hist)):
for w in range(len(marks)):
plotdata = Data(hist[p,w,:,0:3])
plotdata.set_option_colonsep("title","'"+title[p]+" "+marktitle[w]+"'")
plotdata.set_option_colonsep("with","yerrorline")
mp.panel[mpi].adddata(plotdata)
mp.panel[mpi].setyrange(ranges[p][0],ranges[p][1])
mp.panel[mpi].setlogscale(log[p])
del plotdata
from Gnuplot import Gnuplot, Data
from libacepy.multiplot import multiplot
from math import sqrt, atan, pi, acos, cos
parkerwinkel = pi * 0.2
mres = 240
pres = 720
sires = 3600
bintime = 6 * 3600
year = 2007
timeframe = [220, 365]
#Multiplotklasse initialisieren
mp = multiplot(4, 0)
mp.setgv(1)
mp.setxrange(timeframe[0], timeframe[1])
mpi = 0
# Anfang: Magdaten auswerten
mdata = magdata(mres, year, timeframe[0], timeframe[1], "/data/ivar/ace/mag/")
mdata.load()
binres = bintime / mres
mbin = []
for i in range(len(mdata.magbx)):
# Anfang: Winkel ausrechnen
if (mdata.magbx[i][1] < 0):
phi = atan(-mdata.magby[i][1] / mdata.magbx[i][1])
elif (mdata.magbx[i][1] == 0):
if (mdata.magby[i][1] > 0):