Esope source files

   1 : C DOXE      SOURCE    CHAT      05/01/12    22:54:35     5004           
   2 : 
   3 : C=======================================================================
   4 : C=                            D O X E                                  =
   5 : C=                            -------                                  =
   6 : C=                                                                     =
   7 : C=  Fonction :                                                         =
   8 : C=  ----------                                                         =
   9 : C=   Recuperation des coordonnees GLOBALES des noeuds d'un element     =
  10 : C=                                                                     =
  11 : C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  12 : C=  ------------                                                       =
  13 : C=   XCOORD     (E)  Coordonnees de TOUS les POINTS (vient de SMCOORD) =
  14 : C=   IDIM       (E)  Dimension du probleme traite (1 a 3)              =
  15 : C=   NBNN       (E)  Nombre de NOEUDS dans l'element                   =
  16 : C=   IELEM      (E)  Numero (local) de l'element dans le tableau NUM   =
  17 : C=   XE(3,NBNN) (S)  Coordonnees des NBNN noeuds de l'element IELEM    =
  18 : C=                                                                     =
  19 : C=  Remarque :  Quelle que soit la dimension IDIM, les dimensions du   =
  20 : C=  ----------  tableau XE sont toujours XE(3,NBNN).                   =
  21 : C=======================================================================
  22 : 
  23 :       SUBROUTINE DOXE (XCOOR,IDIM,NBNN,NUM,IELEM,XE)
  24 : 
  25 :       IMPLICIT INTEGER(I-N)
  26 :       IMPLICIT REAL*8 (A-H,O-Z)
  27 : 
  28 : -INC CCREEL
  29 : 
  30 :       DIMENSION XCOOR(*),NUM(NBNN,*),XE(3,*)
  31 : 
  32 :       idimp1=IDIM+1
  33 :       IF (IDIM.EQ.1) THEN
  34 :         DO i=1,NBNN
  35 :           j=idimp1*(NUM(i,IELEM)-1)
  36 :           XE(1,i)=XCOOR(j+1)
  37 :           XE(2,i)=XZero
  38 :           XE(3,i)=XZero
  39 :         ENDDO
  40 :       ELSE IF (IDIM.EQ.2) THEN
  41 :         DO i=1,NBNN
  42 :           j=idimp1*(NUM(i,IELEM)-1)
  43 :           XE(1,i)=XCOOR(j+1)
  44 :           XE(2,i)=XCOOR(j+2)
  45 :           XE(3,i)=XZero
  46 :         ENDDO
  47 :       ELSE IF (IDIM.EQ.3) THEN
  48 :         DO i=1,NBNN
  49 :           j=idimp1*(NUM(i,IELEM)-1)
  50 :           XE(1,i)=XCOOR(j+1)
  51 :           XE(2,i)=XCOOR(j+2)
  52 :           XE(3,i)=XCOOR(j+3)
  53 :         ENDDO
  54 :       ENDIF
  55 : 
  56 :       RETURN
  57 :       END
  58 : 
  59 :  
  60 :