def sync_data(self):
from libacepy.swepamdata import sync_h
from libacepy.mag import sync_mag
from libacepy.attitude import sync_attitude
self.syncmag = sync_mag(self.mag, self.iondata)
self.syncswepam = sync_h(self.swepam, self.iondata)
self.syncatt = sync_attitude(self.att, self.iondata)
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()
fe.load()
feoheh = []
feoheh2 = []
feoheh3 = []
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.
sphe6002007 = 0.
nsp3002001 = 0.
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))
x_max.append(max(he2_h[:,1]))
y_max.append(max(fe_o[:,1]))
tmp = get_hist_x_over_y(he2_h[:,1],fe_o[:,1])
while (len(hist) < len(tmp)):
hist.append([])
for i in range(len(tmp)):
def sync_with (self,other): if (self.datatype != "swepam" or type(other) != instance or other.__class__ != "data" or other.datatype != "swepam"): print "data.sync_with : this routine is only to sync swepam with swics data" return self.data = sync_h(self.data,other.data)
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=[]
for i in range(len(iondata.time)):
if (1):
# calculate mag angle