Gibbs2 输出文件与Gnuplot绘图

下载Gibbs2实操包。解压后进入包目录；完整输入、原始数据、提取脚本与高清图随包提供。本专题使用1.0.68固定源码版本、GNU Fortran13.3构建的双精度串行CPU入口，配套Python/Matplotlib与Gnuplot处理结果。每个流程从新的目录开始。

1. 先找对文件#

文件名前缀由ROOT决定，官方例子可能为001_io、001_phases，其他例常为result。表头就是当前文件的单位和列定义，不用记忆另一版本的列号。

2. 实际原生Gnuplot出图#

bash plot-native.sh
ls native-figures
cat native-figures/NATIVE-PLOTS.json

它读取包内实际原生.efit/.aux/.gnu，生成3200×1800 PNG。转换只改绘图终端、线宽与自动y范围，让拟合最低点完整显示；原始数据不改。原生脚本原有EPS/PDF转换命令不在这个PNG流程中执行，不需要额外TeX后处理。

实际原生Si EOS图，含五个输入点及密集拟合曲线。

实际原生C的Vinet图，不是手绘的示意曲线。

3. 提取与组合图#

bash analyze.sh
head -n 20 analysis-reference/SUMMARY.json
ls figures-reference

图、TABLES.json、相界表和独立模型对照均从真实原始文件生成。低温gamma打印缺失保留为空；不会用0或插值假装原程序给出了数字。CSV/JSON里的更多数字也不等于更高材料精度。

从原始结果重新提取、画图#

若想先看配套数据，在另一份新解压目录执行：

bash analyze.sh
ls analysis-reference
ls figures-reference

它读取 reference/ 原始输出，重新提取表格并生成16张高清图。完整重算时，在新目录运行 bash run-all.sh，最终数据在 analysis/、图片在 figures/。同名结果目录已存在时会受到保护；只算一个条件时，不要调用需要全部条件的整包分析。

相关入口：Gibbs2专题、固定版本源码与手册。以本专题对应的集群版本与输入为准；有限体积、温压和频率网格是操作示例，真实材料计算需要检查相应收敛范围。