anchor.diff:       &                         shpb,          shglb,     rowp,     
        &                         colm,          lhsS(1,j), 
        &                         solinc(1,isclr+5), tcorecpscal,
  -     &                         svLS_lhs_S1,   svLS_ls_S1, svls_nfaces)
  +     &                         svLS_lhs_S(isclr),   svLS_ls_S(isclr), svls_nfaces)
                           
                           
                     endif         ! end of scalar type solve
  
positive.code: !  reltol for the GMRES is the stop criterion 
  ! also using Kspaceand maxIters from setup for ACUSIM
  !
           CALL svLS_LS_CREATE(svLS_ls_S1, svLSType, dimKry=Kspace,
       2      relTol=epstol(indx), 
       3      maxItr=maxIters 
       4      )
  
           CALL svLS_COMMU_CREATE(communicator_S1, MPI_COMM_WORLD)
   
                 svLS_nFaces = 1   !not sure about this...should try it with zero
  
              CALL svLS_LHS_CREATE(svLS_lhs_S1, communicator_S1, gnNo, nNo, 
       2         nnz_tot, ltg, colm, rowp, svLS_nFaces)
              
                faIn = 1
  
negative.code:                END IF
                END DO
             
              CALL svLS_BC_CREATE(svLS_lhs_S(isolsc), faIn, facenNo, 
       2         1, BC_TYPE_Dir, gNodes, sV(1,:))
              DEALLOCATE(gNodes)
              DEALLOCATE(sV)
  