Esope source files

   1 : C PROSGA    SOURCE    CHAT      05/01/13    02:35:00     5004           
   2 :       SUBROUTINE PROSGA(ICM,XE,XEC,XX,X,SHPC,SGA,QSI,IG,NBNO)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z) 
   5 :       DIMENSION X(3),XE(3,*),XEC(3,2),XX(3,2),
   6 :      & BJ(3,3),SHPC(6,20,*),SGA(3,5,5),QSI(3,2),DT(3),SIB(3)
   7 :       DJ=0.D0
   8 :       IF(ICM.NE.0) THEN
   9 :       DO 10 I=1,3
  10 :       DT(I)=0.5D0*(XEC(I,2)-XEC(I,1))
  11 :       SIB(I)=DT(I)
  12 :       DJ=DJ+DT(I)*DT(I)
  13 :       DO 10 J=1,3
  14 :   10  BJ(I,J)=0.D0
  15 :       DO 20 I=1,3
  16 :       DO 20 J=1,3
  17 :       DO 20 K=1,NBNO
  18 :   20  BJ(I,J)=BJ(I,J)+XE(I,K)*SHPC(J+1,K,IG)
  19 : C --- LES DRIVEES
  19 :       CALL LINTRI(BJ,DT,X)
  20 :       DJ=SQRT(DJ)
  21 :       DO 30 I=1,3
  22 :   30  SIB(I)=SIB(I)/DJ
  23 :       END IF
  24 :       DO 2000 IS=1,3
  25 : C --- Coordonnees du point d'integration dans le systeme reel (x,y,z)
  26 :       SGA(IS,1,IG)=XX(IS,IG)
  27 : C --- Coordonnees du point d'integration dans le systeme de reference 
  28 : C     (qsi,eta,dze)
  29 :       SGA(IS,2,IG)=QSI(IS,IG)
  30 : C --- Derivee de (qsi,eta,dze) par rapport a miu
  31 :       SGA(IS,3,IG)=X(IS)
  32 : C --- Cosinus du cable avec les axes
  33 :       SGA(IS,4,IG)=SIB(IS)
  34 : C --- Jacobien
  35 :       SGA(IS,5,IG)=DJ
  36 :  2000 CONTINUE
  37 :       END
  38 :  
  39 :