arrays - Fortran エラー: 配列式の形状が一致しません

Question

プログラムでエラーが発生しました。エラーを表示するために、この小さなプログラムと関連するエラーを追加しました。助けていただければ幸いです。

PROGRAM TRY
      IMPLICIT NONE
      REAL(KIND=8) :: A(6,6),B(6,6)
      INTEGER :: I, J
      REAL(KIND=8) :: RHO, H
      REAL(KIND=8) :: C(5,2),D(2,5),E(5,5)
!
      RHO=780
      H=0.5D0
      DO I=1,6
         DO J=1,6
             A(I,J)=1
         ENDDO
      ENDDO 

      B=RHO*H**2.0D0*A
      DO I=1,5
         DO J=1,2
             C(I,J)=I+J
         ENDDO    
      ENDDO

      OPEN(1,FILE='OUTPUT.DAT')
      DO I=1,6
          WRITE(1,20) (B(I,J), J=1,6)
      ENDDO
      DO I=1,5
         DO J=1,2
              D(J,I)=C(I,J)
         ENDDO
      ENDDO
      E=C*D
      DO I=1,5
          WRITE(1,20) (E(I,J), J=1,5)
      ENDDO
      CLOSE(1)

20    FORMAT (999(2x,EN14.4))     

      END PROGRAM TRY

次のエラー メッセージが表示されます。

fortran_worksheet.f(32): error #6366: The shapes of the array expressions do not conform.
      E=C*D
---------^
fortran_worksheet.f(32): error #6366: The shapes of the array expressions do not conform.   [E]
      E=C*D
------^
compilation aborted for fortran_worksheet.f (code 1)