Esope source files

   1 : C CUKSKP    SOURCE    CHAT      07/10/22    21:15:40     5921           
   2 :       SUBROUTINE CUKSKP(IFACE,BKSI)
   3 :       implicit real*8(A-H,O-Z)
   4 :       implicit integer (I-N)
   5 : 
   6 : 
   7 : *      SUBROUTINE CUB_KSI_KP(IFACE,BKSI)
   8 : C                ELEMENT SHB8
   9 : C  DEDIEE AU CALCUL DE KP: POUR LES 4 NOEUDS DE LA FACE PRESSION
  10 : C
  11 : C     ------------------------------------------------------------------
  12 : C
  13 : C        DERIVEE P/R VARIABLES CANONIQUE AUX POINTS D INTEGRATION
  14 : C         POUR LES ELEMENTS CUBES A 8 NOEUDS
  15 : C                                                   H. BUNG     03-93
  16 : C                                        CORRIGE PAR A. LEGAY POUR 4 NOEUDS
  17 : C     ------------------------------------------------------------------
  18 : C         
  19 : C        
  20 : C   ENTREES 
  21 : C     IFACE   : =1 FACE INFERIEURE
  22 : C              =2 FACE SUPERIEURE
  23 : C   SORTIES :  
  24 : C     BKSI(3,8) : LES DERIVEES
  25 : C
  26 : *      IMPLICIT NONE
  27 : C
  28 : C    VARIABLES GLOBALES 
  29 : C
  30 : *      INTEGER IFACE
  31 : *      REAL *8 BKSI(3,4)
  32 :       DIMENSION BKSI(3,4)
  33 : C
  34 : C    VARIABLES LOCALES 
  35 : C
  36 : *      REAL *8 XPOINT(4,3),ZI,UNS8
  37 :       DIMENSION XPOINT(4,3)
  38 :       UNS8=0.125D0
  39 : C
  40 :       XPOINT(1,1)=-1.
  41 :       XPOINT(1,2)=-1.
  42 :       XPOINT(2,1)=1.
  43 :       XPOINT(2,2)=-1.
  44 :       XPOINT(3,1)=1.
  45 :       XPOINT(3,2)=1.
  46 :       XPOINT(4,1)=-1.
  47 :       XPOINT(4,2)=1.
  48 :       IF(IFACE.EQ.1)THEN
  49 :          XPOINT(1,3)=-1.
  50 :          XPOINT(2,3)=-1.
  51 :          XPOINT(3,3)=-1.
  52 :          XPOINT(4,3)=-1.
  53 :          ZI=-1.
  54 :       ENDIF
  55 :       IF(IFACE.EQ.2)THEN
  56 :          XPOINT(1,3)=1.
  57 :          XPOINT(2,3)=1.
  58 :          XPOINT(3,3)=1.
  59 :          XPOINT(4,3)=1.
  60 :          ZI=+1.
  61 :       ENDIF
  62 :       IF(IFACE.NE.1.AND.IFACE.NE.2) THEN
  63 :          WRITE(6,*)'******** FACE NON DEFINIE DANS CUB_KSI_KP *******'
  64 :          STOP
  65 :       ENDIF
  66 : C
  67 :       BKSI(1,1) = UNS8*(1+XPOINT(1,3)*ZI)*XPOINT(1,1)
  68 :       BKSI(1,2) = UNS8*(1+XPOINT(2,3)*ZI)*XPOINT(2,1)
  69 :       BKSI(1,3) = UNS8*(1+XPOINT(3,3)*ZI)*XPOINT(3,1)
  70 :       BKSI(1,4) = UNS8*(1+XPOINT(4,3)*ZI)*XPOINT(4,1)
  71 : C
  72 :       BKSI(2,1) = UNS8*(1+XPOINT(1,3)*ZI)*XPOINT(1,2)
  73 :       BKSI(2,2) = UNS8*(1+XPOINT(2,3)*ZI)*XPOINT(2,2)
  74 :       BKSI(2,3) = UNS8*(1+XPOINT(3,3)*ZI)*XPOINT(3,2)
  75 :       BKSI(2,4) = UNS8*(1+XPOINT(4,3)*ZI)*XPOINT(4,2)
  76 : C
  77 :       BKSI(3,1) = UNS8*XPOINT(1,3)
  78 :       BKSI(3,2) = UNS8*XPOINT(2,3)
  79 :       BKSI(3,3) = UNS8*XPOINT(3,3)
  80 :       BKSI(3,4) = UNS8*XPOINT(4,3)
  81 :       RETURN
  82 :       END
  83 :       
  84 : 
  85 :