Esope source files

   1 : C DRUTRA    SOURCE    CHAT      05/01/12    22:55:27     5004           
   2 : C DRUTRA    SOURCE    JFG        02/07/04               
   3 :       SUBROUTINE DRUTRA(S1,SEQ,BETINSA)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 :       DIMENSION S1(6)
   8 : C
   9 :       SEGMENT BETINSA
  10 :        REAL*8 RT,RC,YOUN,XNU,GFT,GFC,CAR
  11 :        REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
  12 :        INTEGER IFIS,IPLA,IBB,IGAU
  13 :       ENDSEGMENT
  14 : C
  15 :       ALT=1.D0
  16 :       BET=2.D0
  17 : C
  18 :       SX=S1(1)
  19 :       SY=S1(2)
  20 :       SZ=S1(3)
  21 :       SXY=S1(4)
  22 :       SXZ=S1(5)
  23 :       SYZ=S1(6)
  24 :       A = SX-SY
  25 :       B = SY-SZ
  26 :       C = SX-SZ
  27 :       DJ2=(A*A+B*B+C*C+SXY*SXY+SYZ*SYZ+SXZ*SXZ)
  28 :       DJ2=DJ2 * (1.D0/6.D0)
  29 :       DI1=SX+SY+SZ
  30 : C
  31 :       DJ2=3.D0*DJ2
  32 :       DJ2=SQRT(DJ2)
  33 :       SEQ=(DJ2+ALT*DI1)/BET
  34 : C
  35 :       RETURN
  36 :       END
  37 : 
  38 :  
  39 :  
  40 :