C SHROT3    SOURCE    CHAT      05/01/13    03:15:49     5004
      SUBROUTINE SHROT3(QSI,ETA,DZE,IELE,SHP,IRET)
C=======================================================================
C
C     FONCTIONS DE FORME VECTORIELLES POUR L'ELEMENT ROT3 :
C       ELLES SONT LES ROTATIONNELS SURFACIQUES DES FONCTIONS
C       DE LAGRANGE DE L'ELEMENT TRI3 :
C         Wi = Grad(PHIi) x n ou n est la normale exterieure au plan de
C                                l'element
C         Wi = Grad(PHIi) tourne de -pi/2
C         Wi est parallele au cote oppose au sommet i dans le plan de
C                                l'element
C
C=======================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      PARAMETER(XZER=0.D0,UNDEMI=.5D0,UN=1.D0,DEUX=2.D0)
      PARAMETER(TROIS=3.D0,QUATRE=4.D0)
      DIMENSION SHP(6,*)
C=======================================================================
C     QSI ETA DZE COORDONNEES REDUITES DU POINT DE GAUSS
C
C     SHP FONCTIONS DE FORME SHP(1,INO) FONCTIONS DE LAGRANGE DU ND INO
C                            SHP(2,INO) COORDONNEE X DE LA FONCTION DU ND INO
C                            SHP(3,INO) COORDONNEE Y DE LA FONCTION DU ND INO
C                            SHP(4,INO) COORDONNEE Z DE LA FONCTION DU ND INO
C=======================================================================
      IRET=1
C
C    66 INDIQUE QUE L ELEMENT PAS ENCORE IMPLEMENTE
C
C           P1 S2 S3 T3 T4 T6 T7 Q4 Q5 Q8 Q9 R2 R3 C8 C20 P6 P15
      GOTO (66,66,66, 1,66,66,66,66,66,66,66,66,66,66, 66,66, 66,
C        L3 L4 L6 L8 MU T4 T10 P5 P13 ATT SUP RP3 LP6 LP8
     &   66,66,66,66,66,66, 66,66, 66, 66, 66, 66, 66, 66) IELE
C
  66  IRET=0
      GOTO 666
C
C
C     FONCTIONS DE FORME BIDIMENSIONNELLES
C
   1  CONTINUE
C
C     FONCTIONS DE FORME TRIANGLE A 3 NOEUDS
C      NOEUD1 = (0,0)
C      NOEUD2 = (1,0)
C      NOEUD3 = (0,1)
C
      AUX=UN-QSI-ETA
      SHP(1,1)=AUX
      SHP(2,1)=-UN
      SHP(3,1)=UN
      SHP(1,2)=QSI
      SHP(2,2)=XZER
      SHP(3,2)=-UN
      SHP(1,3)=ETA
      SHP(2,3)=UN
      SHP(3,3)=XZER
C
 666  CONTINUE
      RETURN
      END