Sources Esope | Cast3M

drutra
C DRUTRA    SOURCE    CHAT      05/01/12    22:55:27     5004
C DRUTRA    SOURCE    JFG        02/07/04
      SUBROUTINE DRUTRA(S1,SEQ,BETINSA)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION S1(6)
C
      SEGMENT BETINSA
       REAL*8 RT,RC,YOUN,XNU,GFT,GFC,CAR
       REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
       INTEGER IFIS,IPLA,IBB,IGAU
      ENDSEGMENT
C
      ALT=1.D0
      BET=2.D0
C
      SX=S1(1)
      SY=S1(2)
      SZ=S1(3)
      SXY=S1(4)
      SXZ=S1(5)
      SYZ=S1(6)
      A = SX-SY
      B = SY-SZ
      C = SX-SZ
      DJ2=(A*A+B*B+C*C+SXY*SXY+SYZ*SYZ+SXZ*SXZ)
      DJ2=DJ2 * (1.D0/6.D0)
      DI1=SX+SY+SZ
C
      DJ2=3.D0*DJ2
      DJ2=SQRT(DJ2)
      SEQ=(DJ2+ALT*DI1)/BET
C
      RETURN
      END