Esope source files

   1 : C ASKSTB    SOURCE    CHAT      07/10/22    21:15:10     5921           
   2 :       SUBROUTINE ASKSTB(XKSTAB,XK11,XK22,XK33,XK12,XK21,XK13,XK23,
   3 :      !                      XK31,XK32)
   4 :       IMPLICIT REAL*8  (A-H,O-Z)
   5 :       IMPLICIT INTEGER (I-N)
   6 :  
   7 : 
   8 : *      SUBROUTINE ASSE_KSTAB(KSTAB,K11,K22,K33,K12,K21,K13,K23,
   9 : *     !                      K31,K32)
  10 : C ANTOINE LEGAY LE 31/01/01
  11 : *      REAL *8 KSTAB(24,24),K11(8,8),K22(8,8),K33(8,8),K12(8,8)
  12 : *      REAL *8 K21(8,8),K13(8,8),K23(8,8),K31(8,8),K32(8,8)
  13 : *      INTEGER I,J
  14 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC      
  15 : C
  16 : C ASSEMBLAGE DE LA MATRICE DE STABILISATION
  17 : C AVEC LES 3 MATRICES 8*8 XK11 XK22 ET XK33
  18 : C
  19 : 
  20 : 
  21 :       DIMENSION XKSTAB(24,24),XK11(8,8),XK22(8,8),XK33(8,8),XK12(8,8)
  22 :       DIMENSION XK21(8,8),XK13(8,8),XK23(8,8),XK31(8,8),XK32(8,8)
  22 :       CALL ZDANUL(XKSTAB,576)
  23 :       DO I=1,8
  24 :          DO J=1,8
  25 :             XKSTAB(I,J)=XK11(I,J)
  26 :          END DO
  27 :       END DO
  28 :       DO I=9,16
  29 :          DO J=9,16
  30 :             XKSTAB(I,J)=XK22(I-8,J-8)
  31 :          END DO
  32 :       END DO
  33 :       DO I=17,24
  34 :          DO J=17,24
  35 :             XKSTAB(I,J)=XK33(I-16,J-16)
  36 :          END DO
  37 :       END DO
  38 :       DO I=1,8
  39 :          DO J=9,16
  40 :             XKSTAB(I,J)=XK12(I,J-8)
  41 :          END DO
  42 :       END DO
  43 :       DO I=9,16
  44 :          DO J=1,8
  45 :             XKSTAB(I,J)=XK21(I-8,J)
  46 :          END DO
  47 :       END DO
  48 :       DO I=1,8
  49 :          DO J=17,24
  50 :             XKSTAB(I,J)=XK13(I,J-16)
  51 :          END DO
  52 :       END DO
  53 :       DO I=9,16
  54 :          DO J=17,24
  55 :             XKSTAB(I,J)=XK23(I-8,J-16)
  56 :          END DO
  57 :       END DO
  58 :       DO I=17,24
  59 :          DO J=1,8
  60 :             XKSTAB(I,J)=XK31(I-16,J)
  61 :          END DO
  62 :       END DO
  63 :       DO I=17,24
  64 :          DO J=9,16
  65 :             XKSTAB(I,J)=XK32(I-16,J-8)
  66 :          END DO
  67 :       END DO
  68 :       RETURN
  69 :       END
  70 : 
  71 :