anchor.diff:                 svLS_nFaces = 1   !not sure about this...looks like 1 means 0 for array size issues
               END IF
   
  -            CALL svLS_LHS_CREATE(svLS_lhs, communicator, gnNo, nNo, 
  -     2         nnz_tot, ltg, colm, rowp, svLS_nFaces)
  -            
               faIn = 1
               facenNo = 0
               DO i=1, nshg
  
positive.code:                   gNodes(j) = i
                 END IF
                END DO
             
              CALL svLS_BC_CREATE(svLS_lhs_S(isolsc), faIn, facenNo, 
       2         1, BC_TYPE_Dir, gNodes, sV(1,:))
  
negative.code: 
           CALL svLS_COMMU_CREATE(communicator_S(isolsc), MPI_COMM_WORLD)
   
                faIn = 1
                facenNo = 0
                ib=5+isolsc
  