If the solver supports complementarity constraints, which KINTRO.jl does, then you should use;
using JuMP
model = Model()
@variable(model, S_plus >= 0)
@variable(model, S_minus >= 0)
@expression(model, S, S_plus - S_minus)
@constraint(model, S_plus ⟂ S_minus)