铀原子  

 

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)

 

posted on   Stronger_happier  阅读(19)  评论(0编辑  收藏  举报
努力加载评论中...
 
点击右上角即可分享
微信分享提示