@Henrique_Becker
Yes, both sides are scalar
Now the constraint has the following error: Bounds Error: attempet to access 79 elements Array{Int64,1} at index [1.5]
I think the function is correct, but there is an inconsistency in the size of the matrix that I cannot solve.
@constraint(Model1,sum(Q[a,w]*x[a,t] for a in A)-sum(s[a,t] for a in A)<=(RHE[w]+OHE[w])*NE[w]-WH[w,t])