Soner Candas

Such combined optimization can already be realized in this model by setting a price for the CO2 commodity on the "Commodity" sheet

A new numeric, "timestep indexer" `pyomo.Set` can be introduced such as: `m.timestepsnum = range(1,len(timesteps)) ` `m.tnum= pyomo.Set(initialize=m.timestepsnum,ordered=True)` and then the required storage/DSM rules and corresponding pyomo variables such as `m.e_sto_con`...

Looks like you have two (storage?) components with the same indices… can you please share your storage sheet table, in case that is the only part you changed?