Esope source files

   1 : C TSEG2A    SOURCE    CHAT      05/01/13    03:53:27     5004           
   2 :       SUBROUTINE TSEG2A(POINTS,QSIGAU,POIGAU,SHPTOT)                            
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 : ************************************************************************        
   6 : *                                                                               
   7 : *                             T S E G 2 A                                       
   8 : *                             -----------                                       
   9 : *                                                                               
  10 : * FONCTION: ON ACQUIERT:                                                        
  11 : * ---------    - LES COORDONNEES INTRINSEQUES                                   
  12 : *              - LES FONCTIONS DE FORME ET LEURS DERIVEES                       
  13 : *           POUR L'ELEMENT DE THERMIQUE SEG2.                                   
  14 : *                                                                               
  15 : * MODULES UTILISES:                                                             
  16 : * -----------------                                                             
  17 : *                                                                               
  18 : -INC CCOPTIO                                                                    
  19 : *                                                                               
  20 : * PARAMETRES:   (E)=ENTREE   (S)=SORTIE   (+ = CONTENU DANS UN COMMUN)          
  21 : * -----------                                                                   
  22 : *                                                                               
  23 : *     POINTS  (E)  CHAINE DE CARACTERES:                                        
  24 : *                  'NOEUDS' SI ON VEUT LES VALEURS AUX NOEUDS                   
  25 : *                  'GAUSS' POUR LES VALEURS AUX POINTS DE GAUSS.                
  26 : *     QSIGAU  (S) |                                                             
  27 : *     POIGAU  (S)  POIDS D'INTEGRATION,SI ON A 'GAUSS'.                         
  28 : *     SHPTOT  (S)  FONCTIONS DE FORME ET LEURS DERIVEES.                        
  29 : *                                                                               
  30 :       REAL*8 QSIGAU(2),POIGAU(2)                                                
  31 :       REAL*8 SHPTOT(6,2,2)                                                      
  32 :       CHARACTER*(*) POINTS                                                      
  33 : *                                                                               
  34 : * CONSTANTES:                                                                   
  35 : * -----------                                                                   
  36 : *                                                                               
  37 : *                                                                               
  38 :       PARAMETER ( O2=2.D0 )                                                     
  39 :       PARAMETER ( O1=1.D0 )                                                     
  40 :       PARAMETER ( O5=0.5D0 )                                                    
  41 :       PARAMETER ( O577=0.577350269189626D0 )                                    
  42 : *                                                                               
  43 : *                                                                               
  44 : * AUTEUR, DATE DE CREATION:                                                     
  45 : * -------------------------                                                     
  46 : *                                                                               
  47 : *     DENIS ROBERT,LE 03/02/88.                                                 
  48 : *                                                                               
  49 : * LANGAGE:                                                                      
  50 : * --------                                                                      
  51 : *                                                                               
  52 : *             FORTRAN77                                                         
  53 : *                                                                               
  54 : ************************************************************************        
  55 : *                                                                               
  56 :       IF (POINTS(1:5).EQ.'NOEUD') GOTO 66                                       
  57 :       IF (POINTS(1:5).NE.'GAUSS') GOTO 100                                      
  58 : *                                                                               
  59 : *     COORDONNEES DES 2 POINTS DE GAUSS                                         
  60 : *                                                                               
  61 :       NBPGAU=2                                                                  
  62 :       NBNO=2                                                                    
  63 :       QSIGAU(1)=-O577                                                           
  64 :       QSIGAU(2)= O577                                                           
  65 :       POIGAU(1)=O1                                                              
  66 :       POIGAU(2)=POIGAU(1)                                                       
  67 :       GOTO 2                                                                    
  68 : *                                                                               
  69 : *     COORDONNEES DES NOEUDS                                                    
  70 : *                                                                               
  71 :    66 CONTINUE                                                                  
  72 :       NBNO=2                                                                    
  73 :       NBPGAU=NBNO                                                               
  74 :       QSIGAU(1)=-O1                                                             
  75 :       QSIGAU(2)=O1                                                              
  76 :       POIGAU(1)=O1                                                              
  77 :       POIGAU(2)=POIGAU(1)                                                       
  78 : *                                                                               
  79 : *     FONCTIONS DE FORME ET LEURS DERIVEES                                      
  80 : *                                                                               
  81 :     2 DO 1 I=1,NBPGAU                                                           
  82 :       XI=QSIGAU(I)                                                              
  83 :       SHPTOT(1,1,I)=(O1-XI)/O2                                                  
  84 :       SHPTOT(1,2,I)=(O1+XI)/O2                                                  
  85 : *                                                                               
  86 :       SHPTOT(2,1,I)=-O5                                                         
  87 :       SHPTOT(2,2,I)= O5                                                         
  88 :  1    CONTINUE                                                                  
  89 :       RETURN                                                                    
  90 :  100  CONTINUE                                                                  
  90 :       CALL ERREUR(19)                                                           
  91 : *                                                                               
  92 :       END                                                                       
  93 :