def load (self): if (self.datatype == "element"): if (self.source == "swics"): self.data = swicselemdata (self.name,self.timeframe,self.path) self.data.load() elif (self.datatype == "ion"): if (self.source == "swics"): self.data = swicsiondata (self.name,self.timeframe,self.path) self.data.load() if (self.source == "swepam"): self.data = swepamdata (self.res,self.year,self.timeframe[0][0],self.timeframe[0][1],self.path) self.data.load()
def load_swics(self):
from libacepy.swicsiondata import swicsiondata
self.iondata = swicsiondata(self.ionname, self.timeframe,
self.swicspath)
self.iondata.load()
from libacepy.swicsiondata import swicsiondata
from libacepy.swepamdata import swepamdata, sync_h
he2001 = swicsiondata("He2+", [[1., 365.]], "/data/wimmer/ace/1hdata/2001/")
he2004 = swicsiondata("He2+", [[1., 135.], [141., 365.]],
"/data/wimmer/ace/1hdata/2004/")
he2007 = swicsiondata("He2+", [[1., 365.]], "/data/wimmer/ace/1hdata/2007/")
sphe2001 = swepamdata(720., 2001, 1., 365., "/data/ivar/ace/swepam/")
sphe2007 = swepamdata(720., 2007, 1., 365., "/data/ivar/ace/swepam/")
sphe2001.load()
sphe2007.load()
he2001.load()
he2004.load()
he2007.load()
ssphe2001 = sync_h(sphe2001, he2001)
ssphe2007 = sync_h(sphe2007, he2007)
he3002001 = 0.
he6002001 = 0.
he3002004 = 0.
he6002004 = 0.
he3002007 = 0.
he6002007 = 0.
sphe3002001 = 0.
sphe6002001 = 0.
sphe3002007 = 0.
elif (diffhe2o6[i][1] < low and diffhe2o6[i][1] > 1e-20):
low = diffhe2o6[i][1]
plotdata = Data(diffhe2o6)
plotdata.set_option_colonsep("title","'V_{O6} - V_{He2}'")
mp.panel[mpi].adddata(plotdata)
mp.panel[mpi].setyrange(low * 0.7,high* 1.3)
mp.panel[mpi].setlogscale(1)
mp.panel[mpi].setylabel("Vel")
mpi += 1
del plotdata
# Ende: Swicsdaten auswerten
"""
# Anfang: Swics vs Swepam
swics = swicsiondata("He2+", [timeframe],
"/data/wimmer/ace/1hdata/" + str(year) + "/")
swics.load()
swepam = swepamdata(720, year, timeframe[0], timeframe[1],
"/data/ivar/ace/swepam/")
swepam.load()
v_heh = []
dens_heh = []
binres = bintime / 720
vel_low = 1e10
vel_high = 0.
dense_low = 1e10
dense_high = 0.
for i in range(len(swics.time)):
j = i
continue
while (len(blub) <= x_index):
blub.append([])
while (len(blub[x_index]) <= y_index):
blub[x_index].append(0)
blub[x_index][y_index] += 1
return blub
hist = []
x_max = []
y_max= []
for m in range(len(mark)):
print "Mark "+str(m+1)+" von "+str(len(mark))
timeframe = [mark[m][0],mark[m][1]]
he2data = swicsiondata("He2+",[timeframe],"/data/wimmer/ace/1hdata/"+str(year)+"/")
he2data.load()
o6data = swicsiondata("O6+",[timeframe],"/data/wimmer/ace/1hdata/"+str(year)+"/")
o6data.load()
o7data = swicsiondata("O7+",[timeframe],"/data/wimmer/ace/1hdata/"+str(year)+"/")
o7data.load()
pdata = swepamdata(pres,year,timeframe[0],timeframe[1],"/data/ivar/ace/swepam/")
pdata.load()
pdata = sync_h(pdata,he2data)
odata = swicselemdata("O",[timeframe],"/data/wimmer/ace/1hdata/"+str(year)+"/")
odata.load()
fedata = swicselemdata("Fe",[timeframe],"/data/wimmer/ace/1hdata/"+str(year)+"/")
fedata.load()
he2_h = array(densratio(he2data,pdata))
fe_o = array(densratio(fedata,odata))
from libacepy.swicsiondata import swicsiondata
path = "/data/wimmer/ace/1hdata/2007/"
timeframe = [[1., 200.]]
c6 = swicsiondata("C6+", timeframe, path)
c6.load()
n7 = swicsiondata("N7+", timeframe, path)
n7.load()
o8 = swicsiondata("O8+", timeframe, path)
o8.load()
he2 = swicsiondata("He2+", timeframe, path)
he2.load()
co = 0.
con = 0
cn = 0.
cnn = 0
for i in range(len(he2.time)):
if (he2.vel[i] > 600. and he2.vel[i] < 700.):
if (c6.dens[i][0] > 0. and o8.dens[i][0] > 0.):
co += o8.dens[i][0] / c6.dens[i][0]
con += 1
if (c6.dens[i][0] > 0. and n7.dens[i][0] > 0.):
cn += n7.dens[i][0] / c6.dens[i][0]
cnn += 1
co = co / con
cn = cn / cnn
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]) + "/"
datadenssum = []
i = 0
data = swicsprotondata("He1+", jahr[jh][1], pathh1)
data.load()
datah2 = swicsiondata("He2+", jahr[jh][1], pathh2)
datah2.load()
print data
datadenssum.append([])
for swstart in range(300, 800, 100):
print(maxrange / intervalrange)
for vs in range(int(maxrange / intervalrange)):
datadenssum[i].append(0)
for ts in range(len(data.countspec)):
for vs in range(len(data.countspec[ts])):
if (datah2.vel[ts] > swstart
and datah2.vel[ts] < swstart + 100):
interval = int(data.countspec[ts][vs][0] / datah2.vel[ts] /
intervalrange)
#print interval
from libacepy.swicsiondata import swicsiondata
from Gnuplot import Gnuplot, Data
from libacepy.multiplot import multiplot
gp = Gnuplot()
path = "/data/wimmer/ace/1hdata/2007/"
timeframe = [[330, 336]]
ionnames = ["He2+", "O6+", "O7+"]
data = []
datatemp = []
datavel = []
datadens = []
for i in range(len(ionnames)):
data.append(swicsiondata(ionnames[i], timeframe, path))
data[i].load()
datatemp.append([])
datavel.append([])
datadens.append([])
for j in range(len(data[i].time)):
datatemp[i].append([data[i].time[j], data[i].temp[j][0]])
datavel[i].append([data[i].time[j], data[i].vel[j]])
datadens[i].append([data[i].time[j], data[i].dens[j][0]])
if (ionnames.index("O6+") != ValueError
and ionnames.index("O7+") != ValueError):
relationo6o7 = []
for i in range(len(data[ionnames.index("O7+")].time)):
if (data[ionnames.index("O7+")].dens[i][0] > 0.):
relationo6o7.append([
data[ionnames.index("O7+")].time[i],
meanbz = []
meanax = []
meanaz = []
#iondata=swicsiondata(ionlist[0],resswics,year,timeframe,datapath)
#mag.sync_mag(iondata)
#swepam.sync_h(iondata)
#syncatt=sync_attitude(att,iondata)
polarity = zeros([len(timeframe)])
uppervel = 700.
lowervel = 600.
linetype = 0
for ion in ionlist:
iondata = swicsiondata(ion, resswics, year, timeframe, datapath)
#if (ion!="He2+"):
# hedata=swicsiondata("He2+",resswics,year,timeframe,datapath)
#else:
# hedata=iondata
#iondata.load()
swepam = swepamdata(resswepam, year, timeframe, swepampath)
mag = magdata(resmag, year, timeframe, magpath)
att = attdata(mag.year, mag.path)
att.load()
mag.sync_mag(iondata)
swepam.sync_h(iondata)
syncatt = sync_attitude(att, iondata)
dblistin = open("sigb_2007.dat")
dblistin.readline()
datadenssum = []
pdata = []
for i in range(len(ionnames)):
datamin.append(1e10)
datamax.append(0)
for jh in range(len(jahr)):
print "jahr " + str(jahr[jh][0])
path = "/data/wimmer/ace/1hdata/" + str(jahr[jh][0]) + "/"
data.append([])
datadenssum.append([])
pdata.append([])
for i in range(len(ionnames)):
data[jh].append(swicsiondata(ionnames[i], jahr[jh][1], path))
data[jh][i].load()
datadenssum[jh].append([])
for j in range(int(1500 / intervalrange)):
datadenssum[jh][i].append(0)
for j in range(len(data[jh][i].time)):
interval = int(data[jh][i].vel[j] / intervalrange)
if (interval < int(1500 / intervalrange)):
datadenssum[jh][i][interval] += data[jh][i].temp[j][0]
# print "Fuer :" + ionnames[i]
for j in range(len(datadenssum[jh][i])):
datadenssum[jh][i][j] /= len(data[jh][i].time)
if (datadenssum[jh][i][j] > 1.e-24):
summe_rechteck += breite * liste[i][1]
# print "Summe (rechteck):"+str(summe_rechteck)
return ((summe_linfit + summe_rechteck) / 2)
#Ende der Integralfunktion
swbreite = 1.6
gp = Gnuplot()
path = "/data/wimmer/ace/12mdata/2007/"
timeframe = [[100., 150.]]
ionname = "He2+"
data = swicsiondata(ionname, timeframe, path)
data.load()
#Zeittakte durchlaufen
vergleich = []
for takt in range(len(data.densspec)):
#Maximum ermitteln
#maximum = 0
#test = 0
#for i in range(len(data.densspec[takt])):
# if (test < data.densspec[takt][i][1]):
# maximum = i
# test = data.densspec[takt][i][1]
#Ende Maximum Ermitteln
#Index der Schwerpunktsgeschwindigkeit
gp=Gnuplot()
gp2=Gnuplot()
gp3=Gnuplot()
gp4=Gnuplot()
gp5=Gnuplot()
first=1
magdirarr=[]
meanvalf=[]
meanv=[]
meanb=[]
meanbx=[]
meanbz=[]
meanax=[]
meanaz=[]
for ion in ionlist:
iondata=swicsiondata(ion,timeframe,datapath)
iondata.load()
syncmag=sync_mag(mag,iondata)
syncswepam=sync_h(swepam,iondata)
syncatt=sync_attitude(att,iondata)
# build all arrays for histograms!
dvarr=[]
dvarr2=[]
vtharr=[]
vtharr2=[]
temparr=[]
densarr=[]
swepamvth=[]
meandvvsvparr.append([])
for i in range(int(nrxbins)):
#meandvvsvparr[-1].append([float(i*vxbinwidth),0.,0.])
meandvvsvparr[-1].append([10**(-3 + float(i * vxbinwidth)), 0., 0.])
#meandvvsvparr=array(meandvvsvparr)
mass = []
charge = []
ionnr = -1
for ion in ionlist:
ionnr += 1
iondata = []
for ion2 in ion:
#iondata=swicsiondata(ion2,resswics,year,timeframe,datapath)
iondata.append(swicsiondata(ion2, resswics, year, timeframe, datapath))
mass.append(iondata[-1].mass)
charge.append(iondata[-1].charge)
swepam = swepamdata(resswepam, year, timeframe, swepampath)
swepam.sync_h(iondata[0])
mag = magdata(resmag, year, timeframe, magpath)
mag.sync_mag(iondata[0])
dvvsvparr = []
for i in range(nrxbins):
dvvsvparr.append([])
# build all arrays for histograms!
for i in range(len(iondata[0].time)):
p = -1
if (swepam.dens[i] > 0.):
valf = 21.8 * mag.magb[i] / sqrt(swepam.dens[i])
fg = 4.
fk = -1.5
res = 5. # = 1hdata
path = "/data/wimmer/ace/1hdata/2001/"
timeframe = [[85., 86.99], [89., 105.]]
timeframe = [[89.94, 89.97]]
timeframe = [[90.01, 90.04]]
ions = [
"Fe7+", "Fe8+", "Fe9+", "Fe10+", "Fe11+", "Fe12+", "Fe13+", "Fe14+",
"Fe15+", "Fe16+", "Fe17+", "Fe18+", "Fe19+", "Fe20+", "Fe21+", "Fe22+",
"Fe23+", "Fe24+"
]
#ions=["He2+"]
# lade ionen daten
hedata = swicsiondata("He2+", timeframe, path)
hedata.load()
fedata = []
for ion in ions:
fedata.append(swicsiondata(ion, timeframe, path))
fedata[len(fedata) - 1].load()
# erzeugt container für energy pro amu spektren in nativer instrumenteller aufloesung
epaspec = []
for ion in fedata:
epaspec.append([])
for val in ion.countspec[0]:
epaspec[len(epaspec) - 1].append([vel2epa(val[0]), 0., 0., 0., 0.])
# enthält die breiten der E/amu Punkte
mass = []
charge = []
meandvvsvparr = []
dvvsvparr = []
ratioarr = []
for i in range(len(ionlist)):
meandvvsvparr.append([0., 0., 0., 0.])
dvvsvparr.append([])
ratioarr.append([])
p = -1
magarr = []
alfarr = []
masscharge = []
for tf in timeframe:
p += 1
hedata = swicsiondata("He2+", resswics, year, [tf], datapath)
swepam = swepamdata(resswepam, year, [tf], swepampath)
swepam.sync_h(hedata)
mag = magdata(resmag, year, [tf], magpath)
mag.sync_mag(hedata)
ionnr = -1
for ion in ionlist:
ionnr += 1
iondata = swicsiondata(ion[0], resswics, year, [tf], datapath)
meandvvsvparr[ionnr][0] = iondata.mass / iondata.charge
meandvvsvparr[ionnr][3] = iondata.mass / iondata.charge**2
if p == 0:
masscharge.append([iondata.mass, iondata.charge, 0, 0])
for i in range(len(iondata.time)):
if (swepam.dens[i] > 0.):
from libacepy.swicselemdata import swicselemdata
from libacepy.swicsprotondata import swicsprotondata
from libacepy.mag import magdata, sync_mag
from math import cos, acos
from Gnuplot import Gnuplot
swicspath = "/data/wimmer/ace/1hdata/2001/"
ppath = "/data/wimmer/ace/1hdataprotons/2001/"
magpath = "/data/ivar/ace/mag/"
swepampath = "/data/ivar/ace/swepam/"
timeframequietplot = [[4., 9.], [17.4, 25.25], [30.5, 35.], [45.5, 55.],
[60., 63.], [65.5, 70.], [73.25, 80.5], [93., 97.5],
[103., 107.], [119., 126.5], [129., 131.], [134., 135.],
[145.25, 151.1]]
timeframequietplot = [[1., 360.]]
he = swicsiondata("He2+", timeframequietplot, swicspath)
he.load()
h = swepamdata(720, 2001, 1., 362., swepampath)
h.load()
hsync = sync_h(h, he)
"""
hswics=swicsprotondata("H1+",timeframequietplot,ppath)
hswics.load()
"""
mag = magdata(240, 2007, 1., 362., magpath)
mag.load()
magsync = sync_mag(mag, he)
o = swicselemdata("O", timeframequietplot, swicspath)
fe = swicselemdata("Fe", timeframequietplot, swicspath)
o.load()
from libacepy.swicsiondata import swicsiondata
from Gnuplot import Gnuplot, Data
from libacepy.multiplot import multiplot
gp = Gnuplot()
path = "/data/wimmer/ace/1hdata/2007/"
timeframe = [[1., 3.]]
data = swicsiondata("O7+", timeframe, path)
data.load()
dens = []
for i in range(len(data.dens)):
dens.append([data.time[i], data.dens[i][0]])
gp.plot(Data(dens))
