anchor.diff:            implicit none
             include 'fem3Dtet.fd'
   
  -          Real*8  x, y, z, DATA(*), Coef(9, *)
  +          Real*8  x, y, z, DATA(*)
  +
  +c          real*8 Coef(9, *)
  +          real*8 coef(*)
  +
             Integer label, iSYS(*)
   
             Integer i, j
  
positive.code:           Integer label, iSYS(*)
  c === local variables
            integer tet_nodes(4)
            logical calculate_dirac
            integer temp
            real*8 point(3)
            integer i
  c === too bad
  
negative.code: c === ./dirac.f      
            real*8 dirac_value
            integer nearest_node_index
            logical calculate_dirac
  c === ./lininterp.f
            real*8 linterpvalue
  c ==========      
  