Esope source files

   1 : C DRUCO1    SOURCE    CHAT      05/01/12    22:55:14     5004           
   2 : C DRUCO1    SOURCE    JFG        02/07/04               
   3 :       SUBROUTINE DRUCO1(S1,DF,BETINSA)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 :       DIMENSION S1(6),DF(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,IGAU1
  13 :       ENDSEGMENT
  14 : C
  15 :       BETA = 1.2D0
  16 :       ALC=(BETA-1.)/(2.*BETA-1.)
  17 :       BEC=BETA/(2.*BETA-1.)
  18 : C
  19 :       SX=S1(1)
  20 :       SY=S1(2)
  21 :       SZ=S1(3)
  22 :       SXY=S1(4)
  23 :       SXZ=S1(5)
  24 :       SYZ=S1(6)
  25 :       A = SX-SY
  26 :       B = SY-SZ
  27 :       C = SX-SZ
  28 :       DJ2=(A*A+B*B+C*C+SXY*SXY+SYZ*SYZ+SXZ*SXZ)
  29 :       DJ2=DJ2 * (1.D0/6.D0)
  30 :       DJ2=3.D0*DJ2
  31 :       DJ2=SQRT(DJ2)
  32 : C
  33 :       DF(1) = ((2.*SX) - SY -SZ)/2./DJ2 + ALC
  34 :       DF(2) = (-SX +(2.*SY) -SZ)/2./DJ2 + ALC
  35 :       DF(3) = (-SX -SY +(2.*SZ))/2./DJ2 + ALC
  36 :       DF(4) = 6. * SXY/2./DJ2
  37 :       DF(5) = 6. * SXZ/2./DJ2
  38 :       DF(6) = 6. * SYZ/2./DJ2
  39 : C
  40 :       RETURN
  41 :       END
  42 : 
  43 :  
  44 :  
  45 :