anchor.diff:  END DO
   
   DO g = 1, ng
  -    vdum5 = 0.
  +      vdum5 = 0.
       DO h = 1, ng
           DO n = 1, nnod
               vdum5(n) = vdum5(n) + dsigs(n,h,g) * sf(n,h)
  
positive.code:     !Calculate amplitude function
      CALL point(t1, t2, hp, rho, A, beta, amp)
  
      fl = amp * f0
  
      WRITE(ounit,'(I4, F10.3, F10.4, ES15.4, 12F9.2)') step, t2, rho/tbeta, amp, (bpos(n), n = 1, nb)
  
      IF (stime) EXIT
  
negative.code: 
  xrho = Integrate(vdum4 - vdum3) / F2
  
  END SUBROUTINE kinet_par
  
  
  