Stop When Reaching to Zero

[duplicating99]

Hello everyone,

I'm working on repeat a model using Vensim PLE 6.4b (32-bit) from a paper I read, but the model I create does not match the result in that paper. In brief, I need every stocks to be zero when they reached to the zero, since I cannot get negative number in this model. I've tried to use if else then statement like this:

IF THEN ELSE(Iraq Tanks >= 0,-Iraq Tanks Attrition Rate, Iraq Tanks=0),

and I set zero as the minimum for those stocks but it won't work.

Is there any way that I could ask Vensim to avoid getting negative numbers? Thank you!

[tomfid]

Generally the constraint to impose is:

Code: Select all

Stock = INTEG( inflow - outflow, init stock)
init stock = ...
inflow = ...
outflow = MIN( desired outflow, max outflow) 
desired outflow = ...
max outflow = stock/min time to drain stock

Often min time to drain stock is equal to time step for convenience, so you can use the following shortcut:

Code: Select all

outflow = MIN(desired outflow, stock/TIME STEP)

[duplicating99]

Thank you Tom. I searched the forum with some of the key words, but it did not give me what I wanted. However, the website in that post really helps!

Any way, I appreciate your reply!

Generally the constraint to impose is:

Code: Select all

Stock = INTEG( inflow - outflow, init stock)
init stock = ...
inflow = ...
outflow = MIN( desired outflow, max outflow) 
desired outflow = ...
max outflow = stock/min time to drain stock

Often min time to drain stock is equal to time step for convenience, so you can use the following shortcut:

Code: Select all

outflow = MIN(desired outflow, stock/TIME STEP)

[daler6]

Another, perhaps simpler way to do this without adding in a new auxiliary variable (desired outflow) is to modify the equation for the stock from:

Stock = INTEG( inflow - outflow, init stock)

to

Stock = INTEG(max(inflow-outflow,-Stock/TIME STEP), init stock)

As light disadvantage to this is, in your diagram, you now you have an arrow going from TIME STEP directly to your Stock. This is a bit ugly and breaks the (my) rule of not having anything go into your Stock other than flows and initial values.

[Administrator]

Another, perhaps simpler way to do this without adding in a new auxiliary variable (desired outflow) is to modify the equation for the stock from:

Stock = INTEG( inflow - outflow, init stock)

to

Stock = INTEG(max(inflow-outflow,-Stock/TIME STEP), init stock)

As light disadvantage to this is, in your diagram, you now you have an arrow going from TIME STEP directly to your Stock. This is a bit ugly and breaks the (my) rule of not having anything go into your Stock other than flows and initial values.

I'd always advise against this approach. Keep things simple and clear, it means others can easily see what you are doing (and you can as well if you come back to the model after time away).

[tomfid]

An equally simple limit is:

outflow = MIN( desired outflow, stock/TIME STEP )

The assumption here is that the dynamics of the constraint are not of much interest, so that it's OK to assume that the minimum outflow time is equal to TIME STEP. This is usually fine.

This avoids having TIME STEP as an input to the stock on the diagram.

[tomfid]

Related: https://metasd.com/2019/11/the-importance-of-fonfoo/

[daler6]

Thanks for the advice. I guess it is better to have TIME STEP as an input to the flow on the diagram.