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))