Fortran77矩阵格式输出
加油
本次解决了矩阵输出格式的问题。
通过不换输出,模拟矩阵格式。
子程序可以求矩阵的偏应力。
PROGRAM MAIN
REAL XIDENTY(3,3),STRESS(3,3),DEVIATO(3,3)
DO I=1,3
DO J=1,3
XIDENTY(I,J)=0.D0
END DO
XIDENTY(I,I)=1.D0
END DO
C note that the J just exist in the above loop
DO I=1,3
DO J=1,3
STRESS(I,J)=1.D0
END DO
END DO
STRESS(1,1)=1.d0
STRESS(2,2)=0.d0
STRESS(2,2)=0.d0
CALL KDEVIA(STRESS,XIDENTY,DEVIATO)
c***********