编辑代码

program main
real STATEV(78)
real sigma(6)
real s(6,6,2)
INTEGER w
sigma = (/1.0,2.0,3.0,4.0,5.0,6.0/)
DO k = 1, 6
    STATEV(k) = sigma(k)
END DO

DO i = 1, 6
    DO j = 1, 6
        DO w = 1, 2
            s(i,j,w) = 6.0 + (i-1)*6.0*2.0 + (j-1)*2.0 + w
            STATEV(6 + (i-1)*6*2 + (j-1)*2 + w) = s(i,j,w)
        END DO
    END DO
END DO
print*,STATEV(:)
end