python后处理一维epoch数据
import sdf_helper as sh
import matplotlib.pyplot as plt
import numpy as np
from plasmapy.formulary.densities import critical_density as nc
lambda0=1.0e-6
c=3.0e8
omega0=2*np.pi*c/lambda0
n0=0.1*nc(omega0)
for i in range(1,20):
data=sh.getdata('Data/%04d.sdf' % i)
nion=data.Derived_Number_Density_ion.data/n0.value-1
x=np.linspace(-35,35,70*50)
plt.plot(x,nion)
plt.show()
sh.list_variables(data)
