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 __init__(self, year, timeframe):
self.year = year
self.timeframe = timeframe
self.hedata = loadswicsiondata("He2+", 720., self.year, self.timeframe,
"/data/etph/ace/")
self.odata = loadswicsiondata("O6+", 720., self.year, self.timeframe,
"/data/etph/ace/")
self.hdata = swicsdcrtcr("H1+", self.year, self.timeframe,
"/data/etph/ace/pui/12mdata/")
self.swedata = swepamdata(64., self.year, self.timeframe,
"/data/etph/ace/swepam/")
self.swedata.sync_h(self.hedata)
self.mag = magdata(1., self.year, self.timeframe,
"/data/etph/ace/mag/")
self.mag.sync_mag(self.hedata)
self.theta = arccos(cos(self.mag.phi) * cos(self.mag.theta))
self.valf = 21.8 / sqrt(self.swedata.dens) * (self.mag.magb)
self.valf[isinf(self.valf)] = -1.
self.valf[isnan(self.valf)] = -2.
self.effh = zeros([58])
self.effhe = zeros([58])
self.effo = zeros([58])
self.gp = Gnuplot()
self.gp("set xr [-2:3]")
self.gp("set yr [0:1]")
self.gp("set xlabel '(v-v_{sw})/C_{A}'")
self.stdhe = []
self.stdo = []
self.stdh = []
self.skewhe = []
self.skewo = []
self.skewh = []
self.kurthe = []
self.kurto = []
self.kurth = []
self.meanvhe = []
self.meanvo = []
self.meanvh = []
self.maxvhe = []
self.maxvo = []
self.maxvh = []
self.time = []
self.years = []
self.valfs = []
self.magb = []
self.thetas = []
self.sigtheta = []
self.vsw = []
self.dsw = []
self.tsw = []
self.load_eff()
def __init__(self, swicspath, swepampath, magpath, effpath, ionname,
timeframe):
self.swicspath = swicspath
self.swepampath = swepampath
self.magpath = magpath
self.effpath = effpath
self.ionname = ionname
self.timeframe = timeframe
self.swepamres = 720
self.magres = 240
self.year = 2007
swepam = swepamdata(720, 2007, timeframe[0][0],
timeframe[len(timeframe) - 1][1] + 1., swepampath)
swepam.load()
mag = magdata(240, 2007, float(daylist[0]),
float(daylist[len(daylist) - 1]), magpath)
mag.load()
att = attdata(mag.year, mag.path)
att.load()
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()
fe.load()
feoheh = []
def load_swepam(self):
from libacepy.swepamdata import swepamdata
self.swepam = swepamdata(
self.swepamres, self.year, self.timeframe[0][0],
self.timeframe[len(self.timeframe) - 1][1] + 1., self.swepampath)
self.swepam.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.
mp.panel[mpi].adddata(
Data([[timeframe[0], parkerwinkel], [timeframe[1], parkerwinkel]]))
mp.panel[mpi].adddata(
Data([[timeframe[0], pi + parkerwinkel], [timeframe[1],
pi + parkerwinkel]]))
mp.panel[mpi].setyrange(-2, 6)
mp.panel[mpi].setlogscale(0)
mp.panel[mpi].setylabel("B-Winkel")
mpi += 1
del plotdata
del mdata
del mbin
# Ende: Magdaten auswerten
# Anfang: Swepam auswerten
pdata = swepamdata(pres, year, timeframe[0], timeframe[1],
"/data/ivar/ace/swepam/")
pdata.load()
binres = bintime / pres
pvbin = [] #Bins der Geschwindigkeit
ptbin = [] #Bins der Temperatur
pnbin = [] #Bins der Dichte
for i in range(min(len(pdata.vel), len(pdata.temp))):
if (i % binres == 0):
if (i > 0):
pvbin[len(pvbin) - 1][1] /= binres
ptbin[len(ptbin) - 1][1] /= binres
ptbin[len(pnbin) - 1][1] /= binres
pvbin.append([pdata.time[i], pdata.vel[i][1]])
ptbin.append([pdata.time[i], pdata.temp[i][1]])
pnbin.append([pdata.time[i], pdata.dens[i][1]])
timeframe = [[167., 168.]]
timeframe = [[222, 223.]]
timeframe = [[39., 40.]]
timeframe = [[4., 5.]]
timeframe = [[153., 155.]]
timeframe = [[32., 34.]]
#year=2008
#timeframe=[[20.,22.]]
hdata = swicsdcrtcr("H1+", year, timeframe, "/data/etph/ace/pui/12mdata/")
hedata = loadswicsiondata("He2+", 720., year, timeframe, "/data/etph/ace/")
odata = loadswicsiondata("O6+", 720., year, timeframe, "/data/etph/ace/")
fedata = loadswicsiondata("Fe9+", 720., year, timeframe, "/data/etph/ace/")
sidata = loadswicsiondata("C6+", 720., year, timeframe, "/data/etph/ace/")
cdata = loadswicsiondata("Fe8+", 720., year, timeframe, "/data/etph/ace/")
swdata = swepamdata(64., year, timeframe, "/data/etph/ace/swepam/")
magdat = magdata(1., year, timeframe, "/data/etph/ace/mag/")
#magdat.sync_mag(hedata)
#gp=Gnuplot()
def get_tseries(iond):
from numpy import where, mean
timestep = 12. / 86400. # time of single esa-step
d3iondata = []
d3ionv = []
for i in range(len(iond.time)):
ionmax = amax(iond.countspec[i, :, 1])
vtime = 0.
for j in range(58):
wcarr=[]
magtheta=[]
vswarr=[]
timearr=[]
magarr=[]
for year in years:
print year
for doy in range(start,doys[ycount]+1):
#for doy in range(start,stop+1):#doys[ycount]+1):
try:
listin=open(path+"%.4i/%.3i/list.in"%(year,doy),"r")
hesum=0.
hesumeff=0.
hesumeff2=0.
magdoy=magdata(1.,year,[[doy,doy+0.999999]],magpath)
swepamdoy=swepamdata(64.,year,[[doy,doy+0.999999]],swepampath)
for s in listin:
k=s.split()
tmpdat,tmptime,tmpvsw=read_pui(path+"%.4i/%.3i/He1+_%s"%(year,doy,k[0]))
vswepam=mean(swepamdoy.vel[(swepamdoy.time>tmptime)*(swepamdoy.time<tmptime+12./1440.)*(swepamdoy.vel>0.)])
if vswepam>0.:
v=tmpdat[:,0]
#w=tmpdat[:,1]
w=v/vswepam
dcr=tmpdat[:,2]
tcr=tmpdat[:,3]
mag=mean(magdoy.magb[(magdoy.time>tmptime)*(magdoy.time<tmptime+12./1440.)*(magdoy.magb>0.)])
magang=mean(arccos(cos(magdoy.phi[(magdoy.time>tmptime)*(magdoy.time<tmptime+12./1440.)])*cos(magdoy.theta[(magdoy.time>tmptime)*(magdoy.time<tmptime+12./1440.)])))
wcarr.append([w[(w>1.)*(w<2.5)],dcr[(w>1.)*(w<2.5)]+tcr[(w>1.)*(w<2.5)]])
magtheta.append(magang)
vswarr.append(vswepam)
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.):
valf = 21.8 * mag.magb[i] / sqrt(swepam.dens[i])
#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()
sigarr = []
barr = []
i = 0
for s in dblistin:
k = s.split()
from libacepy.funcs import twobimaxmag
ion = "He2+"
ion2 = "C5+"
res = 720.
year = 2008
path = "/data/wimmer/ace/"
#timeframe=array([[60.,90.]])
#for i in range(19):
for i in range(1):
timeframe = array([[float(i) * 20. + 1., float(i + 1) * 20. + 1]])
timeframe = array([[10.25, 10.75]])
hedata = loadswicsiondata(ion, res, year, timeframe, path)
odata = loadswicsiondata(ion2, res, year, timeframe, path)
hdata = swepamdata(64., year, timeframe, path + "swepam/")
hdata.sync_h(hedata)
magdata = loadmagdata(1., year, timeframe, path + "mag/")
magdata.sync_mag(hedata)
att = attdata(year, "/data/ivar/ace/")
att.load()
syncatt = sync_attitude(att, hedata)
thetaarr = arccos(cos(magdata.phi) * cos(magdata.theta))
valf = 21.8 * magdata.magb / sqrt(hdata.dens)
hphi = arctan(hdata.vely / hdata.velx)
htheta = arctan(hdata.velz / sqrt(hdata.velx**2 + hdata.velz**2))
vy = tan(syncatt.phi) * hdata.vel - hdata.vely
vz = tan(syncatt.theta) * hdata.vel - hdata.velz
anx=b
return [anx,anz]
daylist=range(1,153)
timeframe=[[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.,135.],[145.25,151.1]]
#timeframe=[[4.,9.]]
effpath="/data/ivar/ace/efficencies/aspeff/"
datapath="/data/wimmer/ace/12mdata/2007/"
datapath2="/data/zechine/ace/12mprotons/2007/"
swepampath="/data/ivar/ace/swepam/"
magpath="/data/ivar/ace/mag/"
#protondata=swicsprotondata("H1+",timeframe,datapath2)
#protondata.load()
ionlist=["He2+","C6+","C5+","N7+","N6+","O7+","O6+","Ne8+","Mg8+","Mg9+","Si8+","Si9+","S8+","S9+","Ca9+","Ca10+","Fe9+","Fe10+","Fe11+","Fe12+","Fe13+"]
ionlist=["He2+"]
swepam=swepamdata(720,2007,timeframe[0][0],timeframe[len(timeframe)-1][1]+1.,swepampath)
swepam.load()
mag=magdata(240,2007,float(daylist[0]),float(daylist[len(daylist)-1]),magpath)
mag.load()
att=attdata(mag.year,mag.path)
att.load()
pdvhist=[]
pdvhist2=[]
pdvmqhist=[]
ptemphist=[]
pdenshist=[]
pvthhist=[]
pvthhist2=[]
swepamvel=[]
swepamvthhist=[]