indus89.gms : Indus Basin Water Resource Model

Description

This file contains the basic data and definition of the surface water
system. Data is complete for year 1988. Some parameters could be
computed for future years using growth rates provided in this file,
others had to be estimated and entered. Enter the year for which the
setup is desired in set isr (set isr should have only one entry).


Reference

  • Ahmad, M, and Kutcher, G P, Irrigation Planning with Environmental Considerations - A Case Study of Pakistans's Indus Basin. Tech. rep., The World Bank, 1992.

Large Model of Type : LP


Category : GAMS Model library


Main file : indus89.gms

$TITLE Indus Basin Model Revised - IBMR  (INDUS89,SEQ=181)
$ontext
This file contains the basic data and definition of the surface water
system. Data is complete for year 1988. Some parameters could be
computed for future years using growth rates provided in this file,
others had to be estimated and entered. Enter the year for which the
setup is desired in set isr (set isr should have only one entry).


Ahmad, M, and Kutcher, G P, Irrigation Planning with Environmental
Considerations - A Case Study of Pakistans's Indus Basin. Tech. rep.,
The World Bank, 1992.

changes for year 2000 runs
  growth of crop yields set to a maximum of 3 %
  insert this line after growthcy parameter
  growthcy(c,z)$(growthcy(c,z) gt 3) = 3.0 ;

$offtext

$stitle set definition
 sets
  z     agroclimatic zones /
                 nwfp    northwest frontier mixed cropping
                 pmw     punjab wheat-mixed cropping
                 pcw     punjab cotton-wheat
                 psw     punjab sugarcane-wheat
                 prw     punjab rice-wheat
                 scwn    sind cotton-wheat north
                 srwn    sind rice-wheat north
                 scws    sind cotton wheat south
                 srws    sind rice-wheat south  /
  pv        provinces and country /nwfp, punjab, sind, pakistan /
  pv1(pv)   provinces             /nwfp, punjab, sind/
  pv2(pv)   punjab and sind       /punjab, sind/
  pvz(pv,z) province to zone map  /nwfp.nwfp
                                   punjab.(pmw,pcw,psw,prw)
                                   sind.  (scwn,scws,srwn,srws) /
  cq  crop and livestock products
          /basmati,   irri,        cotton,     rab-fod,
           gram,      maize,       mus+rap,    kha-fod,
           sc-mill,   sc-gur,      wheat,      orchard,
           potatoes,  onions,      chilli
           cow-milk,  buff-milk,   meat /
  cc(cq)   consumable commodities /
           basmati,   irri,        gram
           maize  ,   mus+rap,     sc-gur
           wheat  ,   potatoes,    onions
           chilli    /
  c(cq)   crops  /
                 basmati   rice crop
                 irri      rice crop
                 cotton
                 rab-fod   fodder crop
                 gram
                 maize
                 mus+rap
                 kha-fod   fodder crop
                 sc-gur    sugarcane processed at the farm
                 sc-mill   sugarcane for mill
                 wheat
                 orchard
                 potatoes
                 onions
                 chilli    /
  cf(c)   fodder crops     /rab-fod, kha-fod /
  cnf(c)  non-fodder crops

  t       technology      /bullock,   semi-mech /
  s       sequence        /standard   standard sequence
                           la-plant   late planting
                           el-plant   early planting
                           qk-harv    quick harvesting /
  w     water stress level /
                           standard   no stress
                           light      light stress
                           heavy      heavy stress
                           january    water stress in january/
  g     ground water quality types/ fresh, saline/
  gf(g) fresh ground water sub-zone/fresh /
  gs(g) saline ground water sub-zone/saline/
  t1    sub zones by gw quality/ fresh, saline, total /
  r1    resources /
          cca       culturable commanded area of the canal
          ccap      canal capacity at the canal head
          ceff      canal efficiency from barrage to the water course head
          wce-r     water course command efficiency in rabi season
          wce-k     water course command efficiency in kharif season
          flde      field efficiency
          farmpop   farm population in the irrigated
          farmhh    number of agricultural households
          tractors  tractor population in the irrigated area
          tubewells number of tubewells
          twc       existing private tubewell capacity
          bullocks
          cows
          buffalos                         /
  dc(r1)  characteristics of canal command  /
          cca       culturable commanded area of the canal
          ccap      canal capacity at the canal head
          ceff      canal efficiency from barrage to the water course head
          wce-r     water course command efficiency in rabi season
          wce-k     water course command efficiency in kharif season
          flde      field efficiency /
  sa      subareas    /s1*s4/
  wce(dc)       watercourse efficiencies /wce-r, wce-k /
  m1      months and seasons /jan,feb,mar,apr,may,jun,
                              jul,aug,sep,oct,nov,dec,rabi,kharif,annual/
  m(m1)   months      /jan,feb,mar,apr,may,jun
                       jul,aug,sep,oct,nov,dec/
  wcem(wce,m)   mapping from season to months for watercourse efficiencies  /
                       wce-r. (oct,nov,dec,jan,feb,mar)
                       wce-k. (apr,may,jun,jul,aug,sep)/
  sea(m1) seasons     /rabi, kharif /
  seam(sea,m)         mapping from seasons to months/
                      rabi.  (oct,nov,dec,jan,feb,mar)
                      kharif.(apr,may,jun,jul,aug,sep)/
    sea1 /rabi, kharif, annual/
    sea1m(sea1,m) /rabi.  (oct,nov,dec,jan,feb,mar)
                   kharif.(apr,may,jun,jul,aug,sep)/
  ci      crop input outputs    / straw-yld, nitrogen, phosphate , seed /
  p2(ci)                        / nitrogen, phosphate  /
  a       animal types          / cow, bullock, buffalo /
  ai      animals input output  /tdn, dp, labor, cow-milk, buff-milk, meat /
  q(cq)   livestock commodities /cow-milk   milk from cattle cow
                                 buff-milk  milk from buffalo cow
                                 meat       from cows buffaloes and bullocks/
  nt      nutrients for animals /tdn  total digestible nutrients
                                 dp   digestible protein   /
  is      irrigation system scenarios /1980*2000  /
  ps      price scenarios / 87-88  /
*--
* change the set isr to setup data for desired year.
*
      isr(is) irrigation system scenario for this run /1988 /
             ;
 scalar baseyear  base year for crop yields /1988/;

           cnf(c) = yes;  cnf(cf) = no ;
           pvz("pakistan",z) = yes ;
           sea1m("annual",m) = yes;

*  parameters to export data for zone model
*&&z set set1 /tdn,dp,labor,cow-milk,buff-milk,meat,fix-cost/
*&&z     ;
* following parameters are to store data for export to the zone models.
*&&z parameter
*&&z   zone1xxxxx(z,c,t,s,w,m)  bullock requirements (bullock pair hrs per month)
*&&z   zone2xxxxx(z,a,set1)     input output coefficients for livestock
*&&z   zone3xxxxx(z,cq, * )     demand data
*&&z   zone4xxxxx(z,c,ci)       crop input output
*&&z   ;
$stitle crop data

    table land(c,z,t,s,w,m) land occupation by month

                                                                      jan feb mar apr may jun jul aug sep oct nov dec
 (basmati,irri).(pmw,pcw,psw,prw).(bullock,
                        semi-mech).standard.standard                                        1   1   1   1   1  .5
 irri.   scwn.(bullock,semi-mech).standard.standard                                     1   1   1   1   1   1
 irri.   srwn.(bullock,semi-mech).standard.standard                                    .5   1   1   1   1   1  .5
 irri.   scws.(bullock,semi-mech).standard.standard                                .5   1   1   1   1   1   1
 irri.   srws.(bullock,semi-mech).standard.standard                                     1   1   1   1   1   1

 maize.  nwfp.bullock.            standard.standard                                         1   1   1   1   1
 maize. (pcw,psw,prw).bullock.    standard.standard                                        .5   1   1   1   1   1  .5
 maize. (pcw,psw,prw).semi-mech.  standard.standard                                             1   1   1   1   1  .5
 maize.   scwn.      (bullock,semi-mech).standard.standard                                 .5   1   1   1   1  .5
 maize.   scws.      (bullock,semi-mech).standard.standard                             .5   1   1   1   1   1  .5
 maize.  nwfp.semi-mech.standard.standard                                                  .5   1   1   1   1

 mus+rap.nwfp.(bullock,semi-mech).       standard.standard                                      1   1   1   1   1   1
 mus+rap.(pmw,pcw).(bullock,semi-mech).  standard.standard              1   1   1                           1   1   1
 mus+rap. psw.     (bullock,semi-mech).  standard.standard              1   1  .5                      .5   1   1   1
 mus+rap. prw.     (bullock,semi-mech).  standard.standard              1  .5                           1   1   1   1
 mus+rap.(scwn,srwn,srws).(bullock,semi-mech).standard.standard         1   1  .5                       1   1   1   1
 mus+rap.scws.(bullock,semi-mech).       standard.standard              1   1   1                      .5   1   1   1

 (sc-gur,sc-mill).(nwfp,pmw,pcw,psw,prw, scwn,srwn,scws,srws).
                 (bullock,semi-mech).standard.standard                  1   1   1   1   1   1   1   1   1   1   1   1
 kha-fod.(nwfp,pmw,pcw).(bullock,semi-mech).standard.standard                  .5  .5  .5   1   1   1  .5  .5
 kha-fod.psw.(bullock,semi-mech).        standard.standard                     .5  .5  .5   1   1  .5  .5
 kha-fod.prw.(bullock,semi-mech).        standard.standard                     .5  .5  .5   1   1   1  .5  .5
 kha-fod.(scwn).     (bullock,semi-mech).standard.standard                     .5  .5   1   1  .5   1   1 .75
 kha-fod.scws.       (bullock,semi-mech).standard.standard                     .5  .5   1   1  .5   1   1   1
 kha-fod.(srwn,srws).(bullock,semi-mech).standard.standard                     .5  .5   1   1   1   1   1   1
 kha-fod.srws.(bullock,semi-mech).       la-plant. standard                             1   1   1   1   1   1

 rab-fod.nwfp.bullock.                   standard.standard              1   1   1   1   1              .5   1   1   1
 rab-fod.(nwfp,pcw).semi-mech.           standard.standard              1   1   1   1   1             .25   1   1   1
 rab-fod.pmw. bullock.                   standard.standard              1   1   1   1   1                   1   1   1
 rab-fod.pcw. bullock.                   standard.standard              1   1   1   1   1              .5   1   1   1
 rab-fod.psw. bullock.                   standard.standard              1   1   1   1  .5                   1   1   1
 rab-fod.prw. bullock.                   standard.standard              1   1   1   1   1                   1   1   1
 rab-fod.scwn.(bullock,semi-mech).       standard.standard              1   1   1   1                  .5   1   1   1
 rab-fod.srwn.(bullock,semi-mech).       standard.standard              1   1   1 .75                  .5   1   1   1
 rab-fod.scws.(bullock,semi-mech).       standard.standard              1   1   1   1 .25              .5   1   1   1
 rab-fod.srws.(bullock,semi-mech).       standard.standard              1   1   1   1  .5              .5   1   1   1
 rab-fod.pmw. semi-mech.                 standard.standard              1   1   1   1   1                   1   1   1
 rab-fod.psw. semi-mech.                 standard.standard              1   1   1   1  .5                 .25   1   1
 rab-fod.prw. semi-mech.                 standard.standard              1   1   1   1   1                  .5   1   1
 rab-fod.(srwn,srws).(bullock,semi-mech).standard.heavy                 1   1                          .5   1   1   1
 rab-fod.(srwn,srws).(bullock,semi-mech).standard. light                1   1   1                      .5   1   1   1

 cotton. pmw.(bullock,semi-mech).        standard.standard                              1   1   1   1   1   1  .5
 cotton. (pcw,psw,prw).(bullock,semi-mech).standard.standard                            1   1   1   1   1   1  .5
 cotton. scwn.bullock.                   standard.standard                         .5   1   1   1   1   1   1  .5
 cotton. scws.bullock.                   standard.standard                          1   1   1   1   1   1   1
 cotton. srws.bullock.                   standard.standard                         .5   1   1   1   1   1   1
 cotton. scwn.semi-mech.                 standard.standard                              1   1   1   1   1   1  .5
 cotton. scws.semi-mech.                 standard.standard                         .5   1   1   1   1   1   1
 cotton. srws.semi-mech.                 standard.standard                              1   1   1   1   1   1
 cotton. pcw. bullock.  el-plant. standard                                         .5   1   1   1   1   1   1   1  .5
 cotton. pcw. semi-mech.la-plant. standard                                                  1   1   1   1   1   1  .5

 gram.(nwfp,pmw,pcw).bullock.            standard.standard              1   1   1                      .5   1   1   1
 gram.(nwfp,pmw,pcw).semi-mech.          standard.standard              1   1   1                     .25   1   1   1
 gram.(psw,prw).bullock.                 standard.standard              1   1   1                          .5   1   1
 gram.(psw,prw).semi-mech.               standard.standard              1   1   1                         .25   1   1
 gram.(scwn,srwn,scws,srws).(bullock,semi-mech).standard.standard       1   1   1                     .25   1   1   1

 wheat.  (nwfp,pmw,pcw,psw).bullock.la-plant.(standard,light,
                                              heavy,january)            1   1   1   1                           1   1
 wheat.  prw. bullock.  la-plant.(standard,light,heavy,january)         1   1   1   1                               1
 wheat. (scwn,srwn).bullock.la-plant.(standard,light,heavy,january)     1   1   1  .5                           1   1
 wheat. (scws,srws).bullock.la-plant.(standard,light,heavy,january)     1   1   1 .25                           1   1
 wheat. (nwfp,pmw,pcw,psw).bullock.qk-harv.
         (standard,light,heavy,january)                                 1   1   1  .5                      .5   1   1

 wheat.  prw. bullock.  qk-harv.(standard,light,heavy,january)          1   1   1  .5                          .5   1
 wheat.  (scwn,srwn,scws,srws).bullock.qk-harv.(standard,light,
                                                heavy,january)          1   1   1                          .5   1   1

 wheat.  (nwfp,pmw,pcw,psw).bullock.standard.(standard,light,heavy,
                                                       january)         1   1   1   1                          .5   1
 wheat.  prw. bullock.  standard.(standard,light,heavy,january)         1   1   1   1                          .5   1
 wheat. (scwn,srwn).bullock.standard.(standard,light,heavy,january)     1   1   1  .5                      .5   1   1
 wheat. (scws,srws).bullock.standard.(standard,light,heavy,january)     1   1   1 .25                      .5   1   1

 wheat. (nwfp,pmw,pcw,psw).semi-mech.la-plant.(standard,light,heavy)    1   1   1   1                               1
 wheat. (nwfp,pmw,psw).semi-mech.la-plant.january                       1   1   1   1                               1
 wheat.  pcw.semi-mech.la-plant.january                                 1   1   1   1                               1
 wheat.  prw. semi-mech.la-plant.(standard,light,heavy,january)         1   1   1   1                               1
 wheat. (scwn,srwn).semi-mech.la-plant.
         (standard,light,heavy,january)                                 1   1   1  .5                          .5   1
 wheat. (scws,srws).semi-mech.la-plant.
         (standard,light,heavy,january)                                 1   1   1 .25                          .5   1

 wheat.  (nwfp,pmw,psw).semi-mech.qk-harv.
         (standard,light,heavy,january)                                 1   1   1  .5                     .25   1   1
 wheat.  (pcw,prw).semi-mech.qk-harv.
         (standard,light,heavy,january)                                 1   1   1  .5                          .5   1
 wheat.  (scwn,srwn,scws,srws).semi-mech.
          qk-harv.(standard,light,heavy,january)                        1   1   1                               1   1

 wheat.(nwfp,pmw,psw).semi-mech.standard.
         (standard,light,heavy,january)                                 1   1   1   1                          .5   1

 wheat.  (pcw,prw).semi-mech.standard.(standard,light,heavy,january)    1   1   1   1                          .5   1
 wheat.(scwn,srwn).semi-mech.standard.(standard,light,heavy,january)    1   1   1  .5                           1   1
 wheat.(scws,srws).semi-mech.standard.(standard,light,heavy,january)    1   1   1 .25                           1   1
 (orchard). (nwfp,pmw,pcw,psw,prw, scwn,srwn,scws,srws).
                 (bullock,semi-mech).standard.standard                  1   1   1   1   1   1   1   1   1   1   1   1
 potatoes.(scwn,scws,srwn,srws) .semi-mech.standard.standard            1   1                           .5  1   1   1
 potatoes.(nwfp,pmw,pcw,psw,prw).semi-mech.standard.standard            1   1   1   1
 onions. (nwfp,pmw,pcw,psw,prw).semi-mech.standard.standard             1   1   1   1   1
 onions. (scwn,scws,srwn,srws). semi-mech.standard.standard             1                               1   1   1   1
 chilli. (nwfp,pmw,pcw,psw,prw).semi-mech.standard.standard                     .5  1   1   1   1   1   1   1   .5
 chilli. (scwn,scws,srwn,srws). semi-mech.standard.standard             1   1   1  .5                           .5  1
       ;
 set tech(z,c,t,s,w)  technology availability indicator ;
     tech(z,c,t,s,w)$sum(m, land(c,z,t,s,w,m))  = yes;

$eject
    table bullock(c,z,t,s,w,m) bullock power requirements(bullock pair hours per month )

                                             jan  feb  mar  apr  may  jun  jul  aug  sep  oct  nov  dec
 basmati.pmw. bullock.  standard.standard                            22.0 17.2                 2.0 15.6
 basmati.pcw. bullock.  standard.standard                            22.0 17.2                 2.0 15.6
 basmati.psw. bullock.  standard.standard                            22.6 17.3                 2.0 17.4
 basmati.prw. bullock.  standard.standard                            32.5 24.6                 2.0 23.0
 basmati.pmw. semi-mech.standard.standard                                                          13.6
 basmati.pcw. semi-mech.standard.standard                                                          13.6
 basmati.psw. semi-mech.standard.standard                                                          15.4
 basmati.prw. semi-mech.standard.standard                                                          20.0

 rab-fod.nwfp.bullock.  standard.standard    3.0  3.0  3.0  3.0  2.0                 8.0 17.9  2.0  2.0
 rab-fod.pmw. bullock.  standard.standard    3.0  3.0  3.0  3.0  2.0                     17.9  7.6  2.0
 rab-fod.pcw. bullock.  standard.standard    2.0  4.0  4.0  2.0  2.0                10.1 13.0  7.6  2.0
 rab-fod.psw. bullock.  standard.standard    4.0  4.0  4.0  2.0  2.0                     16.5  9.7  2.0
 rab-fod.prw. bullock.  standard.standard    4.0  4.0  4.0  2.0  2.0                     12.7  9.4  2.0
 rab-fod.scwn.bullock.  standard.standard    2.0  2.0  2.0  1.0                      9.6 16.4  2.0  2.0
 rab-fod.srwn.bullock.  standard.standard    2.0  2.0  2.0  1.0                      8.9 14.8  2.0  2.0
 rab-fod.scws.bullock.  standard.standard    2.0  2.0  2.0  2.0  1.0                10.3 18.4  2.0  2.0
 rab-fod.srws.bullock.  standard.standard    2.0  2.0  2.0  2.0  1.0                10.7 18.3  2.0  2.0
 rab-fod.scwn.semi-mech.standard.standard    2.0  2.0  2.0  1.0                                2.0  2.0
 rab-fod.srwn.semi-mech.standard.standard    2.0  2.0  2.0  1.0                                2.0  2.0
 rab-fod.scws.semi-mech.standard.standard    2.0  2.0  2.0  2.0  1.0                           2.0  2.0
 rab-fod.srws.semi-mech.standard.standard    2.0  2.0  2.0  2.0  1.0                           2.0  2.0
 cotton. pmw. bullock.  standard.standard                       14.1  6.0            1.0  1.0  1.0  1.0
 cotton. pcw. bullock.  standard.standard                       17.1 15.2            1.0  1.0  1.0  1.0
 cotton. psw. bullock.  standard.standard                        9.0 13.0            1.0  1.0  1.0  1.0
 cotton. prw. bullock.  standard.standard                       10.3 17.1            1.0  1.0  1.0  1.0
 cotton. scwn.bullock.  standard.standard                   7.4 19.0  5.0            2.0  2.0  4.0
 cotton. scws.bullock.  standard.standard                  26.0       5.0            2.0  4.0
 cotton. srws.bullock.  standard.standard                   9.3 21.4  5.0            2.0  4.0
 cotton. (scwn,scws,srws).semi-mech.
                        standard.standard                             5.0
 cotton. pcw. bullock.  el-plant. standard                 16.0  4.0 12.3            1.0  1.0  1.0  1.0
 gram.   nwfp.bullock.  standard.standard                   7.0                      8.9 14.5
 gram.   pmw. bullock.  standard.standard              7.0                          10.7  5.6
 gram.   pcw. bullock.  standard.standard              7.0                          10.8  5.6
 gram.   psw. bullock.  standard.standard              7.0                               10.2  5.2
 gram.   prw. bullock.  standard.standard              7.0                               10.3  5.3
 gram.   scwn.bullock.  standard.standard              4.0  5.7                      6.4  8.1
 gram.   scws.bullock.  standard.standard              4.0  5.7                      6.4  8.1
 gram.   srwn.bullock.  standard.standard              4.9  6.4                      7.1  9.3
 gram.   srws.bullock.  standard.standard              4.9  6.4                      7.1  9.3

 irri.   pmw. bullock.  standard.standard                            18.9 19.2                 1.5 18.4
 irri.   pcw. bullock.  standard.standard                            18.9 19.2                 1.5 18.4
 irri.   psw. bullock.  standard.standard                            16.3 18.4                 1.5 20.2
 irri.   prw. bullock.  standard.standard                            32.6 24.0                 1.5 22.0
 irri.   scwn.bullock.  standard.standard                       10.5 18.3                 3.5 18.4
 irri.   srwn.bullock.  standard.standard                            10.7 18.6                13.4 10.0
 irri.   scws.bullock.  standard.standard                   9.1 10.3 17.1                 3.6 14.7
 irri.   srws.bullock.  standard.standard                       19.0 17.5                 3.7 10.0 10.0
 irri.   pmw. semi-mech.standard.standard                                                          16.4
 irri.   pcw. semi-mech.standard.standard                                                          16.4
 irri.   psw. semi-mech.standard.standard                                                          18.2
 irri.   prw. semi-mech.standard.standard                                                          20.0
 irri.   scwn.semi-mech.standard.standard                                                     16.4
 irri.   srwn.semi-mech.standard.standard                                                     10.0 10.0
 irri.   scws.semi-mech.standard.standard                                                     12.7
 irri.   srws.semi-mech.standard.standard                                                     10.0 10.0

 maize.  nwfp.bullock.  standard.standard                            41.4                 5.0
 maize.  pcw. bullock.  standard.standard                            10.8  4.5 14.2                 5.0
 maize.  psw. bullock.  standard.standard                            10.2  4.1 13.4                 5.0
 maize.  prw. bullock.  standard.standard                            10.2  4.1 13.4                 5.0
 maize.  scwn.bullock.  standard.standard                            15.1 23.9                 3.0
 maize.  scws.bullock.  standard.standard                       14.3 21.7                      3.5
 maize.  nwfp.semi-mech.standard.standard                             3.5
 maize.  pcw. semi-mech.standard.standard                                       3.0
 maize.  psw. semi-mech.standard.standard                                       3.0
 maize.  prw. semi-mech.standard.standard                                       3.0
 maize.  scwn.semi-mech.standard.standard                                                      3.0
 maize.  scws.semi-mech.standard.standard                                                      3.5

 mus+rap.nwfp.bullock.  standard.standard                                  7.4  8.0 14.5  1.0  1.0  1.0
 mus+rap.pmw. bullock.  standard.standard    1.0  1.0  1.5                               15.6  6.2  1.0
 mus+rap.pcw. bullock.  standard.standard    1.0  1.0  1.0                               10.8 10.1  1.0
 mus+rap.psw. bullock.  standard.standard    1.0  1.0  1.0                          10.2  9.3  1.0  1.0
 mus+rap.prw. bullock.  standard.standard    1.0  1.0                               15.6  1.0  1.0  1.0
 mus+rap.scwn.bullock.  standard.standard    1.0  1.0  1.0                          12.9  8.1  1.0  1.0
 mus+rap.srwn.bullock.  standard.standard    1.0  1.0  1.0                          13.7  7.5  1.0  1.0
 mus+rap.srws.bullock.  standard.standard    1.0  1.0  1.0                          13.7  7.5  1.0  1.0
 mus+rap.scws.bullock.  standard.standard    1.0  1.0  1.0                          12.8  6.9  1.0  1.0

 sc-gur. nwfp.bullock.  standard.standard   14.6 11.3 11.0       2.0  2.0            8.0 11.0 14.0 18.0
 sc-gur. (pmw,pcw,psw,prw).bullock.
                        standard.standard   17.8 11.6  8.2       1.8  1.8                     11.5 13.5
 sc-gur. scwn.bullock.  standard.standard    8.2  6.1 10.9            4.9
 sc-gur. srwn.bullock.  standard.standard    8.5  6.3 11.2            5.6
 sc-gur. scws.bullock.  standard.standard    9.1  7.4 12.1            5.8
 sc-gur. srws.bullock.  standard.standard    9.3  7.2 12.5            5.7
 sc-gur. nwfp.semi-mech.standard.standard    7.0  8.0  3.0       2.0  2.0           10.5  6.0 14.0 11.0
 sc-gur. (pmw,pcw,psw,prw).semi-mech.
                        standard.standard   11.0  6.5  2.5       1.8  1.8                     11.5 11.0
 sc-gur. scwn.semi-mech.standard.standard                             4.9
 sc-gur. srwn.semi-mech.standard.standard                             5.6
 sc-gur. scws.semi-mech.standard.standard                             5.8
 sc-gur. srws.semi-mech.standard.standard                             5.7

 sc-mill.nwfp.bullock.  standard.standard   18.0 18.0 14.0 12.0                     15.9 20.0 15.0 16.0
 sc-mill.(pmw,pcw,psw,prw).bullock.
                        standard.standard   16.8 15.1 15.2 12.0  8.0  1.8                     10.0 12.5
 sc-mill.scwn.bullock.  standard.standard    8.2  6.0 10.9            4.9
 sc-mill.srwn.bullock.  standard.standard    8.5  6.3 11.2            5.6
 sc-mill.scws.bullock.  standard.standard    9.1  7.4 12.1            5.8
 sc-mill.srws.bullock.  standard.standard    9.3  7.2 12.5            5.7
 sc-mill.nwfp.semi-mech.standard.standard    2.0  2.0                                     3.0
 sc-mill.(pmw,pcw,psw,prw).semi-mech.
                        standard.standard         2.5            1.8  1.8
 sc-mill.scwn.semi-mech.standard.standard                             4.9
 sc-mill.srwn.semi-mech.standard.standard                             5.6
 sc-mill.scws.semi-mech.standard.standard                             5.8
 sc-mill.srws.semi-mech.standard.standard                             5.7

 kha-fod.nwfp.bullock.  standard.standard             15.4  1.0  1.0 16.0  1.0 15.0  1.0  0.5
 kha-fod.pmw. bullock.  standard.standard             15.0  1.0  1.0 15.5  1.0 16.0  1.0  0.5
 kha-fod.pcw. bullock.  standard.standard             16.0  1.0  1.0 16.0  1.5 15.0  1.0  0.5
 kha-fod.psw. bullock.  standard.standard             17.0  1.0  1.0 16.5  1.0 16.0  1.0  0.5
 kha-fod.prw. bullock.  standard.standard             16.0  1.0  1.2 15.5  1.5 17.0  1.2  0.5
 kha-fod.scwn.bullock.  standard.standard             10.7  2.0  6.3  2.0  2.0  6.3  2.0  2.0
 kha-fod.scws.bullock.  standard.standard             10.3  2.0  6.4  2.0  2.0  6.4  2.0  2.0
 kha-fod.srwn.bullock.  standard.standard             11.3  2.0  8.7  2.0  2.0  5.8  2.0  2.0
 kha-fod.srws.bullock.  standard.standard             11.3  2.0  8.7  2.0  2.0  5.8  2.0  2.0
 kha-fod.srws.bullock.  la-plant. standard                      10.3  2.0  2.0  6.4  2.0  2.0

 rab-fod.srwn.bullock.  standard. heavy      2.0  2.0                                8.9 14.8  2.0  2.0
 rab-fod.srws.bullock.  standard. heavy      2.0  2.0                               10.7 18.3  2.0  2.0
 rab-fod.srwn.bullock.  standard. light      2.0  2.0  2.0                           8.9 14.8  2.0  2.0
 rab-fod.srws.bullock.  standard. light      2.0  2.0  2.0                          10.7 18.3  2.0  2.0
 rab-fod.srwn.semi-mech.standard. heavy      2.0  2.0                                          2.0  2.0
 rab-fod.srws.semi-mech.standard. heavy      2.0  2.0                                          2.0  2.0
 rab-fod.srwn.semi-mech.standard. light      2.0  2.0  2.0                                     2.0  2.0
 rab-fod.srws.semi-mech.standard. light      2.0  2.0  2.0                                     2.0  2.0

 wheat.  nwfp.bullock.  la-plant. heavy                     4.6  4.6                          38.4
 wheat.  pmw. bullock.  la-plant. heavy                     4.0  4.0                          27.1
 wheat.  pcw. bullock.  la-plant. heavy                     4.8  4.8                          39.2
 wheat.  psw. bullock.  la-plant. heavy                     4.7  4.7                          37.9
 wheat.  prw. bullock.  la-plant. heavy                     4.5  4.5                          17.2 12.9
 wheat.  scwn.bullock.  la-plant. heavy                     5.4                               31.2
 wheat.  srwn.bullock.  la-plant. heavy                     5.5                               33.4
 wheat.  scws.bullock.  la-plant. heavy                2.9  3.0                               32.3
 wheat.  srws.bullock.  la-plant. heavy                3.1  3.1                               31.6
 wheat.  nwfp.bullock.  la-plant. january                   4.6  4.6                          38.4
 wheat.  pmw. bullock.  la-plant. january                   4.0  4.0                          27.1
 wheat.  pcw. bullock.  la-plant. january                   4.8  4.8                          39.2
 wheat.  psw. bullock.  la-plant. january                   4.7  4.7                          37.9
 wheat.  prw. bullock.  la-plant. january                   4.5  4.5                          17.2 12.9
 wheat.  scwn.bullock.  la-plant. january                   5.4                               31.2
 wheat.  srwn.bullock.  la-plant. january                   5.5                               33.4
 wheat.  scws.bullock.  la-plant. january              2.9  3.0                               32.3
 wheat.  srws.bullock.  la-plant. january              3.1  3.1                               31.6
 wheat.  nwfp.bullock.  la-plant. light                     6.0  6.0                          38.4
 wheat.  pmw. bullock.  la-plant. light                     5.2  5.2                          27.1
 wheat.  pcw. bullock.  la-plant. light                     6.3  6.3                          39.2
 wheat.  psw. bullock.  la-plant. light                     6.1  6.1                          37.9
 wheat.  prw. bullock.  la-plant. light                     5.9  5.9                          17.2 12.9
 wheat.  scwn.bullock.  la-plant. light                     7.2                               31.2
 wheat.  srwn.bullock.  la-plant. light                     7.2                               33.4
 wheat.  scws.bullock.  la-plant. light                3.8  3.9                               32.3
 wheat.  srws.bullock.  la-plant. light                4.0  4.1                               31.6
 wheat.  nwfp.bullock.  la-plant. standard                  7.1  7.0                          38.4
 wheat.  pmw. bullock.  la-plant. standard                  6.2  6.1                          27.1
 wheat.  pcw. bullock.  la-plant. standard                  7.4  7.4                          39.2
 wheat.  psw. bullock.  la-plant. standard                  7.2  7.2                          37.9
 wheat.  prw. bullock.  la-plant. standard                  6.9  6.9                          17.2 12.9
 wheat.  scwn.bullock.  la-plant. standard                  8.4                               31.2
 wheat.  srwn.bullock.  la-plant. standard                  8.5                               33.4
 wheat.  scws.bullock.  la-plant. standard             4.5  4.6                               32.3
 wheat.  srws.bullock.  la-plant. standard             4.7  4.8                               31.6
 wheat.  nwfp.bullock.  qk-harv.  heavy                    10.4                          22.3 16.0
 wheat.  pmw. bullock.  qk-harv.  heavy                     9.1                          15.9 11.2
 wheat.  pcw. bullock.  qk-harv.  heavy                    10.9                          18.6 20.6
 wheat.  psw. bullock.  qk-harv.  heavy                    10.7                          18.3 19.6
 wheat.  prw. bullock.  qk-harv.  heavy                    10.3                               30.1
 wheat.  scwn.bullock.  qk-harv.  heavy                3.1  3.1                          17.5 13.7
 wheat.  srwn.bullock.  qk-harv.  heavy                3.2  3.2                          18.3 15.1
 wheat.  scws.bullock.  qk-harv.  heavy                6.7                               13.4 18.9
 wheat.  srws.bullock.  qk-harv.  heavy                7.0                               13.7 17.9
 wheat.  nwfp.bullock.  qk-harv.  january                  10.4                          22.3 16.0
 wheat.  pmw. bullock.  qk-harv.  january                   9.1                          15.9 11.2
 wheat.  pcw. bullock.  qk-harv.  january                  10.9                          18.6 20.6
 wheat.  psw. bullock.  qk-harv.  january                  10.7                          18.3 19.6
 wheat.  prw. bullock.  qk-harv.  january                  10.3                               30.1
 wheat.  scwn.bullock.  qk-harv.  january              3.1  3.1                          17.5 13.7
 wheat.  srwn.bullock.  qk-harv.  january              3.2  3.2                          18.3 15.1
 wheat.  scws.bullock.  qk-harv.  january              6.7                               13.4 18.9
 wheat.  srws.bullock.  qk-harv.  january              7.0                               13.7 17.9
 wheat.  nwfp.bullock.  qk-harv.  light                    13.6                          22.3 16.0
 wheat.  pmw. bullock.  qk-harv.  light                    11.9                          15.9 11.2
 wheat.  pcw. bullock.  qk-harv.  light                    14.3                          18.6 20.6
 wheat.  psw. bullock.  qk-harv.  light                    13.9                          18.3 19.6
 wheat.  prw. bullock.  qk-harv.  light                    13.4                               30.1
 wheat.  scwn.bullock.  qk-harv.  light                4.1  4.0                          17.5 13.7
 wheat.  srwn.bullock.  qk-harv.  light                4.1  4.1                          18.3 15.1
 wheat.  scws.bullock.  qk-harv.  light                8.6                               13.4 18.9
 wheat.  srws.bullock.  qk-harv.  light                9.1                               13.7 17.9
 wheat.  nwfp.bullock.  qk-harv.  standard                 16.0                          22.3 16.0
 wheat.  pmw. bullock.  qk-harv.  standard                 14.0                          15.9 11.2
 wheat.  pcw. bullock.  qk-harv.  standard                 16.8                          18.6 20.6
 wheat.  psw. bullock.  qk-harv.  standard                 16.4                          18.3 19.6
 wheat.  prw. bullock.  qk-harv.  standard                 15.8                               30.1
 wheat.  scwn.bullock.  qk-harv.  standard             4.7  4.8                          17.5 13.7
 wheat.  srwn.bullock.  qk-harv.  standard             4.7  4.8                          18.3 15.1
 wheat.  scws.bullock.  qk-harv.  standard            10.3                               13.4 18.9
 wheat.  srws.bullock.  qk-harv.  standard            10.8                               13.7 17.9
 wheat.  nwfp.bullock.  standard. heavy                     5.2  5.2                     22.3 16.0
 wheat.  pmw. bullock.  standard. heavy                     4.5  4.6                     15.9 11.2
 wheat.  pcw. bullock.  standard. heavy                     5.5  5.5                     18.6 20.6
 wheat.  psw. bullock.  standard. heavy                     5.3  5.4                     18.3 19.6
 wheat.  prw. bullock.  standard. heavy                     5.1  5.2                          30.1
 wheat.  scwn.bullock.  standard. heavy                     6.2                          17.5 13.7
 wheat.  srwn.bullock.  standard. heavy                     6.3                          18.3 15.1
 wheat.  scws.bullock.  standard. heavy                3.3  3.4                          13.4 18.9
 wheat.  srws.bullock.  standard. heavy                3.5  3.5                          13.7 17.9
 wheat.  nwfp.bullock.  standard. january                   5.2  5.2                     22.3 16.0
 wheat.  pmw. bullock.  standard. january                   4.5  4.6                     15.9 11.2
 wheat.  pcw. bullock.  standard. january                   5.5  5.5                     18.6 20.6
 wheat.  psw. bullock.  standard. january                   5.3  5.4                     18.3 19.6
 wheat.  prw. bullock.  standard. january                   5.1  5.2                          30.1
 wheat.  scwn.bullock.  standard. january                   6.2                          17.5 13.7
 wheat.  srwn.bullock.  standard. january                   6.3                          18.3 15.1
 wheat.  scws.bullock.  standard. january              3.3  3.4                          13.4 18.9
 wheat.  srws.bullock.  standard. january              3.5  3.5                          13.7 17.9
 wheat.  nwfp.bullock.  standard. light                     6.8  6.8                     22.3 16.0
 wheat.  pmw. bullock.  standard. light                     6.0  6.0                     15.9 11.2
 wheat.  pcw. bullock.  standard. light                     7.1  7.1                     18.6 20.6
 wheat.  psw. bullock.  standard. light                     6.9  7.0                     18.3 19.6
 wheat.  prw. bullock.  standard. light                     6.7  6.7                          30.1
 wheat.  scwn.bullock.  standard. light                     8.1                          17.5 13.7
 wheat.  srwn.bullock.  standard. light                     8.2                          18.3 15.1
 wheat.  scws.bullock.  standard. light                4.3  4.3                          13.4 18.9
 wheat.  srws.bullock.  standard. light                4.5  4.6                          13.7 17.9
 wheat.  nwfp.bullock.  standard. standard                  8.0  8.0                     22.3 16.0
 wheat.  pmw. bullock.  standard. standard                  7.0  7.0                     15.9 11.2
 wheat.  pcw. bullock.  standard. standard                  8.4  8.4                     18.6 20.6
 wheat.  psw. bullock.  standard. standard                  8.2  8.2                     18.3 19.6
 wheat.  prw. bullock.  standard. standard                  7.9  7.9                          30.1
 wheat.  scwn.bullock.  standard. standard                  9.5                          17.5 13.7
 wheat.  srwn.bullock.  standard. standard                  9.7                          18.3 15.1
 wheat.  scws.bullock.  standard. standard             5.1  5.2                          13.4 18.9
 wheat.  srws.bullock.  standard. standard             5.4  5.4                          13.7 17.9
 wheat.  nwfp.semi-mech.la-plant. heavy                     3.2  3.2
 wheat.  pmw. semi-mech.la-plant. heavy                     3.6  3.6
 wheat.  pcw. semi-mech.la-plant. heavy                     3.4  3.4
 wheat.  psw. semi-mech.la-plant. heavy                     3.2  3.2
 wheat.  prw. semi-mech.la-plant. heavy                     6.3  6.3
 wheat.  scwn.semi-mech.la-plant. heavy                     5.4
 wheat.  srwn.semi-mech.la-plant. heavy                     5.5
 wheat.  scws.semi-mech.la-plant. heavy                2.9  3.0
 wheat.  srws.semi-mech.la-plant. heavy                3.1  3.1
 wheat.  nwfp.semi-mech.la-plant. january                   3.2  3.2
 wheat.  pmw. semi-mech.la-plant. january                   3.6  3.6
 wheat.  pcw. semi-mech.la-plant. january                   3.4  3.4
 wheat.  psw. semi-mech.la-plant. january                   3.2  3.2
 wheat.  prw. semi-mech.la-plant. january                   6.3  6.3
 wheat.  scwn.semi-mech.la-plant. january                   5.4
 wheat.  srwn.semi-mech.la-plant. january                   5.5
 wheat.  scws.semi-mech.la-plant. january              2.9  3.0
 wheat.  srws.semi-mech.la-plant. january              3.1  3.1
 wheat.  nwfp.semi-mech.la-plant. light                     4.1  4.2
 wheat.  pmw. semi-mech.la-plant. light                     4.7  4.7
 wheat.  pcw. semi-mech.la-plant. light                     4.4  4.4
 wheat.  psw. semi-mech.la-plant. light                     4.2  4.2
 wheat.  prw. semi-mech.la-plant. light                     8.2  8.3
 wheat.  scwn.semi-mech.la-plant. light                     7.2
 wheat.  srwn.semi-mech.la-plant. light                     7.2
 wheat.  scws.semi-mech.la-plant. light                3.8  3.9
 wheat.  srws.semi-mech.la-plant. light                4.0  4.1
 wheat.  nwfp.semi-mech.la-plant. standard                  4.9  4.9
 wheat.  pmw. semi-mech.la-plant. standard                  5.5  5.5
 wheat.  pcw. semi-mech.la-plant. standard                  5.2  5.2
 wheat.  psw. semi-mech.la-plant. standard                  4.9  4.9
 wheat.  prw. semi-mech.la-plant. standard                  9.7  9.7
 wheat.  scwn.semi-mech.la-plant. standard                  8.4
 wheat.  srwn.semi-mech.la-plant. standard                  8.5
 wheat.  scws.semi-mech.la-plant. standard             4.5  4.6
 wheat.  srws.semi-mech.la-plant. standard             4.7  4.8
 wheat.  nwfp.semi-mech.qk-harv.  heavy                     7.3
 wheat.  pmw. semi-mech.qk-harv.  heavy                     8.2
 wheat.  pcw. semi-mech.qk-harv.  heavy                     7.7
 wheat.  psw. semi-mech.qk-harv.  heavy                     7.3
 wheat.  prw. semi-mech.qk-harv.  heavy                     7.2
 wheat.  scwn.semi-mech.qk-harv.  heavy                3.1  3.1
 wheat.  srwn.semi-mech.qk-harv.  heavy                3.2  3.2
 wheat.  scws.semi-mech.qk-harv.  heavy                6.7
 wheat.  srws.semi-mech.qk-harv.  heavy                7.0
 wheat.  nwfp.semi-mech.qk-harv.  january                   7.3
 wheat.  pmw. semi-mech.qk-harv.  january                   8.2
 wheat.  pcw. semi-mech.qk-harv.  january                   7.7
 wheat.  psw. semi-mech.qk-harv.  january                   7.3
 wheat.  prw. semi-mech.qk-harv.  january                   7.2
 wheat.  scwn.semi-mech.qk-harv.  january              3.1  3.1
 wheat.  srwn.semi-mech.qk-harv.  january              3.2  3.2
 wheat.  scws.semi-mech.qk-harv.  january              6.7
 wheat.  srws.semi-mech.qk-harv.  january              7.0
 wheat.  nwfp.semi-mech.qk-harv.  light                     9.5
 wheat.  pmw. semi-mech.qk-harv.  light                    10.7
 wheat.  pcw. semi-mech.qk-harv.  light                    10.0
 wheat.  psw. semi-mech.qk-harv.  light                     9.5
 wheat.  prw. semi-mech.qk-harv.  light                     9.4
 wheat.  scwn.semi-mech.qk-harv.  light                4.1  4.0
 wheat.  srwn.semi-mech.qk-harv.  light                4.1  4.1
 wheat.  scws.semi-mech.qk-harv.  light                8.6
 wheat.  srws.semi-mech.qk-harv.  light                9.1
 wheat.  nwfp.semi-mech.qk-harv.  standard                 11.2
 wheat.  pmw. semi-mech.qk-harv.  standard                 12.6
 wheat.  pcw. semi-mech.qk-harv.  standard                 11.8
 wheat.  psw. semi-mech.qk-harv.  standard                 11.2
 wheat.  prw. semi-mech.qk-harv.  standard                 11.1
 wheat.  scwn.semi-mech.qk-harv.  standard             4.7  4.8
 wheat.  srwn.semi-mech.qk-harv.  standard             4.7  4.8
 wheat.  scws.semi-mech.qk-harv.  standard            10.3
 wheat.  srws.semi-mech.qk-harv.  standard            10.8
 wheat.  nwfp.semi-mech.standard. heavy                     3.6  3.7
 wheat.  pmw. semi-mech.standard. heavy                     4.1  4.1
 wheat.  pcw. semi-mech.standard. heavy                     3.8  3.8
 wheat.  psw. semi-mech.standard. heavy                     3.6  3.7
 wheat.  prw. semi-mech.standard. heavy                     3.6  3.6
 wheat.  scwn.semi-mech.standard. heavy                     6.2
 wheat.  srwn.semi-mech.standard. heavy                     6.3
 wheat.  scws.semi-mech.standard. heavy                3.3  3.4
 wheat.  srws.semi-mech.standard. heavy                3.5  3.5
 wheat.  nwfp.semi-mech.standard. january                   3.6  3.7
 wheat.  pmw. semi-mech.standard. january                   4.1  4.1
 wheat.  pcw. semi-mech.standard. january                   3.8  3.8
 wheat.  psw. semi-mech.standard. january                   3.6  3.7
 wheat.  prw. semi-mech.standard. january                   3.6  3.6
 wheat.  scwn.semi-mech.standard. january                   6.2
 wheat.  srwn.semi-mech.standard. january                   6.3
 wheat.  scws.semi-mech.standard. january              3.3  3.4
 wheat.  srws.semi-mech.standard. january              3.5  3.5
 wheat.  nwfp.semi-mech.standard. light                     4.7  4.8
 wheat.  pmw. semi-mech.standard. light                     5.4  5.3
 wheat.  pcw. semi-mech.standard. light                     5.0  5.0
 wheat.  psw. semi-mech.standard. light                     4.7  4.8
 wheat.  prw. semi-mech.standard. light                     4.7  4.7
 wheat.  scwn.semi-mech.standard. light                     8.1
 wheat.  srwn.semi-mech.standard. light                     8.2
 wheat.  scws.semi-mech.standard. light                4.3  4.3
 wheat.  srws.semi-mech.standard. light                4.5  4.6
 wheat.  nwfp.semi-mech.standard. standard                  5.6  5.6
 wheat.  pmw. semi-mech.standard. standard                  6.3  6.3
 wheat.  pcw. semi-mech.standard. standard                  5.9  5.9
 wheat.  psw. semi-mech.standard. standard                  5.6  5.6
 wheat.  prw. semi-mech.standard. standard                  5.5  5.6
 wheat.  scwn.semi-mech.standard. standard                  9.5
 wheat.  srwn.semi-mech.standard. standard                  9.7
 wheat.  scws.semi-mech.standard. standard             5.1  5.2
 wheat.  srws.semi-mech.standard. standard             5.4  5.4
 orchard.(nwfp,pmw,pcw,psw,prw,scwn,scws,
     srwn,srws).bullock.standard.standard    5.   5.   3.5  3.5  3.5  3.5  3.5  3.5  3.5  3.5  3.5  3.5
 orchard.(nwfp,pmw,pcw,psw,prw,scwn,scws,
    srwn,srws).semi-mech.standard.standard   3.   3.   3.   2.   1.                  2.   2.   2.   2.

 chilli.(nwfp,pmw,psw,pcw,prw).semi-mech.
                       standard.standard     16.
 chilli.(scwn,scws,srwn,srws).semi-mech.
                       standard.standard                                                             8.
           ;
*-  convert bullock pair hours to bullock hours
*&&z   zone1xxxxx(z,c,t,s,w,m) = bullock(c,z,t,s,w,m) ;
  bullock(c,z,t,s,w,m) = bullock(c,z,t,s,w,m)*2;
$eject
  table labor(c,z,t,s,w,m) labor requirements for crops(man hours)

                                             jan   feb   mar   apr   may   jun   jul   aug   sep   oct   nov   dec
 basmati.pmw. bullock.  standard. standard                                29.1  88.78 65.9   5.6   5.6  47.9  17.6
 basmati.pcw. bullock.  standard. standard                                29.1  88.78 65.9   5.6   5.6  47.9  17.6
 basmati.psw. bullock.  standard. standard                                26.9  94.9  74.7   6.2   3.1  50.9  19.4
 basmati.prw. bullock.  standard. standard                                44.2  91.8  63.9   8.4   5.6  41.9  25.0
 basmati.pmw. semi-mech.standard. standard                                 8.4  71.6  65.9   5.6   5.6  47.4  16.1
 basmati.pcw. semi-mech.standard. standard                                 8.4  71.6  65.9   5.6   5.6  47.4  16.1
 basmati.psw. semi-mech.standard. standard                                 5.6  74.4  74.7   6.2   3.1  50.4  16.9
 basmati.prw. semi-mech.standard. standard                                11.0  67.7  63.9   8.4   5.6  41.4  22.5

 rab-fod.srwn.bullock.  standard. heavy     27.9  23.4                                       9.1  19.7  19.2  24.8
 rab-fod.srws.bullock.  standard. heavy     27.4  24.8                                      10.7  19.5  16.8  27.9
 rab-fod.srwn.bullock.  standard. light     27.9  23.4  22.5                                 9.1  19.7  19.2  24.8
 rab-fod.srws.bullock.  standard. light     27.4  24.8  19.7                                10.7  19.5  16.8  27.9
 rab-fod.nwfp.bullock.  standard. standard  42.8  42.1  32.6  24.6  10.9                     8.0  14.3  17.8  35.6
 rab-fod.pmw. bullock.  standard. standard  29.3  37.1  36.6  28.1   9.7                          13.35 16.5  34.5
 rab-fod.pcw. bullock.  standard. standard  32.3  41.5  41.8  29.4  18.2                    10.1  15.5  23.5  29.1
 rab-fod.psw. bullock.  standard. standard  31.9  39.8  41.7  28.6  13.4                           8.7  15.7  28.7
 rab-fod.prw. bullock.  standard. standard  31.4  32.0  34.5  23.4  21.9                           4.5  14.4  27.1
 rab-fod.scwn.bullock.  standard. standard  28.0  24.8  24.0  14.3                           9.6  21.5  19.2  24.5
 rab-fod.srwn.bullock.  standard. standard  27.9  23.4  22.5  15.0                           9.1  19.7  19.2  24.8
 rab-fod.scws.bullock.  standard. standard  29.0  23.1  19.8  17.5   8.6                    10.3  19.8  17.5  27.8
 rab-fod.srws.bullock.  standard. standard  27.4  24.8  19.7  18.8   9.3                    10.7  19.5  16.8  27.9
 rab-fod.srwn.semi-mech.standard. heavy     27.9  23.4                                       1.0   8.7  19.2  24.8
 rab-fod.srws.semi-mech.standard. heavy     27.4  24.8                                       1.2   9.3  16.8  27.9
 rab-fod.srwn.semi-mech.standard. light     27.9  23.4  22.5                                 1.0   8.7  19.2  24.8
 rab-fod.srws.semi-mech.standard. light     27.4  24.8  19.7                                 1.2   9.3  16.8  27.9
 rab-fod.nwfp.semi-mech.standard. standard  42.2  43.5  31.2  22.0   9.5                     2.05  8.02 16.29 33.0
 rab-fod.pmw. semi-mech.standard. standard  28.6  36.4  35.6  27.5   8.3                           2.17  7.96 24.08
 rab-fod.pcw. semi-mech.standard. standard  31.1  39.8  40.3  28.7  17.6                     7.5  13.2  21.6  28.4
 rab-fod.psw. semi-mech.standard. standard  30.5  38.4  40.3  27.1  11.9                           7.3  12.7  26.4
 rab-fod.prw. semi-mech.standard. standard  29.6  30.4  29.6  22.5  20.7                           3.4  12.8  25.7
 rab-fod.scwn.semi-mech.standard. standard  26.5  21.4  22.8  12.7                           1.2   8.8  17.5  22.8
 rab-fod.srwn.semi-mech.standard. standard  27.9  23.4  22.5  15.0                           1.0   8.7  19.2  24.8
 rab-fod.scws.semi-mech.standard. standard  29.0  23.1  19.8  17.5   8.6                     1.2   8.9  17.5  27.8
 rab-fod.srws.semi-mech.standard. standard  27.4  24.8  19.7  18.8   9.3                     1.2   9.3  16.8  27.9

 cotton. pcw. bullock.  el-plant. standard                    26.8   6.7  20.5   2.5   7.6  13.2  41.3  57.4  22.1
 cotton. pmw. bullock.  standard. standard                          16.85  9.6         2.9  12.0  50.74 70.5  26.5
 cotton. pcw. bullock.  standard. standard                          35.59 18.4   2.5   7.6  13.2  41.3  57.4  22.1
 cotton. psw. bullock.  standard. standard                          10.2  19.32  2.5   3.9  11.0  40.1  52.9  16.3
 cotton. prw. bullock.  standard. standard                          22.3  20.28        6.7  10.6  40.1  57.1  21.0
 cotton. scwn.bullock.  standard. standard                     7.4  23.5   7.0   2.0   2.0  17.5  39.4  66.2
 cotton. scws.bullock.  standard. standard                    31.2   2.0   7.0   2.0   2.0  38.6  80.4
 cotton. srws.bullock.  standard. standard                     9.5  23.8   7.0   2.0   2.0  28.5  84.1
 cotton. pcw. semi-mech.la-plant. standard                                14.4   2.5   7.6  12.7  40.8  56.9  21.6
 cotton. pmw. semi-mech.standard. standard                          13.36  5.3         2.9  11.5  50.24 70.0  26.0
 cotton. pcw. semi-mech.standard. standard                           4.99  9.42  2.5   7.6  12.7  40.8  56.9  21.6
 cotton. psw. semi-mech.standard. standard                           4.41  9.28  2.5   3.9  10.5  39.6  52.4  15.8
 cotton. prw. semi-mech.standard. standard                           1.45  9.04        6.7  10.1  39.6  56.6  20.5
 cotton. scwn.semi-mech.standard. standard                          13.2                    15.8  37.4  62.2
 cotton. scws.semi-mech.standard. standard                    12.1         7.0   2.0   2.0  36.6  76.4
 cotton. srws.semi-mech.standard. standard                          12.4   7.0   2.0   2.0  28.5  83.2

 gram.   nwfp.bullock.  standard. standard   3.4              23.2                           8.9  18.0
 gram.   pmw. bullock.  standard. standard              21.5                                10.7   9.1   1.9   2.5
 gram.   pcw. bullock.  standard. standard              21.7                                10.8   8.9   0.7   2.5
 gram.   psw. bullock.  standard. standard   2.5        22.3                                      10.2   8.6   0.8
 gram.   prw. bullock.  standard. standard   2.5        22.6                                      10.3   8.5   0.7
 gram.   scwn.bullock.  standard. standard              23.5   9.2                           6.4  12.0   1.0   1.0
 gram.   scws.bullock.  standard. standard              23.5   9.2                           6.4  12.0   1.0   1.0
 gram.   srwn.bullock.  standard. standard              24.8  10.3                           7.1  13.4   1.0   1.0
 gram.   srws.bullock.  standard. standard              24.8  10.3                           7.1  13.4   1.0   1.0
 gram.   nwfp.semi-mech.standard. standard   3.4              20.5                           1.87  6.0
 gram.   pmw. semi-mech.standard. standard              19.0                                 1.63  4.64  0.9   2.5
 gram.   pcw. semi-mech.standard. standard              20.0                                 1.62  4.41  0.7   2.5
 gram.   psw. semi-mech.standard. standard   2.5        21.4                                       1.54  4.39  0.8
 gram.   prw. semi-mech.standard. standard   2.5        21.9                                       1.45  4.19  0.7
 gram.   scwn.semi-mech.standard. standard              21.4   5.5                           0.81  3.4   1.0   1.0
 gram.   scws.semi-mech.standard. standard              21.4   5.5                           0.81  3.4   1.0   1.0
 gram.   srwn.semi-mech.standard. standard              22.7   9.1                           0.96  3.6   1.0   1.0
 gram.   srws.semi-mech.standard. standard              22.7   9.1                           0.96  3.6   1.0   1.0

 irri.   pmw. bullock.  standard. standard                                22.9  122.3 35.9   8.4   6.4  43.9  20.4
 irri.   pcw. bullock.  standard. standard                                22.9  122.3 35.9   8.4   6.4  43.9  20.4
 irri.   psw. bullock.  standard. standard                                22.3  135.7 37.8   6.2   6.2  51.7  22.2
 irri.   prw. bullock.  standard. standard                                122.9 86.7  15.4   8.2   6.8  41.0  24.0
 irri.   scwn.bullock.  standard. standard                          13.6  81.0  13.6  32.0  13.6  39.1  20.4
 irri.   srwn.bullock.  standard. standard                           3.0  12.7  91.4  31.5  12.8  13.1  39.0  14.0
 irri.   scws.bullock.  standard. standard                     9.1  10.3  57.9  13.4  35.3  13.3  41.4  16.7
 irri.   srws.bullock.  standard. standard                          19.0  58.5  13.2  34.5  13.1  40.3  13.5  13.5
 irri.   pmw. semi-mech.standard. standard                                 5.3  105.5 29.9   8.4   6.4  43.4  17.9
 irri.   pcw. semi-mech.standard. standard                                 5.3  105.5 29.9   8.4   6.4  43.4  17.9
 irri.   psw. semi-mech.standard. standard                                 7.9  119.2 28.8   6.2   6.2  51.2  19.7
 irri.   prw. semi-mech.standard. standard                                71.0  80.4  13.7   8.2   6.8  41.5  21.5
 irri.   scwn.semi-mech.standard. standard                           3.3   4.7  13.6  32.0  13.6  36.6  17.9
 irri.   srwn.semi-mech.standard. standard                           3.0   8.7  13.4  31.5  12.8  13.1  36.6  14.0
 irri.   scws.semi-mech.standard. standard                     2.3   2.9  43.5  13.4  35.3  13.3  38.8  14.2
 irri.   srws.semi-mech.standard. standard                           5.0  43.7  13.2  34.5  13.1  37.6  13.5  13.5

 maize.  nwfp.bullock.  standard. standard                                43.4  38.6   0.7   2.7  46.1
 maize.  pcw. bullock.  standard. standard                                10.8   4.5  16.7  44.2   2.5  26.8  27.0
 maize.  psw. bullock.  standard. standard                                10.2   4.1  15.9  42.3   2.0  23.0  23.8
 maize.  prw. bullock.  standard. standard                                10.2   4.1  15.9  42.3   2.0  23.0  23.8
 maize.  scwn.bullock.  standard. standard                                18.65 27.4   3.2   3.2   3.2  57.4
 maize.  scws.bullock.  standard. standard                          17.9  26.8   3.0   3.0   3.0   3.0  59.3
 maize.  nwfp.semi-mech.standard. standard                                14.6  38.6   0.7   2.7  43.8
 maize.  pcw. semi-mech.standard. standard                                       2.41  8.36 44.2   2.5  26.0  26.8
 maize.  psw. semi-mech.standard. standard                                       2.24  8.06 42.3   2.5  22.0  22.5
 maize.  prw. semi-mech.standard. standard                                       2.24  8.06 42.3   2.5  22.0  22.5
 maize.  scwn.semi-mech.standard. standard                                 3.1   5.2   3.2   3.2   3.2  57.4
 maize.  scws.semi-mech.standard. standard                           3.4   5.8   3.0   3.0   3.0   3.0  59.3

 mus+rap.nwfp.bullock.  standard. standard                                       8.9  10.5  18.0  14.6  23.0  23.0
 mus+rap.pmw. bullock.  standard. standard  25.7  15.0  20.5                                      19.3   9.0  23.0
 mus+rap.pcw. bullock.  standard. standard  21.75 20.5   6.0                                      10.8  13.6  13.9
 mus+rap.psw. bullock.  standard. standard  20.5   4.0   3.0                                10.2  12.8  13.2  23.0
 mus+rap.prw. bullock.  standard. standard   3.5   3.0                                      19.1  15.7  25.5  23.0
 mus+rap.scwn.bullock.  standard. standard   1.5   1.5   1.5                                15.1  12.3   1.5   1.5
 mus+rap.srwn.bullock.  standard. standard   1.5   1.5   1.5                                15.2  12.7   1.5   1.5
 mus+rap.srws.bullock.  standard. standard   1.5   1.5   1.5                                15.2  12.7   1.5   1.5
 mus+rap.scws.bullock.  standard. standard   1.5   1.5   1.5                                13.5   7.6   1.5   1.5
 mus+rap.nwfp.semi-mech.standard. standard                                       1.87  3.45  6.0  14.6  23.0  23.0
 mus+rap.pmw. semi-mech.standard. standard  25.7  13.5  18.0                                       3.8   5.5  23.0
 mus+rap.pcw. semi-mech.standard. standard  21.75 20.25  6.0                                       1.62  5.08 13.9
 mus+rap.psw. semi-mech.standard. standard  20.5   3.5   2.5                                 1.54  4.9  13.2  23.0
 mus+rap.prw. semi-mech.standard. standard   4.0   2.0                                       5.94 15.7  25.5  23.0
 mus+rap.scwn.semi-mech.standard. standard   1.0   1.0   1.0                                 3.7   2.8   1.0   1.0
 mus+rap.srwn.semi-mech.standard. standard   1.0   1.0   1.0                                 3.1   2.5   1.0   1.0
 mus+rap.srws.semi-mech.standard. standard   1.0   1.0   1.0                                 3.1   2.5   1.0   1.0
 mus+rap.scws.semi-mech.standard. standard   1.0   1.0   1.0                                 3.3   2.9   1.0   1.0

 sc-gur. nwfp.bullock.  standard. standard  154.6   147 57.4  23.4   5.4   5.74  3.4   3.4  24.0  54.0    216 208.2
 sc-gur.(pmw,pcw,psw,prw).
              bullock.  standard. standard    159 80.4  54.4  22.5   4.75  5.05  3.0   3.0   3.0   1.5  150.3 148.5
 sc-gur. scwn.bullock.  standard. standard  207.6 149.4 53.2  20.8   3.2   8.9   3.2   3.2   3.2   3.2  183.4 195.6
 sc-gur. srwn.bullock.  standard. standard  213.7 145.2 54.4  21.7   3.2   9.0   3.2   3.2   3.2   3.2  187.3   198
 sc-gur. scws.bullock.  standard. standard  207.1 163.5 53.4  31.5   3.2   9.0   3.2   3.2   3.2   3.2  187.4 191.5
 sc-gur. srws.bullock.  standard. standard  208.2 161.4 49.5  37.2   3.2   9.0   3.2   3.2   3.2   3.2  183.2 197.5
 sc-gur. nwfp.semi-mech.standard. standard  143.1 144.1 50.15 23.4   5.4   5.74  3.4   3.4  17.5  54.0    216 201.9
 sc-gur.(pmw,pcw,psw,prw).
              semi-mech.standard. standard  151.2 80.62 48.92 22.5   4.75  5.05  3.0   3.0   3.0   1.5  148.3 142.3
 sc-gur. scwn.semi-mech.standard. standard  201.2 143.4 48.5  16.4   3.2   8.9   3.2   3.2   3.2   3.2  179.2 193.4
 sc-gur. srwn.semi-mech.standard. standard  209.4 139.8 51.4  18.6   3.2   9.0   3.2   3.2   3.2   3.2  184.3 193.7
 sc-gur. scws.semi-mech.standard. standard  201.6 156.5 47.4  28.2   3.2   9.0   3.2   3.2   3.2   3.2  181.5 186.4
 sc-gur. srws.semi-mech.standard. standard  203.5 153.4 44.5  32.6   3.2   9.0   3.2   3.2   3.2   3.2  179.8 191.6

 sc-mill.nwfp.bullock.  standard. standard  120   115   95.0  65.0   5.4   5.74  3.4   3.4  23.4  54.0  107   125
 sc-mill.pmw. bullock.  standard. standard  87.0  85.2  74.0  64.5  44.5   5.05  3.0   3.0   3.0   1.5  74.0  85.5
 sc-mill.pcw. bullock.  standard. standard  90.0  85.0  95.0  72.1  30.0   5.05  3.0   3.0   3.0   1.5  85.0  95.0
 sc-mill.psw. bullock.  standard. standard  88.0  85.0  95.0  70.0  40.0   5.05  3.0   3.0   3.0   1.5  85.0  85.0
 sc-mill.prw. bullock.  standard. standard  90.0  95.0  90.9  65.0  38.5   5.05  3.0   3.0   3.0   1.5  84.0  85.5
 sc-mill.scwn.bullock.  standard. standard  87.5  85.0  92.1  47.5   3.2   8.4  35.4   3.2   3.2   3.2  93.1  87.3
 sc-mill.srwn.bullock.  standard. standard  89.3  87.1  94.2  43.1   3.2   9.2  37.8   3.2   3.2   3.2  94.2  89.1
 sc-mill.scws.bullock.  standard. standard  98.4  95.3  102.4 46.2   3.2   9.8  39.4   3.2   3.2   3.2  92.8  93.7
 sc-mill.srws.bullock.  standard. standard  94.7  96.5  103.4 41.4   3.2   8.9  41.3   3.2   3.2   3.2  91.7  89.8
 sc-mill.nwfp.semi-mech.standard. standard  116   116   90.6  58.0   5.4   5.74  3.4   3.4  22.0  51.0  102   120
 sc-mill.pmw. semi-mech.standard. standard  84.0  81.5  20.0  59.0  59.0  38.5   5.05  3.0   3.0   1.5  69.5  80.5
 sc-mill.pcw. semi-mech.standard. standard  86.5  81.5  90.0  64.0  34.5   5.05  3.0   3.0   3.0   1.5  80.5  90.5
 sc-mill.psw. semi-mech.standard. standard  86.0  82.0  90.5  64.5  35.0   5.05  3.0   3.0   3.0   1.5  81.5  84.0
 sc-mill.prw. semi-mech.standard. standard  87.0  92.5  87.0  62.5  35.0   5.05  3.0   3.0   3.0   1.5  80.0  81.5
 sc-mill.scwn.semi-mech.standard. standard  83.5  81.9  89.2  43.4   3.2   8.4  35.4   3.2   3.2   3.2  89.1  83.8
 sc-mill.srwn.semi-mech.standard. standard  84.1  82.3  90.1  44.5   3.2   9.2  35.7   3.2   3.2   3.2  91.2  84.5
 sc-mill.scws.semi-mech.standard. standard  96.5  92.7  99.4  46.2   3.2   9.8  39.4   3.2   3.2   3.2  89.8  91.1
 sc-mill.srws.semi-mech.standard. standard  92.8  93.7  101.8 39.1   3.2   8.9  40.5   3.2   3.2   3.2  89.5  84.7

 kha-fod.srws.bullock.  la-plant. standard                          13.4   3.0   3.0   8.3   3.0   3.0
 kha-fod.nwfp.bullock.  standard. standard              18.5   2.5   3.0  17.5   2.0  18.0   3.0   1.5
 kha-fod.pmw. bullock.  standard. standard              16.5   3.0   3.5  17.0   4.5  18.5   2.5   0.5
 kha-fod.pcw. bullock.  standard. standard              18.5   2.0   3.0  18.7   4.0  18.0   3.0   1.0
 kha-fod.psw. bullock.  standard. standard              18.5   1.8   2.5  25.5   4.5  27.5   2.5   1.2
 kha-fod.prw. bullock.  standard. standard              17.0   1.5   2.5  20.0   4.0  18.5   3.0   1.5
 kha-fod.scwn.bullock.  standard. standard              10.7   3.0   7.3   3.0   3.0   7.3   3.0   3.0
 kha-fod.scws.bullock.  standard. standard              13.4   3.0   8.2   3.0   3.0   8.3   3.0   3.0
 kha-fod.srwn.bullock.  standard. standard              14.3   3.0  11.4   3.0   3.0   7.5   3.0   3.0
 kha-fod.srws.bullock.  standard. standard              14.3   3.0  11.4   3.0   3.0   7.5   3.0   3.0
 kha-fod.srws.semi-mech.la-plant. standard                           4.0   2.1   0.8   2.3   1.2   1.2
 kha-fod.nwfp.semi-mech.standard. standard               7.0   1.0   2.5   6.0   3.5   7.0   3.0   1.5
 kha-fod.pmw. semi-mech.standard. standard               7.5   1.5   1.0   7.0   1.0   8.0   2.0   1.0
 kha-fod.pcw. semi-mech.standard. standard               6.0   1.3   4.5   6.5   4.2   5.5   2.0   1.5
 kha-fod.psw. semi-mech.standard. standard               5.5   1.2   1.7   7.0   4.5   4.0   2.0   1.5
 kha-fod.prw. semi-mech.standard. standard               6.0   1.5   2.0   8.5   3.8   5.5   2.2   1.5
 kha-fod.scwn.semi-mech.standard. standard               2.1   1.2   2.7   1.2   1.2   2.7   1.2   1.2
 kha-fod.scws.semi-mech.standard. standard               2.0   1.2   2.8   1.2   1.2   2.9   1.2   1.2
 kha-fod.srwn.semi-mech.standard. standard               2.2   1.2   2.1   1.2   1.2   2.3   1.2   1.2
 kha-fod.srws.semi-mech.standard. standard               2.2   1.2   2.1   1.2   1.2   2.3   1.2   1.2

 wheat.  nwfp.bullock.  la-plant. heavy      5.8   3.7   3.7  38.3  14.2                                43.6   3.7
 wheat.  pmw. bullock.  la-plant. heavy      3.4   2.9   1.7  36.3  12.1                                31.9   2.9
 wheat.  pcw. bullock.  la-plant. heavy      4.3   3.9   3.9  33.2  12.3                                45.2   3.9
 wheat.  psw. bullock.  la-plant. heavy      4.4   2.7   2.4  35.9  15.3                                41.4   4.0
 wheat.  prw. bullock.  la-plant. heavy      4.1   2.4   2.4  36.1  15.4                                17.2  19.1
 wheat.  scwn.bullock.  la-plant. heavy      3.2   3.2        24.0                                      34.0   3.2
 wheat.  srwn.bullock.  la-plant. heavy      3.2   3.2        25.4                                      37.5   3.2
 wheat.  scws.bullock.  la-plant. heavy      3.2   3.2  14.0   7.3                                      36.3   3.2
 wheat.  srws.bullock.  la-plant. heavy      3.2   3.2  10.7   7.1                                      36.2   3.2
 wheat.  nwfp.bullock.  la-plant. january    5.8   3.7   3.7  38.3  14.2                                43.6   3.7
 wheat.  pmw. bullock.  la-plant. january    3.4   2.9   1.7  36.3  12.1                                31.9   2.9
 wheat.  pcw. bullock.  la-plant. january    4.3   3.9   3.9  33.2  12.3                                45.2   3.9
 wheat.  psw. bullock.  la-plant. january    4.4   2.7   2.4  35.9  15.3                                41.4   4.0
 wheat.  prw. bullock.  la-plant. january    4.1   2.4   2.4  36.1  15.4                                17.2  19.1
 wheat.  scwn.bullock.  la-plant. january    3.2   3.2        24.0                                      34.0   3.2
 wheat.  srwn.bullock.  la-plant. january    3.2   3.2        25.4                                      37.5   3.2
 wheat.  scws.bullock.  la-plant. january    3.2   3.2  14.0   7.3                                      36.3   3.2
 wheat.  srws.bullock.  la-plant. january    3.2   3.2  10.7   7.1                                      36.2   3.2
 wheat.  nwfp.bullock.  la-plant. light      5.8   3.7   3.7  50.1  18.5                                43.6   3.7
 wheat.  pmw. bullock.  la-plant. light      3.4   2.9   1.7  47.5  15.8                                31.9   2.9
 wheat.  pcw. bullock.  la-plant. light      4.3   3.9   3.9  43.4  16.1                                45.2   3.9
 wheat.  psw. bullock.  la-plant. light      4.4   2.7   2.4  47.0  20.1                                41.4   4.0
 wheat.  prw. bullock.  la-plant. light      4.1   2.4   2.4  47.1  20.1                                17.2  19.1
 wheat.  scwn.bullock.  la-plant. light      3.2   3.2        31.4                                      34.0   3.2
 wheat.  srwn.bullock.  la-plant. light      3.2   3.2        33.2                                      37.5   3.2
 wheat.  scws.bullock.  la-plant. light      3.2   3.2  18.3   9.6                                      36.3   3.2
 wheat.  srws.bullock.  la-plant. light      3.2   3.2  13.9   9.3                                      36.2   3.2
 wheat.  nwfp.bullock.  la-plant. standard   5.8   3.7   3.7  59.0  21.8                                43.6   3.7
 wheat.  pmw. bullock.  la-plant. standard   3.4   2.9   1.7  55.9  18.6                                31.9   2.9
 wheat.  pcw. bullock.  la-plant. standard   4.3   3.9   3.9  51.4  18.9                                45.2   3.9
 wheat.  psw. bullock.  la-plant. standard   4.4   2.7   2.4  55.3  23.6                                41.4   4.0
 wheat.  prw. bullock.  la-plant. standard   4.1   2.4   2.4  55.5  23.7                                17.2  19.1
 wheat.  scwn.bullock.  la-plant. standard   3.2   3.2        37.0                                      34.0   3.2
 wheat.  srwn.bullock.  la-plant. standard   3.2   3.2        39.1                                      37.5   3.2
 wheat.  scws.bullock.  la-plant. standard   3.2   3.2  21.5  11.3                                      36.3   3.2
 wheat.  srws.bullock.  la-plant. standard   3.2   3.2  16.4  10.9                                      36.2   3.2
 wheat.  nwfp.bullock.  qk-harv.  heavy      5.8   3.7   3.7  59.8                                22.3  21.2   3.7
 wheat.  pmw. bullock.  qk-harv.  heavy      3.4   2.9   1.7  55.1                                15.9  16.0   2.9
 wheat.  pcw. bullock.  qk-harv.  heavy      4.3   3.9   3.9  57.3                                18.6  23.9   3.9
 wheat.  psw. bullock.  qk-harv.  heavy      4.4   2.7   2.4  58.4                                18.3  23.1   4.0
 wheat.  prw. bullock.  qk-harv.  heavy      4.1   2.4   2.4  58.6                                      33.9   2.4
 wheat.  scwn.bullock.  qk-harv.  heavy      3.2   3.2  19.2   8.3                                18.9  15.1   3.2
 wheat.  srwn.bullock.  qk-harv.  heavy      3.2   3.2  19.3   9.6                                20.2  17.3   3.2
 wheat.  scws.bullock.  qk-harv.  heavy      3.2   3.2  24.2                                      16.2  20.1   3.2
 wheat.  srws.bullock.  qk-harv.  heavy      3.2   3.2  22.3                                      16.5  19.7   3.2
 wheat.  nwfp.bullock.  qk-harv.  january    5.8   3.7   3.7  59.8                                22.3  21.2   3.7
 wheat.  pmw. bullock.  qk-harv.  january    3.4   2.9   1.7  55.1                                15.9  16.0   2.9
 wheat.  pcw. bullock.  qk-harv.  january    4.3   3.9   3.9  57.3                                18.6  23.9   3.9
 wheat.  psw. bullock.  qk-harv.  january    4.4   2.7   2.4  58.4                                18.3  23.1   4.0
 wheat.  prw. bullock.  qk-harv.  january    4.1   2.4   2.4  58.6                                      33.9   2.4
 wheat.  scwn.bullock.  qk-harv.  january    3.2   3.2  19.2   8.3                                18.9  15.1   3.2
 wheat.  srwn.bullock.  qk-harv.  january    3.2   3.2  19.3   9.6                                20.2  17.3   3.2
 wheat.  scws.bullock.  qk-harv.  january    3.2   3.2  24.2                                      16.2  20.1   3.2
 wheat.  srws.bullock.  qk-harv.  january    3.2   3.2  22.3                                      16.5  19.7   3.2
 wheat.  nwfp.bullock.  qk-harv.  light      5.8   3.7   3.7  78.2                                22.3  21.2   3.7
 wheat.  pmw. bullock.  qk-harv.  light      3.4   2.9   1.7  72.0                                15.9  16.0   2.9
 wheat.  pcw. bullock.  qk-harv.  light      4.3   3.9   3.9  74.9                                18.6  23.9   3.9
 wheat.  psw. bullock.  qk-harv.  light      4.4   2.7   2.4  76.3                                18.3  23.1   4.0
 wheat.  prw. bullock.  qk-harv.  light      4.1   2.4   2.4  76.6                                      33.9   2.4
 wheat.  scwn.bullock.  qk-harv.  light      3.2   3.2  25.1  10.8                                18.9  15.1   3.2
 wheat.  srwn.bullock.  qk-harv.  light      3.2   3.2  25.2  12.6                                20.2  17.3   3.2
 wheat.  scws.bullock.  qk-harv.  light      3.2   3.2  31.7                                      16.2  20.1   3.2
 wheat.  srws.bullock.  qk-harv.  light      3.2   3.2  29.3                                      16.5  19.7   3.2
 wheat.  nwfp.bullock.  qk-harv.  standard   5.8   3.7   3.7  91.9                                22.3  21.2   3.7
 wheat.  pmw. bullock.  qk-harv.  standard   3.4   2.9   1.7  84.8                                15.9  16.0   2.9
 wheat.  pcw. bullock.  qk-harv.  standard   4.3   3.9   3.9  88.1                                18.6  23.9   3.9
 wheat.  psw. bullock.  qk-harv.  standard   4.4   2.7   2.4  89.8                                18.3  23.1   4.0
 wheat.  prw. bullock.  qk-harv.  standard   4.1   2.4   2.4  90.1                                      33.9   2.4
 wheat.  scwn.bullock.  qk-harv.  standard   3.2   3.2  29.5  12.7                                18.9  15.1   3.2
 wheat.  srwn.bullock.  qk-harv.  standard   3.2   3.2  29.7  14.8                                20.2  17.3   3.2
 wheat.  scws.bullock.  qk-harv.  standard   3.2   3.2  37.3                                      16.2  20.1   3.2
 wheat.  srws.bullock.  qk-harv.  standard   3.2   3.2  34.3                                      16.5  19.7   3.2
 wheat.  nwfp.bullock.  standard. heavy      5.8   3.7   3.7  43.6  16.1                          22.3  21.2   3.7
 wheat.  pmw. bullock.  standard. heavy      3.4   2.9   1.7  41.3  13.8                          15.9  16.0   2.9
 wheat.  pcw. bullock.  standard. heavy      4.3   3.9   3.9  41.9  15.5                          18.6  23.9   3.9
 wheat.  psw. bullock.  standard. heavy      4.4   2.7   2.4  40.9  17.5                          18.3  23.1   4.0
 wheat.  prw. bullock.  standard. heavy      4.1   2.4   2.4  41.0  17.6                                33.9   2.4
 wheat.  scwn.bullock.  standard. heavy      3.2   3.2        27.4                                18.9  15.1   3.2
 wheat.  srwn.bullock.  standard. heavy      3.2   3.2        28.9                                20.2  17.3   3.2
 wheat.  scws.bullock.  standard. heavy      3.2   3.2  15.9   8.3                                16.2  20.1   3.2
 wheat.  srws.bullock.  standard. heavy      3.2   3.2  14.2   8.0                                16.5  19.7   3.2
 wheat.  nwfp.bullock.  standard. january    5.8   3.7   3.7  43.6  16.1                          22.3  21.2   3.7
 wheat.  pmw. bullock.  standard. january    3.4   2.9   1.7  41.3  13.8                          15.9  16.0   2.9
 wheat.  pcw. bullock.  standard. january    4.3   3.9   3.9  41.9  15.5                          18.6  23.9   3.9
 wheat.  psw. bullock.  standard. january    4.4   2.7   2.4  40.9  17.5                          18.3  23.1   4.0
 wheat.  prw. bullock.  standard. january    4.1   2.4   2.4  41.0  17.6                                33.9   2.4
 wheat.  scwn.bullock.  standard. january    3.2   3.2        27.4                                18.9  15.1   3.2
 wheat.  srwn.bullock.  standard. january    3.2   3.2        18.9                                20.2  17.3   3.2
 wheat.  scws.bullock.  standard. january    3.2   3.2  15.9   8.3                                16.2  20.1   3.2
 wheat.  srws.bullock.  standard. january    3.2   3.2  14.2   8.0                                16.5  19.7   3.2
 wheat.  nwfp.bullock.  standard. light      5.8   3.7   3.7  57.0  21.1                          22.3  21.2   3.7
 wheat.  pmw. bullock.  standard. light      3.4   2.9   1.7  54.1  18.0                          15.9  16.0   2.9
 wheat.  pcw. bullock.  standard. light      4.3   3.9   3.9  54.7  20.2                          18.6  23.9   3.9
 wheat.  psw. bullock.  standard. light      4.4   2.7   2.4  53.5  22.9                          18.3  23.1   4.0
 wheat.  prw. bullock.  standard. light      4.1   2.4   2.4  53.6  23.0                                33.9   2.4
 wheat.  scwn.bullock.  standard. light      3.2   3.2        35.9                                18.9  15.1   3.2
 wheat.  srwn.bullock.  standard. light      3.2   3.2        37.8                                20.2  17.3   3.2
 wheat.  scws.bullock.  standard. light      3.2   3.2  20.8  10.9                                16.2  20.1   3.2
 wheat.  srws.bullock.  standard. light      3.2   3.2  18.6  10.7                                16.5  19.7   3.2
 wheat.  nwfp.bullock.  standard. standard   5.8   3.7   3.7  67.1  24.8                          22.3  21.2   3.7
 wheat.  pmw. bullock.  standard. standard   3.4   2.9   1.7  63.6  21.2                          15.9  16.0   2.9
 wheat.  pcw. bullock.  standard. standard   4.3   3.9   3.9  64.4  23.8                          18.6  23.9   3.9
 wheat.  psw. bullock.  standard. standard   4.4   2.7   2.4  62.9  26.9                          18.3  23.1   4.0
 wheat.  prw. bullock.  standard. standard   4.1   2.4   2.4  63.1  27.0                                33.9   2.4
 wheat.  scwn.bullock.  standard. standard   3.2   3.2        42.2                                18.9  15.1   3.2
 wheat.  srwn.bullock.  standard. standard   3.2   3.2        44.5                                20.2  17.3   3.2
 wheat.  scws.bullock.  standard. standard   3.2   3.2  24.5  12.8                                16.2  20.1   3.2
 wheat.  srws.bullock.  standard. standard   3.2   3.2  21.9  12.4                                16.5  19.7   3.2
 wheat.  nwfp.semi-mech.la-plant. heavy      5.8   3.7   3.7  37.2  12.4                                13.8   3.7
 wheat.  pmw. semi-mech.la-plant. heavy      3.4   2.9   1.7  33.1  11.1                                 2.3  10.7
 wheat.  pcw. semi-mech.la-plant. heavy      4.3   3.9   3.9  34.3  11.2                                 8.0  11.5
 wheat.  psw. semi-mech.la-plant. heavy      4.4   2.7   2.4  32.4  13.9                                 3.4  15.6
 wheat.  prw. semi-mech.la-plant. heavy      4.1   2.4   2.4  32.3  13.8                                      11.9
 wheat.  scwn.semi-mech.la-plant. heavy      3.2   3.2        22.4                                       3.1   9.1
 wheat.  srwn.semi-mech.la-plant. heavy      3.2   3.2        25.4                                       2.9   9.6
 wheat.  scws.semi-mech.la-plant. heavy      3.2   3.2  14.0   7.3                                       3.2  11.0
 wheat.  srws.semi-mech.la-plant. heavy      3.2   3.2  10.7   7.1                                       3.2  11.0
 wheat.  nwfp.semi-mech.la-plant. january    5.8   3.7   3.7  37.2  12.4                                13.8   3.7
 wheat.  pmw. semi-mech.la-plant. january    3.4   2.9   1.7  33.1  11.1                                 2.3  10.7
 wheat.  pcw. semi-mech.la-plant. january    4.3   3.9   3.9  34.3  11.2                                 8.0  11.5
 wheat.  psw. semi-mech.la-plant. january    4.4   2.7   2.4  32.4  13.9                                 3.4  15.6
 wheat.  prw. semi-mech.la-plant. january    4.1   2.4   2.4  32.3  13.8                                      11.9
 wheat.  scwn.semi-mech.la-plant. january    3.2   3.2        22.9                                       3.1   9.1
 wheat.  srwn.semi-mech.la-plant. january    3.2   3.2        25.4                                       2.9   9.6
 wheat.  scws.semi-mech.la-plant. january    3.2   3.2  14.0   7.3                                       3.2  11.0
 wheat.  srws.semi-mech.la-plant. january    3.2   3.2  10.7   7.1                                       3.2  11.0
 wheat.  nwfp.semi-mech.la-plant. light      5.8   3.7   3.7  48.6  16.2                                13.8   3.7
 wheat.  pmw. semi-mech.la-plant. light      3.4   2.9   1.7  43.3  14.4                                 2.3  10.7
 wheat.  pcw. semi-mech.la-plant. light      4.3   3.9   3.9  44.9  14.6                                 8.0  11.5
 wheat.  psw. semi-mech.la-plant. light      4.4   2.7   2.4  42.4  18.2                                 3.4  15.6
 wheat.  prw. semi-mech.la-plant. light      4.1   2.4   2.4  42.2  18.1                                      11.9
 wheat.  scwn.semi-mech.la-plant. light      3.2   3.2        31.4                                       3.1   9.1
 wheat.  srwn.semi-mech.la-plant. light      3.2   3.2        33.2                                       2.9   9.6
 wheat.  scws.semi-mech.la-plant. light      3.2   3.2  18.3   9.6                                       3.2  11.0
 wheat.  srws.semi-mech.la-plant. light      3.2   3.2  13.9   9.3                                       3.2  11.0
 wheat.  nwfp.semi-mech.la-plant. standard   5.8   3.7   3.7  57.2  19.1                                13.8   3.7
 wheat.  pmw. semi-mech.la-plant. standard   3.4   2.9   1.7  51.0  17.0                                 2.3  10.7
 wheat.  pcw. semi-mech.la-plant. standard   4.3   3.9   3.9  52.8  17.2                                 8.0  11.5
 wheat.  psw. semi-mech.la-plant. standard   4.4   2.7   2.4  49.4  21.4                                 3.4  15.6
 wheat.  prw. semi-mech.la-plant. standard   4.1   2.4   2.4  49.7  21.3                                      11.9
 wheat.  scwn.semi-mech.la-plant. standard   3.2   3.2        37.0                                       3.1   9.1
 wheat.  srwn.semi-mech.la-plant. standard   3.2   3.2        39.1                                       2.9   9.6
 wheat.  scws.semi-mech.la-plant. standard   3.2   3.2  21.5  11.3                                       3.2  11.0
 wheat.  srws.semi-mech.la-plant. standard   3.2   3.2  16.4  10.9                                       3.2  11.0
 wheat.  nwfp.semi-mech.qk-harv.  heavy      5.8   3.7   3.7  56.4                                 4.1   9.7   3.7
 wheat.  pmw. semi-mech.qk-harv.  heavy      3.4   2.9   1.7  50.2                                 2.3   7.8   2.9
 wheat.  pcw. semi-mech.qk-harv.  heavy      4.3   3.9   3.9  44.0                                      15.6   3.9
 wheat.  psw. semi-mech.qk-harv.  heavy      4.4   2.7   2.4  58.4                                 3.4  11.6   4.0
 wheat.  prw. semi-mech.qk-harv.  heavy      4.1   2.4   2.4  52.4                                       9.5   2.4
 wheat.  scwn.semi-mech.qk-harv.  heavy      3.2   3.2  19.2   8.2                                       9.0   3.2
 wheat.  srwn.semi-mech.qk-harv.  heavy      3.2   3.2  19.3   9.6                                       9.3   3.2
 wheat.  scws.semi-mech.qk-harv.  heavy      3.2   3.2  24.2                                            11.0   3.2
 wheat.  srws.semi-mech.qk-harv.  heavy      3.2   3.2  22.3                                            11.0   3.2
 wheat.  nwfp.semi-mech.qk-harv.  january    5.8   3.7   3.7  56.4                                 4.1   9.7   3.7
 wheat.  pmw. semi-mech.qk-harv.  january    3.4   2.9   1.7  50.2                                 2.3   7.8   2.9
 wheat.  pcw. semi-mech.qk-harv.  january    4.3   3.9   3.9  44.0                                      15.6   3.9
 wheat.  psw. semi-mech.qk-harv.  january    4.4   2.7   2.4  58.4                                 3.4  11.6   4.0
 wheat.  prw. semi-mech.qk-harv.  january    4.1   2.4   2.4  52.4                                       9.5   2.4
 wheat.  scwn.semi-mech.qk-harv.  january    3.2   3.2  19.2   8.2                                       9.0   3.2
 wheat.  srwn.semi-mech.qk-harv.  january    3.2   3.2  19.3   9.6                                       9.3   3.2
 wheat.  scws.semi-mech.qk-harv.  january    3.2   3.2  24.2                                            11.0   3.2
 wheat.  srws.semi-mech.qk-harv.  january    3.2   3.2  22.3                                            11.0   3.2
 wheat.  nwfp.semi-mech.qk-harv.  light      5.8   3.7   3.7  73.8                                 4.1   9.7   3.7
 wheat.  pmw. semi-mech.qk-harv.  light      3.4   2.9   1.7  65.7                                 2.3   7.8   2.9
 wheat.  pcw. semi-mech.qk-harv.  light      4.3   3.9   3.9  67.7                                      15.6   3.9
 wheat.  psw. semi-mech.qk-harv.  light      4.4   2.7   2.4  68.9                                 3.4  11.6   4.0
 wheat.  prw. semi-mech.qk-harv.  light      4.1   2.4   2.4  68.6                                       9.5   2.4
 wheat.  scwn.semi-mech.qk-harv.  light      3.2   3.2  25.1  10.8                                       9.0   3.2
 wheat.  srwn.semi-mech.qk-harv.  light      3.2   3.2  25.2  12.6                                       9.3   3.2
 wheat.  scws.semi-mech.qk-harv.  light      3.2   3.2  31.7                                            11.0   3.2
 wheat.  srws.semi-mech.qk-harv.  light      3.2   3.2  29.3                                            11.0   3.2
 wheat.  nwfp.semi-mech.qk-harv.  standard   5.8   3.7   3.7  86.8                                 4.1   9.7   3.7
 wheat.  pmw. semi-mech.qk-harv.  standard   3.4   2.9   1.7  77.3                                 2.3   7.8   2.9
 wheat.  pcw. semi-mech.qk-harv.  standard   4.3   3.9   3.9  79.6                                      15.6   3.9
 wheat.  psw. semi-mech.qk-harv.  standard   4.4   2.7   2.4  81.1                                 3.4  11.6   4.0
 wheat.  prw. semi-mech.qk-harv.  standard   4.1   2.4   2.4  80.7                                       9.5   2.4
 wheat.  scwn.semi-mech.qk-harv.  standard   3.2   3.2  29.5  12.7                                       9.0   3.2
 wheat.  srwn.semi-mech.qk-harv.  standard   3.2   3.2  29.7  14.8                                       9.3   3.2
 wheat.  scws.semi-mech.qk-harv.  standard   3.2   3.2  37.3                                            11.0   3.2
 wheat.  srws.semi-mech.qk-harv.  standard   3.2   3.2  34.3                                            11.0   3.2
 wheat.  nwfp.semi-mech.standard. heavy      5.8   3.7   3.7  42.3  14.1                           4.1   9.7   3.7
 wheat.  pmw. semi-mech.standard. heavy      3.4   2.9   1.7  37.7  12.5                           2.3   7.8   2.9
 wheat.  pcw. semi-mech.standard. heavy      4.3   3.9   3.9  39.1  12.7                                15.6   3.9
 wheat.  psw. semi-mech.standard. heavy      4.4   2.7   2.4  36.9  15.8                           3.4  11.6   4.0
 wheat.  prw. semi-mech.standard. heavy      4.1   2.4   2.4  36.7  15.7                                 9.5   2.4
 wheat.  scwn.semi-mech.standard. heavy      3.2   3.2        27.4                                       9.0   3.2
 wheat.  srwn.semi-mech.standard. heavy      3.2   3.2        28.9                                       9.3   3.2
 wheat.  scws.semi-mech.standard. heavy      3.2   3.2  15.9   8.3                                      11.0   3.2
 wheat.  srws.semi-mech.standard. heavy      3.2   3.2  14.2   8.0                                      11.0   3.2
 wheat.  nwfp.semi-mech.standard. january    5.8   3.7   3.7  42.3  14.1                           4.1   9.7   3.7
 wheat.  pmw. semi-mech.standard. january    3.4   2.9   1.7  37.7  12.5                           2.3   7.8   2.9
 wheat.  pcw. semi-mech.standard. january    4.3   3.9   3.9  39.1  12.7                                15.6   3.9
 wheat.  psw. semi-mech.standard. january    4.4   2.7   2.4  36.9  15.8                           3.4  11.6   4.0
 wheat.  prw. semi-mech.standard. january    4.1   2.4   2.4  36.7  15.7                                 9.5   2.4
 wheat.  scwn.semi-mech.standard. january    3.2   3.2        27.4                                       9.0   3.2
 wheat.  srwn.semi-mech.standard. january    3.2   3.2        28.9                                       9.3   3.2
 wheat.  scws.semi-mech.standard. january    3.2   3.2  15.9   8.3                                      11.0   3.2
 wheat.  srws.semi-mech.standard. january    3.2   3.2  14.2   8.0                                      11.0   3.2
 wheat.  nwfp.semi-mech.standard. light      5.8   3.7   3.7  55.3  18.4                           4.1   9.7   3.7
 wheat.  pmw. semi-mech.standard. light      3.4   2.9   1.7  49.3  16.4                           2.3   7.8   2.9
 wheat.  pcw. semi-mech.standard. light      4.3   3.9   3.9  51.1  16.6                                15.6   3.9
 wheat.  psw. semi-mech.standard. light      4.4   2.7   2.4  48.3  20.7                           3.4  11.6   4.0
 wheat.  prw. semi-mech.standard. light      4.1   2.4   2.4  48.0  20.6                                 9.5   2.4
 wheat.  scwn.semi-mech.standard. light      3.2   3.2        35.9                                       9.0   3.2
 wheat.  srwn.semi-mech.standard. light      3.2   3.2        37.8                                       9.3   3.2
 wheat.  scws.semi-mech.standard. light      3.2   3.2  20.8  10.9                                      11.0   3.2
 wheat.  srws.semi-mech.standard. light      3.2   3.2  18.6  10.7                                      11.0   3.2
 wheat.  nwfp.semi-mech.standard. standard   5.8   3.7   3.7  65.1  21.7                           4.1   9.7   3.7
 wheat.  pmw. semi-mech.standard. standard   3.4   2.9   1.7  58.0  19.3                           2.3   7.8   2.9
 wheat.  pcw. semi-mech.standard. standard   4.3   3.9   3.9  60.1  19.5                                15.6   3.9
 wheat.  psw. semi-mech.standard. standard   4.4   2.7   2.4  56.8  24.3                           3.4  11.6   4.0
 wheat.  prw. semi-mech.standard. standard   4.1   2.4   2.4  56.5  24.2                                 9.5   2.4
 wheat.  scwn.semi-mech.standard. standard   3.2   3.2        42.2                                       9.0   3.2
 wheat.  srwn.semi-mech.standard. standard   3.2   3.2        44.5                                       9.3   3.2
 wheat.  scws.semi-mech.standard. standard   3.2   3.2  24.5  12.8                                      11.0   3.2
 wheat.  srws.semi-mech.standard. standard   3.2   3.2  21.9  12.4                                      11.0   3.2

 orchard.(nwfp,pmw,pcw,psw,prw).
 (bullock,semi-mech).standard.standard        41  41.8  17.3  14.5  14.5   22.   22.    50.  29.   39.  56.5  28.3
 orchard.(scwn,scws,srwn,srws).
 (bullock,semi-mech).standard.standard        31  31.   31.   31.   25.    22.   22.    22.  31.   31.  31.   31.

 potatoes.(scwn,scws,srwn,srws).
               semi-mech.standard.standard   100   100                                        50    50   100    50
 potatoes.(nwfp,pmw,psw,pcw,prw).semi-mech.
                         standard.standard    70    24         150                                               8
 onions.(nwfp,pmw,psw,pcw,prw).semi-mech.
                         standard.standard   160     8    32    32   160
 onions.(scwn,scws,srwn,srws).semi-mech.
                         standard.standard   160                                             160     8    32    32
 chilli.(nwfp,pmw,pcw,psw,prw).semi-mech.
                         standard.standard               100     2    20    10     20    10   20     8     8
 chilli.(scwn,scws,srwn,srws).semi-mech.
                         standard.standard    16    32    28    28                                             100
$eject
     table water(c,z,t,s,w,m) water requirements(acre feet per acre)

                                                              jan  feb  mar  apr  may  jun  jul  aug  sep  oct  nov  dec
 basmati.(pmw,pcw,psw).(bullock,semi-mech).standard.standard                           .37  .84  .93  .93  .47
 basmati.prw.(bullock,semi-mech).standard.standard                                     .56  .84  .93  .93  .56

 rab-fod.(psw,prw).(bullock,semi-mech).standard.standard      .20  .30  .50  .60  .35                      .25  .15  .15
 rab-fod.(scws,srws).(bullock,semi-mech).standard.standard    .27  .37  .53  .64                      .21  .32  .27  .21
 rab-fod.srwn.(bullock,semi-mech).standard.heavy              .20                                     .20  .25  .15  .15
 rab-fod.srws.(bullock,semi-mech).standard.heavy              .27                                     .21  .32  .27  .21
 rab-fod.srwn.(bullock,semi-mech).standard.light              .20  .30                                .20  .25  .15  .15
 rab-fod.srws.(bullock,semi-mech).standard.light              .27  .37                                .21  .32  .27  .21
 rab-fod.nwfp.(bullock,semi-mech).standard.standard           .20  .25  .40  .45  .15                 .20  .25  .15  .15
 rab-fod.pmw.(bullock,semi-mech).standard.standard            .20  .25  .40  .50  .10                      .25  .15  .15
 rab-fod.pcw.(bullock,semi-mech).standard.standard            .20  .30  .45  .55  .25                 .20  .25  .15  .15
 rab-fod.scwn.(bullock,semi-mech).standard.standard           .21  .32  .48  .59                      .21  .27  .16  .16
 rab-fod.srwn.(bullock,semi-mech).standard.standard           .20  .30  .45  .55                      .20  .25  .15  .15

 cotton. pcw. bullock.  el-plant. standard                                   .19  .09  .37  .33  .37  .47  .33  .09
 cotton.(pmw,psw,prw).(bullock,semi-mech).standard.standard                       .19  .23  .28  .37  .37  .23  .09
 cotton.pcw.(bullock,semi-mech).standard.standard                                 .19  .28  .33  .37  .47  .33  .09
 cotton. scwn.bullock.  standard. standard                                   .10  .20  .30  .40  .40  .40  .30
 cotton.scws.(bullock,semi-mech).standard.standard                           .20  .20  .30  .40  .40  .40  .20
 cotton.srws.(bullock,semi-mech).standard.standard                           .20  .20  .30  .40  .40  .40  .15
 cotton. pcw. semi-mech.la-plant. standard                                             .37  .33  .37  .47  .33  .09
 cotton. scwn.semi-mech.standard. standard                                        .25  .35  .40  .40  .40  .30

 gram.nwfp.(bullock,semi-mech).standard.standard              .10                                     .20  .10  .10  .05
 gram.pmw.(bullock,semi-mech).standard.standard               .05  .05                                .20  .10  .10  .15
 gram.pcw.(bullock,semi-mech).standard.standard               .10  .05                                .25  .15  .10  .10
 gram.(psw,prw).(bullock,semi-mech).standard.standard         .10  .10                                     .20  .15  .10
 gram.(scwn,scws).(bullock,semi-mech).standard.standard       .11  .21                                     .21  .21  .11
 gram.(srwn,srws).(bullock,semi-mech).standard.standard       .10  .20                                     .20  .20  .10

 irri.pmw.(bullock,semi-mech).standard.standard                                        .37  .84  .93  .93  .47
 irri.pcw.(bullock,semi-mech).standard.standard                                        .37  .84  .93  .93  .47
 irri.psw.(bullock,semi-mech).standard.standard                                        .33  .79  .89  .89  .47
 irri.prw.(bullock,semi-mech).standard.standard                                        .56  .84  .93  .93  .65
 irri.scwn.(bullock,semi-mech).standard.standard                                       .36  .90 1.01 1.01  .52
 irri.srwn.(bullock,semi-mech).standard.standard                                       .78  .96  .96  .96  .45
 irri.scws.(bullock,semi-mech).standard.standard                                  .20  .85 1.00 1.00  .75
 irri.srws.(bullock,semi-mech).standard.standard                                  .20 1.15 1.15 1.15  .75

 maize.nwfp.(bullock,semi-mech).standard.standard                                      .15  .35  .30  .25  .20
 maize.(pcw,psw,prw).bullock.   standard.standard                                      .10  .10  .30  .30  .40  .20
 maize.(pcw,psw,prw).semi-mech. standard.standard                                           .15  .30  .30  .40  .25
 maize.scwn.(bullock,semi-mech).standard.standard                                      .21  .32  .32  .43  .43  .21
 maize.scws.(bullock,semi-mech).standard.standard                                 .21  .21  .21  .32  .43  .32  .21

 mus+rap.nwfp.(bullock,semi-mech).standard.standard                                         .20  .20  .25  .25  .25
 mus+rap.pmw.(bullock,semi-mech).standard.standard            .20  .20                                     .25  .20  .20
 mus+rap.pcw.(bullock,semi-mech).standard.standard            .20  .15                                     .30  .20  .20
 mus+rap.psw.(bullock,semi-mech).standard.standard            .20  .20                                .20  .20  .20  .20
 mus+rap.prw.(bullock,semi-mech).standard.standard            .20                                     .20  .20  .20  .20
 mus+rap.scwn.(bullock,semi-mech).standard.standard           .21  .21                                .21  .21  .21  .21
 mus+rap.srwn.(bullock,semi-mech).standard.standard           .20  .20                                .20  .20  .20  .20
 mus+rap.srws.(bullock,semi-mech).standard.standard           .20  .20                                .20  .20  .20  .20
 mus+rap.scws.(bullock,semi-mech).standard.standard           .11  .11                                .27  .21  .21  .27

 sc-gur.nwfp.(bullock,semi-mech).standard.standard            .19  .19  .28  .37  .37  .37  .28  .47  .47  .37  .19  .28
 sc-gur. pmw. (bullock,semi-mech).standard.standard           .14  .09  .23  .33  .33  .28  .23  .37  .37  .19  .09  .19
 sc-gur. pcw. (bullock,semi-mech).standard.standard           .09  .14  .23  .37  .37  .33  .37  .47  .47  .28  .37  .19
 sc-gur. psw. (bullock,semi-mech).standard.standard           .09  .09  .23  .33  .33  .33  .37  .47  .47  .23  .23  .09
 sc-gur. prw. (bullock,semi-mech).standard.standard           .14  .19  .28  .37  .37  .28  .33  .37  .37  .28  .19  .19
 sc-gur. scwn.(bullock,semi-mech).standard.standard           .25  .30  .40  .40  .40  .55  .60  .60  .50  .40  .30  .25
 sc-gur. srwn.(bullock,semi-mech).standard.standard           .23  .28  .37  .37  .37  .51  .56  .56  .47  .37  .28  .23
 sc-gur. scws.(bullock,semi-mech).standard.standard           .25  .35  .40  .40  .40  .55  .65  .65  .60  .50  .35  .25
 sc-gur. srws.(bullock,semi-mech).standard.standard           .25  .35  .40  .40  .50  .55  .65  .65  .60  .50  .40  .20

 sc-mill.nwfp.(bullock,semi-mech).standard.standard           .19  .19  .28  .33  .33  .37  .37  .47  .47  .37  .28  .19
 sc-mill.pmw.(bullock,semi-mech).standard.standard            .14  .09  .23  .33  .33  .28  .23  .37  .37  .19  .09  .19
 sc-mill.pcw. (bullock,semi-mech).standard.standard           .09  .14  .23  .37  .37  .33  .37  .47  .47  .28  .37  .19
 sc-mill.psw. (bullock,semi-mech).standard.standard           .09  .09  .23  .33  .33  .33  .37  .47  .47  .23  .23  .09
 sc-mill.prw. (bullock,semi-mech).standard.standard           .14  .19  .28  .37  .37  .28  .33  .37  .37  .28  .19  .19
 sc-mill.scwn.(bullock,semi-mech).standard.standard           .25  .30  .40  .40  .40  .55  .60  .60  .50  .40  .30  .25
 sc-mill.srwn.(bullock,semi-mech).standard.standard           .23  .28  .37  .37  .37  .51  .56  .56  .47  .37  .28  .23
 sc-mill.scws.(bullock,semi-mech).standard.standard           .25  .35  .40  .40  .40  .55  .65  .65  .60  .50  .35  .25
 sc-mill.srws.(bullock,semi-mech).standard.standard           .25  .35  .40  .40  .50  .55  .65  .65  .60  .50  .40  .20

 kha-fod.srws.(bullock,semi-mech).la-plant.standard                               .21  .27  .32  .32  .16  .11
 kha-fod.nwfp.(bullock,semi-mech).standard.standard                          .05  .15  .20  .20  .25  .25  .05
 kha-fod.pmw. (bullock,semi-mech).standard.standard                     .05  .05  .20  .20  .30  .25  .15
 kha-fod.pcw. (bullock,semi-mech).standard.standard                     .10  .15  .25  .30  .40  .30  .20
 kha-fod.psw. (bullock,semi-mech).standard.standard                     .15  .25  .35  .35  .30  .25  .20
 kha-fod.prw. (bullock,semi-mech).standard.standard                     .10  .20  .30  .30  .30  .30  .30
 kha-fod.scwn.(bullock,semi-mech).standard.standard                     .11  .11  .21  .21  .32  .32  .21  .11
 kha-fod.scws.(bullock,semi-mech).standard.standard                     .11  .21  .21  .27  .32  .32  .16  .05
 kha-fod.srwn.(bullock,semi-mech).standard.standard                     .21  .32  .27  .27  .32       .21  .32
 kha-fod.srws.(bullock,semi-mech).standard.standard                     .21  .32  .27  .27  .32       .21  .32

 wheat.  nwfp.bullock.  la-plant. heavy                       .28                                               .25  .25
 wheat.  pmw. bullock.  la-plant. heavy                       .22                                               .28  .16
 wheat.(pcw,psw,srwn).bullock.la-plant.heavy                  .28                                               .28  .19
 wheat.  prw. bullock.  la-plant. heavy                       .28                                               .09  .28
 wheat.(scwn,srws).bullock.la-plant.heavy                     .30                                               .30  .20
 wheat.  scws.bullock.  la-plant. heavy                       .30                                               .30  .35
 wheat.  nwfp.bullock.  la-plant. january                          .28                                          .25  .25
 wheat.  pmw. bullock.  la-plant. january                          .28                                          .28  .16
 wheat.(pcw,psw,srwn).bullock.la-plant.january                     .28                                          .28  .19
 wheat.  prw. bullock.  la-plant. january                          .28                                          .09  .28
 wheat.(scwn,srws).bullock.la-plant.january                        .30                                          .30  .20
 wheat.  scws.bullock.  la-plant. january                          .30                                          .30  .35
 wheat.  nwfp.bullock.  la-plant. light                       .22  .35                                          .25  .25
 wheat.  pmw. bullock.  la-plant. light                       .22  .28                                          .28  .16
 wheat.(pcw,psw).     bullock.la-plant.light                  .23  .37                                          .28  .19
 wheat.  srwn.bullock.  la-plant. light                       .19  .33                                          .28  .19
 wheat.  prw. bullock.  la-plant. light                       .23  .37                                          .09  .28
 wheat.  scwn.bullock.  la-plant. light                       .20  .40                                          .30  .20
 wheat.  scws.bullock.  la-plant. light                       .25  .40                                          .30  .35
 wheat.  srws.bullock.  la-plant. light                       .20  .35                                          .30  .20
 wheat.  nwfp.bullock.  la-plant. standard                    .22  .35  .41                                     .25  .25
 wheat.  pmw. bullock.  la-plant. standard                    .22  .28  .35                                     .28  .16
 wheat.  pcw. bullock.  la-plant. standard                    .23  .37  .47                                     .28  .19
 wheat.  psw. bullock.  la-plant. standard                    .23  .37  .42                                     .28  .19
 wheat.  prw. bullock.  la-plant. standard                    .23  .37  .47                                     .09  .28
 wheat.  scwn.bullock.  la-plant. standard                    .20  .40  .40                                     .30  .20
 wheat.  srwn.bullock.  la-plant. standard                    .19  .33  .37                                     .28  .19
 wheat.  scws.bullock.  la-plant. standard                    .25  .40  .50                                     .30  .35
 wheat.  srws.bullock.  la-plant. standard                    .20  .35  .40                                     .30  .20

 wheat.nwfp.bullock.(qk-harv,standard).heavy                  .28                                          .19  .16  .16
 wheat.pmw. bullock.(qk-harv,standard).heavy                  .22                                          .19  .16  .16
 wheat.pcw. bullock.(qk-harv,standard).heavy                  .28                                          .19  .19  .19
 wheat.psw. bullock.(qk-harv,standard).heavy                  .28                                          .19  .19  .19
 wheat.prw. bullock.(qk-harv,standard).heavy                  .28                                               .19  .19
 wheat.scwn.bullock.(qk-harv,standard).heavy                  .30                                          .20  .20  .20
 wheat.srwn.bullock.(qk-harv,standard).heavy                  .28                                          .19  .19  .19
 wheat.scws.bullock.(qk-harv,standard).heavy                  .30                                          .20  .30  .25
 wheat.srws.bullock.(qk-harv,standard).heavy                  .30                                          .20  .20  .20
 wheat.nwfp.bullock.(qk-harv,standard).january                     .28                                     .19  .16  .16
 wheat.pmw. bullock.(qk-harv,standard).january                     .28                                     .19  .16  .16
 wheat.pcw. bullock.(qk-harv,standard).january                     .28                                     .19  .19  .19
 wheat.psw. bullock.(qk-harv,standard).january                     .28                                     .19  .19  .19
 wheat.prw. bullock.(qk-harv,standard).january                     .28                                          .19  .19
 wheat.scwn.bullock.(qk-harv,standard).january                     .30                                     .20  .20  .20
 wheat.srwn.bullock.(qk-harv,standard).january                     .28                                     .19  .19  .19
 wheat.scws.bullock.(qk-harv,standard).january                     .30                                     .20  .30  .25
 wheat.srws.bullock.(qk-harv,standard).january                     .30                                     .20  .20  .20
 wheat.nwfp.bullock.(qk-harv,standard).light                  .22  .35                                     .19  .16  .16
 wheat.pmw. bullock.(qk-harv,standard).light                  .22  .28                                     .19  .16  .16
 wheat.pcw. bullock.(qk-harv,standard).light                  .23  .37                                     .19  .19  .19
 wheat.psw. bullock.(qk-harv,standard).light                  .23  .37                                     .19  .19  .19
 wheat.prw. bullock.(qk-harv,standard).light                  .23  .37                                          .19  .19
 wheat.scwn.bullock.(qk-harv,standard).light                  .20  .40                                     .20  .20  .20
 wheat.srwn.bullock.(qk-harv,standard).light                  .19  .33                                     .19  .19  .19
 wheat.scws.bullock.(qk-harv,standard).light                  .25  .40                                     .20  .30  .25
 wheat.srws.bullock.(qk-harv,standard).light                  .20  .35                                     .20  .20  .20
 wheat.nwfp.bullock.(qk-harv,standard).standard               .22  .35  .41                                .19  .16  .16
 wheat.pmw. bullock.(qk-harv,standard).standard               .22  .28  .35                                .19  .16  .16
 wheat.pcw. bullock.(qk-harv,standard).standard               .23  .37  .47                                .19  .19  .19
 wheat.psw. bullock.(qk-harv,standard).standard               .23  .37  .42                                .19  .19  .19
 wheat.prw. bullock.(qk-harv,standard).standard               .23  .37  .47                                     .19  .19
 wheat.scwn.bullock.(qk-harv,standard).standard               .20  .40  .40                                .20  .20  .20
 wheat.srwn.bullock.(qk-harv,standard).standard               .19  .33  .37                                .19  .19  .19
 wheat.scws.bullock.(qk-harv,standard).standard               .25  .40  .50                                .20  .30  .25
 wheat.srws.bullock.(qk-harv,standard).standard               .20  .35  .40                                .20  .20  .20

 wheat.nwfp.semi-mech.(qk-harv,standard).heavy                .28                                          .19  .16  .16
 wheat.pmw. semi-mech.(qk-harv,standard).heavy                .22                                          .19  .16  .16
 wheat.pcw. semi-mech.(qk-harv,standard).heavy                .28                                               .28  .19
 wheat.psw. semi-mech.(qk-harv,standard).heavy                .28                                          .19  .19  .19
 wheat.prw. semi-mech.(qk-harv,standard).heavy                .28                                               .19  .19
 wheat.srwn.semi-mech.(qk-harv,standard).heavy                .28                                               .28  .19
 wheat.srws.semi-mech.(qk-harv,standard).heavy                .30                                               .30  .30
 wheat.nwfp.semi-mech.(qk-harv,standard).january                   .28                                     .19  .16  .16
 wheat.pmw. semi-mech.(qk-harv,standard).january                   .28                                     .19  .16  .16
 wheat.pcw. semi-mech.(qk-harv,standard).january                   .28                                          .28  .19
 wheat.psw. semi-mech.(qk-harv,standard).january                   .28                                     .19  .19  .19
 wheat.prw. semi-mech.(qk-harv,standard).january                   .28                                          .19  .19
 wheat.srwn.semi-mech.(qk-harv,standard).january                   .28                                          .28  .19
 wheat.srws.semi-mech.(qk-harv,standard).january                   .30                                          .30  .30
 wheat.nwfp.semi-mech.(qk-harv,standard).light                .22  .35                                     .19  .16  .16
 wheat.pmw. semi-mech.(qk-harv,standard).light                .22  .28                                     .19  .16  .16
 wheat.pcw. semi-mech.(qk-harv,standard).light                .23  .37                                          .28  .19
 wheat.psw. semi-mech.(qk-harv,standard).light                .23  .37                                     .19  .19  .19
 wheat.prw. semi-mech.(qk-harv,standard).light                .23  .37                                          .19  .19
 wheat.srwn.semi-mech.(qk-harv,standard).light                .19  .33                                          .28  .19
 wheat.srws.semi-mech.(qk-harv,standard).light                .20  .35                                          .30  .30
 wheat.nwfp.semi-mech.(qk-harv,standard).standard             .22  .35  .41                                .19  .16  .16
 wheat.pmw. semi-mech.(qk-harv,standard).standard             .22  .28  .35                                .19  .16  .16
 wheat.pcw. semi-mech.(qk-harv,standard).standard             .23  .37  .47                                     .28  .19
 wheat.psw. semi-mech.(qk-harv,standard).standard             .23  .37  .42                                .19  .19  .19
 wheat.prw. semi-mech.(qk-harv,standard).standard             .23  .37  .47                                     .19  .19

 wheat.  scwn.semi-mech.qk-harv.  heavy                       .30                                               .30  .20
 wheat.  scws.semi-mech.qk-harv.  heavy                       .30                                               .35  .30
 wheat.  scwn.semi-mech.qk-harv.  january                          .30                                          .30  .20
 wheat.  scws.semi-mech.qk-harv.  january                          .30                                          .35  .30
 wheat.  scwn.semi-mech.qk-harv.  light                       .20  .40                                          .30  .20
 wheat.  scws.semi-mech.qk-harv.  light                       .25  .40                                          .35  .30
 wheat.  scwn.semi-mech.qk-harv.  standard                    .20  .40  .40                                     .30  .20
 wheat.  srwn.semi-mech.qk-harv.  standard                    .19  .33  .37                                     .28  .19
 wheat.  scws.semi-mech.qk-harv.  standard                    .25  .40  .50                                     .35  .30
 wheat.  srws.semi-mech.qk-harv.  standard                    .20  .35  .40                                     .30  .30
 wheat.  scwn.semi-mech.standard. heavy                       .30                                               .30  .25
 wheat.  scws.semi-mech.standard. heavy                       .30                                               .40  .30
 wheat.  scwn.semi-mech.standard. january                          .30                                          .30  .25
 wheat.  scws.semi-mech.standard. january                          .30                                          .40  .30
 wheat.  scwn.semi-mech.standard. light                       .25  .40                                          .30  .25
 wheat.  scws.semi-mech.standard. light                       .30  .40                                          .30  .30
 wheat.  scwn.semi-mech.standard. standard                    .20  .40  .50                                     .30  .20
 wheat.  srwn.semi-mech.standard. standard                    .19  .33  .47                                     .28  .19
 wheat.  scws.semi-mech.standard. standard                    .25  .45  .55                                     .30  .25
 wheat.  srws.semi-mech.standard. standard                    .20  .40  .55                                     .20  .20

 wheat.  nwfp.semi-mech.la-plant. heavy                       .28                                               .25  .25
 wheat.  pmw. semi-mech.la-plant. heavy                       .22                                               .23  .23
 wheat.  pcw. semi-mech.la-plant. heavy                       .28                                               .23  .23
 wheat.  psw. semi-mech.la-plant. heavy                       .28                                               .28  .19
 wheat.  prw. semi-mech.la-plant. heavy                       .28                                               .09  .28
 wheat.  scwn.semi-mech.la-plant. heavy                       .30                                               .30  .25
 wheat.  srwn.semi-mech.la-plant. heavy                       .28                                               .28  .19
 wheat.  scws.semi-mech.la-plant. heavy                       .30                                               .35  .30
 wheat.  srws.semi-mech.la-plant. heavy                       .30                                               .20  .30
 wheat.  nwfp.semi-mech.la-plant. january                          .28                                          .25  .25
 wheat.  pmw. semi-mech.la-plant. january                          .28                                          .23  .23
 wheat.  pcw. semi-mech.la-plant. january                          .28                                          .23  .23
 wheat.  psw. semi-mech.la-plant. january                          .28                                          .28  .19
 wheat.  prw. semi-mech.la-plant. january                          .28                                          .09  .28
 wheat.  scwn.semi-mech.la-plant. january                          .30                                          .30  .25
 wheat.  srwn.semi-mech.la-plant. january                          .28                                          .28  .19
 wheat.  scws.semi-mech.la-plant. january                          .30                                          .35  .30
 wheat.  srws.semi-mech.la-plant. january                          .30                                          .20  .30
 wheat.  nwfp.semi-mech.la-plant. light                       .22  .35                                          .25  .25
 wheat.  pmw. semi-mech.la-plant. light                       .22  .28                                          .23  .23
 wheat.  pcw. semi-mech.la-plant. light                       .23  .37                                          .23  .23
 wheat.  psw. semi-mech.la-plant. light                       .23  .37                                          .28  .19
 wheat.  prw. semi-mech.la-plant. light                       .23  .37                                          .09  .28
 wheat.  scwn.semi-mech.la-plant. light                       .20  .40                                          .30  .25
 wheat.  srwn.semi-mech.la-plant. light                       .19  .33                                          .28  .19
 wheat.  scws.semi-mech.la-plant. light                       .25  .40                                          .35  .30
 wheat.  srws.semi-mech.la-plant. light                       .20  .35                                          .20  .30
 wheat.  nwfp.semi-mech.la-plant. standard                    .22  .35  .41                                     .25  .25
 wheat.  pmw. semi-mech.la-plant. standard                    .22  .28  .35                                     .23  .23
 wheat.  pcw. semi-mech.la-plant. standard                    .23  .37  .47                                     .23  .23
 wheat.  psw. semi-mech.la-plant. standard                    .23  .37  .42                                     .28  .19
 wheat.  prw. semi-mech.la-plant. standard                    .23  .37  .47                                     .09  .28
 wheat.  scwn.semi-mech.la-plant. standard                    .20  .40  .40                                     .30  .25
 wheat.  srwn.semi-mech.la-plant. standard                    .19  .33  .37                                     .28  .19
 wheat.  scws.semi-mech.la-plant. standard                    .25  .40  .50                                     .35  .30
 wheat.  srws.semi-mech.la-plant. standard                    .20  .35  .40                                     .20  .30

 orchard. nwfp.(bullock,semi-mech).standard.standard          .08  .10  .19  .33  .46  .50  .55  .58  .49  .37  .19  .09
 orchard. pmw. (bullock,semi-mech).standard.standard          .10  .13  .23  .39  .53  .56  .60  .60  .59  .45  .23  .11
 orchard. pcw. (bullock,semi-mech).standard.standard          .11  .14  .30  .44  .54  .59  .57  .61  .58  .47  .23  .12
 orchard. psw. (bullock,semi-mech).standard.standard          .08  .12  .24  .40  .49  .49  .46  .50  .51  .40  .18  .10
 orchard. prw. (bullock,semi-mech).standard.standard          .08  .12  .24  .40  .47  .47  .46  .50  .50  .39  .18  .10
 orchard. scwn.(bullock,semi-mech).standard.standard          .17  .17  .32  .45  .52  .54  .56  .60  .60  .45  .26  .14
 orchard. scws.(bullock,semi-mech).standard.standard          .18  .20  .38  .54  .61  .56  .56  .61  .62  .56  .32  .20
 orchard. srwn.(bullock,semi-mech).standard.standard          .15  .16  .32  .45  .54  .52  .54  .59  .57  .47  .26  .15
 orchard. srws.(bullock,semi-mech).standard.standard          .19  .21  .39  .57  .61  .54  .55  .59  .63  .58  .33  .21

 potatoes.(scwn,scws,srwn,srws).
               semi-mech.standard.standard                                                            .25  .667  .5  .25
 potatoes.(nwfp,pmw,pcw,psw,prw).semi-mech.
                         standard.standard                    .25  .25   .5                                           .5
 onions.(nwfp,pmw,pcw,psw,prw).semi-mech.
                         standard.standard                    .25  .25   .5   .5  .25
 onions.(scwn,scws,srwn,srws).semi-mech.
                         standard.standard                    .25                                     .25  .25   .5   .5
 chilli.(nwfp,pmw,psw,pcw,prw).semi-mech.
                         standard.standard                             .167  .25 .667 .667   .5   .5 .333 .333 .167
 chilli.(scwn,scws,srwn,srws).semi-mech.
                         standard.standard                     .5 .667 .667 .667                                     .25
$eject
   table tractor(c,z,t,s,w,m) tractor requirements (tractor hours per acre)

                                              jan   feb   mar   apr   may   jun   jul   aug   sep   oct   nov   dec
 basmati.pmw. semi-mech.standard.standard                                  2.70  1.70                    1.50  0.50
 basmati.pcw. semi-mech.standard.standard                                  2.70  1.70                    1.50  0.50
 basmati.psw. semi-mech.standard.standard                                  3.30  1.80                    1.50  0.50
 basmati.prw. semi-mech.standard.standard                                  3.10  2.00                    1.50  0.50
 rab-fod.srwn.semi-mech.standard.heavy                                                       1.00  2.80
 rab-fod.srws.semi-mech.standard.heavy                                                       1.50  2.50
 rab-fod.srwn.semi-mech.standard.light                                                       1.00  2.80
 rab-fod.srws.semi-mech.standard.light                                                       1.50  2.50
 rab-fod.nwfp.semi-mech.standard.standard    1.50  1.50  1.50  1.50  1.50                    2.05  1.62  1.50  1.50
 rab-fod.pmw. semi-mech.standard.standard    1.50  1.50  1.50  1.50  1.50                          1.63  0.78  1.50
 rab-fod.pcw. semi-mech.standard.standard    1.50  1.50  1.50  1.50  1.50                    1.26  2.01  2.29  1.50
 rab-fod.psw. semi-mech.standard.standard    3.00  3.00  3.00  1.50  1.50                          1.91  1.23  1.50
 rab-fod.prw. semi-mech.standard.standard    3.00  3.00  3.00  1.50  1.50                          1.49  1.05  1.57
 rab-fod.scwn.semi-mech.standard.standard                                                    1.20  2.60
 rab-fod.srwn.semi-mech.standard.standard                                                    1.00  2.80
 rab-fod.scws.semi-mech.standard.standard                                                    1.50  2.40
 rab-fod.srws.semi-mech.standard.standard                                                    1.50  2.50
 cotton. pcw. semi-mech.la-plant.standard                                  4.70              0.50  0.50  0.50  0.50
 cotton. pmw. semi-mech.standard.standard                            3.06  0.64              0.50  0.50  0.50  0.50
 cotton. pcw. semi-mech.standard.standard                            2.49  2.22              0.50  0.50  0.50  0.50
 cotton. psw. semi-mech.standard.standard                            1.66  1.98              0.50  0.50  0.50  0.50
 cotton. prw. semi-mech.standard.standard                            1.45  1.84              0.50  0.50  0.50  0.50
 cotton. scwn.semi-mech.standard.standard                           12.10                    1.00  1.00  2.00
 cotton. scws.semi-mech.standard.standard                     12.00                          1.00  2.00
 cotton. srws.semi-mech.standard.standard                           12.40                    1.00  2.00
 gram.   nwfp.semi-mech.standard.standard                      3.75                          1.37  2.00
 gram.   pmw. semi-mech.standard.standard                3.75                                1.27  0.78  0.90  2.50
 gram.   pcw. semi-mech.standard.standard                3.75                                1.30  0.79
 gram.   psw. semi-mech.standard.standard                3.75                                      1.25  0.70
 gram.   prw. semi-mech.standard.standard                3.75                                      1.14  0.68
 gram.(scwn,scws).semi-mech.standard.standard                  3.50                          0.81  1.29
 gram.(srwn,srws).semi-mech.standard.standard                  3.80                          0.96  1.39

 irri.   pmw. semi-mech.standard.standard                                  2.60  2.20                    1.00  0.50
 irri.   pcw. semi-mech.standard.standard                                  2.60  2.20                    1.00  0.50
 irri.   psw. semi-mech.standard.standard                                  2.30  1.90                    1.00  0.50
 irri.   prw. semi-mech.standard.standard                                  4.90  1.50                    1.00  0.50
 irri.   scwn.semi-mech.standard.standard                            1.70  2.30                    1.00  0.50
 irri.   srwn.semi-mech.standard.standard                                  4.30                          1.00
 irri.   scws.semi-mech.standard.standard                      2.10  2.30  3.50                    1.00  0.50
 irri.   srws.semi-mech.standard.standard                            4.20  2.50                    1.00

 maize.  nwfp.semi-mech.standard.standard                                  5.40                    2.50
 maize.  pcw. semi-mech.standard.standard                                        1.77  1.58                    2.50
 maize.  psw. semi-mech.standard.standard                                        1.68  1.40                    2.50
 maize.  prw. semi-mech.standard.standard                                        1.68  1.40                    2.50
 maize.  scwn.semi-mech.standard.standard                                  2.60  4.30
 maize.  scws.semi-mech.standard.standard                            2.70  4.50

 mus+rap.nwfp.semi-mech.standard.standard                                        1.37  0.75  2.00  1.00  1.00  1.00
 mus+rap.pmw. semi-mech.standard.standard    1.00  1.00  1.00                                      2.47  1.00  1.00
 mus+rap.pcw. semi-mech.standard.standard    1.00  1.00  1.00                                      1.30  1.26  1.00
 mus+rap.psw. semi-mech.standard.standard    1.00  1.00  1.00                                1.25  1.11  1.00  1.00
 mus+rap.prw. semi-mech.standard.standard    1.00  1.00                                      1.82  1.00  1.00  1.00
 mus+rap.scwn.semi-mech.standard.standard    0.50  0.50  0.50                                2.80  1.90  0.50  0.50
 mus+rap.srwn.semi-mech.standard.standard    0.50  0.50  0.50                                2.40  1.80  0.50  0.50
 mus+rap.srws.semi-mech.standard.standard    0.50  0.50  0.50                                2.40  1.80  0.50  0.50
 mus+rap.scws.semi-mech.standard.standard    0.50  0.50  0.50                                2.40  1.70  0.50  0.50

 sc-gur. nwfp.semi-mech.standard.standard    2.50  0.37  1.25                                                  0.98
 sc-gur.(pmw,pcw,psw,prw).semi-mech.
                        standard.standard    2.41  0.62  0.72                                                  0.25
 sc-gur. scwn.semi-mech.standard.standard    1.70  1.30  2.10
 sc-gur. srwn.semi-mech.standard.standard    1.80  1.20  2.10
 sc-gur. scws.semi-mech.standard.standard    1.90  1.30  2.40
 sc-gur. srws.semi-mech.standard.standard    2.00  1.50  2.30
 sc-mill.nwfp.semi-mech.standard.standard   13.00 13.00 10.00  7.50                          6.00  5.70 12.00 13.00
 sc-mill.(pmw,pcw,psw). semi-mech.
                        standard.standard    6.40  6.60  5.22  4.50  4.00                                5.75  5.75
 sc-mill.prw. semi-mech.standard.standard    6.50  7.00  6.22  4.50  3.00                                5.70  5.75
 sc-mill.scwn.semi-mech.standard.standard    1.10  1.00  2.00
 sc-mill.srwn.semi-mech.standard.standard    1.10  1.00  2.00
 sc-mill.scws.semi-mech.standard.standard    0.40  1.50  2.00
 sc-mill.srws.semi-mech.standard.standard    0.50  1.70  1.80
 kha-fod.srws.semi-mech.la-plant.standard                            1.20  0.50  0.50  0.90  0.50  0.50
 kha-fod.nwfp.semi-mech.standard.standard                3.00  0.50  1.00  0.50  1.00  1.00  1.00  0.50
 kha-fod.pmw. semi-mech.standard.standard                3.50  0.50  2.00  1.00  2.00  1.50  1.00  0.50
 kha-fod.pcw. semi-mech.standard.standard                2.50  0.50  1.80  1.00  2.20  1.50  1.00  0.50
 kha-fod.psw. semi-mech.standard.standard                2.00  0.50  1.50  1.00  2.00  1.20  0.80  0.50
 kha-fod.prw. semi-mech.standard.standard                3.00  0.80  1.50  1.20  2.20  1.50  1.00  0.50
 kha-fod.scwn.semi-mech.standard.standard                1.10  0.50  0.90  0.50  0.50  0.90  0.50  0.50
 kha-fod.scws.semi-mech.standard.standard                1.20  0.50  0.90  0.50  0.50  0.90  0.50  0.50
 kha-fod.srwn.semi-mech.standard.standard                1.30  0.50  0.90  0.50  0.50  0.80  0.50  0.50
 kha-fod.srws.semi-mech.standard.standard                1.30  0.50  0.90  0.50  0.50  0.80  0.50  0.50
 wheat.  nwfp.semi-mech.la-plant.heavy                         0.40  0.40                                8.60

 wheat.  pmw. semi-mech.la-plant.heavy                         0.40  0.50                                2.30  3.00
 wheat.  pcw. semi-mech.la-plant.heavy                         0.40  0.40                                4.00  3.60
 wheat.  prw. semi-mech.la-plant.heavy                         0.80  0.80                                      3.10
 wheat.  scwn.semi-mech.la-plant.heavy                                                                   3.10  3.30
 wheat.  srwn.semi-mech.la-plant.heavy                                                                   2.90  3.00
 wheat. (scws,srws).semi-mech.la-plant.
        (heavy,january,light,standard)                                                                   3.20  4.20
 wheat.  nwfp.semi-mech.la-plant.january                       0.40  0.40                                8.60
 wheat.  pmw. semi-mech.la-plant.january                       0.40  0.50                                2.30  3.00
 wheat.  pcw. semi-mech.la-plant.january                       0.40  0.40                                4.00  3.60
 wheat.  prw. semi-mech.la-plant.january                       0.80  0.80                                      3.10
 wheat.  scwn.semi-mech.la-plant.january                                                                 3.10  3.30
 wheat.  srwn.semi-mech.la-plant.january                                                                 2.90  3.00
 wheat.  nwfp.semi-mech.la-plant.light                         0.50  0.60                                8.60
 wheat.  pmw. semi-mech.la-plant.light                         0.60  0.60                                2.30  3.00
 wheat.  pcw. semi-mech.la-plant.light                         0.50  0.60                                4.00  3.60
 wheat.  prw. semi-mech.la-plant.light                         1.10  1.10                                      3.10
 wheat.  scwn.semi-mech.la-plant.light                                                                   3.10  3.30
 wheat.  srwn.semi-mech.la-plant.light                                                                   2.90  3.00
 wheat.  nwfp.semi-mech.la-plant.standard                      0.60  0.70                                8.60
 wheat.  pmw. semi-mech.la-plant.standard                      0.70  0.70                                2.30  3.00
 wheat.  pcw. semi-mech.la-plant.standard                      0.60  0.70                                4.00  3.60
 wheat.  prw. semi-mech.la-plant.standard                      1.30  1.30                                      3.10
 wheat.  scwn.semi-mech.la-plant.standard                                                                3.10  3.30
 wheat.  srwn.semi-mech.la-plant.standard                                                                2.90  3.00

 wheat.  nwfp.semi-mech.qk-harv. heavy                         1.00                                4.10  4.50
 wheat.  pmw. semi-mech.qk-harv. heavy                         1.00                                2.30  3.00
 wheat.  pcw. semi-mech.qk-harv. heavy                         1.00                                      7.60
 wheat.  psw. semi-mech.qk-harv. heavy                         0.90                                3.40  3.00
 wheat.  prw. semi-mech.qk-harv. heavy                         1.00                                      3.10
 wheat.scwn.semi-mech.(qk-harv,standard).
       (heavy,january,light,standard)                                                                    6.40
 wheat.srwn.semi-mech.(qk-harv,standard).
       (heavy,january,light,standard)                                                                    5.9
 wheat. (scws,srws).semi-mech.
        (qk-harv,standard).
        (heavy,january,light,standard)                                                                   7.40
 wheat.  nwfp.semi-mech.qk-harv. january                       1.00                                4.10  4.50
 wheat.  pmw. semi-mech.qk-harv. january                       1.00                                2.30  3.00
 wheat.  pcw. semi-mech.qk-harv. january                       1.00                                      7.60
 wheat.  psw. semi-mech.qk-harv. january                       0.90                                3.40  3.00
 wheat.  prw. semi-mech.qk-harv. january                       1.00                                      3.10
 wheat.  nwfp.semi-mech.qk-harv. light                         1.30                                4.10  4.50
 wheat.  pmw. semi-mech.qk-harv. light                         1.40                                2.30  3.00
 wheat.  pcw. semi-mech.qk-harv. light                         1.30                                      7.60
 wheat.  psw. semi-mech.qk-harv. light                         1.20                                3.40  3.00
 wheat.  prw. semi-mech.qk-harv. light                         1.30                                      3.10
 wheat.  nwfp.semi-mech.qk-harv. standard                      1.50                                4.10  4.50
 wheat.  pmw. semi-mech.qk-harv. standard                      1.60                                2.30  3.00
 wheat.  pcw. semi-mech.qk-harv. standard                      1.50                                      7.60
 wheat.  psw. semi-mech.qk-harv. standard                      1.40                                3.40  3.00
 wheat.  prw. semi-mech.qk-harv. standard                      1.50                                      3.10
 wheat.  nwfp.semi-mech.standard.heavy                         0.50  0.50                          4.10  4.50
 wheat.  pmw. semi-mech.standard.heavy                         0.50  0.50                          2.30  3.00
 wheat.  pcw. semi-mech.standard.heavy                         0.50  0.50                                7.60
 wheat.  psw. semi-mech.standard.heavy                         0.40  0.50                          3.40  3.00
 wheat.  prw. semi-mech.standard.heavy                         0.50  0.50                                3.10
 wheat.  nwfp.semi-mech.standard.january                       0.50  0.50                          4.10  4.50
 wheat.  pmw. semi-mech.standard.january                       0.50  0.50                          2.30  3.00
 wheat.  pcw. semi-mech.standard.january                       0.50  0.50                                7.60
 wheat.  psw. semi-mech.standard.january                       0.40  0.50                          3.40  3.00
 wheat.  prw. semi-mech.standard.january                       0.50  0.50                                3.10
 wheat.  nwfp.semi-mech.standard.light                         0.60  0.70                          4.10  4.50
 wheat.  pmw. semi-mech.standard.light                         0.70  0.70                          2.30  3.00
 wheat.  pcw. semi-mech.standard.light                         0.60  0.70                                7.60
 wheat.  psw. semi-mech.standard.light                         0.60  0.60                          3.40  3.00
 wheat.  prw. semi-mech.standard.light                         0.60  0.70                                3.10
 wheat.  nwfp.semi-mech.standard.standard                      0.70  0.80                          4.10  4.50
 wheat.  pmw. semi-mech.standard.standard                      0.80  0.80                          2.30  3.00
 wheat.  pcw. semi-mech.standard.standard                      0.70  0.80                                7.60
 wheat.  psw. semi-mech.standard.standard                      0.70  0.70                          3.40  3.00
 wheat.  prw. semi-mech.standard.standard                      0.70  0.80                                3.10
 orchard.(nwfp,pmw,pcw,psw,prw).semi-mech.
                        standard.standard      .1    .2    .1    .1    .1    .1    .1    .1    .1    .2    .3
 orchard.(scwn,scws,srwn,srws).semi-mech.
                        standard.standard      .1    .3    .1    .1    .1    .2    .1    .1    .1    .1    .3
 potatoes.(scwn,scws,srwn,srws).
               semi-mech.standard.standard                                                     2.    2.
 potatoes.(nwfp,pmw,pcw,psw,prw).semi-mech.
                         standard.standard     2.                                                                2.
 onions.(nwfp,pmw,pcw,psw,prw).semi-mech.
                         standard.standard     4.
 onions.(scwn,scws,srwn,srws). semi-mech.
                         standard.standard                                                    4.
 chilli.(nwfp,pmw,pcw,psw,prw).semi-mech.
                         standard.standard                 1.
 chilli.(scwn,scws,srwn,srws). semi-mech.
                         standard.standard                                                                       1.
$eject
  table sylds(c,z,t,s,w,ci) straw yield and seed data
                                                     straw-yld   seed
*                                                   (proportion  (kg)
*                                                    of yield)
 basmati.(pmw,pcw,prw).(bullock,semi-mech).
                    standard.standard                 2.33        6.4
 basmati.psw. (bullock,semi-mech).
                    standard. standard                2.12        6.4

 rab-fod.(srwn,srws).(bullock,semi-mech).
                    standard.(light,heavy)            1           6.0
 rab-fod.(nwfp,pmw,pcw,psw,prw).(bullock,semi-mech).
                    standard.standard                 1           2.0
 rab-fod.(scwn,srwn,scws,srws).(bullock,semi-mech).
                    standard. standard                1           6.0

 cotton. pcw. bullock.            el-plant. standard              9.
 cotton. pcw. semi-mech.          la-plant. standard              9.
 cotton. (pmw,pcw,psw,prw,scwn,scws,srwn,srws).
         (bullock,semi-mech).standard. standard                   9.

 gram.   nwfp.(bullock,semi-mech).standard. standard  1.70       10.
 gram.  (pmw,pcw,psw,prw).
              (bullock,semi-mech).standard. standard  1.52       14.
 gram.  (scwn,scws).(bullock,semi-mech).
                standard. standard                    1.50       12.8
 gram.  (srwn,srws).(bullock,semi-mech).
                standard. standard                    2.00       12.0

 irri.   (pmw,pcw,psw,prw).(bullock,semi-mech).
                     standard.standard                1.80        7.1
 irri.   scwn.(bullock,semi-mech).standard. standard  1.34        4.5
 irri.  (srwn,scws,srws).(bullock,semi-mech).
                standard. standard                    1.32        4.5

 maize.  nwfp.(bullock,semi-mech).standard. standard  2.09       14.2
 maize.  pcw. (bullock,semi-mech).standard. standard  3.00        6.1
 maize. (psw,prw). (bullock,semi-mech).
                standard. standard                    2.50        6.1
 maize. (scwn,scws).
              (bullock,semi-mech).standard. standard  2.80        8.1

 mus+rap.nwfp.(bullock,semi-mech).standard. standard  0.50        1.8
 mus+rap.(pmw,pcw,psw,prw,scwn,srwn,srws,scws).
         (bullock,semi-mech).standard. standard       0.65        2.4

 sc-gur. nwfp.(bullock,semi-mech).standard. standard  1.54       3238
 sc-gur. pmw. (bullock,semi-mech).standard. standard  1.81       3238
 sc-gur. pcw. (bullock,semi-mech).standard. standard  1.66       3238
 sc-gur. psw. (bullock,semi-mech).standard. standard  1.54       3238
 sc-gur. prw. (bullock,semi-mech).standard. standard  2.17       3238
 sc-gur. scwn.(bullock,semi-mech).standard. standard  1.81       3238
 sc-gur. srwn.(bullock,semi-mech).standard. standard  3.08       3238
 sc-gur. scws.(bullock,semi-mech).standard. standard  2.17       3238
 sc-gur. srws.(bullock,semi-mech).standard. standard  2.05       3238

 sc-mill.nwfp.(bullock,semi-mech).standard. standard  0.12       3238
 sc-mill.pmw. (bullock,semi-mech).standard. standard  0.14       3238
 sc-mill.pcw. (bullock,semi-mech).standard. standard  0.13       3238
 sc-mill.psw. (bullock,semi-mech).standard. standard  0.10       3238
 sc-mill.prw. (bullock,semi-mech).standard. standard  0.17       3238
 sc-mill.scwn.(bullock,semi-mech).standard. standard  0.14       3238
 sc-mill.srwn.(bullock,semi-mech).standard. standard  0.25       3238
 sc-mill.scws.(bullock,semi-mech).standard. standard  0.17       3238
 sc-mill.srws.(bullock,semi-mech).standard. standard  0.16       3238

 kha-fod.srws.(bullock,semi-mech).la-plant. standard  1           5.7
 kha-fod.nwfp.(bullock,semi-mech).standard. standard  1          20.2
 kha-fod.(pmw,pcw,psw,prw).
         (bullock,semi-mech).standard. standard       1          27.3
 kha-fod.(scwn,scws,srwn,srws).
         (bullock,semi-mech).standard. standard       1           5.7

 wheat.nwfp.(bullock,semi-mech).
            (standard,la-plant,qk-harv).
            (standard,light,heavy,january)            1.3        40.1
 wheat.pmw. (bullock,semi-mech).
            (standard,la-plant,qk-harv).
            (standard,light,heavy,january)            1.3        34.8
 wheat.(pcw,psw). (bullock,semi-mech).
            (standard,la-plant,qk-harv).
            (standard,light,heavy,january)            1.5        34.8
 wheat.prw. (bullock,semi-mech).
            (standard,la-plant,qk-harv).
            (standard,light,heavy,january)            1.6        34.8
 wheat.(scwn,srwn,scws,srws).(bullock,semi-mech).
            (standard,la-plant,qk-harv).
            (standard,light,heavy,january)            1.5        49.8


 orchard.(nwfp,pmw,pcw,psw,prw,scwn,scws,srwn,srws).
            (bullock,semi-mech).standard.standard                   1

 potatoes.(scwn,scws,srwn,srws).
         semi-mech.standard.standard                             1200
 potatoes.(nwfp,pmw,pcw,psw,prw).
         semi-mech.standard.standard                              600

 onions.(nwfp,pmw,pcw,psw,prw,scwn,scws,srwn,srws).
        semi-mech.standard.standard                                3.

 chilli.(nwfp,pmw,pcw,psw,prw,scwn,scws,srwn,srws).
         semi-mech.standard.standard                               5.
$eject
   table fert(p2,c,z) fertilizer applications  (kg per acre)
                   nwfp   pcw   pmw   prw   psw  scwn  scws  srwn  srws

nitrogen.basmati         26.6  26.6  21.9  23.4
nitrogen.irri            39.4  39.4  23.3  26.8  68.6  61.5  48.9  41.7
nitrogen.cotton    26.7  42.3  30.0  30.0  19.7  55.0  54.9  39.6  39.6
nitrogen.maize     27.0  27.1  27.1  23.6  19.5  42.0  42.0  42.0  42.0
nitrogen.kha-fod   21.0  25.3  18.1  19.2  18.4  49.0  49.0  49.0  49.0
nitrogen.wheat     46.8  40.9  36.5  32.2  33.3  54.9  53.5  29.5  39.8
nitrogen.rab-fod   10.0  25.3  18.1  19.2  18.4  28.0  28.0  28.0  28.0
nitrogen.sc-mill   83.4  44.8  63.2  33.9  33.9  65.1  65.1  65.1  65.1
nitrogen.sc-gur    24.0  19.0  19.0  19.0  19.0  28.0  28.0  28.0  28.0
nitrogen.onions    60.6  48.0  48.0  48.0  48.0  70.7  70.7  70.7  70.7
nitrogen.potatoes  48.6  38.5  38.5  38.5  38.5  56.7  56.7  56.7  56.7
nitrogen.mus+rap   33.6  30.8  30.8  30.8  30.8  45.4  45.4  45.4  45.4
nitrogen.chilli    48.6  38.5  38.5  38.5  38.5  56.7  56.7  56.7  56.7
nitrogen.orchard   60.0  47.5  47.5  47.5  47.5  70.0  70.0  70.0  70.0

phosphate.basmati        13.3  13.3  11.4  10.6   0.0   0.0   0.0   0.0
phosphate.irri            8.0   8.0   6.8   6.4  15.7  16.6  13.5  12.5
phosphate.cotton   26.2  11.3   5.2   7.0   7.0  12.8  14.4  13.0  13.0
phosphate.maize    17.6   8.8   6.9   8.6   8.6   7.7   7.7   7.7   7.7
phosphate.kha-fod  25.9  11.4  11.4  11.4  11.4  15.3  15.3  15.3  15.3
phosphate.wheat    12.3  12.8   8.8  11.4  10.9  15.1  15.1  11.8  17.3
phosphate.rab-fod   5.4   7.6   5.4   5.8   5.6   7.7   7.7   7.7   7.7
phosphate.sc-mill  46.8  14.0  19.8  10.6  10.6  18.5  18.5  18.5  18.5
phosphate.sc-gur   10.4   5.7   5.7   5.7   5.7   7.7   7.7   7.7   7.7
phosphate.onions   26.2   7.0   7.0   7.0   7.0   9.3   9.3   9.3   9.3
phosphate.potatoes 42.0  18.5  18.5  18.5  18.5  24.9  24.9  24.9  24.9
phosphate.mus+rap  25.9  13.3  13.3  13.3  13.3  17.9  17.9  17.9  17.9
phosphate.chilli   26.2  11.6  11.6  11.6  11.6  15.5  15.5  15.5  15.5
phosphate.orchard  13.0   5.7   5.7   5.7   5.7   7.7   7.7   7.7   7.7

 parameter fertgr(c)  fertilizer application growth rate  percent

   /(basmati, irri,   cotton,  rab-fod
     gram   , maize,  mus+rap, kha-fod
     sc-gur , sc-mill,wheat  , orchard)     3,
   (potatoes, onions,chilli)                2.4/
     ;

*&&z   zone4xxxxx(z,c,p2) = fert(p2,c,z) ;


fert(p2,c,z) = fert(p2,c,z)*sum(is$isr(is),
                      (1+fertgr(c)/100)**(ord(is)+1979-baseyear)) ;

$stitle  crop yields
 parameter
  natyield(c) national crop yields 1988 for standard technologies (kgs)
      /  basmati      457
         irri         880
         cotton       695
         rab-fod    15000
         gram         183
         maize        534
         mus+rap      307
         kha-fod    10000
         sc-gur      1270
         sc-mill    15870
         wheat        780
         orchard     3400
         potatoes    3980
         onions      4615
         chilli       567 /


 table yldprpv(c,pv) province yields proportion of national 1987-88
*                    note:  used 3-year average from latest asp

                nwfp   punjab     sind

    wheat      0.926    0.978    1.103
    basmati        0        1        0
    irri           0    0.859    1.001
    cotton         0    1.153    0.861
    sc-mill    1.090    0.919    1.190
    sc-gur     1.090    0.919    1.190
    maize      1.006    1.029    0.408
    mus+rap    0.622    1.192    0.842
    gram       0.911    0.937    1.465
    rab-fod     .400    1.000    1.000
    kha-fod     .770    1.000    1.300
    orchard     .936    1.070        1
    potatoes    .973     .960     .873
    onions     1.172    1.085    0.843
    chilli     0.786    1.167    0.952


  table yldprzs(c,z) zones yields as proportion of province-standard technologies

          nwfp   pcw    pmw    prw    psw   scwn   scws   srwn   srws
  wheat      1  0.96   0.84   1.06   1.13   0.96   1.07   0.82   1.30
  basmati    1   .83   0.00   1.26   0.79   0.00   0.00   0.00   0.00
  irri       1  0.83   0.00   1.26   0.00   1.06   0.91   1.10   0.90
  cotton     1  1.10   0.71   0.70   0.80   1.09   0.96   0.78   0.78
  sc-mill    1  1.04   0.93   0.93   0.98   0.95   1.10   0.64   0.87
  sc-gur     1  1.04   0.93   0.93   0.98   0.95   1.10   0.64   0.87
  maize      1  1.03   0.76   1.08   0.97   1.07   1.06   0.87   1.00
  mus+rap    1  1.03   0.76   1.08   0.97   1.07   1.06   0.87   1.00
  gram       1  1.03   0.76   1.08   0.97   1.07   1.06   0.87   1.00
  rab-fod    1  1.23   0.81   1.17   1.23   1.     1.     1.     1.
  kha-fod    1   .83   0.54   0.82   0.82   1.     1.     1.     1.00
  orchard    1  1.     1.     1.     1.0    1.     1.     1.     1.
  potatoes   1  1.03   0.76   1.08   0.97   1.07   1.06   0.87   1.00
  onions     1  1.03   0.76   1.08   0.97   1.07   1.06   0.87   1.00
  chilli     1  1.03   0.76   1.08   0.97   1.07   1.06   0.87   1.00

  parameter yldprzo(c,s,w) yields as proportion of standard technologies
     /
    rab-fod.standard.heavy        0.70
    rab-fod.standard.light        0.80
    kha-fod.la-plant.standard     1.05

    wheat.   qk-harv .standard                       1.00
    wheat.   la-plant.(heavy,january)                0.57
    wheat.   la-plant.light                          0.75
    wheat.   la-plant.standard                       0.88
    wheat.  (qk-harv,standard).(heavy,january)       0.65
    wheat.  (qk-harv,standard).light                 0.85
      /
 table growthcy(c,z) growth rate of crop yields from 1988 base (percent)

            (nwfp,pmw,pcw,psw,prw,scwn,scws,srwn,srws)
   cotton               5.00
   maize                0.73
   sc-gur               2.29
   sc-mill              2.29
   wheat                0.41
   onions               1.66

   table weedy(z,sea,c) weed yields by crop (tonns per acer)
                                   basmati   irri  cotton  gram  mus+rap  maize  sc-gur  sc-mill  wheat orchard
 (nwfp,pmw,pcw,psw,prw).rabi                                .3      .3            .3      .3       .4     .3
 (scwn,scws,srwn,srws). rabi                                .3      .3            .3      .3       .4     .3
 (nwfp,pmw,pcw,psw,prw).kharif       1        1      1                      1      1       1              1
 (scwn,scws,srwn,srws). kharif                1      1                      1      1       1              1

   table graz(z,sea)  grazing from slack land (tonns per acre )
                                                 rabi     kharif
  prw                                             .15       .5
  pmw                                             .1        .2
 (psw,pcw)                                        .2        .3
 (scwn,scws,srwn,srws)                            .2        .5

 parameter
    yield(c,t,s,w,z) yield by zone crop technology in metric tonns
    growthcyf(c,z)   growth factor for crop yields using growthcy ;

*  growthcy(c,z)$(growthcy(c,z) gt 3) = 3.0 ;
   yield(c,t,"standard","standard",z) =          sum(pv$pvz(pv,z),
                         natyield(c)/1000 *yldprpv(c,pv) *yldprzs(c,z) );

    yield(c,t,s,w,z)$yldprzo(c,s,w) =
                  yield(c,t,"standard","standard",z)*yldprzo(c,s,w) ;

  growthcyf(c,z) =  sum(is$isr(is),
                      (1+growthcy(c,z)/100)**(ord(is)+1979-baseyear));
  yield(c,t,s,w,z) = yield(c,t,s,w,z) * growthcyf(c,z) ;
display baseyear, growthcyf, fert ;
* display yield;
$eject
$ontext
*  report on input data
  set ic crop inputs /land, labor, bullock, water, tractor/
  parameter rep1,rep2 ;
  rep1(z,c,t,s,w,"land",m)    = land(c,z,t,s,w,m)    ;
  rep1(z,c,t,s,w,"bullock",m) = bullock(c,z,t,s,w,m) ;
  rep1(z,c,t,s,w,"labor",m)   = labor(c,z,t,s,w,m)   ;
  rep1(z,c,t,s,w,"water",m)   = water(c,z,t,s,w,m)   ;
  rep1(z,c,t,s,w,"tractor",m) = tractor(c,z,t,s,w,m) ;

  rep1(z,c,t,s,w,ic,"total")  = sum(m, rep1(z,c,t,s,w,ic,m)) ;
  rep2(z,c,t,s,w,ic)          = rep1(z,c,t,s,w,ic,"total")   ;
  rep2(z,c,t,s,w,ci)          = sylds(c,z,t,s,w,ci)          ;
  rep2(z,c,t,s,w,p2)$tech(z,c,t,s,w) = fert(p2,c,z)          ;
  rep2(z,c,t,s,w,"yield")     = yield(c,t,s,w,z)             ;

  options rep1:2:4:1, rep2:2:4:1 ; display rep1, rep2 ;
$offtext
$stitle livestock data
     table iolive(a,z,*) livestock input output coefficients by zones
*     meat and milk yeilds are of 1985.
                                tdn          dp           labor     cow-milk   buff-milk   meat     fix-cost
*                             (metric tons/ season)      man hrs            liters          kgs
*                                                        per month  (- - -  per  year    - - - )
 cow. nwfp                        .68       .061          25.1        359                  13.4       15
 cow. pmw                         .68       .061          25.1        284                  13.4       15
 cow. pcw                         .68       .061          25.1        303                  13.4       15
 cow. psw                         .68       .061          25.1        268                  13.4       15
 cow. prw                         .68       .061          25.1        332                  13.4       15
 cow.scwn                         .68       .046          23.5        460                  10.5       17
 cow.srwn                         .68       .083          23.5        383                  15.4       17
 cow.scws                         .68       .046          23.5        350                  10.5       17
 cow.srws                         .68       .083          23.5        468                  15.4       17

 bullock.(nwfp,pmw,pcw,psw,prw)   .635      .058          15.1                             14.3       50
 bullock.(scwn,srwn,scws,srws)    .635      .058          14.8                             14.1       50

 buffalo.nwfp                    1.04       .095          33.6                    602      17.3       50
 buffalo.pmw                     1.04       .095          33.6                    653      17.3       50
 buffalo.pcw                     1.04       .095          33.6                    652      17.3       50
 buffalo.psw                     1.04       .095          33.6                    643      17.3       50
 buffalo.prw                     1.04       .095          33.6                    590      17.3       50
 buffalo.scwn                    1.04       .073          31.4                    946      16.2       50
 buffalo.srwn                    1.04       .128          31.4                    750      22.5       50
 buffalo.scws                    1.04       .073          31.4                    773      16.2       50
 buffalo.srws                    1.04       .128          31.4                   1018      22.5       50

     table sconv(nt,sea,c)   tdn and dp conversion factor from crop straw

                      basmati  cotton  rab-fod  gram   irri   kha-fod  maize  mus+rap  sc-mill sc-gur  wheat
 tdn.kharif                                       .6             .14             .5                      .6
 dp. kharif                                       .1            .019           .005                    .005
 tdn.rabi                .5               .14            .5               .5               .17    .17
 dp. rabi              .005              .025          .005               .1              .005   .005
     ;
  scalars      repco       reproductive coefficient                            / 2.5  /
               gr          required proportion of green fodder in total fodder / 0.3  /
               growthq     growth rate of milk and meat yields (percent)       / 2.5  /
  parameters   bp(m)       draft power available per bullock(hours per month)
                ;
           bp(m)=96;                  bp("may")=77;             bp("jun")=77;
*&&z    zone2xxxxx(z,a,set1) = iolive(a,z,set1) ;
   iolive(a,z,q) = iolive(a,z,q) * sum(is$isr(is), (1+growthq/100)**(ord(is)+1979 -baseyear) ) ;
    options iolive:3;   display  iolive ;
$stitle canal and agroclimatic zone data
 set cnl  irrigation canals in the indus river irrigation system/
    01-ud            upper dipalpur
    02-cbd           central bari doab canal
    03-ray           raya canal
    04-uc            upper chenab canal
    05-mr            marala ravi canal
    06-sad           sadiqia canal
    07-for           fordwah canal
    08-pak           upper pakpattan+u-bahawal+qaim+u-mailsi canal
    09-ld            lower dipalpur canal
    10-lbd           lower bari doab canal
    11-jha           jhang canal (lcc)
    12-gug           gugera branch canal (lcc)
    13-uj            upper jehlum canal
    14-lj            lower jehlum canal
    15-bah           bahawal canal
    16-mai           lower mailsi+ lower pakpattan canal
    17-sid           sidhnai canal
    18-hav           haveli canal
    19-ran           rangpur canal
    20-pan           panjnad canal
    21-abb           abbasia canal
    22-usw           upper swat canal
    23-lsw           lower swat canal
    24-war           warsak canal
    25-kab           kabul river canal
    26-tha           thal canal
    27-pah           paharpur canal
    28-muz           muzffgarh canal
    29-dgk           dera ghazi khan canal
    31-p+d           pat plus desert canal
    32-beg           begari canal
    33-gho           ghotki canal
    34-nw            north west canal
    35-ric           rice canal
    36-dad           dadu canal
    37-kw            khairpur west canal
    38-ke            khairpur east canal
    39-roh           rohri canal
    41-nar           nara canal
    42-kal           kalri canal
    43-lch           lined channel
    44-ful           fuleli canal
    45-pin           pinyari canal/
 pvcnl(pv,cnl)   province to canals map /
               nwfp.  (22-usw, 23-lsw, 24-war, 25-kab)
               punjab.(01-ud,  02-cbd, 03-ray, 04-uc,  05-mr,  06-sad, 07-for,
                       08-pak, 09-ld,  10-lbd, 11-jha, 12-gug, 13-uj,  14-lj,
                       15-bah, 16-mai, 17-sid, 18-hav, 19-ran, 20-pan, 21-abb,
                       26-tha, 27-pah, 28-muz, 29-dgk)
*       sind canals including baluchistan
               sind.  (31-p+d, 32-beg, 33-gho, 34-nw,  35-ric, 36-dad
                       37-kw,  38-ke,  39-roh, 41-nar, 42-kal, 43-lch
                       44-ful, 45-pin ) /
    gwfg(cnl,sa,g)  subarea identification by the groundwater quality
         ;
$stitle canal command data
  table  comdef(is,dc,cnl) canal command characteristics
*   cca       --   culturable commanded area of the canal( millions of acres)
*   ccap      --   canal capacity at the canal head(millions of acre feet)
*   ceff      --   canal efficiency from barrage to the water course head
*   wce-r     --   water course command efficinecy in rabi season
*   wce-k     --   water course command efficinecy in kharif season
*   flde      --   field efficiency
*
*   note:  ccap for 05-mr, 24-war, 25-kab are set equal to post tarbela average diversion. original capacities were
*                    .166   .044    .047

                  01-ud   02-cbd   03-ray    04-uc    05-mr   06-sad   07-for   08-pak    09-ld   10-lbd    11-jha

 1980.cca          .36      .649     .424    1.017     .158     .969     .426    1.049     .615    1.67      1.168
 1980.ccap         .131     .163     .11      .506     .282     .332     .204     .746     .244     .518      .41
 1980.ceff         .70      .80      .80      .76      .80      .72      .70      .61      .80      .72       .70
 1980.wce-r        .57      .57      .45      .57      .57      .51      .51      .57      .57      .57       .55
 1980.wce-k        .57      .57      .45      .57      .57      .51      .51      .57      .57      .57       .55
 1980.flde         .90      .90      .90      .90      .90      .80      .80      .90      .90      .90       .90

 1988.cca          .36      .649     .424    1.017     .158     .969     .426    1.049     .615    1.67      1.168
 1988.ccap         .131     .163     .11      .506     .282     .332     .204     .746     .244     .518      .41
 1988.ceff         .703     .811     .800     .760     .800     .728     .704     .621     .810     .722      .706
 1988.wce-r        .57      .573     .45      .57      .57      .522     .51      .575     .57      .57       .554
 1988.wce-k        .57      .573     .45      .57      .57      .522     .51      .575     .57      .57       .554
 1988.flde         .90      .90      .90      .90      .90      .80      .80      .90      .90      .90       .90

 1993.cca          .36      .649     .424    1.017     .158     .969     .426    1.049     .615    1.67      1.168
 1993.ccap         .131     .163     .11      .506     .282     .332     .204     .746     .244     .518      .41
 1993.ceff         .705     .815     .808     .765     .800     .733     .708     .625     .810     .724      .713
 1993.wce-r        .58      .585     .46      .58      .58      .543     .53      .59      .58      .58       .564
 1993.wce-k        .58      .585     .46      .58      .58      .543     .53      .59      .58      .58       .564
 1993.flde         .90      .90      .90      .90      .90      .80      .80      .90      .90      .90       .90

 2000.cca          .36      .649     .424    1.017     .158     .969     .426    1.049     .615    1.67      1.168
 2000.ccap         .131     .163     .11      .506     .282     .332     .204     .746     .244     .518      .41
 2000.ceff         .705     .815     .808     .765     .800     .733     .708     .625     .810     .724      .713
 2000.wce-r        .59      .595     .47      .59      .59      .553     .54      .60      .59      .59       .574
 2000.wce-k        .59      .595     .47      .59      .59      .553     .54      .60      .59      .59       .574
 2000.flde         .90      .90      .90      .90      .90      .80      .80      .90      .90      .90       .90

          +      12-gug    13-uj    14-lj   15-bah   16-mai   17-sid   18-hav   19-ran   20-pan   21-abb    22-usw

 1980.cca         1.866     .544    1.50      .605     .996     .869     .179     .344    1.348     .154      .279
 1980.ccap         .433     .345     .464     .346     .362     .268     .068     .176     .671     .08       .117
 1980.ceff         .74      .80      .64      .77      .77      .72      .77      .70      .72      .64       .75
 1980.wce-r        .55      .51      .51      .57      .57      .57      .57      .46      .57      .57       .52
 1980.wce-k        .55      .51      .51      .57      .57      .57      .57      .46      .57      .57       .52
 1980.flde         .90      .90      .90      .90      .90      .90      .85      .80      .90      .90       .90

 1988.cca         1.866     .544    1.50      .605     .996     .869     .179     .344    1.348     .154      .279
 1988.ccap         .433     .345     .464     .346     .362     .268     .068     .176     .773     .08       .117
 1988.ceff         .749     .800     .642     .782     .780     .721     .773     .700     .728     .640      .750
 1988.wce-r        .553     .51      .51      .57      .57      .57      .57      .46      .59      .59       .53
 1988.wce-k        .553     .51      .51      .57      .57      .57      .57      .46      .59      .59       .53
 1988.flde         .90      .90      .90      .90      .90      .90      .85      .80      .90      .90       .90

 1993.cca         1.866     .544    1.500     .605     .996     .869     .179     .344    1.348     .154      .279
 1993.ccap         .433     .345     .464     .346     .362     .268     .068     .176     .773     .08       .132
 1993.ceff         .751     .808     .644     .783     .783     .727     .776     .706     .73      .641      .757
 1993.wce-r        .565     .52      .52      .58      .58      .58      .58      .48      .59      .59       .53
 1993.wce-k        .565     .52      .52      .58      .58      .58      .58      .48      .59      .59       .53
 1993.flde         .90      .90      .90      .90      .90      .90      .85      .80      .90      .90       .90

 2000.cca         1.866     .544    1.573     .605     .996     .869     .179     .344    1.348     .154      .326
 2000.ccap         .433     .345     .505     .346     .362     .268     .068     .176     .773     .08       .162
 2000.ceff         .751     .808     .644     .783     .783     .727     .776     .706     .73      .641      .757
 2000.wce-r        .575     .53      .53      .59      .59      .59      .59      .49      .60      .60       .53
 2000.wce-k        .575     .53      .53      .59      .59      .59      .59      .49      .60      .60       .53
 2000.flde         .90      .90      .90      .90      .90      .90      .85      .80      .90      .90       .90

          +      23-lsw   24-war   25-kab   26-tha   27-pah   28-muz   29-dgk   31-p+d   32-beg   33-gho     34-nw

 1980.cca          .182     .119     .048    1.641     .104     .809     .909    1.075    1.002     .858     1.215
 1980.ccap         .06      .049     .047     .577     .037     .495     .529     .799    1.155     .648      .566
 1980.ceff         .75      .80      .72      .65      .73      .70      .70      .83      .82      .76       .80
 1980.wce-r        .52      .52      .52      .48      .51      .46      .54      .55      .55      .45       .55
 1980.wce-k        .52      .52      .52      .48      .51      .46      .54      .60      .65      .45       .60
 1980.flde         .90      .90      .90      .80      .85      .80      .85      .85      .85      .85       .85

 1988.cca          .182     .119     .048    1.641     .570     .809     .909    1.075    1.002     .858     1.215
 1988.ccap         .120     .049     .047     .577     .330     .495     .529     .799    1.155     .648      .566
 1988.ceff         .750     .800     .720     .663     .730     .720     .703     .830     .820     .763      .800
 1988.wce-r        .57      .57      .53      .48      .51      .46      .54      .55      .55      .45       .55
 1988.wce-k        .57      .57      .53      .48      .51      .46      .54      .60      .65      .45       .60
 1988.flde         .90      .90      .90      .80      .85      .80      .85      .85      .85      .85       .85

 1993.cca          .182     .119     .048    1.641     .570     .809     .909    1.075    1.002     .858     1.215
 1993.ccap         .120     .049     .047     .577     .330     .495     .529     .799    1.155     .648      .566
 1993.ceff         .751     .80      .723     .671     .73      .725     .707     .831     .82      .765      .80
 1993.wce-r        .57      .60      .53      .49      .56      .47      .55      .56      .56      .47       .56
 1993.wce-k        .57      .60      .53      .49      .56      .47      .55      .61      .66      .47       .61
 1993.flde         .90      .90      .90      .80      .85      .80      .85      .85      .85      .85       .85

 2000.cca          .182     .119     .048    1.641     .570     .809    1.232    1.215    1.002     .858     1.215
 2000.ccap         .120     .049     .047     .577     .330     .495     .731     .900    1.155     .648      .566
 2000.ceff         .751     .80      .723     .671     .73      .725     .707     .831     .82      .765      .80
 2000.wce-r        .58      .61      .54      .50      .57      .48      .56      .594     .57      .48       .57
 2000.wce-k        .58      .61      .54      .50      .57      .48      .56      .65      .67      .48       .62
 2000.flde         .90      .90      .90      .80      .85      .80      .85      .85      .85      .85       .85

          +      35-ric   36-dad    37-kw    38-ke   39-roh   41-nar   42-kal   43-lch   44-ful   45-pin

 1980.cca          .519     .584     .417     .373    2.561    2.176     .592     .502     .923     .758
 1980.ccap         .829     .319     .157     .207     .981     .873     .546     .205     .894     .831
 1980.ceff         .85      .80      .74      .73      .80      .80      .80      .80      .80      .82
 1980.wce-r        .55      .55      .45      .45      .45      .45      .55      .55      .55      .55
 1980.wce-k        .75      .60      .45      .45      .45      .45      .60      .61      .65      .62
 1980.flde         .85      .85      .85      .85      .85      .85      .85      .85      .85      .85

 1988.cca          .519     .584     .417     .373    2.561    2.176     .592     .502     .923     .758
 1988.ccap         .829     .319     .157     .214     .981     .873     .546     .205     .894     .831
 1988.ceff         .850     .800     .748     .736     .807     .811     .800     .800     .800     .820
 1988.wce-r        .55      .55      .45      .455     .456     .45      .55      .55      .55      .55
 1988.wce-k        .75      .60      .45      .455     .456     .45      .60      .61      .65      .62
 1988.flde         .85      .85      .85      .85      .85      .85      .85      .85      .85      .85

 1993.cca          .519     .584     .417     .373    2.561    2.176     .592     .502     .923     .758
 1993.ccap         .829     .319     .157     .214     .981    1.041     .546     .205     .894     .831
 1993.ceff         .854     .804     .75      .738     .81      .816     .80      .80      .80      .821
 1993.wce-r        .56      .56      .46      .465     .47      .469     .56      .56      .56      .56
 1993.wce-k        .74      .61      .46      .465     .47      .469     .61      .62      .65      .63
 1993.flde         .85      .85      .85      .85      .85      .85      .85      .85      .85      .85

 2000.cca          .519     .584     .417     .373    2.561    2.176     .592     .502     .923     .758
 2000.ccap         .829     .319     .157     .214     .981    1.041     .546     .205     .894     .831
 2000.ceff         .854     .804     .75      .738     .811     .816     .80      .80      .80      .821
 2000.wce-r        .57      .57      .47      .475     .49      .479     .57      .57      .57      .57
 2000.wce-k        .75      .62      .47      .475     .49      .479     .61      .62      .65      .63
 2000.flde         .85      .85      .85      .85      .85      .85      .85      .85      .85      .85

* 1988 parameters are derived considering the following projects;
* command water management project  50% completion
* scarp vi, scarp mardan, khairpur tile, fourth drainage, ofwmi and ii
* irrigation system rehablitation phase i, chasma right bank.
*


      table subdef(sa,cnl) sub-area definition(proportion of cca)  by canals

                  01-ud   02-cbd   03-ray    04-uc    05-mr   06-sad   07-for   08-pak    09-ld   10-lbd    11-jha
  s1                1       .5       1         1        1       1        1        1         1       .5        .32
  s2                        .5                                                                      .5        .49
  s3                                                                                                          .19

          +      12-gug    13-uj    14-lj   15-bah   16-mai   17-sid   18-hav   19-ran   20-pan   21-abb    22-usw
  s1               .26       1        .64     .8       .65      1        1        1        .7       1         1
  s2               .53                .36     .2       .35                                 .3
  s3               .21

          +      23-lsw   24-war   25-kab   26-tha   27-pah   28-muz   29-dgk   31-p+d   32-beg   33-gho     34-nw
  s1               1        1        1        .35      1        .25      1        1        .5       .5         1
  s2                                          .17               .75                        .5       .5
  s3                                          .30
  s4                                          .18
          +      35-ric   36-dad    37-kw    38-ke   39-roh   41-nar   42-kal   43-lch   44-ful   45-pin
  s1               1        1         1        1       .39      .2       1        1        1        1
  s2                                                   .20      .8
  s3                                                   .16
  s4                                                   .25
                  ;
$eject
  set    zsa(z,cnl,sa)   canal-subarea to agroclimatic zone mapping /

          nwfp. (22-usw.  s1,             23-lsw.  s1,      24-war.  s1,      25-kab.  s1)
          pmw.  (26-tha. (s1,s2,s3,s4),   27-pah.  s1,      28-muz.  s1)
          pcw.  (01-ud.   s1,             02-cbd.  s2,      06-sad.  s1,      07-for.  s1,      08-pak. s1
                 09-ld.   s1,             10-lbd. (s1,s2),  15-bah. (s1,s2),  16-mai. (s1,s2),  17-sid. s1
                 19-ran.  s1,             20-pan. (s1,s2),  21-abb.  s1,      28-muz.  s2,      29-dgk. s1)
          psw.  (11-jha. (s2,s3),         12-gug. (s2,s3),  13-uj.   s1,      14-lj.  (s1,s2),  18-hav. s1)
          prw.  (02-cbd.  s1,             03-ray. s1,       04-uc.   s1,      05-mr.   s1,      11-jha. s1, 12-gug.s1)

          scwn. (33-gho. (s1,s2),         37-kw. s1,        38-ke. s1,        39-roh. (s1,s2),  41-nar. s1)
          srwn. (31-p+d. s1,              32-beg. (s1,s2),  34-nw. s1,        35-ric.  s1,      36-dad. s1)
          scws. (39-roh. (s3,s4),         41-nar.  s2   )
          srws. (42-kal. s1,              43-lch. s1,       44-ful. s1,       45-pin. s1 )  /

    gwf(cnl,sa)    subareas with fresh ground water /
          01-ud.  s1,    02-cbd. (s1,s2),03-ray. s1,  04-uc.  s1,     05-mr.   s1
          07-for. s1,    08-pak.  s1,    09-ld.  s1,  10-lbd.(s1,s2), 11-jha. (s1,s2)
          12-gug.(s1,s2),13-uj.   s1,    14-lj.  s1,  15-bah. s1,     16-mai.  s1
          17-sid. s1,    19-ran.  s1,    20-pan. s1,  21-abb. s1,     22-usw.  s1
          23-lsw. s1,    24-war.  s1,    25-kab. s1,  26-tha.(s1,s3), 27-pah.  s1
          28-muz. s1,    29-dgk.  s1,    31-p+d. s1,  32-beg. s1,     33-gho.  s1
          37-kw.  s1,    39-roh. (s1,s3) /
    ;


    gwfg(cnl,sa,"saline")$subdef(sa,cnl)  = yes;
    gwfg(cnl,sa,"saline")$gwf(cnl,sa)     = no ;   gwfg(cnl,sa,"fresh")$gwf(cnl,sa) = yes ;
  parameter  carea(cnl,*)     cca classified by groundwater quality for each canal
               ;
  loop(isr,
  carea(cnl,g)        = sum(sa$gwfg(cnl,sa,g), subdef(sa,cnl)*comdef(isr,"cca",cnl)  ) );
  carea(cnl,"total")  = sum(g,  carea(cnl,g) );
 display carea  ;
$ontext
*- report and check on acz definition

 set zsa1(cnl,sa,z);
 parameter rep3
           rep4
             ;
  rep3(is,cnl,sa)       = subdef(sa,cnl)*comdef(is,"cca",cnl) ;
  rep3(is,cnl,g)        = sum(sa$gwfg(cnl,sa,g), subdef(sa,cnl)*comdef(is,"cca",cnl)  ) ;
  rep3(is,cnl,"total")  = sum(g, rep3(is,cnl,g) )     ;
  rep3(is,"total",t1)   = sum(cnl, rep3(is,cnl,t1))   ;
  rep3(is,cnl,"cca")    = comdef(is,"cca",cnl)        ;
  rep3(is,"total","cca")= sum(cnl, comdef(is,"cca",cnl)) ;

  rep4(is,z,g)       = sum((cnl,sa)$(zsa(z,cnl,sa)$gwfg(cnl,sa,g)), comdef(is,"cca",cnl)*subdef(sa,cnl) );
  rep4(is,z,"total") = sum((cnl,sa)$zsa(z,cnl,sa), comdef(is,"cca",cnl)*subdef(sa,cnl) );
  rep4(is,"total",t1)= sum(z, rep4(is,z,t1)    )   ;
  zsa1(cnl,sa,z)$zsa(z,cnl,sa) = yes ;
  display zsa1, gwf, gwfg, rep3, rep4 ;
$offtext
*-
$stitle climatic data

   table evap(cnl,m) pan evaporation(feet)
                   jan      feb      mar      apr      may     jun     jul     aug     sep     oct     nov     dec
 01-ud             .244     .314     .545     .786    .982   1.040    .913    .814    .729    .571    .357    .244
 02-cbd            .200     .314     .529     .771    .968    .998    .829    .743    .686    .543    .329    .229
 03-ray            .200     .314     .514     .771    .951    .968    .812    .743    .671    .543    .326    .214
 04-uc             .200     .283     .514     .771    .940    .998    .829    .743    .671    .543    .326    .214
 05-mr             .200     .286     .529     .769    .942    .998    .829    .743    .671    .543    .326    .214
 06-sad            .286     .414     .648     .814   1.139   1.254   1.012    .940    .743    .629    .402    .286
 07-for            .271     .400     .614     .786   1.083   1.196    .968    .893    .714    .564    .386    .286
 08-pak            .271     .343     .588     .829   1.012   1.056    .970    .885    .771    .614    .400    .288
 09-ld             .243     .312     .586     .786    .969   1.056    .955    .871    .771    .586    .371    .286
 10-lbd            .243     .314     .557     .786    .982   1.040    .925    .938    .745    .581    .360    .243
 11-jha            .229     .300     .514     .729    .990   1.012    .870    .757    .700    .486    .343    .217
 12-gug            .214     .312     .526     .769    .965   1.008    .826    .740    .683    .555    .326    .226
 13-uj             .229     .300     .514     .729    .967   1.012    .870    .760    .700    .571    .343    .214
 14-lj             .257     .314     .557     .788    .982   1.026    .930    .843    .743    .600    .371    .260
 15-bah            .271     .360     .657     .860   1.069   1.098   1.011    .899    .771    .657    .400    .271
 16-mai            .271     .340     .629     .842   1.042   1.056    .940    .869    .757    .643    .400    .271
 17-sid            .257     .343     .588     .842   1.026   1.026    .940    .860    .757    .643    .414    .271
 18-hav            .271     .329     .571     .800    .982   1.025    .940    .857    .757    .600    .386    .271
 19-ran            .244     .314     .557     .786    .982   1.040    .913    .814    .729    .571    .357    .244
 20-pan            .200     .314     .529     .771    .968   1.112   1.026    .913    .800    .543    .329    .457
 21-abb            .271     .338     .671     .858   1.085   1.098   1.026    .913    .788    .671    .400    .288
 (22-usw,23-lsw,
  24-war,25-kab)   .214     .257     .429     .643    .913   1.025    .969    .870    .657    .529    .329    .200
 26-tha            .260     .331     .517     .760   1.057   1.143   1.057    .899    .721    .631    .402    .245
 27-pah            .257     .329     .502     .729   1.094   1.225   1.154    .943    .829    .671    .357    .243
 28-muz            .271     .340     .600     .870   1.056   1.039    .968    .885    .771    .657    .414    .271
 29-dgk            .286     .357     .657     .885   1.082   1.098   1.011    .913    .798    .671    .414    .271
 31-p+d            .371     .457     .758     .956   1.142   1.126   1.026    .942    .830    .715    .486    .343
 (32-beg,34-nw,
 35-ric)           .386     .474     .700     .870   1.068    .949    .943    .949    .757    .657    .457    .343
 33-gho            .429     .429     .700     .871   1.040   1.098    .998    .900    .800    .629    .457    .329
 36-dad            .357     .429     .700     .956   1.169   1.185   1.027    .940    .786    .700    .429    .357
 37-kw             .386     .414     .674     .900   1.083   1.112    .998    .899    .786    .643    .429    .343
 38-ke             .429     .443     .745     .926   1.152   1.154   1.025    .955    .829    .700    .486    .371
 39-roh            .371     .414     .700     .940   1.154   1.139    .982    .926    .757    .698    .457    .371
 41-nar            .457     .500     .843    1.054   1.212   1.155    .999    .914    .829    .786    .557    .443
 42-kal            .529     .557     .857     .982   1.094    .970    .814    .757    .786    .800    .614    .529
 (43-lch,44-ful)   .471     .529     .870    1.095   1.210   1.111    .985    .886    .843    .814    .586    .471
 45-pin            .557     .586     .899    1.044   1.140   1.012    .843    .800    .829    .843    .643    .557

   table rain(cnl,m) rain (inches)
                       jan      feb     mar     apr     may     jun      jul     aug     sep     oct      nov     dec
 01-ud                .770     .520    .660    .470    .450    1.180    4.480   4.070   1.530    .160    .220    .360
 02-cbd               .878     .700    .712    .510    .506    1.364    4.868   4.502   1.798    .192    .172    .404
 03-ray              1.180    1.120   1.040    .740    .710    1.820    5.610   6.450   2.460    .230    .140    .490
 04-uc               1.141    1.089   1.013    .724    .693    1.779    5.444   6.278   2.396    .224    .138    .478
 05-mr               1.180    1.120   1.040    .740    .710    1.820    5.610   6.450   2.460    .230    .140    .490
 (06-sad,07-for)      .500     .520    .520    .360    .310     .720    3.500   3.000    .860    .090    .110    .240
 08-pak               .505     .520    .523    .362    .313     .729    3.520   3.021    .873    .091    .112    .242
 09-ld                .689     .520    .618    .437    .408    1.042    4.186   3.749   1.329    .139    .187    .324
 10-lbd               .564     .522    .543    .380    .360     .778    3.545   3.130   1.018    .111    .122    .278
 11-jha               .627     .686    .662    .509    .465    1.239    3.286   4.035   1.564    .139    .112    .315
 12-gug               .548     .624    .608    .476    .430    1.156    2.954   3.690   1.436    .126    .108    .290
 13-uj               1.210    1.170   1.240    .850    .690    1.500    5.320   5.610   2.030    .230    .160    .480
 14-lj                .920     .889    .942    .646    .524    1.140    4.043   4.264   1.543    .175    .122    .365
 15-bah               .230     .290    .380    .250    .170     .260    2.520   1.930    .190    .020    .120    .130
 16-mai               .402     .434    .463    .316    .270     .553    3.020   2.530    .618    .068    .107    .208
 17-sid               .420     .435    .445    .305    .325     .495    2.570   2.265    .665    .085    .080    .245
 18-hav               .390     .500    .500    .410    .360     .990    2.290   3.000   1.180    .100    .100    .240
 19-ran               .370     .380    .400    .270    .330     .550    2.010   1.820    .540    .080    .060    .240
 20-pan               .220     .252    .296    .184    .140     .184    1.716   1.444    .152    .010    .108    .120
 21-abb               .221     .263    .326    .205    .143     .206    1.998   1.615    .163    .011    .102    .121
 (22-usw,23-lsw,
  24-war,25-kab)     1.440    1.530   2.440   1.760    .770     .310    1.260   2.030    .810    .230    .310    .670
 26-tha               .610     .717    .882    .621    .448     .809    2.946   2.769    .944    .132    .130    .298
 27-pah               .450     .670    .960    .690    .390     .610    2.290   1.900    .630    .110    .150    .240
 28-muz               .375     .400    .435    .298    .335     .555    2.028   1.825    .548    .083    .068    .240
 29-dgk               .260     .330    .270    .220    .180     .380     .740   1.230    .310    .030    .070    .210
 (31-p+d,32-beg)      .199     .282    .204    .179    .100     .257     .447    .905    .187   0.000    .097    .179
 33-gho               .210     .230    .250    .140    .120     .140    1.230   1.160    .130   0.000    .110    .110
 34-nw                .171     .336    .263    .171    .100     .242     .612    .994    .330   0.000    .035    .173
 35-ric               .173     .341    .273    .173    .100     .251     .706   1.020    .346   0.000    .032    .163
 36-dad               .170     .335    .270    .170    .100     .275    1.330   1.170    .400   0.000    .035    .130
 37-kw                .190     .240    .200    .140    .100     .250    1.310   1.210    .270   0.000    .090    .100
 38-ke                .150     .200    .200    .100    .100     .250    1.350   1.450    .350   0.000    .050    .100
 39-roh               .138     .235    .188    .143    .150     .368    2.297   1.876    .564    .034    .051    .095
 41-nar               .091     .151    .119    .103    .159     .397    3.003   2.097    .658    .071    .056    .060
 (42-kal,43-lch,
  44-ful,45-pin)      .140     .300    .120    .060    .140     .720    4.020   2.080    .700    .060    .060    .120

$stitle historic canal diversions(maf)

 table divpost(cnl,m1) average(1976-77 to 1987-88) canal diversions(maf)
*      source : indus basin irrigation system
*               historic rivers and canals discharge data,
*               water resources management, directorate. wapda lahore

           apr   may   jun    jul   aug   sep   oct   nov   dec   jan   feb   mar   rabi kharif   annual
01-ud    .0238 .1028 .1091  .1155 .1183 .1263 .0730 .0236 .0072 .0100 .0186 .0343  .1668  .5958   0.7626
02-cbd   .1102 .1339 .1372  .1278 .1151 .1202 .1284 .1294 .1113 .0500 .0958 .1010  .6159  .7443   1.3602
03-ray   .0025 .0370 .0884  .0902 .0820 .0855 .0473 .0025 .0020 .0005 .0024 .0007  .0553  .3856   0.4409
04-uc    .0984 .1829 .2771  .3155 .2713 .2577 .1928 .0822 .0343 .0428 .0631 .0566  .4719 1.4029   1.8747
05-mr    .0163 .0964 .1218  .2160 .2366 .1665 .0402 .0122 .0107 .0111 .0176 .0103  .1020  .8536   0.9556
06-sad   .2686 .2865 .2838  .2762 .2677 .2693 .2538 .2851 .2719 .1106 .2329 .2796 1.4339 1.6522   3.0861
07-for   .0340 .1224 .1481  .1433 .1415 .1481 .1044 .0102 .0091 .0061 .0137 .0137  .1572  .7375   0.8946
08-pak   .2183 .3408 .3542  .3536 .3881 .3739 .2978 .2190 .2399 .1134 .1833 .2470 1.3003 2.0289   3.3292
09-ld    .0479 .1923 .2140  .2295 .2320 .2278 .1535 .0451 .0258 .0256 .0479 .0295  .3273 1.1436   1.4709
10-lbd   .4055 .4822 .4717  .4564 .4611 .4590 .4244 .4269 .3860 .1539 .3521 .4069 2.1503 2.7360   4.8862
11-jha   .2838 .3425 .3423  .3177 .3010 .3391 .3284 .3025 .2970 .1279 .1941 .2728 1.5227 1.9265   3.4492
12-gug   .3186 .3845 .3843  .3567 .3379 .3807 .3687 .3396 .3334 .1436 .2179 .3063 1.7094 2.1626   3.8720
13-uj    .0791 .1180 .1285  .1486 .1294 .1305 .1117 .0884 .0732 .0461 .0457 .0462  .4114  .7341   1.1454
14-lj    .2421 .3128 .3155  .3099 .2902 .2901 .2987 .2509 .2151 .1084 .2138 .2118 1.2988 1.7607   3.0594
15-bah   .1205 .2033 .2293  .2448 .2282 .2359 .2008 .1238 .1212 .0958 .0827 .1475  .7718 1.2619   2.0337
16-mai   .1441 .2961 .3286  .3163 .3375 .3478 .2875 .1336 .1392 .1309 .0897 .1643  .9452 1.7704   2.7156
17-sid   .1469 .2226 .2334  .2216 .2264 .2332 .2203 .1631 .1593 .0975 .1124 .1632  .9159 1.2841   2.2000
18-hav   .0314 .0455 .0474  .0447 .0464 .0513 .0472 .0252 .0282 .0224 .0255 .0318  .1803  .2666   0.4469
19-ran   .0314 .0809 .0904  .0823 .0689 .0817 .0663 .0174 .0147 .0206 .0076 .0242  .1507  .4356   0.5863
20-pan   .2684 .4531 .5359  .5554 .5569 .5577 .5056 .2296 .1913 .1355 .1759 .2235 1.4613 2.9275   4.3888
21-abb   .0408 .0529 .0602  .0599 .0604 .0602 .0567 .0375 .0328 .0202 .0280 .0357  .2109  .3343   0.5453
22-usw   .1025 .1131 .1101  .1094 .1028 .1048 .1052 .0876 .0802 .0631 .0654 .0869  .4884  .6426   1.1310
23-lsw   .0413 .0518 .0514  .0461 .0401 .0455 .0487 .0388 .0372 .0051 .0239 .0308  .1845  .2762   0.4607
24-war   .0377 .0441 .0436  .0385 .0307 .0412 .0454 .0436 .0406 .0021 .0145 .0263  .1725  .2358   0.4083
25-kab   .0339 .0434 .0466  .0440 .0366 .0425 .0434 .0392 .0328 .0017 .0177 .0341  .1688  .2470   0.4158
26-tha   .3777 .4206 .3864  .4027 .3991 .4096 .4016 .3934 .3680 .1064 .3175 .3615 1.9483 2.3962   4.3445
27-pah   .0269 .0340 .0355  .0363 .0307 .0340 .0348 .0301 .0315 .0045 .0259 .0288  .1556  .1976   0.3532
28-muz   .1041 .3154 .3886  .3937 .3812 .3749 .2665 .1312 .0745 .0707 .0920 .1057  .7407 1.9578   2.6985
29-dgk   .1874 .3408 .4056  .4062 .3455 .3699 .2960 .1787 .1091 .0935 .1275 .1580  .9629 2.0555   3.0183
31-p+d   .0020 .0812 .5064  .5786 .4625 .4169 .2224 .0942 .0700 .1628 .0665 .0632  .6791 2.0475   2.7266
32-beg   .0017 .0798 .7547 1.0590 .6712 .5522 .1528 .0023 .0122 .1344 .0059 .0038  .3115 3.1186   3.4300
33-gho   .0134 .2218 .4264  .4093 .4541 .4266 .3847 .2007 .1010 .2016 .0720 .1217 1.0817 1.9517   3.0334
34-nw    .0917 .1070 .3522  .4994 .4373 .3499 .3011 .2055 .1948 .0973 .2070 .1844 1.1900 1.8375   3.0275
35-ric   .0000 .1199 .5937  .5945 .4428 .3722 .1454 .0000 .0000 .0000 .0000 .0000  .1454 2.1232   2.2686
36-dad   .0775 .0820 .1964  .2978 .2643 .2269 .1894 .1350 .1378 .0620 .1494 .1400  .8137 1.1449   1.9587
37-kw    .0904 .1003 .1086  .1123 .1112 .1283 .1187 .1064 .1073 .0342 .0908 .0938  .5512  .6511   1.2024
38-ke    .1140 .1357 .1464  .1564 .1537 .1504 .1609 .1456 .1456 .0442 .1185 .1186  .7334  .8567   1.5901
39-roh   .6860 .7981 .7821  .8123 .8474 .8517 .8026 .7187 .7404 .2608 .7123 .7525 3.9873 4.7776   8.7649
41-nar   .6064 .6819 .6671  .6932 .6901 .6967 .6916 .6246 .6078 .2227 .5419 .5689 3.2576 4.0355   7.2931
42-kal   .1342 .1644 .2888  .4046 .3561 .3739 .2843 .1374 .0777 .1522 .0963 .0861  .8340 1.7220   2.5559
43-lch   .1155 .1476 .1645  .1671 .1421 .1451 .1332 .0954 .0612 .0618 .0752 .0780  .5048  .8819   1.3867
44-ful   .1021 .3667 .6712  .6460 .4545 .4078 .2991 .1375 .0732 .1594 .1005 .0591  .8287 2.6483   3.4770
45-pin   .0464 .1899 .4162  .5242 .4391 .3334 .2264 .1068 .0469 .1041 .0611 .0349  .5802 1.9493   2.5295

$stitle govt. tubewell pumpage and depth to water table
   table gwt(cnl,m)  public tunewell pumpage(kaf)
                 jan  feb  mar  apr  may  jun  jul  aug  sep  oct  nov  dec
  01-ud            0    0    1    0    1    0    1    1    1    0    0    0
  02-cbd           2    2    3    2    2    2    6    7    8    2    2    2
  03-ray          15   29   41   45   40   50   30   17   48   43   19    7
  04-uc           39   49   49   30   49   49   55   58   58   16   36   15
  05-mr            3    6   20   17    9   11    8    6   12    9    5    2
  09-ld            1    1    1    1    1    1    0    0    0    1    1    1
  10-lbd           2    2    2    2    2    2    2    3    3    1    2    2
  11-jha          38   42   49   57   57   57   50   54   60   57   48   32
  12-gug          50   62   69   78   79   79   69   73   82   85   67   34
  13-uj           67   88  104  119  141  121  114  114  140  155  117    5
  14-lj           51   65   88   70   81   72   87   81   94   91   70   77
  15-bah           0    0    1    0    1    1    2    2    2    0    1    1
  18-hav           6    6    8    6    7    7    6    6    6    5    7    7
  19-ran          55   66   82   84   68   52   36   57   59   52   51   43
  20-pan           6    6    8    6    6    6    6    6    5    5    6    6
  21-abb           1    1    2    1    2    2    1    2    1    1    2    2
  28-muz         114  114  159  147  109   89   80   90   94   45   79  119
  32-beg           2    3    4    2    2    1    0    1    1    2    2    2
  37-kw           18   20   24   19   16   20   18   18   16   10   11   11
  39-roh           6    7    8    6    5    7    6   48   42    3    4    4


   table dep1(cnl,is, * )     depth to water table (feet)
*     1980 depth is from ibm and irrgation supporting report of rap.
*     year 1993 and 2000 depths are estimated considering the drainage
*     projects under construction or planning.
                               1980.depth    (1988,1993,2000).depth
  01-ud                          17.0             17
  02-cbd                         10.5             10
 (03-ray,05-mr)                  15.0             15
  04-uc                          15.1             15
  06-sad                          5.5              7
 (07-for,22-usw,23-lsw ,
  24-war,25-kab)                  8.0              8
  08-pak                         19.9             15
  09-ld                          17.9             17
  10-lbd                         18.1             18
  11-jha                         13.6             14
  12-gug                         14.8             14
  13-uj                          10.0             10
  14-lj                           8.1              9
 (15-bah,19-ran,28-muz)           9.0              9
 (16-mai,17-sid)                 19.0             19
  18-hav                         10.1             10
  20-pan                          8.4              9
  21-abb                          8.2              9
  26-tha                          8.7              8
  27-pah                         10.0             12
  29-dgk                          7.8              8
  (31-p+d,32-beg,36-dad,39-roh)                    7
  (34-nw,35-ric,42-kal,43-lch,
            44-ful,45-pin)                         6
  33-gho                          5.5              7
  37-kw                           5.6              7
  38-ke                           5.4              7
  41-nar                          4.9              7

   table dep2(is,cnl,m)   depth to water table (feet)
                                      jan  feb  mar  apr  may  jun  jul  aug  sep  oct  nov  dec
  1980. 31-p+d                        6.9  6.9  6.9  6.9  6.9  4.5  2.5  2.5  2.5  3.5  4.5  6.9
  1980. 32-beg                        6.6  6.6  6.6  6.6  6.6  4.5  2.5  2.5  2.5  3.5  4.5  6.6
  1980.(34-nw,35-ric)                 4.6  4.6  4.6  4.6  4.6  4.5  2.5  2.5  2.5  3.5  4.5  4.6
  1980. 36-dad                        5.6  5.6  5.6  5.6  5.6  4.5  2.5  2.5  2.5  3.5  4.5  5.6
  1980. 39-roh                        5.9  5.9  5.9  5.9  5.9  5.8  5.8  5.8  5.8  5.9  5.9  5.9
  1980.(42-kal,43-lch,44-ful,45-pin)  5.0  5.0  5.0  5.0  5.0  3.5  2.0  2.0  2.5  4.0  5.0  5.0

parameters
   depth(cnl,m)      depth to groundwater (feet)
   efr(cnl,m)        effective rainfall in feet
   eqevap(cnl,m)     evaporation from the equaifer   (feet)
   subirr(cnl,m)     water supplied by capillary action from the aquifer
   subirrfac(z)      maximum sub-irrigation in saline areas as proportion of crop req.(net of rain)
                         /(nwfp, pmw,pcw,psw,prw,srwn,scws) .4,
                                                       scwn .3,
                                                       srws .1/
                 ;
 scalars drc  run-off portion of rainfall /.15/
         the1 portion of equaifer evaporation used by crops /0.6/;

*  sub-irrigation and effective rain calculations
 loop(isr,
    depth(cnl,m) = dep1(cnl,isr,"depth");
    depth(cnl,m)$(dep2(isr,cnl,m) ne 0 ) = dep2(isr,cnl,m)  ;
    efr(cnl,m)   = (1.0 - drc - (1-comdef(isr,"flde",cnl)))*
                                           rain(cnl,m)/12.0
      ) ;
    eqevap(cnl,m) = min(1., 10.637/depth(cnl,m)**2.558)*evap(cnl,m)  ;
    subirr(cnl,m) = eqevap(cnl,m) * the1 ;

   options depth:1;   display depth ;

$stitle network nodes and arc definition

 set n  nodes of the indus river system/
          sulem-b    sulemanki barrage                 - sutlej river
          islam-b    islam barrage                     - sutlej river
          panjnad-b  panjnad barrage                   - panjnad river
          ravi-i     ravi inflow at madhopur           - ravi river
          balloki-b  balloki barrage                   - ravi river
          sidhnai-b  sidhani barrage                   - ravi river
          marala-b   marrala barrage                   - chenab river
          khanki-b   khanki barrage                    - chenab river
          qadira-b   qadirabad barrage                 - chenab river
          trimmu-b   trimmu barrage                    - chenab river
          mangla-r   mangla reservoir                  - jehlum river
          rasul-b    rasul barrage                     - jhelum river
          tarbela-r  tarbela reservoir                 - indus river
          amanda-h   amandara head works               - swat river
          munda-h    munda head works                  - swat river
          warsak-d   warsak reservoir                  - kabul river
          k-s-jct    kabul and swat river jct          - kabul river
          kalabagh-r kalabagh reservoir                - indus river
          chasma-r   chasma reservoir                  - indus river
          taunsa-b   taunsa barrage                    - indus river
          gudu-b     gudu barrage                      - indus river
          sukkur-b   sukkur barrage                    - indus river
          nara-jct   nara junction                     - nara complex
          nara-head  head works for irri. diversion    - nara complex
          chotiari-r chotiari reservoir                - nara complex
          kotri-b    kotri barrage                     - indus river
          a-sea      arabian sea
          a1         diversion point for uj uj link and r.p.c
          a2         div. for brbd link uc link and uc int.
          a3         div. for mr cross link
          a4         mr cross tail and brbd link
          a5         div. for cbd and ud canals
          a6         diversion to lcc feeder
          a7         div. for gugera and jhang canals
          a8         divesion point for maili+l pakpattan canals
          a9         smb link to l-bahawl canal
          a10        diversion point for lj
          /
     i  system inflows /
          swat     swat river at chakdara
          kabul    kabul river at warsak
          indus    indus river at tarbela
          haro     haro river at gariall
          soan     soan river at dhok pathan
          jehlum   jehlum river at mangls
          chenab   chenab river at marrala
          ravi     ravi river below madhopur
          sutlej   sutlej river below ferozpur /
$eject
     nc(n,cnl) node to canal map/
          sulem-b.   (06-sad,07-for,08-pak),     a8.         16-mai
          a9.         15-bah               ,     panjnad-b. (20-pan,21-abb)

          balloki-b. (10-lbd,09-ld)        ,     sidhnai-b.  17-sid

          marala-b.   05-mr                ,     a2.        (04-uc,03-ray)
          a5.        (02-cbd,01-ud)        ,     a7.        (11-jha,12-gug)
          trimmu-b.  (19-ran,18-hav)
          a1.         13-uj                ,     a10.        14-lj
          amanda-h.   22-usw               ,     munda-h.    23-lsw
          warsak-d.   24-war               ,     warsak-d.   25-kab
          chasma-r.  (26-tha,27-pah)
          taunsa-b.  (28-muz,29-dgk)       ,     gudu-b.    (31-p+d,32-beg,33-gho)
          sukkur-b.  (34-nw,35-ric,36-dad,37-kw,38-ke,39-roh)
          nara-head.  41-nar
          kotri-b.   (42-kal,43-lch,44-ful,45-pin)/
             ;
$eject
 alias (n,n1);
 set nn(n,n1) water flow system node to node /
*--  sutlej ravi system
          sulem-b.         balloki-b
          islam-b.         sulem-b
          a9.             (a8,       islam-b)
          a8.              sidhnai-b
          panjnad-b.      (trimmu-b,sidhnai-b, islam-b, taunsa-b)
*--  ravi chenab system
          ravi-i.          a3
          a3.              marala-b
          a4.             (a3,a2)
          a5.              a4
          a2.              marala-b
          balloki-b.      (ravi-i,   a2,        a6)
          a6.              qadira-b
          a7.             (a6,       khanki-b)
          sidhnai-b.      (trimmu-b, balloki-b)
*--  chenab jehlum system
          khanki-b.       (a1,       marala-b)
          qadira-b.       (rasul-b,  khanki-b)
          trimmu-b.       (qadira-b, rasul-b, chasma-r)
          a1.              mangla-r
          rasul-b.        (mangla-r, a10)
          a10.            (rasul-b,  a1)
*--  kabul swat system
          munda-h.         amanda-h
          k-s-jct.        (munda-h,  warsak-d)
*--  indus river
          kalabagh-r.     (k-s-jct,  tarbela-r)
          chasma-r.        kalabagh-r
          taunsa-b.        chasma-r
          gudu-b.         (panjnad-b,taunsa-b)
          sukkur-b.        gudu-b
          nara-jct.        sukkur-b
          chotiari-r.      nara-jct
          nara-head.      (nara-jct,  chotiari-r)
          kotri-b.         sukkur-b
          a-sea.           kotri-b /

     ni(n,i) node to rim station inflow map /
          amanda-h.    swat    ,  warsak-d.    kabul,
          tarbela-r.   indus   ,  kalabagh-r.  haro
          kalabagh-r.  soan    ,  mangla-r.    jehlum
          marala-b.    chenab  ,  ravi-i.      ravi
          sulem-b.     sutlej  /

     nb(n)    ;  nb(n) = yes; nb("a-sea")=no;
$eject
*-    flow capacities

 parameter
    ncap(n,n1) node to node transfer capacity (maf)/
         chasma-r.   trimmu-b  1.307,     taunsa-b.  panjnad-b   .724
         sukkur-b.   nara-jct  1.190,     nara-jct.  nara-head  1.190
         nara-jct.   chotiari-r .416,     chotiari-r.nara-head   .135

         mangla-r.  a1         .685,      a1.        a10         .199
         rasul-b.   a10        .328,      a1.        khanki-b    .422
         rasul-b.   qadira-b  1.232,

         marala-b.  a3         .985,      a3.        a4          .199
         a3.        ravi-i     .938,      marala-b.  a2          .995
         a2.        balloki-b  .698,      a2.        a4          .310
         a4.        a5         .310

         khanki-b.  a7         .684,      qadira-b.  a6         1.146
         a6.        a7         .270,      a6.        balloki-b   .808

         balloki-b. sulem-b   1.172,      trimmu-b.  sidhnai-b  1.03
         sidhnai-b. a8         .662,      a8.        a9          .640
         islam-b.   a9         .299 /

    lloss(n,n1) link canal loss factors /
         chasma-r.   trimmu-b   .096,     taunsa-b.  panjnad-b   .0615
         sukkur-b.   nara-jct   .055,     nara-jct.  nara-head   .0001
         nara-jct.  chotiari-r  .001,     chotiari-r. nara-head  .0001

         a1.        a10        .004 ,    a1.        khanki-b   .099
         rasul-b.   qadira-b   .039 ,    marala-b.  a3         .0705
         a3.        a4         .0075,    a2.        balloki-b  .0705
         a2.        a4         .097 ,    a4.        a5         .007

         khanki-b.  a7        .0435 ,    a6.        balloki-b .087

         balloki-b. sulem-b   .0525,      trimmu-b.  sidhnai-b .075
         sidhnai-b. a8        .0495 /
    lceff(n,n1)  link canal efficiency from head to tail
           ;
    lceff(n,n1)$lloss(n,n1) = 1 - lloss(n,n1) ;

* display n,nc,nn, ni;
* display ncap, lloss, lceff  ;
$stitle flow routing coefficints
 set  cd /c,d/
 table rivercd(n,cd)   coefficients for river routing
                            c        d
   rasul-b               .02732     .93766
   trimmu-b              .15348
   khanki-b                        2.2113
   balloki-b
   sidhnai-b             .14063    1.0024
   islam-b               .12294
   panjnad-b             .08077
   gudu-b                .10947
   kalabagh-r            .04856
   chasma-r              .04526
   taunsa-b              .07205
   sukkur-b              .02218
   kotri-b               .17054

table riverb(n,n1)  coefficients for river routing

            mangla-r  rasul-b  qadira-b marala-b ravi-i  balloki-b
rasul-b      1.01841
trimmu-b              .82803   .9068
khanki-b                                 .93361
balloki-b                                        1.2181
sidhnai-b                                                 .70555

  +         sulem-b     trimmu-b sidhnai-b islam-b taunsa-b panjnad-b
islam-b      .9086
panjnad-b               .81033    1.13626  .80359
gudu-b                                             .94003   1.0

+           k-s-jct tarbela-r  kalabagh-r  chasma-r gudu-b sukkur-b


kalabagh-r  1.04621 1.08686
chasma-r                       .9849
taunsa-b                                   .9497
sukkur-b                                            .9908
kotri-b                                                    .70625

   ;

 riverb(n,n1)$(riverb(n,n1) eq 0 )    = 1;
 rivercd(n,"d")$(rivercd(n,"d") eq 0) = 1;


$stitle rim station and tributory inflows
  set s58/50, 80 /

  table infl5080(s58,i,m1) system inflows measured atthe rim stations (maf)

*     50 - 50 % probability flows and 80 are 80 % probability flows.
*
* source:   water resouce management directorate, wapda publications and
*           rap irrigation and drainage supporting report
*           historic data used: indus  oct 1936 - march 1988
*                               jehlum april 1922 - march 1988
*                               chenab april 1922 - march 1988
*                      ravi and sutlej april 1966 - march 1988
*                     all others       april 1966 - march 1976

                          apr    may    jun    jul    aug    sep kharif
50.indus                1.931  3.984  7.191 16.148 17.437  6.196 52.887
50.jehlum               3.826  4.355  2.993  3.429  2.236  0.846 17.685
50.chenab               1.361  2.281  4.036  5.279  5.097  2.540 20.594
50.ravi                 0.226  0.161  0.038  1.089  1.732  0.841  4.088
50.sutlej               0.031  0.045  0.052  0.483  1.003  1.198  2.810

50.haro                 0.016  0.013  0.014  0.048  0.204  0.073  0.369
50.soan                 0.012  0.007  0.090  0.225  0.426  0.068  0.829
50.swat                 0.483  0.636  1.049  0.863  0.642  0.346  4.020
50.kabul                1.512  1.883  3.467  3.490  2.619  0.975 13.945

80.indus                1.871  3.760  9.634 13.135 10.524  5.234 44.159
80.jehlum               1.526  2.242  2.594  2.632  3.140  1.623 13.757
80.chenab               1.120  2.064  3.133  4.515  3.537  2.844 17.213
80.ravi                 0.049  0.020  0.011  0.696  1.346  0.313  2.434
80.sutlej               0.000  0.001  0.000  0.145  0.291  0.004  0.441

80.haro                 0.014  0.019  0.011  0.020  0.101  0.022  0.187
80.soan                 0.006  0.004  0.020  0.252  0.145  0.021  0.448
80.swat                 0.323  0.574  0.848  0.579  0.463  0.448  3.235
80.kabul                0.963  1.737  2.238  2.157  2.303  1.534 10.934

 +                        oct    nov    dec    jan    feb    mar  rabi
50.indus                2.297  1.279  1.188  0.952  0.973  1.468  8.157
50.jehlum               0.988  0.521  0.412  0.460  0.621  1.417  4.419
50.chenab               0.728  0.487  0.407  0.480  0.711  1.079  3.891
50.ravi                 0.188  0.134  0.131  0.090  0.026  0.455  1.023
50.sutlej               0.129  0.009  0.004  0.060  0.000  0.000  0.202

50.haro                 0.037  0.020  0.016  0.019  0.019  0.051  0.163
50.soan                 0.010  0.005  0.020  0.022  0.028  0.043  0.128
50.swat                 0.195  0.122  0.088  0.079  0.068  0.140  0.692
50.kabul                0.523  0.405  0.354  0.379  0.440  0.687  2.787

80.indus                2.133  1.281  1.060  0.898  0.835  1.025  7.232
80.jehlum               0.405  0.270  0.309  0.306  0.526  1.514  3.331
80.chenab               0.516  0.325  0.285  0.292  0.804  0.711  2.933
80.ravi                 0.076  0.048  0.074  0.135  0.112  0.135  0.580
80.sutlej               0.004  0.002  0.004  0.005  0.002  0.000  0.017

80.haro                 0.024  0.016  0.018  0.017  0.018  0.018  0.112
80.soan                 0.012  0.010  0.009  0.010  0.019  0.031  0.091
80.swat                 0.103  0.070  0.061  0.064  0.085  0.212  0.597
80.kabul                0.273  0.420  0.354  0.331  0.429  0.359  2.167


  table tri(s58,n1,n,m1) tributary inflows (maf)
*      source: irrigation drainage and flood management supporting report of rap. planning division wapda.
                          apr    may    jun    jul    aug    sep kharif
50.marala-b. khanki-b   0.004  0.006  0.023  0.120  0.240  0.056  0.449
50.mangla-r. rasul-b    0.120  0.180  0.150  0.170  0.320  0.210  1.150
50.chasma-r. taunsa-b   0.051  0.013  0.021  0.048  0.021  0.012  0.166
50.tarbela-r.kalabagh-r 0.041  0.087  0.182  0.350  0.602  0.180  1.442
50.balloki-b.sidhnai-b  0.000  0.001  0.009  0.096  0.206  0.041  0.353
50.amanda-h. munda-h    0.036  0.337  0.637  0.854  0.544  0.062  2.470
50.munda-h.  k-s-jct    0.015  0.013  0.005  0.008  0.007  0.007  0.054
50.warsak-d. k-s-jct    0.126  0.102  0.059  0.070  0.056  0.056  0.469

80.marala-b. khanki-b   0.004  0.005  0.006  0.136  0.131  0.011  0.293
80.mangla-r. rasul-b    0.120  0.180  0.150  0.170  0.320  0.210  1.150
80.chasma-r. taunsa-b   0.021  0.014  0.006  0.047  0.008  0.003  0.098
80.tarbela-r.kalabagh-r 0.052  0.094  0.106  0.126  0.416  0.079  0.873
80.balloki-b.sidhnai-b  0.000  0.000  0.006  0.075  0.164  0.027  0.271
80.amanda-h. munda-h    0.121  0.410  0.161  0.132  0.000  0.019  0.844
80.munda-h.  k-s-jct    0.006  0.009  0.003  0.005  0.005  0.004  0.031
80.warsak-d. k-s-jct    0.050  0.077  0.027  0.039  0.036  0.032  0.261

 +                        oct    nov    dec    jan    feb    mar   rabi
50.marala-b. khanki-b   0.017  0.007  0.008  0.007  0.006  0.007  0.052
50.mangla-r. rasul-b    0.060  0.000  0.000  0.000  0.000  0.000  0.060
50.chasma-r. taunsa-b   0.004  0.005  0.010  0.006  0.015  0.026  0.065
50.tarbela-r.kalabagh-r 0.100  0.064  0.038  0.055  0.052  0.130  0.439
50.balloki-b.sidhnai-b  0.002  0.003  0.008  0.005  0.007  0.006  0.030
50.amanda-h. munda-h    0.133  0.153  0.150  0.119  0.119  0.221  0.894
50.munda-h.  k-s-jct    0.003  0.003  0.004  0.004  0.004  0.009  0.028
50.warsak-d. k-s-jct    0.024  0.027  0.035  0.035  0.038  0.070  0.229

80.marala-b. khanki-b   0.007  0.006  0.006  0.003  0.011  0.006  0.039
80.mangla-r. rasul-b    0.060  0.000  0.000  0.000  0.000  0.000  0.060
80.chasma-r. taunsa-b   0.009  0.006  0.010  0.006  0.008  0.011  0.050
80.tarbela-r.kalabagh-r 0.088  0.039  0.041  0.040  0.042  0.045  0.295
80.balloki-b.sidhnai-b  0.000  0.001  0.001  0.000  0.000  0.000  0.002
80.amanda-h. munda-h    0.090  0.075  0.045  0.047  0.019  0.121  0.397
80.munda-h.  k-s-jct    0.005  0.004  0.003  0.003  0.000  0.005  0.022
80.warsak-d. k-s-jct    0.020  0.016  0.020  0.020  0.027  0.065  0.169
     ;
  parameter inflow(i,m)  inflows for this run (maf)
            trib(n1,n,m) tributary inflows for this run (maf) ;
 inflow(i,m) = infl5080("50",i,m);  trib(n1,n,m) = tri("50",n1,n,m);

$stitle  reservoir characteristics
  table   rrcap(n,is)    live storage capacity of reservoirs (maf)
*
*
                            1980       1988     2000
           mangla-r         5.053      4.881    4.617
           tarbela-r        9.197      8.861    8.357
           kalabagh-r                           6.100
           chasma-r         0.571       .435     .231
           chotiari-r                           1.100

* storage capacities are provided by planning division of wapda.
* 1988 capacities are from recent hydrographic surveys.
*

 table rulelo(n,m)  lower rule curve
*                   end of month contents as percent of live capacity

            jan  feb  mar  apr  may  jun  jul  aug  sep  oct  nov  dec
 mangla-r    16    6    0    7   21   48   74   95   83   55   35   24
 tarbela-r   43   30   14    4    0   10   46  100  100   64   59   53
 kalabagh-r  55   49   34   28    0    0   24   70  100   96   78   61
 chasma-r     0    0    0    0    0    0    0    0    0    0    0    0
 chotiari-r   0    0    0    0    0    0    0    0    0    0    0    0

     table ruleup(n,m) upper rule curve
*                   end of month contents as percent of live capacity

                 jan  feb  mar  apr  may  jun  jul  aug  sep  oct  nov  dec
 mangla-r         39   28   19   31   42   59   87  100  100   90   68   49
 tarbela-r        69   47   37   26   13   48   73  100  100  100   87   77
 kalabagh-r      100  100  100  100  100  100  100  100  100  100  100  100
 chasma-r        100  100  100  100  100  100  100  100  100  100  100  100
 chotiari-r      100  100  100  100  100  100  100  100  100  100  100  100

     table  revapl(n,m)  evaporation losses from reservoirs (kaf)

             jan  feb  mar  apr  may  jun    jul    aug  sep  oct  nov  dec
 mangla-r     1    5    4    5    12   17   -1.9   -2.5    4   22   10    5
 tarbela-r    1    5    7    9    16   28     10      1   12   29   12    5
 kalabagh-r   1    4    6    9    12    5      9      1    9   24    9    5
 chasma-r     3    3    1                      6      3    7   15    6    3

$eject
set  pow /r-ele  reservoir elevation (feet from spd)
          p-cap  installed capacity of the power house at r-ele
          g-cap  generation capability (kwh per af) at r-ele
          r-cap  gross reservoir capacity (maf) at r-ele/
     pn(n) nodes with power house /mangla-r,tarbela-r/
     v   /1*27/;

    table powerchar(n,*,v)    power generation chractersitics of hrdro stations

                     1      2     3     4     5     6     7     8     9    10    11    12    13
 tarbela-r.r-ele    1300  1310  1320  1330  1340  1350  1360  1370  1380  1390  1400  1410  1420
 tarbela-r.p-cap     480   552   630   700   780   832   894   942   992  1052  1114  1184  1266
 tarbela-r.g-cap     136   166   183   192   202   209   216   224   233   237   247   257   263
 tarbela-r.r-cap   1.313 1.439 1.571 1.707 1.848 1.995 2.236 2.486 2.745 3.014 3.297 3.631 3.977

 mangla-r .r-ele    1040  1050  1055  1060  1065  1070  1075  1080  1085  1090  1095  1100  1105
 mangla-r .p-cap     384   424   440   464   488   504   528   544   568   584   608   632   648
 mangla-r .g-cap     141   155   164   173   183   192   198   205   209   212   220   228   233
 mangla-r .r-cap   0.417 0.509 0.555 0.601 0.648 0.696 0.744 0.792 0.883 0.973 1.064 1.155 1.271



       +              14    15    16    17    18    19    20    21    22    23    24    25     26     27
 tarbela-r.r-ele     1430  1440  1450  1460  1470  1480  1490  1500  1510  1520  1530  1540   1550   1550
 tarbela-r.p-cap     1308  1366  1416  1476  1530  1596  1678  1730  1792  1848  1908  1936   1952   1952
 tarbela-r.g-cap      269   281   288   295   303   310   318   336   346   356   356   368    376    376
 tarbela-r.r-cap    4.334 4.705 5.095 5.520 5.960 6.412 6.882 7.369 7.865 8.381 8.917 9.478 10.063 20.000

 mangla-r .r-ele     1110  1115  1120  1130  1135  1140  1145  1155  1160  1165  1170  1180   1202   1202
 mangla-r .p-cap      672   696   712   760   776   800   816   856   880   896   920   960    960    960
 mangla-r .g-cap      237   247   252   263   263   269   275   281   288   295   303   318    356    356
 mangla-r .r-cap    1.386 1.501 1.617 1.933 2.091 2.249 2.453 2.862 3.066 3.316 3.566 4.067  5.364 20.000
  ;
* segment 27 would be used if the
* dam is raised to create additional storage on top of the existing.

* select capacity and adjust the  capacity curve for sediment/dead storage
 parameter rcap(n)  live capacity of resrvoirs (maf);
    loop(isr,
          rcap(n)           = rrcap(n,isr) );

powerchar(pn,"r-cap",v) = max(0, powerchar(pn,"r-cap",v)-
                                (powerchar(pn,"r-cap","26") -rcap(pn)));
powerchar(pn,"r-cap","1")= 0;

* $ontext
*-  report on the inflows etc.

   parameter rep7
             rep8
               ;
 rep7(i,m) = inflow(i,m);
 rep7(i,sea1)=sum(m$sea1m(sea1,m), inflow(i,m));
 rep7("total",m1)        = sum(i, rep7(i,m1));
 display " system inflows at rim stations (maf) ", rep7;

*
 rep8(n,n1,m)=trib(n,n1,m);
 rep8(n,n1,sea1) = sum(m$sea1m(sea1,m), trib(n,n1,m));
 display " tributory inflow in (maf)", rep8;
 rep8(n,n1,m1) = 0 ;
 rep8(pn,v,pow) = powerchar(pn,pow,v) ;
 option rep8:3:1:1 ;
 display  rcap
         "r-ele  reservoir elevation (feet from spd)"
         "p-cap  installed capacity of the power house at r-ele"
         "g-cap  generation capability (kwh per af) at r-ele"
         "r-cap  live  reservoir capacity (maf) at r-ele ",  rep8 ;
*-  -  -  -
* $offtext

$stitle  prices
 set p3 /financial,economic, export,import/

  table prices(ps,cq,p3)  1988  prices

                  financial  economic    export  import
 87-88.basmati        6.         6.       4.8
 87-88.irri           2.9        2.4      2.3
 87-88.cotton         4.5        6.       4.5
 87-88.gram           3.9        3.8
 87-88.maize          1.6         .6
 87-88.mus+rap        3.75       3.4
 87-88.sc-gur         3.0        3.0
 87-88.sc-mill         .3         .21             .5
 87-88.wheat          2.0        2.2      2.0
 87-88.orchard        4.8        4.3
 87-88.potatoes       2.6        2.6
 87-88.onions         2.9        2.9      2.8
 87-88.chilli         10.       10.
 87-88.cow-milk        5.4      2.46
 87-88.buff-milk       5.4      3.24
 87-88.meat           13.8      14.4


 table finsdwtpr(c,ps,*)  prices of seed(rs per kg) and water(rs per acre )

              87-88.seed  87-88.water 87-88.miscc
basmati        5.          32.           75.
irri           2.5         32.           75.
cotton         5.5         34.          500.
rab-fod        40.         11.
gram           5.7         16.
maize          4.          14.
mus+rap        4.3         22.
kha-fod        2.5         14.
sc-gur          .2         64.
sc-mill         .2         64.
wheat          2.9         22.
orchard                    84.
potatoes       2.6         16.         1650.
onions         2.9         16.         1020.
chilli         10.         42.          525.

$eject
 table ecnsdwtpr(c,ps,*)  prices of seed(rs per kg) and water(rs per acre )

           87-88.seed 87-88.water 87-88.miscc
basmati     7.          16.          75.
irri        3.7         16.          75.
cotton      6.7         14.4        500.
rab-fod    36.           6.4
gram        5.1          8.
maize       7.           9.6
mus+rap     5.9         12.
kha-fod     2.25         6.4
sc-gur       .43        33.
sc-mill      .43        33.
wheat       3.6         10.
orchard                 64.
potatoes    2.6         12.        1650.
onions      2.9         12.        1020.
chilli     10           32.         525.

*        wages, fertilizer and other input prices
set p1 /  nitrogen, phosphate , protein, twinvt, trinvt,twopc,tropc/
    p11 /financial, economic /

 table pri1(ps,p11,p1)  fertilizer tubewell tractor and protein prices

*       fertilizer and protein prices are in rs/kg, twinvt and trinvt are
*       annualized cost for a tubewell and tractor (rupees)
*       twopc and tropc are cost of tubewell water(rs/acre foot) and
*       cost of tractor (rs/tractor hour)
*
                 nitrogen phosphate  protein twinvt trinvt twopc tropc
 87-88.financial  5.8       7.0        9.    10000   25000   225   45.
 87-88.economic   8.3      10.9        9.    10000   25000   170   50.

     table wageps(ps,p11,m)    wage rates rs per man hour

                jan   feb   mar  apr may  jun  jul  aug  sep  oct nov dec
87-88.financial 3.    3.    3.   6.  6.   3.   3.   3.   3.   6.  6.  3.
87-88.economic  2.75  2.75  2.75 5.  5.   2.75 2.75 2.75 2.75 5.  5.  2.75


* miscellaneous parametres
  scalars    lstd           standard labor limit ( hours per month )                    /  200     /
             trcap          tractor capacity in tractor hours per month                 /  250     /
             twcap          nameplate capacity of the private tubewell( af per month)   /   59.5041/
             ntwucap        effective capacity of new tubewells(af per month)
             twefac         factor to convert wc losses to from private tubewell losses /    0.5   /
             labfac         factor to convert wage to the reservation wage              /    0.5   /
               ;
   parameter twutil(*) effective capacity of tubewells(proportion of name plate capacity)
              / existing   .60
                new        .75 /;
$stitle demand and consumption data
table totprod(z,cq ) total production  1988 (000's tons)

*   sc-mill is total cane production (including that used for gur)
*   sc-gur is refined portion of total cane - same as consumption

       wheat  basmati   irri maize sc-mill sc-gur cotton mus+rap   gram  chilli potatoes onions

nwfp   161.8      0.0    0.0  42.1  2473.7  49.5     0.0    14.8   22.4     4.8    29.7   25.2
pcw   4636.8    182.8  175.7  23.0 10885.2 217.7  3318.4   103.8  132.4    38.7   199.5  191.5
pmw    747.1      0.0    0.0   0.0  1364.8  27.3   153.6     7.5   50.6     7.1    36.8   35.1
prw   1720.1    636.0  279.6   5.6  1065.7  21.3    46.9    16.6   43.0    15.2    78.5   75.0
psw   2327.3    141.2    0.0  27.2  7148.6 143.0   397.0    14.6   36.4    21.6   110.8  106.2
scwn   830.1      0.0   18.1   0.0  3030.7  60.6   485.5    17.4   16.7     5.9    33.7   27.8
scws   847.7      0.0  157.3   0.0  4947.9  99.0   437.4     6.8    6.2     7.0    39.9   32.4
srwn   414.2      0.0  920.3   0.0   289.5   5.8     6.1    27.8   23.1     7.2    41.2   33.6
srws   112.7      0.0  257.3   0.0  1832.0  36.6    18.9     3.4   10.1     2.2    12.3   10.1

  table farmcons(z,cq ) on-farm consumption 1988 (000's tons)

       wheat  basmati   irri  maize  sc-gur mus+rap gram chilli potatoes onions

nwfp   110.0                   18.9    49.5    6.7  10.1   2.2    13.4   11.3
pcw   3013.9     27.4   21.1   10.3   217.7   46.7  59.6  17.4    89.8   86.2
pmw    485.6                           27.3    3.4  22.8   3.2    16.6   15.8
prw   1118.1     95.4   33.6           21.3    7.5  19.4   6.8    35.3   33.7
psw   1512.7     21.2          12.3   143.0    6.6  16.4   9.7    49.9   47.8
scwn   697.3            11.8           60.6    7.8   7.5   2.7    15.2   12.5
scws   712.1           102.2           99.0    3.1   2.8   3.1    18.0   14.6
srwn   347.9           171.1            5.8   12.5  10.4   2.9    16.8   13.5
srws    94.6            27.1           36.6    1.5   4.5   1.3     7.4    6.1

 table demand(*,cq ) market demand by zone (000 tons or million liters)

*   basmati & irri adjusted for expected exports
*   sc-mill adjusted for expected imports of refined sugar

      wheat basmati   irri maize sc-mill cotton mus+rap  gram chilli potatoes onions cow-milk buff-milk

nwfp   51.8                 23.1  2337.6    0.0     8.1  12.3    2.6    16.3   13.9   167       163
pcw  1622.9   155.4   49.3       10286.5 3318.4    57.1  72.8   21.3   109.7  105.3   190       267
pmw   261.5                       1289.8  153.6     4.1  27.8    3.9    20.2   19.3   709      2298
prw   602.0   290.6   32.2        1007.1            9.1  23.7    8.3    43.2   41.2   260      1591
psw   814.5   120.2         15.0  6755.5  397.0     8.0  20.0   11.9    61.0   58.4   180      1142
scwn  132.8            6.3        2864.0  485.5     9.6   9.2    3.3    18.5   15.3   201       374
scws   66.3           55.0        4675.8  437.4     3.8   3.4    3.8    22.0   17.8   110       216
srwn  18.03          171.2         273.5           15.3  12.7    3.9    22.7   18.5   187       469
srws   28.0           27.1        1731.3            1.8   5.5    1.2     6.7    5.5   206       366

   ;

* adujust the demands for cow and buffalo population in the irrigated area
*&&z    zone3xxxxx(z,cq, "demand") = demand(z,cq ) ;
 scalar
  cowf adjustment factor for cows population in the irrigated areas/.5/
  buff adjustment factor for buffloes pop.   in the irrigated areas/.8/
     ;


   demand(z,"cow-milk")  = demand(z,"cow-milk" )*cowf;
   demand(z,"buff-milk") = demand(z,"buff-milk")*buff;


  parameter elast(cq ) elasticity of demand for crop and livestock comodities
    /
 (basmati,irri)            -.4
 cotton                    -.9
 gram                      -.22
 maize                     -.4
 mus+rap                   -.35
 sc-mill                   -.7
 wheat                     -.23
 (potatoes,onions,chilli)  -.7
 (cow-milk,buff-milk)      -.44
      /

    growthrd(cq )  growth rate of reference demand (percent)/
           basmati   4.5,  irri       4.0,  cotton    5.0,
           gram      4.7,  wheat      4.1,  mus+rap   3.4,
           sc-mill   5.5,
           cow-milk  6.3,  buff-milk  6.3,  meat      6.6 /
       ;

    table consratio(z,g) proportion of consumption by growundwater type

                   fresh
        nwfp       1
        pmw         .85
        pcw         .9
        psw         .83
        prw        1
        scwn        .55
        srwn        .33
        scws        .34
        srws       0
         ;
    parameter natexp(cq) national exports (000 tons)
              /basmati  250, irri 1100, cotton 700, onions 10 /
              explimit(z,cq) export limits by zone
    scalar    explimitgr     growth rate of export limits(percent) /5 /
               ;
    table exppv(pv,cq) provincial exports as proportion of national

                       basmati     irri   cotton   onions
    nwfp                                             0.14
    punjab                1.00     0.29     0.84     0.33
    sind                           0.71     0.16     0.53

    table expzo(z,cq)  zonal exports as proportion of provincial

                       basmati     irri   cotton   onions
    nwfp                                             1
    pcw                            0.33     1.00     0.52
    prw                   1.00     0.67              0.22
    pmw                                              0.09
    psw                                              0.17

    scwn                                    0.74     0.41
    srwn                           0.74              0.30
    scws                                    0.26     0.22
    srws                           0.26              0.08
      ;

   growthrd(cq)$(growthrd(cq) eq 0) = 3.0 ;
   consratio(z,"saline") = 1 - consratio(z,"fresh") ;
   explimit(z,cq) = natexp(cq)*
                         sum(pv$pvz(pv,z), exppv(pv,cq))*expzo(z,cq) ;
   explimit(z,cq) = sum(is$isr(is),(1+explimitgr/100)**
                          (ord(is)+1979-baseyear)*explimit(z,cq)) ;

*&&z    zone3xxxxx(z,cq,"totprod")  = totprod(z,cq)  ;
*&&z    zone3xxxxx(z,cq,"farmcons") = farmcons(z,cq) ;
*&&z    zone3xxxxx(z,cq,"elast")    = elast(cq)      ;
*&&z    zone3xxxxx(z,cq,"growthrd") = growthrd(cq)   ;
*&&z    zone3xxxxx(z,cq,"explimit") = explimit(z,cq) ;


$stitle data transformation to zones
   set sr1(dc) /cca, ccap /
   alias (g,g1)
  parameter
      zwt(z,cnl,sa)     weighting factor to map rain evap and efficiencies to zones
      eqevapz(z,m)      evaporation from the equaifer by acz (feet)
      subirrz(z,m)      subirrigation by acz                 (feet)
      efrz(z,m)         effective rain  by acz               (feet)
      resource(z,g,r1)  endowments by acz and groundwater quality
      cneff(cnl)        canal efficiency from canal head to the watercourse head
      wceff(cnl,m)      watercourse command delivery efficiency
      tweff(cnl,m)      delivery efficiency from private tubewell to the root zone
      cneffz(z)         weighted canal delivery efficiency from canal head to watercourse head
      tweffz(z,m)       weighted private tubewell delivery efficiency by zone
      wceffz(z,m)       weighted water course command delivery efficiency by zone
      fleffz(z)         weighted field efficiency by zone
      canalwz(z,g,m)    canal water availablility at the canal head(maf)
      canalwrtz(z,g,m)  canal water availablility at the root zone(maf)
      gwtsa(cnl,sa,m)   government tubewell pumpage by canal and subarea(kaf)
      gwt1(z,g,m)       public tubewell pumpage at the root zone (kaf)
      ratiofs(z,g)      fresh and saline cca as a proportion off total
           ;
  loop(isr,
  resource(z,g,sr1)            = sum((cnl,sa)$(zsa(z,cnl,sa)$gwfg(cnl,sa,g)), comdef(isr,sr1,cnl) *subdef(sa,cnl) );
  zwt(z,cnl,sa)$zsa(z,cnl,sa)  =  comdef(isr,"cca",cnl)*subdef(sa,cnl)/sum(g,resource(z,g,"cca") );
  cneff(cnl)       = comdef(isr,"ceff",cnl)   ;
  wceff(cnl,m)     = sum( wce$wcem(wce,m), comdef(isr,wce,cnl)  );
  tweff(cnl,m)     = (1- (1-wceff(cnl,m)/comdef(isr,"flde",cnl))*twefac)*comdef(isr,"flde",cnl);
  cneffz(z)        = sum( (cnl,sa),   comdef(isr,"ceff",cnl)*zwt(z,cnl,sa)    )            ;
  fleffz(z)        = sum( (cnl,sa),   comdef(isr,"flde",cnl)*zwt(z,cnl,sa)    )            ;
  canalwrtz(z,g,m) = sum( (cnl,sa)$(zsa(z,cnl,sa)$gwfg(cnl,sa,g)), divpost(cnl,m)*subdef(sa,cnl)*
                                                                   comdef(isr,"ceff",cnl)*wceff(cnl,m) );
  gwtsa(cnl,sa,m)$carea(cnl,"fresh")=subdef(sa,cnl)$gwf(cnl,sa)*comdef(isr,"cca",cnl)/carea(cnl,"fresh")*gwt(cnl,m)
   );

  ratiofs(z,g)     = resource(z,g,"cca")/ sum(g1, resource(z,g1,"cca") ) ;
  canalwz(z,g,m)   = sum( (cnl,sa)$(zsa(z,cnl,sa)$gwfg(cnl,sa,g)), divpost(cnl,m)*subdef(sa,cnl) );
  eqevapz(z,m)     = sum( (cnl,sa),      eqevap(cnl,m)*zwt(z,cnl,sa)   )       ;
  subirrz(z,m)     = sum( (cnl,sa),      subirr(cnl,m)*zwt(z,cnl,sa)   )       ;
  efrz(z,m)        = sum( (cnl,sa),      efr(cnl,m)*zwt(z,cnl,sa)      )       ;
  tweffz(z,m)      = sum( (cnl,sa),      tweff(cnl,m)*zwt(z,cnl,sa)     )      ;
  wceffz(z,m)      = sum( (cnl,sa),      wceff(cnl,m)*zwt(z,cnl,sa)     )      ;
  gwt1(z,"fresh",m) = sum( (cnl,sa)$zsa(z,cnl,sa), gwtsa(cnl,sa,m)*wceff(cnl,m)  ) ;

 display   wceff, tweff, gwtsa, gwt1 ;
$eject
$ontext
*- report water components
    set crt /canal, canal-rt,govt-tw, gwt-rt  /
        ras /rain, subirr, total /
 parameter rep5(z,*,*)
           rep51
           rep6
            ;
     rep5(z,"rain",m)      = efrz(z,m)*12    ;
     rep5(z,"subirr",m)    = subirrz(z,m)*12 ; rep5(z,"total",m) = sum(ras, rep5(z,ras,m) ) ;
     rep5(z,ras,sea1)      = sum(m$sea1m(sea1,m), rep5(z,ras,m)) ;

     rep51(z,"cnldeleff",m)= cneffz(z)*wceffz(z,m) ;
     rep51(z,"wcdeleff",m) = wceffz(z,m) ;
     rep51(z,"twdeleff",m) = tweffz(z,m) ;

     rep6(z,g,"canal",m)          = canalwz(z,g,m)             ;
     rep6(z,"total","canal",m)    = sum(g,canalwz(z,g,m) )     ;
     rep6(z,g,"canal-rt",m)       = canalwrtz(z,g,m)           ;
     rep6(z,"total","canal-rt",m) = sum(g, canalwrtz(z,g,m))   ;
     rep6(z,"fresh","govt-tw",m)  = sum( (cnl,sa)$zsa(z,cnl,sa), gwtsa(cnl,sa,m) )/1000;
     rep6(z,"total","govt-tw",m)  = rep6(z,"fresh","govt-tw",m);
     rep6(z,g,"gwt-rt",m)         = gwt1(z,g,m)/1000;
     rep6(z,"total","gwt-rt",m)   = gwt1(z,"fresh",m)/1000;
     rep6(z,t1,crt,sea1)          = sum(m$sea1m(sea1,m), rep6(z,t1,crt,m));
     display zwt  ;
     display "rain       effective rain(inches) "
             "subirr     water available to the crops by cappilary action from the gw(inches) "
             "total      total water available from rain and subirr(inches) ", rep5 ;
     display "cnldeleff  weighted canal delivery efficiency from canal head to the root zone "
             "wcdeleff   weighted watercourse command  eff. from canal head to the root zone "
             "twdeleff   weighted perivate tubewell eff. from canal head to the root zone ", rep51 ;
     display "canal      canal water available at the canal head  (maf)    "
             "cannal-rt  canal water available at the root zone (maf)      "
             "govt-tw    government tubewell pumpage at the watercourse head(maf)"
             "gwt-rt     government tubewell pumpage at the root zone(maf)", rep6 ;
 display eqevapz ;
$offtext
*--
$stitle resource stocks
 set ftt(r1)/farmpop,farmhh,tractors,tubewells,bullocks,cows,buffalos/;

       table res88(r1 ,z) available resources 1988

*  farmpop -- thosand workers, farmhh -- thousand households
*  bullocks, cows and buffalos are in thousands
*  tubewells are 1 cfs tubewells

            nwfp     pmw     pcw     psw     prw    scwn    srwn   scws    srws

farmpop      371     478    2589    1441    1016     484     589    567     176
farmhh       174     226    1123     586     513     211     253    275      88
tractors    7465   12044   86817   45777   44564    8177    4154   8839    2702
tubewells   2638   12265  109658   37952   49897    9573    4863  10348    3163
bullocks     151     332    1978     952     368     399     486    444     143
cows         505     736    2549    1055     590     442     570    630     256
buffalos     319     487    4146    2911    2279     584     587    557     249


 table croparea(z,c) cropped area 1988 (000's acres)

*     sc-mill area includes gur area

      wheat basmati irri maize sc-mill cotton mus+rap gram chilli potatoes onions orchard rab-fod kha-fod
nwfp    224       0    0    78     143             52  132      9      8      6    27       47       14
pcw    6324     396  277    41     722   3820     275  746     55     51     40   300     1174     1332
pmw    1168      25   29           100    270      27  389     14     13     10    93      158      328
prw    2141    1189  272     9      79     78      42  233     21     19     15    56      506      374
psw    2688     391   42    51     495    614      41  218     33     30     24   242      675      698
scwn   1004       0   19           169    743      63   58     10      9      7    49      131       72
scws    921       0  196           238    758      25   22     12     11      8   103      137       59
srwn    590       0  951            24     13     124  100     15     14     11     7      117        4
srws    101       0  325           111     41      13   38      4      4      3    41       30        7

    table growthres(r1 ,z) growth rate of farm population tractors and tubwells (percent)

               (nwfp,pmw,pcw,psw,prw)   (scwn,scws,srwn,srws)
  farmpop
  farmhh
  tractors          3                      3
  tubewells         1
  bullocks          2.35                   2.35
  cows              1.21                   1.21
  buffalos          1.21                   1.21

  parameter
    orcharea(z)   area under orchards by zone (thousand acres )
    orchgrowth(z) growth rate of orchard area
                /(nwfp,pmw,pcw,psw,prw) 6.43
                 (scwn,scws,srwn,srws)  2.591/
                        ;
  demand(z,cq) = demand(z,cq)* sum(is$isr(is),
                     (1+growthrd(cq)/100)**(ord(is)+1979-baseyear) );
  orcharea(z)  = croparea(z,"orchard")* sum(is$isr(is),
                     (1+orchgrowth(z)/100)**(ord(is)+1979-baseyear) );
  resource(z,g,ftt) = res88(ftt,z)*ratiofs(z,g)*sum(is$isr(is),
                    (1+growthres(ftt,z)/100)**(ord(is)+1979-baseyear) );

 ntwucap  = twutil("new")*twcap ;
 resource(z,"fresh","tubewells")$resource(z,"fresh","cca") =
                                   sum(g, resource(z,g,"tubewells"));
 resource(z,"saline","tubewells")= 0 ;
 resource(z,"fresh","twc")$resource(z,"fresh","cca") =
                             resource(z,"fresh","tubewells")*twcap*
                             twutil("existing") /1000  ;

  parameter  scmillcap(z)  sugarcane mill capacity (thousand tonns per year )
            /nwfp  6937,
             pmw   5750,    pcw  8550,    psw  12000,    prw  3600
             scwn  6000,    scws 8370,    srwn 3800,     srws 3730 /;
 option totprod:1, farmcons:1, demand:1, consratio:2 ;
 display resource, totprod, farmcons,consratio, demand,explimit ;
 set cnl1(cnl)      canals excluding nwfp canals      ;
     cnl1(cnl) = yes;   cnl1(cnl)$pvcnl("nwfp",cnl) = no ;
 parameter  postt    average canal diversions by season
            protarb  diversions as proportion of total(punjab and sind) ost tarbela
            ;
   postt(cnl,sea) = sum(m$seam(sea,m),   divpost(cnl,m));
   postt(pv2,sea)  = sum(cnl$pvcnl(pv2,cnl), postt(cnl,sea));
   protarb(cnl1,sea) = 0.999 * postt(cnl1,sea)/ ( postt("punjab",sea) + postt("sind",sea)) ;
   protarb(pv2,sea)   = 0.999 * postt(pv2,sea)  / ( postt("punjab",sea) + postt("sind",sea)) ;
   option postt:5, protarb:5 ;
 display cnl1, postt, protarb ;
$stitle model setup
 sets
  psr(ps)  price scenario for the model(financial prices) /87-88  /
  psr1(ps) price scenario for report(economic prices)     /87-88 /
  z1(z)    zone selection for this run
         / nwfp,pmw,pcw,psw,prw,scwn,scws,srwn,srws /
*
*  total comodities are 18,
*        with endogenous prices= 13, fixed prices=2, fodder=2
*        consumption only = 1

  cn(cq)     comodities  endogenous prices
            /basmati,   irri,       cotton,
             gram,      maize,      mus+rap,
             sc-mill,   wheat,      potatoes,
             onions,    chilli,
             cow-milk,  buff-milk  /
  ccn(cq)    crop comodities with endogenous prices
            /basmati,   irri,       cotton,
             gram,      maize,      mus+rap,
             sc-mill,   wheat,      potatoes,
             onions,    chilli /
  qn(cq)     livestock comodities endogenous prices
            /cow-milk, buff-milk/
  ncn(cq)    crops with fixed prices excluding fodder
            /orchard, meat /
  ce(cq)     exportable comodities
            / basmati, irri, cotton,onions,wheat /
  cm(cq)     comodities which could be imported/sc-mill/

  ex(z,g)          to check fresh or saline area within a zone
  techc(z,cq)      comodities by zones
         ;
  ex(z1,g)$resource(z1,g,"cca") = yes ;
  display  cq,cn, ccn, qn, ncn,ce,cm,ex;
  display  fert;

        table tec(c,t,s,w,z)  crop technology disabled for 1988 run
                                                                     scwn  srwn  scws  srws
  wheat. (bullock,semi-mech).la-plant.(heavy,january)                  1     1     1     1
  wheat. (bullock,semi-mech).(qk-harv,standard). heavy                       1     1     1
  wheat. (bullock,semi-mech).(qk-harv,standard).january                      1           1
  wheat. (bullock,semi-mech).(qk-harv,standard,la-plant).light                           1
     ;
 scalars big big number used for artifical production /4000/
         pawat  big number for artificial water       /999999/
         pafod  big number for artifical fodder       /1000/
 parameter
    divnwfp(m)  monthy diversion to the nwfp zone  (maf)
              /jan .158, feb .187, mar .270, apr .184,
               may .232, jun .325, jul .333, aug .333,
               sep .332, oct .200, nov .063, dec .133 /

    rval(n)  value of water stored in the reservoirs
                 /tarbela-r 1, kalabagh-r .8, mangla-r .6
                  chasma-r .2, chotiari-r .2, a-sea    .1 /
   fsalep(cq)   financial sale price for crop and livestock comodities (rs per kg or per liter)
   pp           financial purchase price of protein (rs per kgs)
   misc(*)      financial miscellenious prices
   seedp(c)     financial seed price                (rs per kgs)
   wage(m)      financial wage rates                (rs per man-hour)
   miscct(c)    financial water charges and miscillenious costs(rs per acre)

   esalep(cq)   economic sale price for crop and livestock comodities  (rs per kg or per liter)
   epp          economic price of protein concentrate    (rs per kgs)
   emisc(*)     economic miscellenious prices
   eseedp(c)    economic seed price
   ewage(m)     economic wage rate     (rs per man-hour)
   emiscct(c)   economic water charges and miscillenious costs(rs per acre)

   importp(cq)  import prices for the  scenario
   exportp(cq)  export prices for the  scenario

   wnr(c,z,t,s,w,m) water requirements net of rain
               ;

  loop (psr,
     fsalep(cq)    =  prices(psr,cq,"financial")      ;
     pp            =  pri1(psr,"financial","protein");

     misc(p1)      =  pri1(psr,"financial",p1)    ;
     seedp(c)      =  finsdwtpr(c,psr,"seed")     ;
     wage(m)       =  wageps(psr,"financial",m)   ;
     miscct(c)     =  finsdwtpr(c,psr,"water")+finsdwtpr(c,psr,"miscc");
     importp(cq)   =  prices(psr,cq,"import") ;
     exportp(cq)   =  prices(psr,cq,"export")
  );
  loop (psr1,
     esalep(cq)    =  prices(psr1,cq ,"economic")       ;
     epp           =  pri1(psr1,"economic","protein") ;

     emisc(p1)     =  pri1(psr1,"economic",p1)        ;
     eseedp(c)     =  ecnsdwtpr(c,psr1,"seed")        ;
     ewage(m)      =  wageps(psr1,"economic",m)       ;
     emiscct(c)    =  ecnsdwtpr(c,psr1,"water")+ecnsdwtpr(c,psr1,"miscc")
  );
$ontext
* economic prices are set to finacial
     esalep(cq) =  fsalep(cq);
     emisc(p1)  =  misc(p1)   ;
     eseedp(c)  =  seedp(c)   ;
     epp        =  pp         ;
     ewage(m)   =  wage(m)    ;
$offtext
* if we know that water is positive runs 10 times as fast
*
*  wnr(c,z1,t,s,w,m)$water(c,z1,t,s,w,m) = max(0, (water(c,z1,t,s,w,m)-efrz(z1,m))  ) ;
   wnr(c,z1,t,s,w,m)    = max(0, (water(c,z1,t,s,w,m)-efrz(z1,m))  ) ;
*- set coton in prw to no
   tec("cotton",t,s,w,"prw") = 1   ;
   tec("maize",t,s,w,"prw")  = 1   ;

   tech(z1,c,t,s,w)$tec(c,t,s,w,z1) = no ;
   techc(z1,c)$sum( (t,s,w), tech(z1,c,t,s,w) ) = yes   ;
   techc(z1,cf)=no;  techc(z1,q) = yes ;
*  display tech;
   display   techc ;

*  note: inflows from ravi and sutlej river are set to zero here.
*  inflow("ravi",m) = 0 ;  inflow("sutlej",m) = 0 ;

 scalars
    tolcnl  allowed deviation from proportional allocation by canal     /0.0 /
    tolpr   allowed deviation from proportional allocation by province / 0.0 /
    tolnwfp nwfp diversion tolerance / 0 /
          ;
*-
  parameter beta(cq,z1)  gradient comodities demand curve
            alpha(cq,z1) demand curve intecept
            ;
  scalar    betaf         beta factor /.5 /  ;
  beta(cn,z1 )$demand(z1 ,cn) = fsalep(cn) / demand(z1 ,cn) /  elast(cn) ;
  alpha(cn,z1 )               = fsalep(cn)  - beta(cn,z1 )*demand(z1 ,cn)   ;

*-- linearization of the demand function.

  set   p   grid points for linearization / 1*20 /
 parameter
      pmax(cq,z1)     maximum price for segments
      pmin(cq,z1)     minimum price for segments
      qmax(cq,z1)     max national consumption
      qmin(cq,z1)     min national consumption
      incr(cq,z1)     increment
      ws(cq,z1,p)     welfare segments   (million rupees)
      rs(cq,z1,p)     revenue definition (million rupees)
      qs(cq,z1,p)     quantity definition(thousand tons or million liters)
      endpr(cq,z1,p)  price (rupees per kgs or liter )
              ;
 pmin(cn,z1) = 0.5*fsalep(cn)  ;
 pmax(cn,z1) = min(alpha(cn,z1), 2*fsalep(cn) ) ;
* pmax(ce,z1) = min(alpha(ce,z1 ),2*fsalep(ce) ) ;
* pmin(ce,z1) =     fsalep(ce)  ;
 qmin(cn,z1)$beta(cn,z1 ) = (pmax(cn,z1 )-alpha(cn,z1 ))/beta(cn,z1 );
 qmax(cn,z1)$beta(cn,z1 ) = (pmin(cn,z1 )-alpha(cn,z1 ))/beta(cn,z1 );
 incr(cn,z1) = (qmax(cn,z1 )-qmin(cn,z1 ))/(card(p)-1);

 qs(cn,z1,p) = qmin(cn,z1) + incr(cn,z1)*(ord(p)-1);
 ws(cn,z1,p) = alpha(cn,z1)*qs(cn,z1,p) + betaf*beta(cn,z1)*sqr(qs(cn,z1,p)) ;
 rs(cn,z1,p) = alpha(cn,z1)*qs(cn,z1,p) + beta(cn,z1)*sqr(qs(cn,z1,p));
 endpr(cn,z1,p)= alpha(cn,z1)           + beta(cn,z1)*qs(cn,z1,p);
 display pmax, pmin, qmax, qmin,incr, qs, ws, rs, endpr ;
*--
 display alpha , beta ;
$stitle    equations and variables
 variables   cps   consumer plus producers surplus (million rupees)
 positive variables
      acost(z,g)          farm cost in           (million rupees)
      ppc(z,g,sea)        purchases of protein concentrates                 (thousand metric tons)
      x(z,g,c,t,s,w)      cropped area by technology                        (thousand acres)
      animal(z,g,a)       production of livestock type a                    (thousands)
      prodt(z,g,cq)       production (crop commodities 000 metric tons livestock comm mill. kgs or liters)
      proda(z,g,cq )      artificial supply
      import(z ,cq )      import of comodities (crop comm. 000 m. tons livestock mill. kgs or liters)
      export(z ,cq )      export of comodities                               (000 metric tonns)
      consump(z,g,cq )    on farm consumption                                (000 metric tonns)
      familyl(z,g,m)      family labor used                                  (million man hours)
      hiredl(z,g,m)       hired labor used                                   (million man hours)
      itw(z)              investment in increased private tubewell capacity  (kaf per month)
      tw(z,m)             private tubewell water used  by month m            (kaf)
      itr(z,g)            investment in increased tractor capacity           (000 tractor-hrs per month)
      ts(z,g,m)           private tractor services use by month              (thousand hrs)
      f(n,n1,m)           flow to node n from node n1                        (maf)
      rcont(n,m)          end of the month resrvoir contents                 (maf)
      canaldiv(cnl,m)     canal diversion at the canal head                  (maf)
      cnldivsea(cnl,sea)  canal diversion by season                          (maf)
      prsea(pv2,sea)      canal diversion by province (sind and punjab)      (maf)
      tcdivsea(sea)       total canal diversion in sind and punjab by season (maf)
      wdivrz(z1,g,m)      surface water diversion at the root zone           (kaf)
      slkland(z,g,m)      slack land                                         (thousand acres)
      slkwater(z,g,m)     slack water at the root zone                       (kaf)
      artfod(z1,g,sea)    artificial fodder supply  equaivalent of rab-fod   (000 tonns)
      artwater(z,g,m)     water from imaginary source at the root zone       (kaf)
      artwaternd(n,m)     water from imaginary source at nodes               (maf)
      nat(cq,z,p)         provincial demand linearized
      natn(cq,z)          provincial demand non-linear

 equations
    objz   objective function for the zone model linear version      (million rupees)
    objzn  objective function for the zone model non-linear version  (million rupees)
    objn   objective function for the indus model linear version     (million rupees)
    objnn  objective function for the indus model non-linear version (million rupees)
    cost(z,g)           annual farm cost                             (million rupees)
    conv(cq,z)          convex combination for aggregate consumption
    demnat(cq,z)        provincial demand balance linear (000 tons or million liters)
    demnatn(cq,z)       zonal demand balance non-linear  (000 tons or million liters)
    ccombal(z,g,c)      commodity balances for crops     (000 tons)
    qcombal(z,g,q)      livestock comodity balances      (000 tons or m liters)
    consbal(z,g,cq)     consumption balance              (000 tons or m liters)
    laborc(z,g,m)       monthly labor constraint         (million man hours)
    fodder(z,g,sea)     seasonal maintenance of fodder supplies     (000 metric tons)
    protein(z,g,sea)    protein requirements of livestock by season (000 metric tons)
    grnfdr(z,g,sea)     green fodder requirements                   (000 metric tons)
    bdraft(z,g,m)       bullock draft power constraint              (million bullock hours)
    brepco(z,g)         bullock reproduction constraint
    bullockc(z1)        bullock population constraint               (000 bullocks)
    tdraft(z,g,m)       tractor draft power balance                 (000 tractor hours)
    trcapc(z,m)         tractor capacity constraint                 (000 tractor hours)
    twcapc(z,m)         tubewell capacity constraint                (kaf)
    landc(z,g,m)        land constraint                             (000 acres)
    orchareac(z)        orchard area constraint                     (000 acres)
    scmillc(z)          sugar cane to mill constraint               (000 acres)
    waterbaln(z,g,m)    water balance at the root zone              (kaf)
    watalcz(z,g,m)      surface water by zone                       (kaf)
    subirrc(z,g,m)      subirrigation constraint                    (kaf)
    nbal(n,m)           water balance at a node                     (maf)
    watalcsea(cnl,sea)  water allocations by season                 (maf)
    divsea(sea)         total canal diversions in sind and punjab   (maf)
    divcnlsea(cnl,sea)  canal diversion by season                   (maf)
    watalcpro(pv,sea)   water allocation by province                (maf)
    prseaw(pv,sea)      diversions by province and season           (maf)
    nwfpalc(m)          water allocations to the nwfp acz           (maf)
           ;
$eject
objz..
 cps  =e=
   sum(z1, sum(g$ex(z1,g), sum(ncn, fsalep(ncn)*prodt(z1,g,ncn)   )
                         - acost(z1,g) - sum(sea, artfod(z1,g,sea))*pafod
                         - sum(m, artwater(z1,g,m))*pawat
                         - sum(cq$techc(z1,cq), proda(z1,g,cq )*big)) )
   - sum(z1,  sum(cm$techc(z1,cm),      import(z1,cm)*importp(cm))    )
   + sum(z1,  sum(ce$techc(z1,ce),      export(z1,ce)*exportp(ce))    )
   + sum(z1,  sum( (cn,p)$techc(z1,cn), nat(cn,z1,p)*ws(cn,z1 ,p))    ) ;

objzn..
 cps  =e=
   sum(z1, sum(g$ex(z1,g), sum(ncn, fsalep(ncn)*prodt(z1,g,ncn)   )
                         - acost(z1,g) - sum(sea, artfod(z1,g,sea))*pafod
                         - sum(m, artwater(z1,g,m))*pawat
                         - sum(cq$techc(z1,cq), proda(z1,g,cq )*big)) )
   - sum(z1,  sum(cm$techc(z1,cm),      import(z1,cm)*importp(cm))    )
   + sum(z1,  sum(ce$techc(z1,ce),      export(z1,ce)*exportp(ce))    )
   + sum(z1,  sum(cn$techc(z1,cn),  alpha(cn,z1)*natn(cn,z1)
                            + betaf*beta(cn,z1)*sqr(natn(cn,z1))  )   ) ;

objn..
 cps  =e=
   sum(z1, sum(g$ex(z1,g), sum(ncn, fsalep(ncn)*prodt(z1,g,ncn)   )
                         - acost(z1,g) - sum(sea, artfod(z1,g,sea))*pafod
                         - sum(m, artwater(z1,g,m))*pawat
                         - sum(cq$techc(z1,cq), proda(z1,g,cq )*big)) )
   - sum(z1,  sum(cm$techc(z1,cm),      import(z1,cm) *importp(cm))   )
   + sum(z1,  sum(ce$techc(z1,ce),      export(z1,ce) *exportp(ce))   )
   + sum(z1,  sum( (cn,p)$techc(z1,cn), nat(cn,z1,p)*ws(cn,z1,p))     )
   + sum((n,m), -artwaternd(n,m)*pawat + rval(n)*rcont(n,m)$rcap(n)
                               + rval("a-sea")*f("a-sea","kotri-b",m) ) ;

 objnn..
 cps  =e=
   sum(z1, sum(g$ex(z1,g), sum(ncn, fsalep(ncn)*prodt(z1,g,ncn)   )
                         - acost(z1,g) - sum(sea, artfod(z1,g,sea))*pafod
                         - sum(m, artwater(z1,g,m))*pawat
                         - sum(cq$techc(z1,cq), proda(z1,g,cq )*big)) )
   - sum(z1,  sum(cm$techc(z1,cm),      import(z1,cm)*importp(cm))    )
   + sum(z1,  sum(ce$techc(z1,ce),      export(z1,ce)*exportp(ce))    )
   + sum(z1,  sum(cn$techc(z1,cn),  alpha(cn,z1)*natn(cn,z1)
                       + betaf*beta(cn,z1)*sqr(natn(cn,z1))  )        )
   + sum((n,m), -artwaternd(n,m)*pawat + rval(n)*rcont(n,m)$rcap(n)
                               + rval("a-sea")*f("a-sea","kotri-b",m) ) ;


  cost(z1,g)$ex(z1,g)..   acost(z1,g) =e= (sum((c,t,s,w)$tech(z1,c,t,s,w), (sum(p2, fert(p2,c,z1)*misc(p2))+
                                                       miscct(c)+seedp(c)*sylds(c,z1,t,s,w,"seed") )*x(z1,g,c,t,s,w) )+
                                           sum(m, misc("twopc")*tw(z1,m)$gf(g) + misc("tropc")*ts(z1,g,m) ) +
                                           misc("twinvt")*itw(z1)$gf(g)  + misc("trinvt")*itr(z1,g) +
                                           sum(a, iolive(a,z1,"fix-cost")*animal(z1,g,a))  )/1000+
                                           sum(sea, pp*ppc(z1,g,sea)) +
                                           sum(m, (familyl(z1,g,m)*labfac + hiredl(z1,g,m))*wage(m) ) ;
  conv(cn,z1)$techc(z1,cn)..     sum(p, nat(cn,z1,p)) =l= 1;

  demnat(cq,z1 )$techc(z1,cq)..
     sum(g$ex(z1,g),
           prodt(z1,g,cq) - consump(z1,g,cq)$cc(cq)+proda(z1,g,cq) )-
     export(z1,cq)$ce(cq) + import(z1,cq)$cm(cq)
                     =g= sum(p, nat(cq,z1,p)*qs(cq,z1 ,p))$cn(cq);

  demnatn(cq,z1 )$techc(z1,cq)..
     sum(g$ex(z1,g),
           prodt(z1,g,cq) - consump(z1,g,cq)$cc(cq)+proda(z1,g,cq) )-
     export(z1,cq)$ce(cq) + import(z1,cq)$cm(cq)
                             =g= natn(cq,z1)$cn(cq) ;

  ccombal(z1,g,c)$(cnf(c)$ex(z1,g)).. sum((t,s,w)$tech(z1,c,t,s,w), yield(c,t,s,w,z1)*x(z1,g,c,t,s,w))=e= prodt(z1,g,c);

  qcombal(z1,g,q)$ex(z1,g)..          sum(a, iolive(a,z1,q)*animal(z1,g,a) ) /1000                   =e= prodt(z1,g,q);

  consbal(z1,g,cc)$(ex(z1,g)$techc(z1,cc) )..
           prodt(z1,g,cc) + proda(z1,g,cc)  =g= consump(z1,g,cc) ;

  laborc(z1,g,m)$ex(z1,g)..    (sum( (c,t,s,w)$tech(z1,c,t,s,w), labor(c,z1,t,s,w,m)*x(z1,g,c,t,s,w) ) +
                                sum( a,  iolive(a,z1,"labor")*animal(z1,g,a)))/1000 =l= familyl(z1,g,m)+hiredl(z1,g,m);

  fodder(z1,g,sea)$ex(z1,g)..  sum(a, iolive(a,z1,"tdn")*animal(z1,g,a) )  =l= sum( (c,t,s,w)$tech(z1,c,t,s,w), (
                               yield(c,t,s,w,z1)*sylds(c,z1,t,s,w,"straw-yld")*sconv("tdn",sea,c) + weedy(z1,sea,c)*
                               sconv("tdn","rabi","rab-fod")    ) * x(z1,g,c,t,s,w)  )  +
                               sum(m$seam(sea,m), slkland(z1,g,m))*graz(z1,sea)*sconv("tdn","rabi","rab-fod") +
                               artfod(z1,g,sea)*sconv("tdn","rabi","rab-fod");

  protein(z1,g,sea)$ex(z1,g).. sum(a, iolive(a,z1,"dp")*animal(z1,g,a) ) =l= ppc(z1,g,sea)+
                               sum( (c,t,s,w)$tech(z1,c,t,s,w), (yield(c,t,s,w,z1)*sylds(c,z1,t,s,w,"straw-yld")*
                               sconv("dp",sea,c) + weedy(z1,sea,c)*sconv("dp","rabi","rab-fod") )*x(z1,g,c,t,s,w) )+
                               sum(m$seam(sea,m), slkland(z1,g,m))*graz(z1,sea)*sconv("dp","rabi","rab-fod") +
                               artfod(z1,g,sea)*sconv("dp","rabi","rab-fod");

  grnfdr(z1,g,sea)$ex(z1,g)..  gr*sum(a, iolive(a,z1,"tdn")*animal(z1,g,a) )  =l=  sum( (cf,t,s,w)$tech(z1,cf,t,s,w),
                               yield(cf,t,s,w,z1)*sconv("tdn",sea,cf)*x(z1,g,cf,t,s,w)  ) +
                               sum( (c,t,s,w)$tech(z1,c,t,s,w),
                               weedy(z1,sea,c)*sconv("tdn","rabi","rab-fod")*x(z1,g,c,t,s,w)  ) +
                               artfod(z1,g,sea)*sconv("tdn","rabi","rab-fod");

  bdraft(z1,g,m)$ex(z1,g)..    sum((c,t,s,w)$tech(z1,c,t,s,w), bullock(c,z1,t,s,w,m)*x(z1,g,c,t,s,w) )/1000 =l=
                                                                                 bp(m)*animal(z1,g,"bullock")/1000 ;

  brepco(z1,g)$ex(z1,g)..      animal(z1,g,"bullock")   =l=   repco*animal(z1,g,"cow")    ;

  bullockc(z1)..               sum(g$ex(z1,g), animal(z1,g,"bullock") )  =l=  res88("bullocks",z1)   ;

  tdraft(z1,g,m)$ex(z1,g)..    sum((c,t,s,w)$tech(z1,c,t,s,w), tractor(c,z1,t,s,w,m)*x(z1,g,c,t,s,w) ) =e= ts(z1,g,m) ;

  trcapc(z1,m)..   sum(g$ex(z1,g), ts(z1,g,m) ) =l= sum(g$ex(z1,g),
               (resource(z1,g,"tractors")/1000 +itr(z1,g))*trcap )  ;

  twcapc(z1,m)$ex(z1,"fresh")..

        tw(z1,m) =l=  resource(z1,"fresh","twc") + ntwucap*itw(z1) ;

  landc(z1,g,m)$ex(z1,g)..   sum((c,t,s,w)$tech(z1,c,t,s,w), land(c,z1,t,s,w,m)*x(z1,g,c,t,s,w)) + slkland(z1,g,m)=e=
                                                                                           resource(z1,g,"cca")*1000 ;
  orchareac(z1)..            sum( (g,t,s,w)$ex(z1,g), x(z1,g,"orchard",t,s,w)$tech(z1,"orchard",t,s,w))=l=orcharea(z1);

  scmillc(z1)..   sum(g$ex(z1,g),  prodt(z1,g,"sc-mill") ) =l= scmillcap(z1)   ;


  waterbaln(z1,g,m)$ex(z1,g)..
     sum((c,t,s,w)$tech(z1,c,t,s,w),
         max( (wnr(c,z1,t,s,w,m)-subirrz(z1,m)*land(c,z1,t,s,w,m)) , 0.0)*
              x(z1,g,c,t,s,w))+ slkwater(z1,g,m) =e=
              tweffz(z1,m)*tw(z1,m)$gf(g)  + gwt1(z1,g,m) +
              artwater(z1,g,m) + wdivrz(z1,g,m) ;

   watalcz(z1,g,m)$ex(z1,g)..  wdivrz(z1,g,m) =e= sum( (cnl,sa)$(zsa(z1,cnl,sa)$gwfg(cnl,sa,g) ),
                                                       cneff(cnl)*wceff(cnl,m)*canaldiv(cnl,m)*subdef(sa,cnl)*1000 ) ;

   divcnlsea(cnl,sea)..
      cnldivsea(cnl,sea) =e= sum(m$seam(sea,m), canaldiv(cnl,m));
   prseaw(pv2,sea)..
      prsea(pv2,sea)      =e= sum(cnl1$pvcnl(pv2,cnl1),
                                   cnldivsea(cnl1,sea)) ;
   divsea(sea)..
      tcdivsea(sea)      =e= sum(pv2,   prsea(pv2,sea) )          ;
   watalcsea(cnl1,sea)..
      protarb(cnl1,sea)*(1-tolcnl)*tcdivsea(sea)   =l=
                               cnldivsea(cnl1,sea)                ;
   watalcpro(pv2,sea)..
      protarb(pv2,sea)*(1-tolpr)*tcdivsea(sea)      =l=
                                              prsea(pv2,sea)      ;

   nwfpalc(m)..  sum(cnl$pvcnl("nwfp",cnl), canaldiv(cnl,m))=g=
                                               (1-tolnwfp)*divnwfp(m) ;

  subirrc(z1,g,m)$(ex(z1,g)$gs(g)).. wdivrz(z1,g,m) =g=
                                 (1-subirrfac(z1) )*sum((c,t,s,w)$tech(z1,c,t,s,w), wnr(c,z1,t,s,w,m)*x(z1,g,c,t,s,w)) ;

 nbal(n,m)$nb(n)..
   sum(i$ni(n,i), inflow(i,m)) +
   sum(n1, rivercd(n,"d")*trib(n1,n,m)+rivercd(n,"c")*trib(n1,n,m--1) )+
   sum(n1$nn(n,n1),
        f(n,n1,m)*lceff(n1,n)$lceff(n1,n) +
        (riverb(n,n1)*f(n,n1,m)+
                 rivercd(n,"c")*f(n,n1,m--1)  )$(lceff(n1,n) eq 0) )-
   sum(n1$nn(n1,n),
        f(n1,n,m)) + ( rcont(n,m--1) - rcont(n,m) - revapl(n,m)/1000)$rcap(n) -
   sum(cnl$nc(n,cnl), canaldiv(cnl,m)) + artwaternd(n,m) =e= 0 ;

*- bounds for the network link canal capacities.
   artwaternd.lo(n,m) = 0   ;
   f.up(n,n1,m)       = inf ;
   f.up(n,n1,m)$( ncap(n1,n) ne 0 ) = ncap(n1,n);

*-  bounds on family labor
   familyl.up(z1,g,m) =  resource(z1,g,"farmpop")*lstd/1000 ;

* consumption bounds
   consump.fx(z1,g,cq) = farmcons(z1,cq)*consratio(z1,g) ;
   export.up(z1,ce)    = explimit(z1,ce) ;
   itr.fx(z1,g)=0; itw.fx(z1)  =0;

*-
 model wsisz agroclimatic zones model linear objective /
       objz,      cost,      conv,      demnat,    ccombal, qcombal,
       consbal,   laborc,    fodder,    protein,   grnfdr,  bdraft,
       brepco,    tdraft,    trcapc,    twcapc,    landc,   orchareac,
       scmillc,   waterbaln, watalcz,   subirrc               /

 model wsiszn agroclimatic zones model non-linear objective /
       objzn,     cost,                 demnatn,   ccombal, qcombal,
       consbal,   laborc,    fodder,    protein,   grnfdr,  bdraft,
       brepco,    tdraft,    trcapc,    twcapc,    landc,   orchareac,
       scmillc,   waterbaln, watalcz,   subirrc               /
 model wsisn  ibmr model with water network  linear /
       objn,      cost,      conv,      demnat,    ccombal,   qcombal,
       consbal,   laborc,    fodder,    protein,   grnfdr,    bdraft,
       brepco,    tdraft,    trcapc,    twcapc,    landc,     orchareac,
                  waterbaln, watalcz,   subirrc,
       divcnlsea, prseaw,    divsea,    watalcsea, watalcpro, nwfpalc,
       nbal /

 model wsisnn ibmr model with water network  non-linear /
       objnn,     cost,                 demnatn,   ccombal,   qcombal,
       consbal,   laborc,    fodder,    protein,   grnfdr,    bdraft,
       brepco,    tdraft,    trcapc,    twcapc,    landc,     orchareac,
                  waterbaln, watalcz,   subirrc,
       divcnlsea, prseaw,    divsea,    watalcsea, watalcpro, nwfpalc,
       nbal /

           ;
*-
* irrigation canal capicity bounds
  canaldiv.up(cnl,m) = sum(isr, comdef(isr,"ccap",cnl) )  ;
  canaldiv.lo(cnl,m) = divpost(cnl,m) ;

* reservoir operating rule bounds
*
   rcont.lo(n,m) = rulelo(n,m)*rcap(n)/100;
   rcont.up(n,m) = ruleup(n,m)*rcap(n)/100 ;
*
   tolnwfp    = 1 ;
*- river flow tests
   trib("chasma-r","taunsa-b",m)    = 0 ;
   trib("tarbela-r","kalabagh-r",m) = 0 ;
   inflow("haro",m) = 0;  inflow("soan",m) = 0 ;

option limrow  =      0
       limcol  =      0;

   solve wsisn  maximizing cps using  lp     ;