anchor.diff:        Subroutine Det(qw1,qw2,qw3,qw4,qw5,qw6,qw7,qw8,qw9,qw)
         implicit none
   
  -      Real*8 qw1,qw2,qw3,qw4,qw5,qw6
  +      Real*8 qw1, qw2, qw3, qw4, qw5, qw6
         Real*8 qw7,qw8,qw9,qw
   
         qw=(qw1*qw5*qw9+qw7*qw2*qw6+qw3*qw4*qw8-
  
positive.code:       iSYS(1) = 3
        iSYS(2) = 3
  
        if (label.EQ.1 .or. label.EQ.2 .or.
       & label.EQ.3) then
  
        Coef(1, 1) = DATA(1)
        Coef(2, 2) = DATA(4)
        Coef(3, 3) = DATA(6)
  
        Coef(1, 2) = DATA(2)
        Coef(2, 1) = DATA(2)
  
        Coef(1, 3) = DATA(3)
        Coef(3, 1) = DATA(3)
  
        Coef(3, 2) = DATA(5)
        Coef(2, 3) = DATA(5)
  
        Else
  
        Coef(1, 1) = DATA(label+7)
        Coef(2, 2) = DATA(label+7)
        Coef(3, 3) = DATA(label+7)
  
        Coef(1, 2) = 0D0
        Coef(2, 1) = 0D0
  
        Coef(1, 3) = 0D0
        Coef(3, 1) = 0D0
  
        Coef(3, 2) = 0D0
        Coef(2, 3) = 0D0
        End if
  
  
  
        Ddiff = TENSOR_SYMMETRIC
  
        Return
  
negative.code:             End if
  
  
          End if
        End do
  
  
  