anchor.diff:  
            CALL svLS_COMMU_CREATE(communicator_S(isolsc), MPI_COMM_WORLD)
    
  -               svLS_nFaces = 1   !not sure about this...should try it with zero
  -
  -            CALL svLS_LHS_CREATE(svLS_lhs_S(isolsc), communicator_S(isolsc), gnNo, nNo, 
  -     2         nnz_tot, ltg, colm, rowp, svLS_nFaces)
  -            
                 faIn = 1
                 facenNo = 0
                 ib=5+isolsc
  
positive.code:                   IF (.NOT.BTEST(iBC(i),5)) sV(3,j) = 1D0
                 END IF
              END DO
              CALL svLS_BC_CREATE(svLS_lhs, faIn, facenNo, 
       2         nsd, BC_TYPE_Dir, gNodes, sV)
              DEALLOCATE(gNodes)
  
negative.code:                   gNodes(j) = i
                 END IF
                END DO
  
              svLS_nFaces = 1   !not sure about this...should try it with zero
              CALL svLS_LHS_CREATE(svLS_lhs_S(isolsc), communicator_S(isolsc), gnNo, nNo,
       2         nnz_tot, gNodes, colm, rowp, svLS_nFaces)
             
              CALL svLS_BC_CREATE(svLS_lhs_S(isolsc), faIn, facenNo, 
       2         1, BC_TYPE_Dir, gNodes, sV(1,:))
  