anchor.diff:            iSYS(1) = 1
             iSYS(2) = 1
   
  -          coef(1, 1) = -thermo_coefficient
  +          coef(1, 1) = thermo_coefficient
   
             Ddiff = TENSOR_SCALAR
   
  
positive.code: 
                    ibc = Dbc(x, y, z, lbF(k), DATAFEM, iSYS, eBC)
                     
                    label = lbF(k)
  
                    Call fem3Dtri(XYP(1, k), XYP(1, l), XYP(1, m),
  
negative.code:      &                fixed_elements,                ! IEV
       &                work_buffer,                   ! iW
       &                work_buffer,                   ! iW
       &                "mesh.ani")   
  c     &                "../../../../data/cube.ani")   
          write (*,*) "loadMani has finished!"
  
          do i = 1, total_boundary_faces
  