= 开发手册 [http://www.lanseyilin.cn:8080/trac/wiki/code_page_Eng Development Manual]
== Matlab Code
ISSDE分析文件
分类：输入设置，输入分析，输出分析，脚本

输入设置：
{{{
get_input_file.m
get_n_T_for_read_data.m
get_psirz_psidr_psidz.m
get_ripple_deltarz_deltar_deltaz.m
get_xyz_vxyz.m
manipulateGfile.m
readGfile.m
MHD_input_setting:
    R_z_from_psi.m
    poloidal_theta_calculation.m
MHD_trace_input_setting
    R_phi_z_each_psi.m
    R_phi_z_pitch_energy_each_psi.m
    get_cartesian_coordinate_from_guiding_center_each_psi_mhd.m
}}}

输入分析：
{{{
NBI_position_distribution.m
magnetic_line_psi.m
n_T_plot.m
particle_orbit_type.m
psi_distribution_fun.m
q_profile.m
ripple_delta_engineering.m
}}}

理论模拟分析：
{{{
magnetic_island_width_calculation.m
}}}

输出分析：
{{{
poincare_map_plot.m
magnetic_trace.m
}}}

脚本：
{{{
run_diff_E_alpha_mhd_phase_space.sh
}}}

== C++ Code
=== ISSDE程序流程图
[[Image(htdocs:flow_chart_english.png,50%,title=ISSDE与拟牛顿法流程图)]]

=== ISSDE程序框架
[[Image(htdocs:ISSDE_different_scenarios.png,50%,title=ISSDE程序框架)]]

=== ISSDE归一化
程序采用的归一化单位是电子质量[[math(m_{e}=9.10938356\times10^{-31}\rm{kg})]]，电荷量大小[[math(q=1.602176565\times10^{-19}\rm{C})]]，光速[[math(c=3.0\times10^{8}\rm{m/s})]]，以及真空介电常数[[math(\varepsilon_0=8.854187817\times 10^{-12})]]。然后对其他量进行归一化。使用的归一化物理量是[[math(m_{e},e,c,\varepsilon_{0})]]，结果如下表
||  变量  ||  符号  ||  单位  ||
||  质量  ||  [[math(m)]]  ||  [[math(m_e)]]  ||
||  电荷  ||  [[math(q)]]  ||  [[math(|e|)]]  ||  
||  速度  ||  [[math(\bf{v})]]  ||  [[math(c)]]  || 
||  位置  ||  [[math(\bf{x})]]  ||  [[math(\frac{e^2}{m_e c^2 \varepsilon_0})]]  ||
||  时间  ||  [[math(t)]]  ||  [[math(\frac{e^2}{m_e c^3 \varepsilon_0})]]  ||
||  电场  ||  [[math(\bf{E})]]  ||  [[math(\frac{m_e ^2 e^4 \varepsilon_0}{|e|^3})]]  ||
||  磁场  ||  [[math(\bf{B})]]  ||  [[math(\frac{m_e ^2 e^3 \varepsilon_0}{|e|^3})]]  ||
||  磁面函数  ||  [[math(\psi)]]  ||  [[math(\frac{|e|}{c\varepsilon_0})]]  ||
||  磁面函数的导数  ||  [[math(\partial\psi/\partial R)]]  ||  [[math(m_ec)]]  ||


== 版本库分支图

[http://www.lanseyilin.cn:8080/trac/wiki 前页]