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=[]