C PROSGA SOURCE CHAT 05/01/13 02:35:00 5004 SUBROUTINE PROSGA(ICM,XE,XEC,XX,X,SHPC,SGA,QSI,IG,NBNO) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) DIMENSION X(3),XE(3,*),XEC(3,2),XX(3,2), & BJ(3,3),SHPC(6,20,*),SGA(3,5,5),QSI(3,2),DT(3),SIB(3) DJ=0.D0 IF(ICM.NE.0) THEN DO 10 I=1,3 DT(I)=0.5D0*(XEC(I,2)-XEC(I,1)) SIB(I)=DT(I) DJ=DJ+DT(I)*DT(I) DO 10 J=1,3 10 BJ(I,J)=0.D0 DO 20 I=1,3 DO 20 J=1,3 DO 20 K=1,NBNO 20 BJ(I,J)=BJ(I,J)+XE(I,K)*SHPC(J+1,K,IG) C --- LES DRIVEES CALL LINTRI(BJ,DT,X) DJ=SQRT(DJ) DO 30 I=1,3 30 SIB(I)=SIB(I)/DJ END IF DO 2000 IS=1,3 C --- Coordonnees du point d'integration dans le systeme reel (x,y,z) SGA(IS,1,IG)=XX(IS,IG) C --- Coordonnees du point d'integration dans le systeme de reference C (qsi,eta,dze) SGA(IS,2,IG)=QSI(IS,IG) C --- Derivee de (qsi,eta,dze) par rapport a miu SGA(IS,3,IG)=X(IS) C --- Cosinus du cable avec les axes SGA(IS,4,IG)=SIB(IS) C --- Jacobien SGA(IS,5,IG)=DJ 2000 CONTINUE END