This was done using the convention that for a well w in W, W.ctrls[w]->current < 0 means that it should be treated as shut. This allows us to keep the matrix structure related to wells constant, without needing to remake it when wells shut or open.