anchor.diff:  !  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,
  +         CALL svLS_LS_CREATE(svLS_ls_S(isolsc), svLSType, dimKry=Kspace,
        2      relTol=epstol(indx), 
        3      maxItr=maxIters 
        4      )
   
  -         CALL svLS_COMMU_CREATE(communicator_S1, MPI_COMM_WORLD)
  +         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_S1, communicator_S1, gnNo, nNo, 
  +            CALL svLS_LHS_CREATE(svLS_lhs_S(isolsc), communicator_S(isolsc), gnNo, nNo, 
        2         nnz_tot, ltg, colm, rowp, svLS_nFaces)
               
                 faIn = 1
  
positive.code:      &                         shpb,          shglb,     rowp,     
       &                         colm,          lhsS(1,j), 
       &                         solinc(1,isclr+5), tcorecpscal,
       &                         svLS_lhs_S1,   svLS_ls_S1, svls_nfaces)
                          
                          
                    endif         ! end of scalar type solve
  
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)
  