这是使用 Eddingtons 公式生成 IC 的代码
如何编译?
只需运行“make”即可。 确保安装/配置了“gsl”(GNU 科学库)。
如何运行代码
要编译和运行代码,请执行以下操作:
make
./EddIso Output.bin 1000000
这将创建一个包含 1000000 个粒子的 Gadget2 文件 (Output.bin)。 请注意,输出文件中省略了具有电位的块。
参数
这里列出了一些可以修改的参数:
密度分布由以下几行决定(在 EddIso.c 中):
//DENSITY PROFILE: (inner slope, outer slope, transition...)
#define ALPHA (1.0)
#define BETA (4.0)
#define GAMMA (1.0)
(根据给定的参数,它是 Hernquist 配置文件)
密度
评论0
最新资源