Di = max[(POBi * (PMIN * 0.5 * (REMi + IEMi)) - Gti), 0]