具体内容请参考: https://www.tabirstrees.top/2022/02/12/Eorupainternalstructure/
parameter.py
中更改参数(一般是星球的质量分数);
ode_model.ipynb
是计算的主程序(迭代)以及绘制最后的图表;
eos_no_tem.py
是定义计算函数的地方,求解微分方程的位置,具体的计算方式可以看PPT中的描述;
core.txt
,curst.txt
,low_mantle.txt
,mantle.txt
分别是定义各个圈层的热力学参数的文件
欢迎讨论!!!(*^_^*)