C DOXE SOURCE CHAT 05/01/12 22:54:35 5004 C======================================================================= C= D O X E = C= ------- = C= = C= Fonction : = C= ---------- = C= Recuperation des coordonnees GLOBALES des noeuds d'un element = C= = C= Parametres : (E)=Entree (S)=Sortie = C= ------------ = C= XCOORD (E) Coordonnees de TOUS les POINTS (vient de SMCOORD) = C= IDIM (E) Dimension du probleme traite (1 a 3) = C= NBNN (E) Nombre de NOEUDS dans l'element = C= IELEM (E) Numero (local) de l'element dans le tableau NUM = C= XE(3,NBNN) (S) Coordonnees des NBNN noeuds de l'element IELEM = C= = C= Remarque : Quelle que soit la dimension IDIM, les dimensions du = C= ---------- tableau XE sont toujours XE(3,NBNN). = C======================================================================= SUBROUTINE DOXE (XCOOR,IDIM,NBNN,NUM,IELEM,XE) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) -INC CCREEL DIMENSION XCOOR(*),NUM(NBNN,*),XE(3,*) idimp1=IDIM+1 IF (IDIM.EQ.1) THEN DO i=1,NBNN j=idimp1*(NUM(i,IELEM)-1) XE(1,i)=XCOOR(j+1) XE(2,i)=XZero XE(3,i)=XZero ENDDO ELSE IF (IDIM.EQ.2) THEN DO i=1,NBNN j=idimp1*(NUM(i,IELEM)-1) XE(1,i)=XCOOR(j+1) XE(2,i)=XCOOR(j+2) XE(3,i)=XZero ENDDO ELSE IF (IDIM.EQ.3) THEN DO i=1,NBNN j=idimp1*(NUM(i,IELEM)-1) XE(1,i)=XCOOR(j+1) XE(2,i)=XCOOR(j+2) XE(3,i)=XCOOR(j+3) ENDDO ENDIF RETURN END