Sources Esope | Cast3M

drucom
C DRUCOM    SOURCE    CHAT      05/01/12    22:55:19     5004
C DRUCOM    SOURCE    JFG        02/07/04
      SUBROUTINE DRUCOM(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
      BETA = 1.2D0
      ALC=(BETA-1.)/(2.*BETA-1.)
      BEC=BETA/(2.*BETA-1.)
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+ALC*DI1)/BEC
C
      RETURN
      END