Esope source files

   1 : C CUBKSI    SOURCE    CHAT      07/10/22    21:15:34     5921           
   2 :       SUBROUTINE CUBKSI(NPINT,XXG,BKSI)
   3 :       implicit real*8(A-H,O-Z)
   4 :       implicit integer (I-N)
   5 :  
   6 : *      SUBROUTINE CUB_KSI(NPINT,XXG,BKSI)
   7 : *
   8 : *     ------------------------------------------------------------------
   9 : *
  10 : *         DERIVEE P/R VARIABLES CANONIQUE AUX POINTS D INTEGRATION
  11 : *          POUR LES ELEMENTS CUBES A 8 NOEUDS
  12 : *                                                    H. BUNG     03-93
  13 : *     ------------------------------------------------------------------
  14 : *
  15 : *
  16 : *   ENTREES 
  17 : *     NPINT   : NBRE DE  PTS D INTEGRATION
  18 : *     XG      : COOR. CANONIQUES DES  PTS D INTEGRATION
  19 : *
  20 : *   SORTIES :  
  21 : *     BKSI(3,8,NPINT) : LES DERIVEES
  22 : *
  23 : *      IMPLICIT NONE
  24 : *
  25 : *---    VARIABLES GLOBALES 
  26 : *
  27 : *      INTEGER NPINT
  28 : *      REAL *8 XG(3,5),XXG(5),GS(8,4)
  29 : *      REAL *8 BKSI(3,8,5)
  30 :        dimension XG(3,5),XXG(*),GS(8,4),BKSI(3,8,*)
  31 : 
  32 : *
  33 : *---    VARIABLES LOCALES 
  34 : *
  35 : *      REAL *8 AUX
  36 : *      REAL *8 UN,UNS8
  37 : *      INTEGER IP,I
  38 :        dimension XP(3),XM(3)
  39 : 
  40 :       UN=1.D0
  41 :       UNS8=0.125D0
  42 : *
  43 : C MODIF ANTOINE:
  44 : C ON REMPLIT LE TABLEAU XG
  45 : C
  46 :       DO IP=1,NPINT
  47 :       XG(1,IP)=0.d0
  48 :       XG(2,IP)=0.d0
  49 :       XG(3,IP)=XXG(IP)
  50 :       END DO
  51 : *
  52 : *
  53 :       DO IP=1,NPINT
  54 : *-
  55 :         DO  I=1,3
  56 :           XP(I) = UN + XG(I,IP)
  57 :           XM(I) = UN - XG(I,IP)
  58 :         END DO
  59 : *-
  60 :         AUX= -UNS8*XM(3)
  61 :         BKSI(1,1,IP)= AUX*XM(2)
  62 :         BKSI(1,2,IP)=-BKSI(1,1,IP)
  63 :         BKSI(1,3,IP)=-AUX*XP(2)
  64 :         BKSI(1,4,IP)=-BKSI(1,3,IP)
  65 :         BKSI(2,1,IP)= AUX*XM(1)
  66 :         BKSI(2,2,IP)= AUX*XP(1)
  67 :         BKSI(2,3,IP)=-BKSI(2,2,IP)
  68 :         BKSI(2,4,IP)=-BKSI(2,1,IP)
  69 :         AUX= -UNS8*XP(3)
  70 :         BKSI(1,5,IP)= AUX*XM(2)
  71 :         BKSI(1,6,IP)=-BKSI(1,5,IP)
  72 :         BKSI(1,7,IP)=-AUX*XP(2)
  73 :         BKSI(1,8,IP)=-BKSI(1,7,IP)
  74 :         BKSI(2,5,IP)= AUX*XM(1)
  75 :         BKSI(2,6,IP)= AUX*XP(1)
  76 :         BKSI(2,7,IP)=-BKSI(2,6,IP)
  77 :         BKSI(2,8,IP)=-BKSI(2,5,IP)
  78 :         AUX= -UNS8*XM(2)
  79 :         BKSI(3,1,IP)= AUX*XM(1)
  80 :         BKSI(3,2,IP)= AUX*XP(1)
  81 :         BKSI(3,5,IP)=-BKSI(3,1,IP)
  82 :         BKSI(3,6,IP)=-BKSI(3,2,IP)
  83 :         AUX= -UNS8*XP(2)
  84 :         BKSI(3,3,IP)= AUX*XP(1)
  85 :         BKSI(3,4,IP)= AUX*XM(1)
  86 :         BKSI(3,7,IP)=-BKSI(3,3,IP)
  87 :         BKSI(3,8,IP)=-BKSI(3,4,IP)
  88 : *
  89 :       END DO
  90 : *
  91 :       END
  92 :       
  93 : 
  94 :