Esope source files

   1 : C MKREDI    SOURCE    AM        11/05/26    21:16:01     6982           
   2 :       subroutine MKREDI(da,XMKel)
   3 : C     Matrice de rigiditÃ© Ã©lastique redimensionnÃ©e
   4 : 
   5 :         IMPLICIT INTEGER(I-N)
   6 :         IMPLICIT REAL*8(A-H,O-Z)
   7 :         DIMENSION da(20), XMKel(5,5), XMKelinv(5,5)
   8 : *
   9 :         XKelz=da(4)
  10 :         XKelh=da(5)
  11 :         XKelry=da(6)
  12 :         XKelrx=da(7)
  13 : *
  14 :         do i=1,5
  15 :            do j=1,5
  16 :               XMKel(i,j)=0.
  17 :               XMKelinv(i,j)=0.
  18 :            enddo
  19 :         enddo
  20 : *
  21 :         XMKel(1,1)=XKelz
  22 :         XMKel(2,2)=XKelh
  23 :         XMKel(3,3)=XKelry
  24 :         XMKel(4,4)=XKelh
  25 :         XMKel(5,5)=XKelrx
  26 :         XMKelinv(1,1)=1./(XKelz)
  27 :         XMKelinv(2,2)=1./(XKelh)
  28 :         XMKelinv(3,3)=1./(XKelry)
  29 :         XMKelinv(4,4)=1./(XKelh)
  30 :         XMKelinv(5,5)=1./(XKelrx)
  31 : *
  32 :         return
  33 :         end
  34 :