数组表达式的形状不符合 Fortran 90

Question

我正在尝试实现与代码最后一行相同的目标，但您可以使用您可以在我的代码中看到的 do 循环。问题是，当我尝试使用类似命令形成循环时

CURV_KVEC%KNOTS(knot)=(/0.d0/)

我收到错误

The shapes of the array expressions do not conform. 

任何帮助将不胜感激！谢谢。

cdegree=2 
count=5
do knot=0,cdegree+count
if (knot.le.cdegree) then
    CURV_KVEC%KNOTS(knot)=0.d0
    test(knot+1)=0.d0
elseif (knot.ge.count) then
    test(knot+1)=1.d0
    CURV_KVEC%KNOTS(knot)=1.d0   
elseif(nvec.eq.0) then
    num=num+1
    do kn=1,cdegree
       CURV_KVEC%KNOTS(knot+kn-1)=0.5d0  !dble(num)/(dble(count-cdegree-1)/dble(cdegree)+1.d0)
        test(knot+kn)=0.5d0  !dble(num)/(dble(count-cdegree-1)/dble(cdegree)+1.d0)
        nvec=cdegree-1
    enddo
elseif(nvec.ne.0) then 
    nvec=nvec-1
endif
write(14,*)knot,num,nvec,test(knot+1)
enddo    
! CURV_KVEC%KNOTS(0:CURV_KVEC%LENGTH) = (/0.D0,0.D0,0.D0,0.5d0,0.5d0,1.d0,1.d0,1.d0/)

Answer 1

CURV_KVEC%KNOTS(knot)是一个数组元素，一个标量。 (/0.d0/)表示长度为 1 的数组。这些不是一回事...

简单写作

CURV_KVEC%KNOTS(knot) = 0.d0

会做这个伎俩。

Answer 2

等于两个数组，左右两侧的数组必须具有相同的形状。 在这里，对于一维数组，这仅仅意味着它们具有相同的长度。 由于左侧的运行时确定长度为 CURV_KVEC%LENGTH + 1，因此不能在右侧使用固定长度的常量数组。

如果您希望 LHS 上数组的所有元素都具有相同的值，则可以按照 @Alexander Vogt 的建议使用右侧的缩放器，该值将被广播到数组。 我将示例代码解释为您想要不同的值。 如果您想要相同的值，这是最佳解决方案。