gfortran数组与逻辑标量的if语句问题及编译错误排查

现有代码的错误点 #

1. 括号不匹配：if语句中all(的左括号没有对应闭合的右括号，正确写法需要在all(...)后补上右括号，即if (all(...)) then。
2. 无效的冒号表达式：0.0:0.0是Fortran的范围表达式（用于生成序列），但这里你只是要和标量0.0比较，完全不需要冒号，直接写0.0即可。x(:) > 0.0会自动生成和x同维度的逻辑数组，供all()函数判断所有元素是否满足条件。

修正后的全局判断版本 #

如果你确实需要判断所有x元素都大于0才计算sigmay，否则设为0，修正后的代码如下：

if (all(x(:) > 0.0)) then 
  sigmay = ay * abs(x)**by * x
else
  sigmay = 0.0
end if

逐元素处理的优化方案 #

如果你实际需求是对x的每个元素单独判断（即某个元素大于0时计算对应sigmay值，否则设为0），用Fortran的where语句或merge函数更高效，也符合数组运算的风格：

方案1：where语句 #

sigmay = 0.0  ! 先初始化所有元素为0
where(x > 0.0)  ! 对满足条件的元素赋值
  sigmay = ay * abs(x)**by * x
end where

方案2：merge函数 #

sigmay = merge(ay * abs(x)**by * x, 0.0, x > 0.0)

这两种方式都避免了逻辑值直接参与实数运算的问题，同时实现逐元素处理，比循环更简洁高效。

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