anchor.diff:  END DO
   
   NATOM = N - 1
  -WRITE(6,325) NATOM
  +WRITE(6,*) NATOM
   
   EPS = 10.0**(-NEPS)
   S1 = 0.0
  
positive.code:     E3 = E2 * E(K)
      S4 = S4 + E3
  
      DO L= 1, NATOM
          EL2 = E(L)**2
          IF (K .LT. L) THEN
              D2 = 0.0
  
negative.code:             B = 1.0 - FACT
  
              ! Calculation of the diagonal part of the Oseen tensor only
              DO I=1,IP
                  V=A+(DIS(I)*DIS(I)*B/D2)
                  T(K,L,I)=V*6.0/(8.0*DISTAN)
                  T(L,K,I)=T(K,L,I)
              END DO
          END IF
  
      END DO
  END DO
  