|
| | Dispersion.dirname = sys.argv[1] |
| | Dispersion.ptnoninteractive = int(os.environ.get('PTNOINTERACTIVE', '0')) |
| bool | Dispersion.have_tqdm = True |
| x | Dispersion.tqdm = lambda x |
| str | Dispersion.do_windowing = "both" |
| list | Dispersion.timesteps = [] |
| | Dispersion.parts = filename.split("/")[-1].split(".") |
| | Dispersion.tsize = len(timesteps) |
| | Dispersion.f = analysator.vlsvfile.VlsvReader(dirname+"/bulk."+"{:07d}".format(t)+".vlsv") |
| | Dispersion.fg_b = f.read_fsgrid_variable("fg_b") |
| | Dispersion.B0vec = np.array([np.average(fg_b[:,0]), np.average(fg_b[:,1]), np.average(fg_b[:,2])]) |
| | Dispersion.B0 = np.sqrt(np.sum(B0vec**2)) |
| | Dispersion.config = f.get_config() |
| | Dispersion.dt = f.read_parameter("dt") |
| | Dispersion.dtout = float(config["io"]["system_write_t_interval"][0]) |
| | Dispersion.xmin = f.read_parameter("xmin") |
| | Dispersion.xmax = f.read_parameter("xmax") |
| tuple | Dispersion.dx = (xmax-xmin)/xsize |
| | Dispersion.ni = float(config["proton_Dispersion"]["rho"][0]) |
| | Dispersion.Ti = float(config["proton_Dispersion"]["Temperature"][0]) |
| | Dispersion.ne |
| | Dispersion.Te |
| | Dispersion.Wci = SI.e * B0 / SI.mp |
| | Dispersion.Wce = SI.e * B0 / SI.me |
| | Dispersion.wpi = np.sqrt(ni * SI.e**2 / SI.mp / SI.eps0) |
| | Dispersion.wpe = np.sqrt(ne * SI.e**2 / SI.me / SI.eps0) |
| | Dispersion.vthi = np.sqrt(2.*SI.kB * Ti / SI.mp) |
| | Dispersion.vA = B0 / np.sqrt(SI.mu0 * (SI.me*ne + SI.mp*ni)) |
| | Dispersion.vthe = np.sqrt(2.*SI.kB * Te / SI.me) |
| | Dispersion.di = SI.c / wpi |
| | Dispersion.de = SI.c / wpe |
| | Dispersion.ri = vthi / Wci |
| | Dispersion.re = vthe / Wce |
| | Dispersion.lD = vthe / wpe |
| | Dispersion.B = np.zeros( (len(timesteps), xsize, 5) , dtype=complex) |
| list | Dispersion.t = timesteps[i] |
| | Dispersion.spatial_window = np.hamming(xsize) |
| | Dispersion.temporal_window = np.hamming(tsize) |
| | Dispersion.window = np.outer(spatial_window, temporal_window).T |
| list | Dispersion.componentnames = ["x","y","z", "left", "right"] |
| | Dispersion.total |
| | Dispersion.X = np.linspace(xmin, xmax, xsize) |
| | Dispersion.T = np.linspace(timesteps[0], timesteps[-1], len(timesteps)) |
| | Dispersion.vmax = np.amax(abs(B[2:,:,c])) |
| | Dispersion.im = plt.pcolormesh(X/ri, T*Wci, np.real(B[:,:,c]), shading="gouraud", vmin=-vmax, vmax=vmax) |
| | Dispersion.label |
| | Dispersion.kB = np.fft.fftshift(np.fft.fft2(B[:,:,c]*window)) |
| int | Dispersion.w = 2.*np.pi*np.fft.fftshift(np.fft.fftfreq(tsize, d=dtout)) |
| int | Dispersion.kx = 2.*np.pi*np.fft.fftshift(np.fft.fftfreq(xsize, d=dx)) |
| int | Dispersion.kleft = 1./SI.c * np.sqrt(w**2 - wpe**2/(1.+Wce/w) - wpi**2/(1.-Wci/w)) |
| int | Dispersion.kright = 1./SI.c * np.sqrt(w**2 - wpe**2/(1.-Wce/w) - wpi**2/(1.+Wci/w)) |
| int | Dispersion.powerB = kB.real**2 + kB.imag**2 |
| int | Dispersion.vmin = vmax - 7 |
| | Dispersion.color |
| | Dispersion.linestyle |
| | Dispersion.linewidth |
| | Dispersion.sB = np.fft.fftshift(np.fft.fft(B[:,:,c]*window, axis=1), axes=1) |
| int | Dispersion.spowerB = sB.real**2 + sB.imag**2 |
| int | Dispersion.mask = 0. |
1.16.1