能否在Fortran的data语句中自动提供双精度常量？

这确实是Fortran里很容易踩的单/双精度字面量坑——你遇到的情况是：不带d0的实数字面量默认是单精度（real*4），把它们赋值给双精度数组时，会先把单精度值转成双精度，但单精度的有效位数有限，所以才会出现那些看起来像“随机尾数位”的精度损失。

关于自动给数值添加d0，有几个实用方案：

1. 编译器选项强制默认实数字面量为双精度 #

很多Fortran编译器支持全局调整默认实数字面量的类型，比如：

• GCC/gfortran：使用-fdefault-real-8选项，这个选项会把默认的real类型设为8字节双精度，同时不带后缀的实数字面量也会被当作双精度处理。这样你原来的data语句不用改，就能直接得到正确的双精度值。
• Intel Fortran：可以用-r8选项，效果类似，强制默认实型为双精度。

⚠️ 注意：这个选项是全局生效的，如果你代码里还有需要单精度的部分，可能会受影响，需要提前确认。

2. 利用预处理宏自动添加后缀 #

如果不想全局修改编译器选项，可以开启Fortran的预处理功能，定义一个宏来自动给数字加上d0：
首先，编译时加上预处理选项（比如gfortran的-cpp），然后在代码开头添加：

#define D(x) x##d0

之后在data语句里就可以这样写：

data xsolar/D(12.00), D(10.93), D(1.05), D(1.38), D(2.70)/

预处理阶段会自动把D(12.00)替换成12.00d0，实现自动添加后缀的效果。

3. 用脚本批量修改现有代码 #

如果你的代码里有大量类似的data语句需要修改，可以写个简单的脚本批量处理：
比如用Python或者sed命令，匹配data语句里的实数字面量，自动添加d0。举个sed的例子（仅针对简单场景，复杂代码可能需要调整正则）：

sed -E 's/(data [a-zA-Z0-9]+\/)([0-9]+\.[0-9]+)(,|\/)/\1\2d0\3/g' your_code.f90 > modified_code.f90

这个命令会找到data语句里的小数，在后面加上d0。

另外，从代码规范的角度，更推荐的长期方案是使用kind参数来明确类型，比如：

integer, parameter :: dp = kind(1.d0)
real(dp) :: xsolar(5)
data xsolar/12.00_dp, 10.93_dp, 1.05_dp, 1.38_dp, 2.70_dp/

这样不仅更清晰，也能避免单双精度混淆的问题，而且有些IDE（比如Visual Studio Code配Fortran插件）会自动提示或者补全后缀，减少手动输入的麻烦。

内容的提问来源于stack exchange，提问作者Evgenii