anchor.diff:            Integer label, iSYS(*)
   c === local variables
             integer tet_nodes(4)
  -          logical calculate_dirac
  -          integer temp
  +          integer swap_temporary
             real*8 point(3)
             integer i
   c === too bad
  
positive.code:               if (tet_nodes(i) .eq. nearest_node_index) then
                    calculate_dirac = .true.
                    
                    temp = tet_nodes(i)
                    tet_nodes(i) = tet_nodes(4)
                    tet_nodes(4) = temp
                    goto 10
                end if
            end do
  
negative.code: c === ./dirac.f      
            real*8 dirac_value
            integer nearest_node_index
            logical calculate_dirac
  c === ./lininterp.f
            real*8 linterpvalue
  c ==========      
  