Esope source files

   1 : C HYSHPT    SOURCE    AM        08/02/14    21:24:21     6050           
   2 :       SUBROUTINE HYSHPT(NEF,NBDDL,IPINTE,IPTHYB)
   3 : C-----------------------------------------------------------------------
   4 : C Calcul les fonctions de bases de la formulation hybride pour l'ef NEF
   5 : C aux points contenus dans le segment d'integration IPINTE.
   6 : C-----------------------------------------------------------------------
   7 : C
   8 : C---------------------------
   9 : C Parametres Entree/Sortie :
  10 : C---------------------------
  11 : C
  12 : C E/  NEF     : Numero de l'element fini dans NOMTP (cf common CCHAMP)
  13 : C E/  NBDDL   : Nombre de noeuds portant des degres de liberte
  14 : C E/  IPINTE  : Pointeur vers le segment d'integration des fonctions
  15 : C               de forme de l'element fini NEF
  16 : C  /S IPTHYB  : Pointeur du segment HYBSTO contenant les fonctions de
  17 : C               bases en formulation hybride
  18 : C
  19 : C----------------------
  20 : C Variables en COMMON :
  21 : C----------------------
  22 : C
  23 : C E/  IDIM    : Dimension de l'espace   (cf CCOPTIO)
  24 : C E/  NOMTP   : Noms des elements finis (cf CCHAMP)
  25 : C  /S MOTERR  : Chaine de char pour erreur (cf CCOPTIO)
  26 : C
  27 : C----------------------
  28 : C Tableaux de travail :
  29 : C----------------------
  30 : C
  31 : C     SHP     : Contient les fonctions de base en un point
  32 : C
  33 : C-----------------------------------------------------------------------
  34 : C
  35 : C Langage : ESOPE + FORTRAN77
  36 : C
  37 : C Auteurs : F.DABBENE 08/93
  38 : C
  39 : C-----------------------------------------------------------------------
  40 :       IMPLICIT INTEGER(I-N)
  41 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  42 : C
  43 : -INC CCHAMP
  44 : -INC SMINTE                                                                     
  45 : -INC CCOPTIO                                                                    
  46 : C
  47 :       SEGMENT HYBSTO
  48 :          REAL*8 HYBASE(NDIM,NBDDL,NBPP)
  49 :       ENDSEGMENT                                                                
  50 : C
  51 :       SEGMENT SHXX
  52 :          REAL*8 SHY(NDIM,NBDDL)
  53 :       ENDSEGMENT
  54 : C                                                                               
  55 : C= Initialisations
  56 : C                                                                               
  57 :       MINTE = IPINTE
  58 :       SEGACT MINTE
  59 :       NBPP  = SHPTOT(/3)
  60 :       NDIM  = IDIM * (IDIM+1)
  61 :       SEGINI HYBSTO , SHXX
  62 : *
  63 : *= Calcul des fonctions de base en chaque point du segment d'integration
  64 : *
  65 :       DO 20 II=1,NBPP
  66 :          XX = QSIGAU(II)
  67 :          YY = ETAGAU(II)
  68 :          ZZ = DZEGAU(II)
  68 :          CALL HYSHP(XX,YY,ZZ,NEF,NBDDL,NDIM,SHY,IRET)
  69 :          IF (IRET.NE.0) GOTO 100
  70 :          DO 10 ID=1,NDIM
  71 :             DO 10 NO=1,NBDDL
  72 :                HYBASE(ID,NO,II) = SHY(ID,NO)
  73 :   10     CONTINUE
  74 :   20  CONTINUE
  75 : C
  76 : C= Desactivation et menage
  77 : C
  78 :       IPTHYB = HYBSTO
  79 :       SEGDES MINTE , HYBSTO
  80 :       SEGSUP SHXX
  81 :       RETURN
  82 : *
  83 : *= ERREUR : Fonctions de forme non encore implementees
  84 : *
  85 :  100  CONTINUE
  86 :       MOTERR(1:4) = NOMTP(NEF)
  86 :       CALL ERREUR(68)                                                           
  87 :       SEGSUP SHXX , HYBSTO
  88 :       SEGDES MINTE
  89 :       IPTHYB = 0
  90 :       RETURN                                                                    
  91 :       END                                                                       
  92 :  
  93 :  
  94 :