anchor.diff:       &                    face_nodes,
        &                    face_material_labels) 
   
  -      write (*,*) "GMVscalarTet has finished!"
  +        write (*,*) "GMVscalarTet has finished!"
   
   c ===   I may be wrong about that      
         end program thermoconductivity
  
positive.code: C ======================================================================
  
        Integer Function Dreact(x, y, z, label, DATA, iSYS, Coef)
        implicit none
        include 'fem3Dtet.fd'
  
        Real*8  x, y, z, DATA(*), Coef(*)
        Integer label, iSYS(*)
  
        iSYS(1) = 1
        iSYS(2) = 1
  c     A == 0
        Coef(1) = 0D0
        Dreact = TENSOR_SCALAR
  
        Return
        End
  
  C ======================================================================
  
negative.code:      &                 sle_right_hand_side,                              ! b
       &                 sle_solution)                                     ! x
  
          write (*,*) "LSolver has finished!"
  
  c ===   write results to VTK file
          call GMVscalarTet(sle_solution, 
  