Jorge Gonzalez
0c27b98e2e
Now n_e is given as the density at the boundary and n_i at the boundary is calculated once Z is known. This aims to eliminate the iterative process.
20 lines
515 B
Python
20 lines
515 B
Python
import matplotlib.pyplot as plt
|
|
import glob
|
|
import numpy as np
|
|
import readBC
|
|
|
|
|
|
paths = ['../2025-04-10_08.40.09/']
|
|
time, n_i, u_i, T_i, Zinj, Z_Tne = readBC.read(paths[0] + 'bc.csv')
|
|
|
|
fig, ax = plt.subplots()
|
|
|
|
n_e = n_i * Zinj
|
|
|
|
plt.plot(time, n_e / n_e[0], label = f"$n_e$ ($\\times {n_e[0] * 1e-6:.0e} \\; cm^{{-3}})$")
|
|
plt.plot(time, T_i / T_i[0], label = f"$T \\; (\\times{T_i[0]:.1f} \\; eV)$")
|
|
plt.plot(time, Zinj, label = "Zinj")
|
|
plt.plot(time, Z_Tne, label = "Z_Tne")
|
|
plt.semilogy()
|
|
plt.legend()
|
|
plt.show()
|