egypt1

# egypt1 - # = # EGYPT: GAMS Egyptian Fertilizer Model # = #...

This preview shows pages 1–3. Sign up to view the full content.

# ================================================= # EGYPT: GAMS Egyptian Fertilizer Model # ================================================= # # ORDERED PAIRS VERSION # Source: "On the Development of a General Algebraic Modeling # System in a Strategic Planning Environment" by Johannes Bisschop # and Alexander Meeraus # -------------------------------------- # Sets # -------------------------------------- # set center; # Locations from which final product may be shipped s set port within center; # Locations at which imports can be received s set plant within center; # Locations of plants s set region; # Demand regions s set unit; # Productive units s set proc; # Processes s set nutr; # Nutrients s set c_final; # Final products (fertilizers) s set c_inter; # Intermediate products s set c_ship within c_inter; # Intermediates for shipment s set c_raw; # Domestic raw materials and miscellaneous inputs s set commod := c_final union c_inter union c_raw; # All commodities # -------------------------------------- # Parameters # -------------------------------------- # param cf75 {region,c_final} >= 0; # Consumption of fertilizer 1974-75 (1000 tpy) param fn {c_final,nutr} >= 0; # Nutrient content of fertilizers param cn75 {r in region, n in nutr} := sum {c in c_final} cf75[r,c] * fn[c,n]; # Consumption of nutrients 1974-75 (1000 tpy) param road {region,center} >= 0; # Road distances param rail_half {plant,plant} >= 0;

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
param rail {p1 in plant, p2 in plant} := rail_half[p1,p2] + rail_half[p2,p1]; p # Interplant rail distances (kms) param impd_barg {plant} >= 0; param impd_road {plant} >= 0; # Import distances (kms) by barge and road param tran_final {pl in plant, r in region} := if road[r,pl] > 0 then (.5 + .0144 * road[r,pl]) else 0; param tran_import {r in region, po in port} := if road[r,po] > 0 then (.5 + .0144 * road[r,po]) else 0; param tran_inter {p1 in plant, p2 in plant} :=
This is the end of the preview. Sign up to access the rest of the document.

## egypt1 - # = # EGYPT: GAMS Egyptian Fertilizer Model # = #...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online