Esope source files

   1 : C DISTRR    SOURCE    CHAT      05/01/12    22:51:39     5004           
   2 : 
   3 : C=======================================================================
   4 : C=                          D I S T R R                                =
   5 : C=                          -----------                                =
   6 : C=  Fonction :                                                         =
   7 : C=  ----------                                                         =
   8 : C=  Calcul de la coordonnee radiale en un point d'un element fini (en  =
   9 : C=  general un point de Gauss) pour les modes axisymetriques 1D et 2D, =
  10 : C=  et le mode spherique 1D.                                           =
  11 : C=                                                                     =
  12 : C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  13 : C=  ------------                                                       =
  14 : C=   XEL(3,NBNO)  (E)  Coordonnees GLOBALES des noeuds de l element    =
  15 : C=   SHP(6,NBNO)  (E)  Valeurs des fonctions de forme (SHP(1,*)) et de =
  16 : C=                     leurs derivees (SHP(2 a 6,*)) au point considere=
  17 : C=   NBNO         (E)  Nombre de noeuds de l element fini              =
  18 : C=   RGau         (S)  Coordonnee radiale au point considere (Gauss)   =
  19 : C=======================================================================
  20 : 
  21 :       SUBROUTINE DISTRR (XEL,SHP,NBNO,RGau)
  22 : 
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8 (A-H,O-Z)
  25 : 
  26 : -INC CCREEL
  27 : 
  28 :       DIMENSION XEL(3,*),SHP(6,*)
  29 : 
  30 :       RGau=XZero
  31 :       DO i=1,NBNO
  32 :         RGau=RGau+SHP(1,i)*XEL(1,i)
  33 :       ENDDO
  34 : 
  35 :       RETURN
  36 :       END
  37 : 
  38 :  
  39 :