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:       Real*8   C(3, 3), G(3), XYP(3, 4)
        Real*8   x, y, z, eBC(1)
  
        Integer  i,j,k,l,m, ir, ic, label, ibc
  
        Integer  iref(5), ip(4)
        DATA     iref/1,2,3,4,1/
  
  
  C ======================================================================
  
negative.code:      & -1D0,XY3(1),XY3(2),-1D0,D3)
  
        if (D.NE.0) then
          A = D1 / D
          B = D2 / D
          C = D3 / D
          in = 1
        else
          stop 'determinant=0'
        end if
  
        Else if (t1.EQ.0) then
  