anchor.diff:            implicit none
             Include 'assemble.fd'
   
  -          Real*8  x, y, z, DATA(*), Coef(9, *)
  +          Real*8  x, y, z, DATA(*), Coef(*)
             Integer label, iSYS(*)
   
             iSYS(1) = 1
             iSYS(2) = 1
   
  -          Dbc = BC_NEUMANN
  +c          Dbc = BC_NEUMANN
  +          Dbc = BC_DIRICHLET
   
   c === \int\limits_{\partial\Omega} g_N d S = 1.
   c === where \Omega = [0, 1]^3 
   c === and g_N = const
   
  -          Coef(1, 1) = -1d0 / 6d0
  +c          Coef(1, 1) = -1d0 / 6d0
  +          if (label.eq.1)  Then
  +              Coef(1) = -1d0 / 3.14159265359
  +              Dbc = BC_NEUMANN
  +          end if
  +          if (label.eq.2)  Then
  +              Coef(1) = 10d0
  +              Dbc = BC_DIRICHLET
  +          end if
  +          if (label.eq.0) Then
  +              Coef(1) = 1d0
  +              write(*,*) "label = 0 in Dbc"
  +              Dbc = BC_DIRICHLET
  +          end if
   
             Return
         End
  
positive.code:               end if
            End do
  
            Return
        End
  
  
negative.code: 
                    ibc = Dbc(x, y, z, lbF(k), DATAFEM, iSYS, eBC)
                     
              If(ibc .EQ. BC_NEUMANN) Then
                    label = lbF(k)
  
                    Call fem3Dtri(XYP(1, k), XYP(1, l), XYP(1, m),
  