* AIMMS/GAMS file October 2003; by Olga Kiuila

* $TITLE A GENERAL EQUILIBRIUM MODEL FOR EESTONIA


$offsymlist offsymxref
 option limrow=23
 option limcol=2
*option solprint=off
option decimals=4


$INCLUDE "NEWIO971.INC"

* A file is included with all the sets used in this model, and with the 
* I-O-table for Estonia in 1997.

*-------------------------------------------------------------------------

 SCALAR TIME  NUMBER OF YEARS FROM INITIAL POINT  /0/;

*--------------------------------------------------------------------------


***********************************************************************************
* abbreviations: {q}- quantity, {p}-price, {v}- value, {c}-coefficient, FUELe=HEAT*
***********************************************************************************


* INPUT-OUTPUT-COEFFICIENTS (DERIVED FROM IO)


* Parameter IO will be used for scaling the problem {v}

 PARAMETER IO(S,S1);
           IO(S,S1) = SA(S,S1)/1;

* Parameter INOUT is the i-o-value for the intermediate demand{v}

 PARAMETER INOUT(I,J);
           INOUT(I,J) = IO(I,J);

* Parameter AT is the i-o-coefficients for the intermediate demand{c}.

 PARAMETER AT(I,J) ;
           AT(I,J) = INOUT(I,J)/IO("GROSSPROD",J);

* Parameter AM is the i-o-coefficient for the import demand.{c}

 PARAMETER AM(J) ;
           AM(J) = IO("IMPORT",J)/IO("GROSSPROD",J);

* Parameter TTAX is the i-o-coefficient for the payments to the government.{c}

 PARAMETER TTAX(J);
           TTAX(J) = IO("TAXES",J) / IO("GROSSPROD",J) ;

* Parameter PROF is the i-o-coefficient for the catagory profits.{c}

 PARAMETER PROF(J);
           PROF(J) = IO("PROFIT",J) / IO("GROSSPROD" ,J);

* Parameter A is the i-o-coefficient for the input of production factors.{c}

 PARAMETER A(J) ;
           A(J) = 1 - SUM(I,AT(I,J)) - AM(J) - TTAX(J) - PROF(J) ;


* Parameter EIO is the i-o-coefficient for the input of industrial
* products in the energy sectors.{c}

 PARAMETER EIO(I,EN) ;
           EIO(I,EN) = IO(I,EN) / IO ("GROSSPROD",EN);

* Parameter ENTAX is the i-o-coefficient for the payments to the government.{c}

 PARAMETER ENTAX(EN);
           ENTAX(EN) = IO("TAXES",EN) / IO("GROSSPROD",EN) ;

* Parameter PROFEN is the i-o-coefficient for profits in the energy sectors.{c}

 PARAMETER PROFEN(EN) ;
           PROFEN(EN) = IO("PROFIT",EN)/IO("GROSSPROD",EN);

* Parameter ME is the i-o-coefficient for the imports in the energy sectors.{c}

 PARAMETER ME(EN) ;
           ME(EN) = IO("IMPORT",EN)/IO("GROSSPROD",EN);

 PARAMETER AE(EN1);
           AE(EN1)=1 - SUM(I,EIO(I,EN1)) - ME(EN1) - ENTAX(EN1) - PROFEN(EN1);



*--------------------------------------------------------------------------
* PARAMETERS IN CES-PRODUCTION-FUNCTION

* FACTOR SHARES

* Factor shares are defined as the quantity of the factor added, divided
* by the quantity of all factors taken into account until then.
* A nested CES-function is considered, so there will be one factor of
* production added in each step.
* QF is a combination of COALe and FUEL.
* RFe is a combination of QF and FUELe
* VE is a combination of QF and ELECTRICITY.
* UK is a combination of VE and CAPITAL.
* YL is a combination of UK and LABOURH.
* LL is a combination of YL and LABOURL.


* The following 12 parameters are coefficients{c}:

 PARAMETER DELQF(J) ;
       DELQF(J) = IO("COALe",J) / (IO("COALe",J)+IO("FUEL",J));

 PARAMETER DELRFe (J);
        DELRFe(J) = IO("FUELe",J) / (IO("FUELe",J)+IO("COALe",J)+IO("FUEL",J));

 PARAMETER DELVE(J) ;
       DELVE(J) = IO("ELECe",J) / (IO("ELECe",J)+IO("COALe",J)+IO("FUEL",J)+IO("FUELe",J));

 PARAMETER DELUK(J) ;
       DELUK(J) = IO("CAPITAL",J) / (IO("CAPITAL",J)+IO("ELECe",J)+
                  IO("COALe",J)+IO("FUEL",J)+ IO("FUELe",J));

 PARAMETER DELYL(J) ;
       DELYL(J) = IO("LABOURH",J) / (IO("LABOURH",J)+IO("CAPITAL",J)+
                  IO("ELECe",J)+IO("COALe",J)+IO("FUEL",J)+ IO("FUELe",J));

 PARAMETER DELLL(J) ;
       DELLL(J) = IO("LABOURL",J) / (IO("LABOURL",J)+IO("LABOURH",J)+IO("CAPITAL",J)+
                  IO("ELECe",J)+IO("COALe",J)+IO("FUEL",J)+ IO("FUELe",J));


 PARAMETER EDELQF(EN) ;
       EDELQF(EN) = IO("COALe",EN) / (IO("COALe",EN)+IO("FUEL",EN));

 PARAMETER EDELRFe (EN);
        EDELRFe(EN) = IO("FUELe",EN) / (IO("FUELe",EN)+IO("COALe",EN)+IO("FUEL",EN));

 PARAMETER EDELVE(EN) ;
       EDELVE(EN) = IO("ELECe",EN) / (IO("ELECe",EN)+IO("COALe",EN)+IO("FUEL",EN)+IO("FUELe",EN));

 PARAMETER EDELUK(EN) ;
       EDELUK(EN) = IO("CAPITAL",EN) / (IO("CAPITAL",EN)+IO("ELECe",EN)+
                  IO("COALe",EN)+IO("FUEL",EN)+ IO("FUELe",EN));

 PARAMETER EDELYL(EN) ;
       EDELYL(EN) = IO("LABOURH",EN) / (IO("LABOURH",EN)+IO("CAPITAL",EN)+
                  IO("ELECe",EN)+IO("COALe",EN)+IO("FUEL",EN)+ IO("FUELe",EN));

 PARAMETER EDELLL(EN) ;
       EDELLL(EN) = IO("LABOURL",EN) / (IO("LABOURL",EN)+IO("LABOURH",EN)+IO("CAPITAL",EN)+
                  IO("ELECe",EN)+IO("COALe",EN)+IO("FUEL",EN)+ IO("FUELe",EN));


*---------------------------------------------------------------------------

* SUBSTITUTION ELASTICITIES

* These parameters give the elasticity of substitution between two
* succesive factors of production. 

TABLE SIGMA(SIG,J)   SUBSTITUTION ELASTICITIES
          METALe   MINERALe  WOODPAP  ELENGEE   LIGHT    FOOD   CHEMICe 
    QF     0.25      0.30     0.30     0.50     0.50     0.50    0.40
    RFe    0.25      0.30     0.30     0.50     0.00     0.50    0.40  
    VE     0.40      0.70     0.35     0.55     0.50     0.70    0.40
    UK     0.80      0.10     0.10     0.95     0.95     0.80    0.80
    YL     0.60      0.75     0.85     0.80     0.80     0.90    0.60
    LL     0.30      0.35     0.45     0.30     0.30     0.40    0.30
+
          CONSTR    AGRICUL  TRANSPe SERVICES PUBLSERV  PUBLSERVe
    QF     0.50      0.50     0.05     0.05     0.25     0.25  
    RFe    0.50      0.50     0.05     0.05     0.60     0.52  
    VE     0.70      0.70     0.50     0.50     0.52     0.52   
    UK     0.80      0.80     0.10     0.55     0.22     0.22    
    YL     0.90      0.90     0.80     0.80     0.70     0.70    
    LL     0.40      0.40     0.30     0.30     0.20     0.20            ;


TABLE SIGME(SIG,EN)   SUBSTITUTION ELASTICITIES
          COALe     FUEL     FUELe    ELECe 
    QF     0.25      0.25     0.25     0.25     
    RFe    0.25      0.25     0.25     0.25    
    VE     0.40      0.40     0.40     0.40     
    UK     0.80      0.80     0.80     0.80    
    YL     0.60      0.60     0.60     0.60     
    LL     0.30      0.30     0.30     0.30             ;


*-----------------------------------------------------------------------
* VAT rate for coal, fuel and electricity {c}

 PARAMETER ALFAVATE(EN)
           / COALe  1.0
             FUEL   1.18
             FUELe  1.18
             ELECe  1.18  /

 SCALAR ALFAVATM;
        ALFAVATM = 1.0 ;
           
 PARAMETER ALFAVAT(I)
     / METALe    1.18
       MINERALe  1.18
       WOODPAP   1.18
       ELENGEE   1.18
       LIGHT     1.18
       FOOD      1.18
       CHEMICe   1.18
       CONSTR    1.18
       AGRICUL   1.18
       TRANSPe   1.18
       SERVICES  1.18
       PUBLSERV  1.0
       PUBLSERVe 1.18  /

 SCALAR RR;
        RR = 1.0;

*-----------------------------------------------------------------------

* ENDOWMENTS

* The endowments of the input factors are considered to grow by a
* certain growth rate annually.

* Growth rates of input factors {c}

* Capital growth rate
  SCALAR CAPGRO;
         CAPGRO=0.00;

* Labour growth rate low educated (labor force)
  SCALAR LABLGRO;
         LABLGRO=0.00;

* Labour growth rate high educated
  SCALAR LABHGRO;
         LABHGRO=0.0;


* Capital endowment.{q}
  SCALAR K ;
     K =((1+CAPGRO)**TIME)*IO("CAPITAL","TOTDEL");

*-----------------------------------------------------------------------

PARAMETER EMP0(I)
* 1997 [thous.people] {q}
  / METALe        0.4
    MINERALe      3.8
    WOODPAP      38.0
    ELENGEE      23.5
    LIGHT        28.5
    FOOD         27.7
    CHEMICe      10.4
    CONSTR       44.8
    AGRICUL      56.8
    TRANSPe      39.2
    SERVICES    186.1
    PUBLSERV    130.2
    PUBLSERVe     8.5  /

PARAMETER EMPE0(EN)
* 1997 energy sectors {q}
  / COALe          6.8
    FUEL           0.0
    FUELe          4.1
    ELECe          7.6  /
*{q}
PARAMETER EMPTOT0;
          EMPTOT0=SUM(I,EMP0(I))+SUM(EN,EMPE0(EN));

*{c}
PARAMETER SHARELOW(I)
  / METALe      0.71
    MINERALe    0.63
    WOODPAP     0.78
    ELENGEE     0.67
    LIGHT       0.81
    FOOD        0.70
    CHEMICe     0.68
    CONSTR      0.72
    AGRICUL     0.84
    TRANSPe     0.74
    SERVICES    0.48
    PUBLSERV    0.28
    PUBLSERVe   0.65  /


PARAMETER SHARELOWE(EN)
  / COALe        0.95
    FUEL         0.00
    FUELe        0.67
    ELECe        0.67  /

*The following 6 parameters are quantities {q}:

PARAMETER EMPL0(I);
          EMPL0(I)=EMP0(I)*SHARELOW(I);

PARAMETER EMPLE0(EN);
          EMPLE0(EN)=EMPE0(EN)*SHARELOWE(EN);

PARAMETER EMPH0(I);
          EMPH0(I)=EMP0(I)*(1-SHARELOW(I));

PARAMETER EMPHE0(EN);
          EMPHE0(EN)=EMPE0(EN)*(1-SHARELOWE(EN));

PARAMETER EMPLTOT0;
          EMPLTOT0=SUM(I,EMPL0(I))+SUM(EN,EMPLE0(EN));

PARAMETER EMPHTOT0;
          EMPHTOT0=SUM(I,EMPH0(I))+SUM(EN,EMPHE0(EN));

SCALAR UNL0;
       UNL0=50.4;

SCALAR UNH0;
       UNH0=14.4;


* SECTORAL GROSS WAGE IN BASE YEAR [tys.zl/os.] {p}

 PARAMETER  WAGEL0(I);
            WAGEL0(I)=IO("LABOURL",I)/EMPL0(I);

 PARAMETER  WAGELE0(EN);
            WAGELE0(EN)=IO("LABOURL",EN)/EMPLE0(EN);

 PARAMETER  WAGEH0(I);
            WAGEH0(I)=IO("LABOURH",I)/EMPH0(I);

 PARAMETER  WAGEHE0(EN);
            WAGEHE0(EN)=IO("LABOURH",EN)/EMPHE0(EN);

 SCALAR WAGELL0; 
        WAGELL0=IO("LABOURL","TOTDEL")/(SUM(I,EMPL0(I))+SUM(EN,EMPLE0(EN)));  
 SCALAR WAGELH0; 
        WAGELH0=IO("LABOURH","TOTDEL")/(SUM(I,EMPH0(I))+SUM(EN,EMPHE0(EN)));

***********

PARAMETER TRWAL0(I)
* 1997 [bln EKK] {q}
  / METALe       0.00906
    MINERALe     1.24885
    WOODPAP      2.59606
    ELENGEE      2.02329
    LIGHT        0.88759
    FOOD         3.00759
    CHEMICe      1.98923
    CONSTR       1.43206
    AGRICUL      2.00134
    TRANSPe      1.91914
    SERVICES    23.03222
    PUBLSERV     1.26314
    PUBLSERVe    3.22325  /

PARAMETER TRWALE0(EN)
* 1997 energy sectors {q}
  / COALe         0.57696
    FUEL          0.01244
    FUELe         1.34485
    ELECe         2.45898  /

* SECTORAL CAPITAL RATE IN BASE YEAR [capital units] {p}

 PARAMETER  RATE0(I);
            RATE0(I)=TRWAL0(I)/IO("CAPITAL",I);

 PARAMETER  RATEE0(EN);
            RATEE0(EN)=TRWALE0(EN)/IO("CAPITAL",EN);


 SCALAR RATTOT0; 
        RATTOT0=(SUM(I,TRWAL0(I))+SUM(EN,TRWALE0(EN)))/IO("CAPITAL","TOTDEL");

*-----------------------------------------------------------------------

* Exogenous current account surplus.

*{c}
 SCALAR SBARGRO;
        SBARGRO=0.00;

*{v}
 SCALAR SBAR ;
       SBAR =(IO("GROSSPROD","EXPORT")-IO("IMPORT","TOTDEL"))*(1+SBARGRO)**TIME;

* Exogenous demand for output of sectors.{q}
 PARAMETER EXDEM(I);
           EXDEM(I) = IO(I,"STOCK") + IO(I,"EXTDEM");

* Exogenous demand for coal.{q}
  SCALAR EXDEMC;
         EXDEMC = IO("COALe","STOCK") + IO("COALe","EXTDEM");

* Exogenous demand for electricity{q}
  SCALAR EXDEME;
         EXDEME = IO("ELECe","STOCK") + IO("ELECe","EXTDEM");

* Exogenous demand for fuels{q}
  SCALAR EXDEMF;
         EXDEMF = IO("FUEL","STOCK") + IO("FUEL","EXTDEM");

  SCALAR EXDEMFe;
         EXDEMFe = IO("FUELe","STOCK") + IO("FUELe","EXTDEM");

* Exogenous demand for imports.{q}
   SCALAR EXDEMM;
          EXDEMM = IO("IMPORT","STOCK") + IO("IMPORT","EXTDEM");


*----------------------------------------------------------------

* BENCHMARK VALUES
* All benchmark values are based on the I-O-table 

* Benchmark production level.{q}
 PARAMETER X0(I) ;
           X0(I) =  IO(I,"TOTDEL");

 PARAMETER XE0(EN) ;
           XE0(EN) =  IO(EN,"TOTDEL");

* Benchmark government consumption level.{q}
 PARAMETER D0GV(I) ;
           D0GV(I) =  IO(I,"CONGOV");

* Benchmark household consumption level.{q}
 PARAMETER D0HR(I) ;
           D0HR(I) =  IO(I,"CONHR");

 PARAMETER D0HP(I) ;
           D0HP(I) =  IO(I,"CONHP");

* Benchmark export level for non-price-taking sectors.{q}
 PARAMETER ZN0(I);
           ZN0(I)$N(I) = IO(I,"EXPORT");

PARAMETER EZN0(EN);
           EZN0(EN)= IO(EN,"EXPORT");

* Benchmark government consumption of import.{q}
 SCALAR DM0GV ;
        DM0GV = IO("IMPORT","CONGOV");

* Benchmark household consumption of import.{q}
 SCALAR DM0HR ;
        DM0HR =  IO("IMPORT","CONHR");

 SCALAR DM0HP ;
        DM0HP =  IO("IMPORT","CONHP");

* Benchmark government consumption of coal.{q}
 SCALAR DC0GV ;
        DC0GV =  IO("COALe","CONGOV");

* Benchmark household consumption of coal.{q}
 SCALAR DC0HR ;
        DC0HR =  IO("COALe","CONHR");

 SCALAR DC0HP ;
        DC0HP =  IO("COALe","CONHP");

* Benchmark government consumption of electricity.{q}
 SCALAR DE0GV;
        DE0GV = IO("ELECe","CONGOV");

* Benchmark household consumption of electricity.{q}
 SCALAR DE0HR;
        DE0HR =  IO("ELECe","CONHR");

 SCALAR DE0HP;
        DE0HP =  IO("ELECe","CONHP");

* Benchmark government consumption of fuel.{q}
 SCALAR DF0GV ;
        DF0GV = IO("FUEL","CONGOV");

 SCALAR DFe0GV ;
        DFe0GV = IO("FUELe","CONGOV");

* Benchmark household consumption of fuel.{q}
 SCALAR DF0HR ;
        DF0HR =  IO("FUEL","CONHR");

 SCALAR DFe0HR ;
        DFe0HR =  IO("FUELe","CONHR");

 SCALAR DF0HP ;
        DF0HP =  IO("FUEL","CONHP");

 SCALAR DFe0HP ;
        DFe0HP =  IO("FUELe","CONHP");


* Benchmark allocation of capital.{q}
 PARAMETER K0(J) ;
           K0(J) =  IO("CAPITAL",J);

* Benchmark allocation of labour.{q}
 PARAMETER LHIGH0(J) ;
           LHIGH0(J) =  IO("LABOURH",J);

 PARAMETER LLOW0(J) ;
           LLOW0(J) =  IO("LABOURL",J);

* Benchmark allocation of coal.{q}
 PARAMETER C0(J) ;
           C0(J) =  IO("COALe",J);

 PARAMETER CN0(EN) ;
           CN0(EN) =  IO("COALe",EN);

* Benchmark allocation of electricity.{q}
 PARAMETER E0(J) ;
           E0(J) =  IO("ELECe",J);

 PARAMETER EN0(EN) ;
           EN0(EN) =  IO("ELECe",EN);

* Benchmark allocation of fuel.{q}
 PARAMETER F0(J) ;
           F0(J) =  IO("FUEL",J);

 PARAMETER FN0(EN) ;
           FN0(EN) =  IO("FUEL",EN);

 PARAMETER Fe0(J) ;
           Fe0(J) =  IO("FUELe",J);

 PARAMETER FeN0(EN) ;
           FeN0(EN) =  IO("FUELe",EN);

* Benchmark input coefficient for capital.{c}
 PARAMETER KX0(J);
           KX0(J) = IO("CAPITAL",J) / IO("GROSSPROD",J) ;

 PARAMETER KXE0(EN);
           KXE0(EN) = IO("CAPITAL",EN) / IO("GROSSPROD",EN) ;

* Benchmark input coefficient for labour.{c}
 PARAMETER LHIGHX0(J);
           LHIGHX0(J) = IO("LABOURH",J) / IO("GROSSPROD",J) ;

 PARAMETER LHIGHXE0(EN);
           LHIGHXE0(EN) = IO("LABOURH",EN) / IO("GROSSPROD",EN) ;

 PARAMETER LLOWX0(J);
           LLOWX0(J) = IO("LABOURL",J) / IO("GROSSPROD",J) ;

 PARAMETER LLOWXE0(EN);
           LLOWXE0(EN) = IO("LABOURL",EN) / IO("GROSSPROD",EN) ;

* Benchmark input coefficient for electricity.{c}
 PARAMETER EX0(J) ;
           EX0(J) = IO("ELECe",J) / IO("GROSSPROD",J) ;

  PARAMETER EXE0(EN) ;
            EXE0(EN) = IO("ELECe",EN) / IO("GROSSPROD",EN) ;


* Benchmark input coefficient for fuel.{c}
 PARAMETER FX0(J) ;
           FX0(J) = IO("FUEL",J) / IO("GROSSPROD",J) ;

 PARAMETER FXE0(EN) ;
           FXE0(EN) = IO("FUEL",EN) / IO("GROSSPROD",EN) ;


 PARAMETER FeX0(J) ;
           FeX0(J) = IO("FUELe",J) / IO("GROSSPROD",J) ;

 PARAMETER FeXE0(EN) ;
           FeXE0(EN) = IO("FUELe",EN) / IO("GROSSPROD",EN) ;

* Benchmark input coefficient for coal.{c}
 PARAMETER CX0(J) ;
           CX0(J) = IO("COALe",J) / IO("GROSSPROD",J) ;

 PARAMETER CXE0(EN) ;
           CXE0(EN) = IO("COALe",EN) / IO("GROSSPROD",EN) ;

 PARAMETER CAPENS0(EN);
        CAPENS0(EN)=IO("CAPITAL",EN);

* Benchmark import in energy sectors.{q}
 SCALAR IMPEN0;
        IMPEN0 = SUM(EN, IO("IMPORT",EN));

 PARAMETER IMPENS0(EN) ;
        IMPENS0(EN) = IO("IMPORT",EN);

* Benchmark demand for goods in energy sectors.{q}
 PARAMETER ED0(I);
           ED0(I)=SUM(EN, IO(I,EN));

 PARAMETER EDS0(I,EN);
           EDS0(I,EN)= IO(I,EN);

*----------------------------------------------------------------------
* PARAMETERS IN LINEAR EXPENDITURE SYSTEM

* EXPENDITURE SHARES

* The marginal expenditure shares are defined equal to the base
* year average expenditure shares.


* Government expenditure share for coal {c}

  SCALAR BETCGV ;
        BETCGV = IO("COALe","CONGOV")/(IO("GROSSPROD","CONGOV")-
                 IO("TAXES","CONGOV")) ;

* Household expenditure share for coal {c}

  SCALAR BETCHR;
      BETCHR = (IO("COALe","CONHR")+IO("COALe","CONHR")*(ALFAVATE("COALe")-1))
   /(IO("GROSSPROD","CONHR")) ;

  SCALAR BETCHP ;
       BETCHP = (IO("COALe","CONHP")+IO("COALe","CONHP")*(ALFAVATE("COALe")-1))
   /(IO("GROSSPROD","CONHP")) ;

* Government expenditure share for fuels {c}

 SCALAR BETFGV ;
        BETFGV = IO("FUEL", "CONGOV")/( IO("GROSSPROD","CONGOV")-
                 IO("TAXES","CONGOV")) ;

 SCALAR BETFeGV ;
        BETFeGV =IO("FUELe", "CONGOV")/( IO("GROSSPROD","CONGOV")-
                 IO("TAXES","CONGOV")) ;

* Household expenditure share for fuels {c}

 SCALAR BETFHR ;
        BETFHR = (IO("FUEL","CONHR")+IO("FUEL","CONHR")*(ALFAVATE("FUEL")-1))
   /(IO("GROSSPROD","CONHR")) ;

 SCALAR BETFeHR ;
        BETFeHR = (IO("FUELe","CONHR")+IO("FUELe","CONHR")*(ALFAVATE("FUELe")-1))
   /(IO("GROSSPROD","CONHR")) ;

 SCALAR BETFHP ;
        BETFHP = (IO("FUEL","CONHP")+IO("FUEL","CONHP")*(ALFAVATE("FUEL")-1))
   /(IO("GROSSPROD","CONHP")) ;

 SCALAR BETFeHP ;
        BETFeHP = (IO("FUELe","CONHP")+IO("FUELe","CONHP")*(ALFAVATE("FUELe")-1))
   /(IO("GROSSPROD","CONHP")) ;

* Government expenditure share for electricity {c}

 SCALAR BETEGV ;
        BETEGV = IO("ELECe","CONGOV")/(IO("GROSSPROD","CONGOV") -
                 IO("TAXES","CONGOV")) ;

* Household expenditure share for electricity {c}

 SCALAR BETEHR ;
        BETEHR = (IO("ELECe","CONHR")+IO("ELECe","CONHR")*(ALFAVATE("ELECe")-1))
   /(IO("GROSSPROD","CONHR")) ;

 SCALAR BETEHP ;
        BETEHP = (IO("ELECe","CONHP")+IO("ELECe","CONHP")*(ALFAVATE("ELECe")-1))
   /(IO("GROSSPROD","CONHP")) ;

* Government expenditure share for imports {c}

 SCALAR BETMGV ;
        BETMGV = IO("IMPORT","CONGOV")/(IO("GROSSPROD","CONGOV")-
                 IO("TAXES","CONGOV")) ;

* Household expenditure share for imports {c}

 SCALAR BETMHR ;
        BETMHR = IO("IMPORT","CONHR")/(IO("GROSSPROD","CONHR")) ;

 SCALAR BETMHP ;
        BETMHP = IO("IMPORT","CONHP")/(IO("GROSSPROD","CONHP")) ;

* Government marginal expenditure shares for industrial goods {c}

 PARAMETER BETGV(I) ;
           BETGV(I) = IO(I,"CONGOV")/(IO("GROSSPROD","CONGOV") -
                      IO("TAXES","CONGOV")) ;

* Household marginal expenditure shares for industrial goods {c}

 PARAMETER BETHR(I) ;
           BETHR(I) = (IO(I,"CONHR")+IO(I,"CONHR")*(ALFAVAT(I)-1))
   /(IO("GROSSPROD","CONHR")) ;

 PARAMETER BETHP(I) ;
           BETHP(I) = (IO(I,"CONHP")+IO(I,"CONHP")*(ALFAVAT(I)-1))
   /(IO("GROSSPROD","CONHP")) ;


* LES-constants
* The constant part of the LES-function is equal to 75, 50, 25 percent of the
* base year consumption level

*The following 16 parameters are  quantities {q}:

 SCALAR DCGVBAR ;
        DCGVBAR =  0.75 * ( IO("COALe","CONGOV"));

 SCALAR DCHRBAR ;
        DCHRBAR =  0.5 * (IO("COALe","CONHR"));

 SCALAR DCHPBAR ;
        DCHPBAR =  0.25 * (IO("COALe","CONHP"));

 SCALAR DFGVBAR ;
        DFGVBAR =  0.75 * (IO("FUEL","CONGOV"));

 SCALAR DFHRBAR ;
        DFHRBAR =  0.5 * (IO("FUEL","CONHR"));

 SCALAR DFHPBAR ;
        DFHPBAR =  0.25 * (IO("FUEL","CONHP"));

 SCALAR DFeGVBAR ;
        DFeGVBAR =  0.75 * (IO("FUELe","CONGOV"));

 SCALAR DFeHPBAR ;
        DFeHPBAR =  0.25 * (IO("FUELe","CONHP"));

 SCALAR DFeHRBAR ;
        DFeHRBAR =  0.5 * (IO("FUELe","CONHR"));

 SCALAR DEGVBAR ;
        DEGVBAR =  0.75 * (IO("ELECe","CONGOV"));

 SCALAR DEHRBAR ;
        DEHRBAR =  0.5 * (IO("ELECe","CONHR"));

 SCALAR DEHPBAR ;
        DEHPBAR =  0.25 * (IO("ELECe","CONHP"));

 SCALAR DMGVBAR ;
        DMGVBAR =  0.75 * ( IO("IMPORT","CONGOV"));

 SCALAR DMHRBAR ;
        DMHRBAR =  0.5 * (IO("IMPORT","CONHR"));

 SCALAR DMHPBAR ;
        DMHPBAR =  0.25 * (IO("IMPORT","CONHP"));

 PARAMETER DGVBAR(I) ;
           DGVBAR(I) =  0.75 * (IO(I,"CONGOV"));

 PARAMETER DHRBAR(I) ;
           DHRBAR(I) =  0.5 * (IO(I,"CONHR"));

 PARAMETER DHPBAR(I) ;
           DHPBAR(I) =  0.25 * (IO(I,"CONHP"));

*---------------------------------------------------------------------------

* PRICES

* For the non-price-taking sectors the world-market price is exogenously
* The domestic price is determined endogenously.

  PARAMETER PWN(I)
     / 
       METALe    1
       MINERALe  1
       WOODPAP   1  
       ELENGEE   1
       LIGHT     1
       FOOD      1
       CHEMICe   1
       CONSTR    1
       AGRICUL   1
       TRANSPe   1
       SERVICES  1
       PUBLSERV  1
       PUBLSERVe 1  /

  PARAMETER PWEN(EN)

     / COALe     1.0
       FUEL      1.0
       FUELe     1.0
       ELECe     1.0  /



* Estonia can be considered as a price-taker in the fuel-sector, so an exogenous
* (world-market) price is introduced.

* The world market price for fuel will grow by a certain rate annually.{c}

 SCALAR PFGRO   ANNUAL GROWTH RATE FOR FUEL PRICE (IN %)  /0.00/;

*The world market proce for fuel {p}
 SCALAR PFX;
        PFX = (1+PFGRO)**TIME; 

*{c}
 SCALAR PLLGRO   ANNUAL GROWTH RATE FOR LABOUR PRICE LOW EDUCATED (IN %) /0.00/;

 SCALAR PLHGRO  ANNUAL GROWTH RATE FOR LABOUR PRICE HIGH EDUCATED (IN %) /0.00/;

 SCALAR LLFGRO   ANNUAL GROWTH RATE FOR EMPLOYMENT OF LOW EDUCATED (IN %) /0.00/;

 SCALAR LHFGRO  ANNUAL GROWTH RATE FOR EMPLOYMENT OF HIGH EDUCATED (IN %) /0.00/;


* BUTTON FOR REGIME SWITCHES LABOUR MARKET

 SCALAR LSWITCH FIXED ENDOWMENT WHEN   /1/;

 SCALAR HSWITCH FIXED ENDOWMENT WHEN NO  /0/;


 PARAMETER GAMMAHTAX(I)
* [0.14=14%] {c}

     / METALe    .14
       MINERALe  .14
       WOODPAP   .14
       ELENGEE   .14
       LIGHT     .14
       FOOD      .14
       CHEMICe   .14
       CONSTR    .14
       AGRICUL   .14
       TRANSPe   .14
       SERVICES  .14
       PUBLSERV  .14
       PUBLSERVe .14 /

 PARAMETER GAMMAEHTAX(EN)
     / COALe     .14
       FUEL      .14
       FUELe     .14
       ELECe     .14  /

 PARAMETER GAMMAHSOC(I)
     / METALe    .21
       MINERALe  .21
       WOODPAP   .21
       ELENGEE   .21
       LIGHT     .21
       FOOD      .21
       CHEMICe   .21
       CONSTR    .21
       AGRICUL   .21
       TRANSPe   .21
       SERVICES  .21
       PUBLSERV  .21
       PUBLSERVe .21 /

 PARAMETER GAMMAEHSOC(EN)
     / COALe     .21
       FUEL      .21
       FUELe     .21
       ELECe     .21  /

 PARAMETER GAMMALTAX(I)
     / METALe    .14
       MINERALe  .14
       WOODPAP   .14
       ELENGEE   .14
       LIGHT     .14
       FOOD      .14
       CHEMICe   .14
       CONSTR    .14
       AGRICUL   .14
       TRANSPe   .14
       SERVICES  .14
       PUBLSERV  .14
       PUBLSERVe .14 /

 PARAMETER GAMMAELTAX(EN)
     / COALe     .14
       FUEL      .14
       FUELe     .14
       ELECe     .14  /

 PARAMETER GAMMALSOC(I)
     / METALe    .21
       MINERALe  .21
       WOODPAP   .21
       ELENGEE   .21
       LIGHT     .21
       FOOD      .21
       CHEMICe   .21
       CONSTR    .21
       AGRICUL   .21
       TRANSPe   .21
       SERVICES  .21
       PUBLSERV  .21
       PUBLSERVe .21 /

 PARAMETER GAMMAELSOC(EN)
     / COALe     .21
       FUEL      .21
       FUELe     .21
       ELECe     .21  /

*{c}
 SCALAR PMGRO   ANNUAL GROWTH RATE FOR IMPORT PRICE (IN %)  /0.00/;

* Prices of import are considered exogenously{p}
 SCALAR PM;
        PM = (1+PMGRO)**TIME;


*----------------------------------------------------------------------------- 

* EMISSION COEFFICIENTS (these values depends on i/o table that we use)

* Coal emission coefficients [million t / billion EKK] {c}

* Producing sectors
 TABLE CEMS(EN,J)   COALe EMISSION COEFFICIENTS

       METALe   MINERALe  WOODPAP   ELENGEE    LIGHT     FOOD     CHEMICe

COALe  0.00000  101.81784   0.00000   0.00000   2.77555   8.50629   0.00000

 +

       CONSTR    AGRICUL   TRANSPe   SERVICES  PUBLSERV  PUBLSERVe

COALe   2.14381   0.82553   6.62624   1.53989    2.53712     0.00000          
                                                                               ;

* Households
 PARAMETER CEMH;
           CEMH=8.797353;
     

* Energy sectors
 TABLE     CEME(EN,EN1)
            COALe     FUEL     FUELe      ELECe

COALe       46.46076   0.000     8.55163     0.00000    
                                                        ;

*  emission coefficients
 TABLE NEMS(I,J)     EMISSION COEFFICIENTS

          METALe MINERALe  WOODPAP   ELENGEE   LIGHT     FOOD     CHEMICe
METALe    0.0
MINERALe  0.0
WOODPAP   0.0     0.09851   0.18537   0.07213   0.02429   0.03407  0.0
ELENGEE   0.0
LIGHT     0.0
FOOD      0.0
CHEMICe   0.20906 2.81153   0.23518   0.12810   0.05581   0.49054  0.81456
CONSTR    0.0
AGRICUL   0.0
TRANSPe   0.0
SERVICES  0.0
PUBLSERV  0.0
PUBLSERVe 0.0     1.29242   0.25145   0.22838   0.27482   0.63285  4.43028    

 +

          CONSTR    AGRICUL   TRANSPe   SERVICES  PUBLSERV  PUBLSERVe
METALe    0.0
MINERALe  0.0
WOODPAP   0.00761    0.43090   0.07057   0.07020   0.28702   0.0
ELENGEE   0.0
LIGHT     0.0
FOOD      0.0
CHEMICe   0.20350    0.32086   0.52943   0.20526   0.14370   0.00551
CONSTR    0.0
AGRICUL   0.0
TRANSPe   0.0
SERVICES  0.0
PUBLSERV  0.0
PUBLSERVe 0.14840    0.48930   0.11919   0.24825   0.08155   0.00183    
       
                                                                                ;

* Households
 PARAMETER NEMH(I)  / 
METALe    0.0
MINERALe  0.0
WOODPAP   3.13852     
ELENGEE   0.0
LIGHT     0.0
FOOD      0.0
CHEMICe   0.57345
CONSTR    0.0
AGRICUL   0.0
TRANSPe   0.0
SERVICES  0.0
PUBLSERV  0.0
PUBLSERVe 0.23061          
                                        /;

* Energy sectors
 TABLE NEME(I,EN)
         COALe       FUEL      FUELe     ELECe
METALe    0.0
MINERALe  0.0
WOODPAP   1.25529     0.00046   9.83767     0.0
ELENGEE   0.0
LIGHT     0.0
FOOD      0.0
CHEMICe   2.36897     0.19153   2.59172   225.24436
CONSTR    0.0
AGRICUL   0.0
TRANSPe   0.0
SERVICES  0.0
PUBLSERV  0.0
PUBLSERVe 0.0         0.0       1.30113     7.44530     ;




* Process emission coefficients
  PARAMETER PEMS(J)

     / METALe    0.0
*       MINERALe  0.00008
       WOODPAP   0.0
       ELENGEE   0.0
       LIGHT     0.0
       FOOD      0.0
       CHEMICe   0.0
       CONSTR    0.0
       AGRICUL   0.0
       TRANSPe   0.0
       SERVICES  0.0
       PUBLSERV  0.0
       PUBLSERVe 0.0   /


    
* Energy sectors
 PARAMETER PEME(EN)
     / COALe     0.0
       FUEL      0.0
       FUELe     0.0
       ELECe     0.0  /              
                                                ;

                                                     
* Emission limit for CO2 [million t] {q}
 PARAMETER EMLIM     TOTAL EMISSION LIMITS;
           EMLIM=10E+3;       
                        

 
* Total quantity of emission in the base case for 2012.
* This parameter is used to calculate the relative quantity of cleaning [million t] {q}

 PARAMETER REF   REFERENCE EMISSION LEVEL;
           REF=0.0246043 ;      
                       


* The elasticity of the export demand will be used for the non-price taking
* sectors only. {c}

 PARAMETER EPSI(I)       PRICE ELASTICITY OF EXPORT DEMAND
     / 
       METALe    -2.00 
       MINERALe  -2.00
       WOODPAP   -2.00
       ELENGEE   -1.50
       LIGHT     -1.00
       FOOD      -0.25
       CHEMICe   -1.50
       CONSTR    -0.50
       AGRICUL   -0.25
       TRANSPe   -0.25
       SERVICES  -0.50
       PUBLSERV  -0.50
       PUBLSERVe -0.50
                       /

 SCALAR EPSIX;
        EPSIX=1.00;


 PARAMETER EPSIEN(EN)       PRICE ELASTICITY OF EXPORT DEMAND
     / 
       COALe     -2.0
       FUEL      -2.0
       FUELe     -0.50
       ELECe     -1.0  /


PARAMETER SUB(I)    SUBSIDIES PAID BY THE GOVERNMENT PER SECTOR
*[%] {c}
    
     / METALe    0.0
       MINERALe  0.0
       WOODPAP   0.0
       ELENGEE   0.0
       LIGHT     0.0
       FOOD      0.0
       CHEMICe   0.0
       CONSTR    0.0
       AGRICUL   0.0
       TRANSPe   0.0
       SERVICES  0.0
       PUBLSERV  0.0
       PUBLSERVe 0.0  /

SCALAR SUBC    SUBSIDIES PAID BY THE GOVERNMENT TO THE COALe SECTOR   /0.0/ ; 

SCALAR SUBF    SUBSIDIES PAID BY THE GOVERNMENT TO THE FUEL            /0.0/ ;

SCALAR SUBFe   SUBSIDIES PAID BY THE GOVERNMENT TO THE FUELe          /0.0/ ;

SCALAR SUBE    SUBSIDIES PAID BY THE GOVERNMENT TO THE ELECTRICITY SECTOR /0.0/; 

SCALAR SUBLL   SUBSIDIES PAID ON LOW EDUCATED LABOUR (REDUCTION OF TAXES) /0.0/; 

SCALAR SUBLH   SUBSIDIES PAID ON HIGH EDUCATED LABOUR (REDUCTION OFTAXES) /0.0/; 

SCALAR SUBK    SUBSIDIES PAID ON CAPITAL                               / 0.0 / ; 

SCALAR SUBHR   SUBSIDIES PAID TO CONSUMERS                              /0.0/;

SCALAR SUBHP   SUBSIDIES PAID TO CONSUMERS                              /0.0/;


* Annual growth rates for (capital and labour) productivity growth and energy efficiency {c}

 PARAMETER LAMBDA(I)        ANNUAL RATE OF TOTAL PRODUCTIVITY GROWTH;
           LAMBDA(I)=0.0;

PARAMETER LAMBDE(EN)        ANNUAL RATE OF TOTAL PRODUCTIVITY GROWTH;
           LAMBDE(EN)=0.0;

* Annual growth rates for (capital and labour) productivity growth and energy efficiency {c}

 PARAMETER LAMBDALL(I)     ANNUAL RATE OF LABOUR PRODUCTIVITY GROWTH
*[%] 
     / METALe    0.0
       MINERALe  0.0
       WOODPAP   0.0
       ELENGEE   0.0
       LIGHT     0.0
       FOOD      0.0
       CHEMICe   0.0
       CONSTR    0.0
       AGRICUL   0.0
       TRANSPe   0.0
       SERVICES  0.0
       PUBLSERV  0.0
       PUBLSERVe 0.0   /

PARAMETER LAMBDALLE(EN)

     / COALe     0.0
       FUEL      0.0
       FUELe     0.0
       ELECe     0.0  /

 PARAMETER LAMBDALH(I)     ANNUAL RATE OF LABOUR PRODUCTIVITY GROWTH
*[%] 
     / METALe    0.0
       MINERALe  0.0
       WOODPAP   0.0
       ELENGEE   0.0
       LIGHT     0.0
       FOOD      0.0
       CHEMICe   0.0
       CONSTR    0.0
       AGRICUL   0.0
       TRANSPe   0.0
       SERVICES  0.0
       PUBLSERV  0.0
       PUBLSERVe 0.0   /

PARAMETER LAMBDALHE(EN)

     / COALe     0.0
       FUEL      0.0
       FUELe     0.0
       ELECe     0.0 /

PARAMETER LAMBDAK(I)    ANNUAL RATE OF TOTAL CAPITAL PRODUCTIVITY GROWTH

     / METALe    0.0
       MINERALe  0.0
       WOODPAP   0.0
       ELENGEE   0.0
       LIGHT     0.0
       FOOD      0.0
       CHEMICe   0.0
       CONSTR    0.0
       AGRICUL   0.0
       TRANSPe   0.0
       SERVICES  0.0
       PUBLSERV  0.0
       PUBLSERVe 0.0   /

PARAMETER LAMBDAKE(EN)

     / COALe     0.0
       FUEL      0.0
       FUELe     0.0
       ELECe     0.0  /


 TABLE ENEFF(I,EN)      ANNUAL RATE OF ENERGY EFFICIENCY GROWTH PER SECTOR AND FUEL TYPE
*[%] {c}
                 COALe    FUEL    FUELe     ELECe        
       METALe    0        0        0         0
       MINERALe  0        0        0         0
       WOODPAP   0        0  
       ELENGEE   0        0        0         0
       LIGHT     0        0        0         0
       FOOD      0        0        0         0
       CHEMICe   0        0        0         0
       CONSTR    0        0        0         0
       AGRICUL   0        0        0         0
       TRANSPe   0        0        0         0
       SERVICES  0        0        0         0
       PUBLSERV  0        0        0         0
       PUBLSERVe 0        0        0         0   ;

 TABLE ENEFFE (EN,EN1)

                 COALe    FUEL    FUELe     ELECe        
       COALe     0        0        0         0
       FUEL      0        0        0         0
       FUELe     0        0        0         0
       ELECe     0        0        0         0     ;

PARAMETER ENEFFX (EN)    ANNUAL RATE OF ENERGY EFFICIENCY GROWTH PER FUEL TYPE FOR CONSUMERS
                        
    /  COALe     0       
       FUEL      0       
       FUELe     0        
       ELECe     0  /;

 PARAMETER EXPGRO(I)     ANNUAL RATE OF EXPORT DEMAND GROWTH
* [%] {c}

     / METALe    0.0
       MINERALe  0.0
       WOODPAP   0.0
       ELENGEE   0.0
       LIGHT     0.0
       FOOD      0.0
       CHEMICe   0.0
       CONSTR    0.0
       AGRICUL   0.0
       TRANSPe   0.0
       SERVICES  0.0
       PUBLSERV  0.0
       PUBLSERVe 0.0   /


 PARAMETER EXPGROEN(EN)     ANNUAL RATE OF EXPORT DEMAND GROWTH
* [%]
     / COALe    0.0
       FUEL     0.0
       FUELe    0.0
       ELECe    0.0   /



PARAMETER EMISTAX   EMISSION TAX;
* [bln EKK/mln t] {p}

         EMISTAX=0.0;

PARAMETER CTAX;
* excise tax: P(1+T) value tax [%] {c}
*CTAX=0, FTAX=0, FeTAX=0.10 (for 1997), but it is already has been included in the i/o table,
*i.e. 0.10 for Fe in the model will be 0.0 (scale shift). It does not have influence on the results 
*because it is included only in the equations where it does not have influence on marginal values
*or it has additive conversion (It is for KFTAX, PCC, PCT, PFC, PFT, PFec, PFeT)

        CTAX=0.0;

PARAMETER FTAX;
        FTAX=0.0;

PARAMETER FeTAX;
        FeTAX=0.0;

PARAMETER ETAX;
        ETAX=0.0;

PARAMETER NTAX(I)
     / METALe    0.0
       MINERALe  0.0
       WOODPAP   0.0
       ELENGEE   0.0
       LIGHT     0.0
       FOOD      0.0
       CHEMICe   0.0
       CONSTR    0.0
       AGRICUL   0.0
       TRANSPe   0.0
       SERVICES  0.0
       PUBLSERV  0.0
       PUBLSERVe 0.0   /;



* P+T quantity tax [EKK/MJ] has been transformed into [%] P+T
PARAMETER NCTAX;
        NCTAX=0.0;

PARAMETER NFTAX;
        NFTAX=0.0;

PARAMETER NFeTAX;
        NFeTAX=0.0;

PARAMETER NETAX;
        NETAX=0.0;

PARAMETER NNTAX(I)
     / METALe    0.0
       MINERALe  0.0
       WOODPAP   0.0
       ELENGEE   0.0
       LIGHT     0.0
       FOOD      0.0
       CHEMICe   0.0
       CONSTR    0.0
       AGRICUL   0.0
       TRANSPe   0.0
       SERVICES  0.0
       PUBLSERV  0.0
       PUBLSERVe 0.0   /;


*--------------------------------------------------------------------------

 VARIABLES

ABAT    UNIT ABATEMENT COST [bln EKK\mln t] {c}
CLEAN   CLEANING OF EMISSIONS {q}
CLEANSEC(J)
CLEANSECEN(EN)

CT(I)       COALe DEMAND IN (NON-ENERGY) SECTORS {q}
CTE(EN)

DCHR        HOUSEHOLD CONSUMPTION DEMAND FOR COALe {q}
DCHP           HHHH {q}
DCGV        GOVERNMENT CONSUMPTION DEMAND FOR COALe{q}
DEHR        HOUSEHOLD CONSUMPTION DEMAND FOR ELECTRICTY{q}
DEHP           VV {q}
DEGV        GOVERNMENT CONSUMPTION DEMAND FOR ELECTRICITY{q}
DFHR        HOUSEHOLD CONSUMPTION DEMAND FOR FUEL{q}
DFHP           V {q}
DFeHR         V {q}
DFeHP           V{q}
DFGV        GOVERNMENT CONSUMPTION DEMAND FOR FUEL{q}
DFeGV           V{q}
DGV(I)      GOVERNMENT CONSUMPTION DEMAND FOR GOODS{q}
DHR(I)      HOUSEHOLD CONSUMPTION DEMAND FOR GOODS{q}
DHP(I)          V{q}
DMHR        HOUSEHOLD CONSUMPTION DEMAND FOR IMPORT{q}
DMHP            V{q}
DMGV        GOVERNMENT CONSUMPTION DEMAND FOR IMPORT{q}

DPQDPFT(J) PARTIAL OF PQ WRT PF(I) share of PF in (shadow) price PQ for production sectors {c}
DPQDPCT(J)  PARTIAL OF PQ WRT PC(I)
DPRDPFeT(J) PARTIAL OF PR WRT PFe(I)
DPRDPQT(J)  PARTIAL OF PR WRT PQ(I)
DPUDPKT(J)  PARTIAL OF PU WRT PK(I)
DPUDPVT(J)  PARTIAL OF PU WRT PV(I)
DPVDPET(J)  PARTIAL OF PV WRT PE(I)
DPVDPRT(J)  PARTIAL OF PV WRT PR(I)
DPYDPLH(J)  PARTIAL OF PY WRT PLH(I)
DPYDPPU(J)  PARTIAL OF PY WRT PPU(I)
DPLLPLL(J)  PARTIAL OF PL WRT PLL(I)
DPLLPYT(J)  PARTIAL OF PL WRT PYT(I)


DPQDPFTE(EN) PARTIAL OF PQ WRT PF(I)  {c}
DPQDPCTE(EN)  PARTIAL OF PQ WRT PC(I)
DPRDPFeTE(EN) PARTIAL OF PR WRT PFe(I)
DPRDPQTE(EN)  PARTIAL OF PR WRT PQ(I)
DPUDPKTE(EN)  PARTIAL OF PU WRT PK(I)
DPUDPVTE(EN)  PARTIAL OF PU WRT PV(I)
DPVDPETE(EN)  PARTIAL OF PV WRT PE(I)
DPVDPRTE(EN)  PARTIAL OF PV WRT PR(I)
DPYDPLHE(EN)  PARTIAL OF PY WRT PLH(I)
DPYDPPUE(EN)  PARTIAL OF PY WRT PPU(I)
DPLLPLLE(EN)  PARTIAL OF PL WRT PLL(I)
DPLLPYTE(EN)  PARTIAL OF PL WRT PYT(I)


DPXDPCT(J)  PARTIAL OF PX WRT PC(I)  {c}
DPXDPET(J)  PARTIAL OF PX WRT PE(I)
DPXDPFeT(J) PARTIAL OF PX WRT PFe(I)
DPXDPFT(J)  PARTIAL OF PX WRT PF(I) 
DPXDPKT(J)  PARTIAL OF PX WRT PK(I)
DPXDPYT(J)  PARTIAL OF PX WRT PYT(I)
DPXDPLL(J)  PARTIAL OF PX WRT PLL(I)

DPXDPCTE(EN)  PARTIAL OF PX WRT PC(I)  {c}
DPXDPETE(EN)  PARTIAL OF PX WRT PE(I)
DPXDPFeTE(EN) PARTIAL OF PX WRT PFe(I)
DPXDPFTE(EN)  PARTIAL OF PX WRT PF(I) 
DPXDPKTE(EN)  PARTIAL OF PX WRT PK(I)
DPXDPYTE(EN)  PARTIAL OF PX WRT PYT(I)
DPXDPLLE(EN)  PARTIAL OF PX WRT PLL(I)


DUM1        OBJECTIVE VARIABLE{q}
DUME        UNDER UTILIZATION OF ELECTRICITY{q}
ED(I)       TOTAL DEMAND FOR NON-ENERGY GOODS IN ENERGY SECTORS{q}
EDS(I,EN)   DEMAND FOR NON-ENERGY GOODS PER ENERGY SECTOR{q}
EMRED   SHARE OF EMISSION REDUCTION {c}
EMIS(J)  EMISSION PER SECTOR {q}
HOUSEMIS HOUSEHOLD EMISSION
ENEMIS(EN) EMISSION PER ENERGY SECTOR {q}
ET(I)       ELECTRICITY DEMAND IN (NON-ENERGY) SECTORS {q}
ETE(EN)

EXGOV       government expenditures {v}
EXH         households expenditures
EXHR
EXHP

F           TOTAL DEMAND FOR FUEL {q}
C
E
Fe
FT(I)       FUELe DEMAND IN (NON-ENERGY) SECTORS {q}
FTE(EN)
FeT(I)         V {q}
FeTE(EN)
FUND    YELD EMISSION TAX {v}
IMPEN       TOTAL IMPORT IN ENERGY SECTORS {q} 
IMPENS(EN)  IMPORT PER ENERGY SECTOR {q}
IMPORTS     VALUE OF TOTAL IMPORTS{v}
EXPORTS     VALUE OF TOTAL EXPORTS{v}

IN          DISPOSABLE INCOME {v}
INHR         
INHP   
INGV        

KFTAX       TAX RE-ALLOCATION {q}
KPFTAX      GOVERNMENT INCOME FROM FTAXes {v}
KT(I)       CAPITAL INPUT IN (NON-ENERGY) SECTORS{q}
KTE(EN)
RATE(I)
RATEE(EN)
RATTOT
TRWAL(I)
TRWALE(EN)
TRWALTOT

LLT(I)      LOW EDUCATED LABOUR DEMAND IN (NON-ENERGY) SECTORS{q}
LLTE(EN)
LHT(I)      HIGH EDUCATED LABOUR DEMAND IN (NON-ENERGY) SECTORS{q}
LHTE(EN)
M           TOTAL IMPORT {q}
MCT(I)      MARGINAL COST IN PRODUCTION {p}
MCEN(EN)
MTOT        TOTAL IMPORT INCLUDING FUEL {q}

PC          PRICE (PRODUCER) OF COALe (excluding subsidies) {p}
PCC         PRICE OF COALe FOR CONSUMERS (including subsidies and tax)
PCT(I)      PRICE OF COALe FOR SECTORS (including subsidies and tax) - INTERMIDIATE DEMAND
PCEN(EN)
PE          PRICE (PRODUCER) OF ELECTRICITY (excluding subsidies)
PEC
PET(J)
PEEN(EN)
PEM     PRICE OF EMISSION PERMITS
PF          PRICE (PRODUCER) OF FUELS (excluding subsidies, including excise tax)
PFe         PRICE (PRODUCER) OF FUELeS (excluding subsidies)
PFC         PRICE OF FUELS FOR CONSUMERS (including subsidies and tax))
PFeC        PRICE OF FUELeS FOR CONSUMERS (including subsidiesand tax)
PFT(I)      PRICE OF FUELS FOR SECTORS (including subsidies and tax) - INTERMEDIATE DEMAND
PFEN(EN)
PFeT(I)     PRICE OF FUELeS FOR SECTORS (including subsidies and tax) - INTERMEDIATE DEMAND
PFeEN(EN)
PK          PRICE OF CAPITAL (excluding subsidies)
PLHIGH      PRICE OF HIGH EDUCATED LABOUR INCLUDING SUBSIDY (=marginal costs - subsidy)
PLHIGHX     PRICE OF HIGH EDUCATED LABOUR EXCLUDING SUBSIDY (=marginal costs)
PLLOW       PRICE OF LOW EDUCATED LABOUR INCLUDING SUBSIDY (=marginal costs - subsidy)
PLLOWX      PRICE OF LOW EDUCATED LABOUR IEXLUDING SUBSIDY (=marginal costs)

LH          HIGH EDUCATED LABOUR ENDOWMENT{q}
LABSUB      TOTAL SUBSIDIES PAID BY THE GOVERNMENT ON LABOUR COSTS{v}
SUBTOT      TOTAL SUBSIDIES (quantitive) PAID BY THE GOVERNMENT{v}
LL          LOW EDUCATED LABOUR ENDOWMENT{q}

LABHTAX(I)  TAXES IN LABOUR COSTS HIGH EDUCATED {v}
LABHSOC(I)  SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS HIGH EDUCATED{v}
LABLTAX(I)  TAXES IN LABOUR COSTS LOW EDUCATED{v}
LABLSOC(I)   SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS LOW EDUCATED{v}
LABHTAXE(EN) TAXES IN LABOUR COSTS HIGH EDUCATED{v}
LABHSOCE(EN) SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS HIGH EDUCATED{v}
LABLTAXE(EN) TAXES IN LABOUR COSTS LOW EDUCATED{v}
LABLSOCE(EN) SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS LOW EDUCATED{v}
LABTAXTOT   TAXES IN LABOUR COSTS{v}
LABSOCTOT   SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS{v}

UNL         UNEMPLOYMENT{q}
UNH         UNEMPLOYMENT{q}

PN(I)       PRICE OF OUTPUT IN NON-PRICE-TAKING SECTORS{p}
PNC(I)
PNT(I,J)
PNEN(I,EN)
PQT(J)      IMPLICIT PRICE OF Q(I)
PRT(J)      IMPLICIT PRICE OF R(I)
PVT(J)      IMPLICIT PRICE OF V(I)
PPU(J)      IMPLICIT PRICE OF U(I)
PYT(J)      IMPLICIT PRICE OF Y(I)
PLL(J)      IMPLICIT PRICE OF L(I)

PQTE(EN)      IMPLICIT PRICE OF Q(I)
PRTE(EN)      IMPLICIT PRICE OF R(I)
PVTE(EN)      IMPLICIT PRICE OF V(I)
PPUE(EN)      IMPLICIT PRICE OF U(I)
PYTE(EN)      IMPLICIT PRICE OF Y(I)
PLLE(EN)      IMPLICIT PRICE OF L(I)

TOTEM   TOTAL EMISSION{q}
X(I)        OUTPUT IN PRODUCING SECTORS{q}
XGRO(I)     OUTPUT GROWTH FROM BASE YEAR PER SECTOR{c}
XEE(EN)     GROSS TOTAL DEMAND OF ENERGY CARRIERS (WITHOUT ENERGY SECTORS){q}
XEEGRO(EN)  OUTPUTGROWTH FROM BASE YEAR IN ENERGY SECTORS{c}
XPGRO       value of XGRO
XPEEGRO     value of XEEGRO


GDP         Gross Domestic Product - expenditures method
GDP0        GDP IN BASE YEAR
PKB         GDP - income method #1
PKB0
PKKB        gdp- income method #2
PKKB0

ZN(I)       NET EXPORT OF NON-PRICE-TAKERS {q}
ZNE
ZNC
ZNF
ZNFe


DTGVBAR     exogenous input LES government {v}
DTHRBAR     exogenous input LES households {v}
DTHPBAR       V  {v}

VAT         total VAT income for the government {v}
VATR        VAT income for the government from rich household consumption
VATP        VAT income for the government from poor households consumption

VATN(I)     VAT PAID ON non-PRICE TAKING SECTORS USE
VATC        VAT PAID ON COALe USE
VATF        VAT PAID ON FUEL USE
VATFe       VAT PAID ON FUELe USE
VATE         VAT PAID ON ELECTRICITY USE
VATM         VAT PAID ON IMPORTS USE
VATNR(I)     VAT PAID ON non-PRICE TAKING SECTORS USE BY RICH HOUSEHOLDS
VATCR        VAT PAID ON COALe USE BY RICH HOUSEHOLDS
VATFR        VAT PAID ON FUEL USE BY RICH HOUSEHOLDS
VATFeR       VAT PAID ON FUELe USE BY RICH HOUSEHOLDS
VATER        VAT PAID ON ELECTRICITY USE BY RICH HOUSEHOLDS
VATMR        VAT PAID ON IMPORTS USE BY RICH HOUSEHOLDS
VATNP(I)     VAT PAID ON non-PRICE TAKING SECTORS USEBY POOR HOUSEHOLDS
VATCP        VAT PAID ON COALe USEBY POOR HOUSEHOLDS
VATFP        VAT PAID ON FUEL USE BY POOR HOUSEHOLDS
VATFeP       VAT PAID ON FUELe USE BY POOR HOUSEHOLDS
VATEP        VAT PAID ON ELECTRICITY USE BY POOR HOUSEHOLDS
VATMP        VAT PAID ON IMPORTS USE BY POOR HOUSEHOLDS

PCON(I)     CONSUMER PRICE PER SECTOR INCLUDING VAT {p} 
PCONC       CONSUMER PRICE PER ENERGY SECTOR INCLUDING VAT
PCONF       CONSUMER PRICE PER ENERGY SECTOR INCLUDING VAT
PCONFe
PCONE       CONSUMER PRICE PER ENERGY SECTOR INCLUDING VAT

*====== Employment =====

EMPL(I)     LOW EDUCATED EMPLOYMENT PER SECTOR{q}
EMPLE(EN)   LOW EDUCATED EMPLOYMENT PER ENERGY SECTOR
EMPH(I)     HIGH EDUCATED EMPLOYMENT PER SECTOR
EMPHE(EN)   HIGH EDUCATED EMPLOYMENT PER ENERGY SECTOR
EMPLTOT     TOTAL LOW EDUCATED EMPLOYMENT
EMPHTOT     TOTAL HIGH EDUCATED EMPLOYMENT
EMPTOT      TOTAL EMPLOYMENT
EMPLBAR     TOTAL LOW EDUCATED EMPLOYMENT (POTENTIAL)
EMPHBAR     TOTAL HIGH EDUCATED EMPLOYMENT (POTENTIAL)
EMPBAR      TOTAL EMPLOYMENT (POTENTIAL)

EMPHGRO(I)    GROWTH OF HIGH EDUCATED EMPLOYMENT PER SECTOR{c}
EMPHEGRO(EN)  GROWTH OF HIGH EDUCATED EMPLOYMENT PER ENERGY SECTOR
EMPLGRO(I)    GROWTH OF LOW EDUCATED EMPLOYMENT PER SECTOR
EMPLEGRO(EN)  GROWTH OF LOW EDUCATED EMPLOYMENT PER ENERGY SECTOR

WAGEL(I)        LOW EDUCATED WAGE RATE IN PRODUCTION SECTORS{p}
WAGELE(EN)       LOW EDUCATED WAGE RATE IN ENERGY SECTOR
WAGEH(I)
WAGEHE(EN)
WAGELL
WAGELH

POSITIVE VARIABLES   X, XEE, DUME, ED,PT, PC, CLEAN, ABAT, PEM, DHP,DHR,VATR, VATP, 
                     PLL, PYT, PPU, PVT, PQT, PRT, PLLE, PYTE, PPUE, PVTE, PQTE, PRTE,
                     PE, IMPENS, TOTEM;

*----------------------------------------------------------------------------

 EQUATIONS

*======= DEMAND ======

 DEFDTGVBAR      LES government
 DEFDTHRBAR      LES households
 DEFDTHPBAR

 DEFGOVTOT       DEFINITION OF TOTAL GOVERNMENT DEMAND
 CCONSGV         GOVERNMENT CONSUMPTION DEMAND FOR COALe
 CCONSHR         HOUSEHOLD CONSUMPTION DEMAND FOR COALe
 CCONSHP
 CDEMAND(I)      SECTOR DEMAND FOR COALe
CDEMANDE(EN)
 CEN             TOTAL DEMAND FOR COALe FROM NON-ENERGY SECTORS,

 CTOT            DEFINITION OF TOTAL COALe DEMAND
 ECONSHR         HOUSEHOLD CONSUMPTION DEMAND FOR COALe
 ECONSHP
 ECONSGV         GOVERNMENT CONSUMPTION DEMAND FOR ELECTRICITY
 EDEMAND(I)      SECTOR DEMAND FOR ELECTRICITY
EDEMANDE(EN)
 EEN             TOTAL DEMAND FOR ELECTRICITY FROM NON-ENERGY SECTORS,

 ENDEMEN(EN)     DEFINITION OF TOTAL ENERGY DEMAND IN ENERGY SECTORS
 ENDEMSEC(I,EN)  DEFINITION OF NON-ENERGY DEMAND IN ENERGY SECTORS
 ENDEMTOT(I)     TOTAL NON-ENERGY DEMAND IN ENERGY SECTORS
 ENSEC(EN,EN1)   DEFINITION OF ENERGY DEMAND IN ENERGY SECTORS
 ENSECC(EN,EN1) 
 ENSECF(EN,EN1) 
 ENSECFe(EN,EN1) 
 ENSECE(EN,EN1) 
 ETOT            DEFINITION OF TOTAL ELECTRICITY DEMAND
 FCONSHR         HOUSEHOLD CONSUMPTION DEMAND FOR FUELS
 FCONSHP
 FeCONSHR
 FeCONSHP
 FCONSGV         CONSUMPTION DEMAND FOR FUELS
 FeCONSGV
 FDEMAND(I)      SECTOR DEMAND FOR FUELS
FDEMANDE(EN)
 FeDEMAND(I)
FeDEMANDE(EN)

 FEN             TOTAL DEMAND FOR FUELe FROM NON-ENERGY SECTORS,

 FeEN

 FTOT            DEFINITION OF TOTAL FUELe DEMAND
 FeTOT

 KDEMAND(I)      SECTOR DEMAND FOR CAPITAL
KDEMANDE(EN)
 LLDEMAND(I)     SECTOR DEMAND FOR LOW EDUCATED LABOUR
LLDEMANDE(EN)
 LHDEMAND(I)     SECTOR DEMAND FOR HIGH EDUCATED LABOUR
LHDEMANDE(EN)
 MCONSHR         HOUSEHOLD CONSUMPTION DEMAND FOR IMPORTS
 MCONSHP
 MCONSGV         GOVERNEMENT CONSUMPTION DEMAND FOR IMPORTS
 NCONSHR(I)      HOUSEHOLD CONSUMPTION DEMAND FOR GOODS OF NON-PRICE-TAKING SECTORS
 NCONSHP(I)
 NCONSGV(I)      GOVERNMENT CONSUMPTION DEMAND FOR GOODS OF NON-PRICE-TAKING SECTORS
 TCONSHR(I)      HOUSEHOLD CONSUMPTION DEMAND FOR GOODS OF PRICE-TAKING SECTORS
 TCONSHP(I)
 TCONSGV(I)      GOVERNMENT CONSUMPTION DEMAND FOR GOODS OF PRICE-TAKING SECTORS

 UTILITYGV
 UTILITYHR
 UTILITYHP
 INDUTILHP
 INDUTILHR

*======= TRADE ======

 COMPL(I)        COMPLEMENTARY CONDITION FOR TRADABLES
 DEFIMPORT       DEFINITION OF TOTAL IMPORTS
 DEFIMPTOT       DEFINITION OF TOTAL IMPORT (INCLUDING FUEL)
 DEFIMPORTS      DEFINITION OF VALUE OF TOTAL IMPORTS
 DEFEXPORTS      DEFINITION OF VALUE OF TOTAL EXPORTS
 DEFEXPN(I)      DEFINITION OF THE EXPORTS OF THE NON-PRICE TAKING SECTORS 
DEFEXPNE
DEFEXPNC
DEFEXPNF
DEFEXPNFe 
DEFEXPC

*======= PRICES AND COSTS ======

 DEFMCE          DEFINITION OF MARGINAL COST IN ELECTRICITY SECTOR
DEFEX 
DEFEXH
 DEFEXHR
 DEFEXHP
 DEFEXGOV

DEFINGV
 DEFIN           DEFINITION OF TOTAL INCOME (HOUSEHOLDS)    
 DEFINH  
 DEFINHR         DEFINITION OF THE INCOME OF HOUSEHOLDS
 DEFINHP
 DEFPCC          DEFINITION OF CONSUMER PRICE OF COALe
 DEFPCT          DEFINITION OF SECTOR PRICE OF COALe
DEFPCE
 DEFPE           PRICE EQUATION OF ELECTRICITY
DEFPEC
DEFPET(J)
DEFPEE(EN)
 DEFPFC          DEFINITION OF CONSUMER PRICE OF FUELS
 DEFPFeC
 DEFPFT         DEFINITION OF SECTOR PRICE OF FUELS
DEFPFE
 DEFPFeT
DEFPFeE
 PRLMCT(J)       PRODUCER PRICE - MARGINAL COST INEQUALITY CONDITION
 PRZMCT(I)       PRODUCER PRICE - MARGINAL COST INEQUALITY CONDITION
PRCMCE
PRFMCE
PRFeMCE
PREMCE
 DEFMCT(J)       DEFINITION OF MARGINAL COST IN PRODUCING SECTORS
DEFMCEN(EN)
DEFPNE(I,EN)
DEFPNT(I,J)
DEFPNC(I)

*======= DEFINITIONS ======

 DEFPPU(J)       DEFINITION OF PPU
 DEFPQT(J)       DEFINITION OF PQT
 DEFPRT(J)       DEFINITION OF PRT
 DEFPVT(J)       DEFINITION OF PVT
 DEFPYT(J)       DEFINITION OF PYT
 DEFPLT(J)       DEFINITION OF PLT

 DEFPPUE(EN)       DEFINITION OF PPU
 DEFPQTE(EN)       DEFINITION OF PQT
DEFPQTEC(EN)
DEFPQTEF(EN)
 DEFPRTE(EN)       DEFINITION OF PRT
DEFPRTEF(EN)
 DEFPVTE(EN)       DEFINITION OF PVT
DEFPVTEE(EN)
 DEFPYTE(EN)       DEFINITION OF PYT
 DEFPLTE(EN)       DEFINITION OF PLT

 DEFDPQDPFT(J)   DEFINITION OF DPQDPFT
 DEFDPQDPCT(J)   DEFINITION OF DPQDPCT
 DEFDPRDPFeT(J)  DEFINITION OF DPRDPFeT
 DEFDPUDPKT(J)   DEFINITION OF DPUDPKT
 DEFDPUDPVT(J)   DEFINITION OF DPUDPVT
 DEFDPVDPET(J)   DEFINITION OF DPVDPET
 DEFDPRDPQT(J)   DEFINITION OF DPRDPQT
 DEFDPVDPRT(J)   DEFINITION OF DPVDPRT
 DEFDPYDPPU(J)   DEFINITION OF DPYDPPU
 DEFDPLLPLL(J)   DEFINITION OF DPLLPLL
 DEFDPLLPYT(J)   DEFINITION OF DPLLPYT
 DEFDPYDPLH(J)   DEFINITION OF DPYDPLH

 DEFDPQDPFTE(EN)   DEFINITION OF DPQDPFT
DEFDPQDPFTEF(EN)
 DEFDPQDPCTE(EN)   DEFINITION OF DPQDPCT
DEFDPQDPCTEC(EN)
 DEFDPRDPFeTE(EN)  DEFINITION OF DPRDPFeT
DEFDPRDPFeTEF(EN)
 DEFDPUDPKTE(EN)   DEFINITION OF DPUDPKT
 DEFDPUDPVTE(EN)   DEFINITION OF DPUDPVT
 DEFDPVDPETE(EN)   DEFINITION OF DPVDPET
DEFDPVDPETEE(EN)
 DEFDPRDPQTE(EN)   DEFINITION OF DPRDPQT
 DEFDPVDPRTE(EN)   DEFINITION OF DPVDPRT
 DEFDPYDPPUE(EN)   DEFINITION OF DPYDPPU
 DEFDPLLPLLE(EN)   DEFINITION OF DPLLPLL
 DEFDPLLPYTE(EN)   DEFINITION OF DPLLPYT
 DEFDPYDPLHE(EN)   DEFINITION OF DPYDPLH

 DEFDPXDPCT(J)   DEFINITION OF DPXDPCT
 DEFDPXDPFeT(J)  DEFINITION OD DPXDPFeT
 DEFDPXDPET(J)   DEFINITION OF DPXDPET
 DEFDPXDPFT(J)  DEFINITION OF DPXDPFT
 DEFDPXDPKT(J)   DEFINITION OF DPXDPKT
 DEFDPXDPYT(J)   DEFINITION OF DPXDPYT
 DEFDPXDPLL(J)   DEFINITION OF DPXDPLL

 DEFDPXDPCTE(EN)   DEFINITION OF DPXDPCT
 DEFDPXDPFeTE(EN)  DEFINITION OD DPXDPFeT
 DEFDPXDPETE(EN)   DEFINITION OF DPXDPET
 DEFDPXDPFTE(EN)  DEFINITION OF DPXDPFT
 DEFDPXDPKTE(EN)   DEFINITION OF DPXDPKT
 DEFDPXDPYTE(EN)   DEFINITION OF DPXDPYT
 DEFDPXDPLLE(EN)   DEFINITION OF DPXDPLL



*======= MARKET CLEARING ======

 CMARKET         MARKET CLEARING FOR COALe
 EMARKET         MARKET CLEARING FOR ELECTRICITY
 FMARKET         MARKET CLEARING FOR FUELS
 FeMARKET

 KMARKET         MARKET CLEARING FOR CAPITAL
 LLMARKET        MARKET CLEARING FOR LOW EDUCATED LABOUR
 LHMARKET        MARKET CLEARING FOR HIGH EDUCATED LABOUR
 NMARKET(I)      MARKET CLEARING FOR NON-PRICE-TAKING SECTORS
 TMARKET(I)      MARKET CLEARING FOR PRICE-TAKING SECTORS
 DEFXGRO(I)      DEFINITION OF OUTPUT GROWTH FROM BASE YEAR PER SECTOR
 DEFXEEGRO(EN)   DEFINITION OF OUTPUT GROWTH FROM BASE YEAR IN ENERGY SECTORS
 DEFXPGRO(I)
 DEFXPEEGRO(EN)

*======= INPUT ======

 CAPITSEC(EN)    DEFINITION OF CAPITAL INPUT IN ENERGY SECTORS
 CAPITTOT        TOTAL CAPITAL INPUT IN ENERGY SECTORS

 DEFKFTAX        CAPITAL RE-ALLOCATION (FROM ENERGY TAX)
 DEFKPFTAX      
DEFRATE(I)
DEFRATEE(EN)
DEFRATTOT
DEFTRWAL(I)
DEFTRWALE(EN)
DEFTRWALTOT
 
 IMPORSEC(EN)    DEFINITION OF IMPORT IN ENERGY SECTORS
 IMPOREN
 LABOURLSEC(EN)  DEFINITION OF LOW EDUCATED LABOUR INPUT IN ENERGY SECTORS
 LABOURLTOT      TOTAL LOW EDUCATED LABOUR INPUT IN ENERGY SECTORS
 LABOURHSEC(EN)  DEFINITION OF HIGH EDUCATED LABOUR INPUT IN ENERGY SECTORS
 LABOURHTOT      TOTAL HIGH EDUCATED LABOUR INPUT IN ENERGY SECTORS

*======= ENVIRONMENT ======

 COMPCLEAN   COMPLEMENTARY CONDITION FOR EMISSION CLEANING
 COMPLEM     COMPLEMENTARY CONDITION FOR THE PERMITS MARKET
 DEFE            QUANTITY EQUATION OF ELECTRICITY
 DEFEMIS     DEFINITION OF TOTAL EMISSIONS
 DEFEMRED    DEFINITION OF EMISSION REDUCTION       
 DEFENEMIS(EN) DEFINITIOM OF EMISSIONS PER ENERGY SECTOR   
 EMCOSTCO2       ABATEMENT COST FOR CO2
 EMISCALC(J)  CALCULATION OF EMISSION PER SECTOR
 EMISHOUS  CALCULATION OF HOUSEHOLD EMISSION 
 FUNDCALC    CALCULATION OF EMISSION TAX YIELD
 PERMARKET   MARKET CLEARING FOR EMISSION PERMITS
 PROFCLEAN   ZERO PROFIT CONDITION FOR EMISSION CLEANING

 COMPCLEANJ(J)
 COMPCLEANE(EN)
 DEFCLEANSEC
 WARUNEK(J)
 WARUNEKE(EN)

*======= Value Added Tax system ======

 DEFVAT           TOTAL VAT income for the government

 DEFVATN(I)       VAT PER SECTOR
 DEFVATC          VAT COAL SECTOR
 DEFVATF          VAT FUEL SECTOR
 DEFVATFe
 DEFVATE          VAT ELECTRICITY AND HEAT SECTOR
 DEFVATM          VAT ON IMPORTS

 DEFVATNR(I)       VAT PER SECTOR
 DEFVATCR          VAT COAL SECTOR
 DEFVATFR          VAT FUEL SECTOR
 DEFVATFeR
 DEFVATER          VAT ELECTRICITY AND HEAT SECTOR
 DEFVATMR          VAT ON IMPORTS

 DEFVATNP(I)       VAT PER SECTOR
 DEFVATCP          VAT COAL SECTOR
 DEFVATFP          VAT FUEL SECTOR
 DEFVATFeP
 DEFVATEP          VAT ELECTRICITY AND HEAT SECTOR
 DEFVATMP          VAT ON IMPORTS
 DEFVATR
 DEFVATP

*======= Prices ======

DEFPCONT(I)     CONSUMER PRICE INCLUDING VAT
DEFPCONN(I)
DEFPCONC        CONSUMER PRICE INCLUDING VAT COAL
DEFPCONF        CONSUMER PRICE INCLUDING VAT FUEL
DEFPCONFe
DEFPCONE        CONSUMER PRICE INCLUDING VAT ELECTRICITY
DEFPF           DEFINITION OF FUEL PRICE

DEFPLL          PRICE OF LOW EDUCATED LABOUR
DEFPLLE

DEFPLLOW        PRICE OF LOW EDUCATED LABOUR
DEFPLLOWX       PRICE OF LOW EDUCATED LABOUR
DEFLABSUB       TOTAL SUBSIDIES PAID BY THE GOVERNMENT ON LABOUR COSTS
DEFSUBTOT       TOTAL SUBSIDIES PAID BY THE GOVERNMENT
DEFLL           DEMAND FOR LOW EDUCATED LABOUR
DEFUNL          NUMBER OF UNEMPLOYED
DEFLLOW(I)
DEFPLHIGH       PRICE OF HIGH EDUCATED LABOUR
DEFPLHIGHX      PRICE OF HIGH EDUCATED LABOUR EXCLUDING SUBSIDY
DEFLH           DEMAND FOR HIGH EDUCATED LABOUR
DEFUNH          NUMBER OF UNEMPLOYED HIGH EDUCATED
DEFLHIGH(I)
DEFHIGH
DEFLABHTAX(I)   TAXES IN LABOUR COSTS HIGH EDUCATED
DEFLABHSOC(I)   SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS HIGH EDUC.
DEFLABLTAX(I)   TAXES IN LABOUR COSTS LOW EDUCATED
DEFLABLSOC(I)   SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS LOW EDUC.
 DEFLABHTAXE(EN) TAXES IN LABOUR COSTS HIGH EDUCATED
DEFLABHSOCE(EN) SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS HIGH EDUC.
DEFLABLTAXE(EN) TAXES IN LABOUR COSTS LOW EDUCATED
DEFLABLSOCE(EN) SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS LOW EDUC.
DEFLABTAXTOT    TAXES IN LABOUR COSTS
DEFLABSOCTOT    SOCIAL SECURITY CONTRIBUTIONS IN LABOUR COSTS

DEFWAGEL(I)
DEFWAGELE(EN)
DEFWAGEH(I)
DEFWAGEHE(EN)
DEFWAGELL
DEFWAGELH

*====== Employment =====

DEFEMPL         LOW EDUCATED EMPLOYMENT PRODUCTION SECTORS
DEFEMPLE        LOW EDUCATED EMPLOYMENT ENERGY SECTORS
DEFEMPH         HIGH EDUCATED EMPLOYMENT PRODUCTION SECTORS
DEFEMPHE        HIGH EDUCATED EMPLOYMENT ENERGY SECTORS
DEFEMPLT        TOTAL EMPLOYMENT PRODUCTION SECTORS
DEFEMPHT        TOTAL EMPLOYMENT ENERGY SECTORS
DEFEMPT         TOTAL EMPLOYMENT
DEFEMPLB
DEFEMPHB
DEFEMPB
DEFEMPHG
DEFEMPHEG
DEFEMPLG
DEFEMPLEG

*======= MISCELLANEOUS ======

ACCOUNT         EQUILIBRIUM ON INTERNATIONAL MARKET
DEFGNP          DEFINITION OF GNP
DEFGDP          DEFINITION OF GDP - expenditures method
DEFPKB          
DEFPKKB
DEFVA           TOTAL VALUE ADDED
DUMELEC         UNDER UTILIZATION OF ELECTRICITY
DEFGNP0         DEFINITION OF GNP0 (BASE YEAR LEVEL)
DEFGDP0         DEFINITION OF GDP0 (BASE YEAR LEVEL)
DEFPKB0       
DEFPKKB0

;



* Definition of the implicit prices of composites of 2 "factors of production"
* CES-production functions


 DEFPLL(J)..     PLL(J)=E=(DELLL(J)*PLLOW**(1-SIGMA("LL",J))+
                          (1-DELLL(J))*PYT(J)**
                          (1-SIGMA("LL",J)))**(1/(1-SIGMA("LL",J)));


 DEFPYT(J)..     PYT(J)=E=(DELYL(J)*PLHIGH**(1-SIGMA("YL",J))+
                          (1-DELYL(J))*PPU(J)**
                          (1-SIGMA("YL",J)))**(1/(1-SIGMA("YL",J)));

 DEFPPU(J)..     PPU(J)=E=(DELUK(J)*(PK/(1+SUBK))**(1-SIGMA("UK",J))+(1-
                          DELUK(J))*PVT(J)**
                          (1-SIGMA("UK",J)))**(1/(1-SIGMA("UK",J)));

 DEFPVT(J)..     PVT(J)=E=(DELVE(J)*PET(J)**(1-SIGMA("VE",J))+(1-DELVE(J))*PRT(J)**
                          (1-SIGMA("VE",J)))**(1/(1-SIGMA("VE",J)));

 DEFPRT(J)..     PRT(J)=E=(DELRFe(J)*PFeT(J)**(1-SIGMA("RFe",J))+(1-DELRFe(J))*PQT(J)**
                          (1-SIGMA("RFe",J)))**(1/(1-SIGMA("RFe",J)));

 DEFPQT(J)..     PQT(J)=E=(DELQF(J)*PCT(J)**(1-SIGMA("QF",J))+(1-DELQF(J))*PFT(J)**
                          (1-SIGMA("QF",J)))**(1/(1-SIGMA("QF",J)));




 DEFPLLE(EN)..     PLLE(EN)=E=(EDELLL(EN)*PLLOW**(1-SIGME("LL",EN))+
                          (1-EDELLL(EN))*PYTE(EN)**
                          (1-SIGME("LL",EN)))**(1/(1-SIGME("LL",EN)));

 DEFPYTE(EN)..     PYTE(EN)=E=(EDELYL(EN)*PLHIGH**(1-SIGME("YL",EN))+
                          (1-EDELYL(EN))*PPUE(EN)**
                          (1-SIGME("YL",EN)))**(1/(1-SIGME("YL",EN)));

 DEFPPUE(EN)..     PPUE(EN)=E=(EDELUK(EN)*(PK/(1+SUBK))**(1-SIGME("UK",EN))+(1-
                          EDELUK(EN))*PVTE(EN)**
                          (1-SIGME("UK",EN)))**(1/(1-SIGME("UK",EN)));

 DEFPVTE(EN)..     PVTE(EN)=E=(EDELVE(EN)*PEEN(EN)**(1-SIGME("VE",EN))+(1-EDELVE(EN))*PRTE(EN)**
                          (1-SIGME("VE",EN)))**(1/(1-SIGME("VE",EN)));

 DEFPRTE(EN)..    PRTE(EN)=E=(EDELRFe(EN)*PFeEN(EN)**(1-SIGME("RFe",EN))+(1-EDELRFe(EN))*
                                PQTE(EN)**(1-SIGME("RFe",EN)))**(1/(1-SIGME("RFe",EN)));

 DEFPQTE(EN)..    PQTE(EN)=E=(EDELQF(EN)*PCEN(EN)**(1-SIGME("QF",EN))+(1-EDELQF(EN))*
                                PFEN(EN)**(1-SIGME("QF",EN)))**(1/(1-SIGME("QF",EN)));



* Definition of marginal (and average) cost of production
* CEMS, FEMS,K,E  has been already included in PLL


 DEFMCT(J)..     MCT(J)=E=A(J)*PLL(J)/EXP(LAMBDA(J)*TIME)+
                          SUM(I$N(I),PN(I)*AT(I,J)) +
                          PROF(J)
                          +TTAX(J)+ PM*AM(J)+(PEM+EMISTAX)*PEMS(J);

* Definition of marginal cost in the electricity sector


 DEFMCEN(EN)..     MCEN(EN)=E=AE(EN)*PLLE(EN)/EXP(LAMBDE(EN)*TIME)+
                          SUM(I$N(I), PN(I)*EIO(I,EN))+
                          PROFEN(EN)+
                          ENTAX(EN) 
                          +PM*ME(EN)$FU(EN)+PFX*ME(EN)$FL(EN)+
                          (PEM+EMISTAX)*PEME(EN);


* Profit maximization

 PRZMCT(I)$N(I)..   PN(I)-MCT(I)/(1+SUB(I))=E=0;



 PRCMCE..  PC-MCEN("COALe")/(1+SUBC)=E=0;
 PRFMCE..  PFX-MCEN("FUEL")/(1+SUBF)=E=0;
 PRFeMCE.. PFe-MCEN("FUELe")/(1+SUBFe)=E=0;
 PREMCE..  PE-MCEN("ELECe")/(1+SUBE)=E=0;


* Derivatives of the implicit prices with respect to their "elements"
* used for calculation of marginal costs

 DEFDPLLPLL(J)..    DPLLPLL(J)=E=DELLL(J)*(PLL(J)**SIGMA("LL",J))*(PLLOW**(-SIGMA("LL",J)));

 DEFDPLLPYT(J)..    DPLLPYT(J)=E=(1-DELLL(J))*(PLL(J)**SIGMA("LL",J))*(PYT(J)**(-SIGMA("LL",J)));


 DEFDPYDPLH(J)..    DPYDPLH(J)=E=DELYL(J)*(PYT(J)**SIGMA("YL",J))
                                       *(PLHIGH**(-SIGMA("YL",J)));

 DEFDPYDPPU(J)..    DPYDPPU(J)=E=(1-DELYL(J))*(PYT(J)**SIGMA("YL",J))
                                 *(PPU(J)**(-SIGMA("YL",J)));


 DEFDPUDPKT(J)..    DPUDPKT(J)=E=DELUK(J)*(PPU(J)**SIGMA("UK",J))
                                *(PK**(-SIGMA("UK",J)));

 DEFDPUDPVT(J)..    DPUDPVT(J)=E=(1-DELUK(J))*(PPU(J)**SIGMA("UK",J))
                                 *(PVT(J)** (-SIGMA("UK",J)));


 DEFDPVDPET(J)..    DPVDPET(J)=E=DELVE(J)*(PVT(J)**SIGMA("VE",J))
                                 *(PET(J)**(-SIGMA("VE",J)));

 DEFDPVDPRT(J)..    DPVDPRT(J)=E=(1-DELVE(J))*(PVT(J)**SIGMA("VE",J))
                                 *(PRT(J)**(-SIGMA("VE",J)));

 
 DEFDPRDPQT(J)..    DPRDPQT(J)=E=(1-DELRFe(J))*(PRT(J)**SIGMA("RFe",J))
                                  *(PQT(J)**(-SIGMA("RFe",J)));

 DEFDPRDPFeT(J)..   DPRDPFeT(J)=E=DELRFe(J)*(PRT(J)**SIGMA("RFe",J))
                                  *(PFeT(J)**(-SIGMA("RFe",J)));


 DEFDPQDPFT(J)..    DPQDPFT(J)=E=(1-DELQF(J))*(PQT(J)**SIGMA("QF",J))
                                 *(PFT(J)**(-SIGMA("QF",J)));


 DEFDPQDPCT(J)..    DPQDPCT(J)=E=DELQF(J)*(PQT(J)**SIGMA("QF",J))
                                 *(PCT(J)**(-SIGMA("QF",J)));





 DEFDPLLPLLE(EN)..    DPLLPLLE(EN)=E=EDELLL(EN)*(PLLE(EN)**SIGME("LL",EN))
                                        *(PLLOW**(-SIGME("LL",EN)));

 DEFDPLLPYTE(EN)..    DPLLPYTE(EN)=E=(1-EDELLL(EN))*(PLLE(EN)**SIGME("LL",EN))
                                        *(PYTE(EN)**(-SIGME("LL",EN)));


 DEFDPYDPLHE(EN)..    DPYDPLHE(EN)=E=EDELYL(EN)*(PYTE(EN)**SIGME("YL",EN))
                                       *(PLHIGH**(-SIGME("YL",EN)));

 DEFDPYDPPUE(EN)..    DPYDPPUE(EN)=E=(1-EDELYL(EN))*(PYTE(EN)**SIGME("YL",EN))
                                 *(PPUE(EN)**(-SIGME("YL",EN)));


 DEFDPUDPKTE(EN)..     DPUDPKTE(EN)=E=EDELUK(EN)*(PPUE(EN)**SIGME("UK",EN))
                                *(PK**(-SIGME("UK",EN)));

 DEFDPUDPVTE(EN)..    DPUDPVTE(EN)=E=(1-EDELUK(EN))*(PPUE(EN)**SIGME("UK",EN))
                                 *(PVTE(EN)** (-SIGME("UK",EN)));


 DEFDPVDPETE(EN)..    DPVDPETE(EN)=E=EDELVE(EN)*(PVTE(EN)**SIGME("VE",EN))
                                 *(PEEN(EN)**(-SIGME("VE",EN)));

 DEFDPVDPRTE(EN)..            DPVDPRTE(EN)=E=(1-EDELVE(EN))*(PVTE(EN)**SIGME("VE",EN))
                                 *(PRTE(EN)**(-SIGME("VE",EN)));

 
 DEFDPRDPFeTE(EN)..   DPRDPFeTE(EN)=E=EDELRFe(EN)*(PRTE(EN)**SIGME("RFe",EN))
                                  *(PFeEN(EN)**(-SIGME("RFe",EN)));

 DEFDPRDPQTE(EN)..            DPRDPQTE(EN)=E=(1-EDELRFe(EN))*(PRTE(EN)**SIGME("RFe",EN))
                                  *(PQTE(EN)**(-SIGME("RFe",EN)));


 DEFDPQDPFTE(EN)..     DPQDPFTE(EN)=E=(1-EDELQF(EN))*(PQTE(EN)**SIGME("QF",EN))
                                 *(PFEN(EN)**(-SIGME("QF",EN)));


 DEFDPQDPCTE(EN)..    DPQDPCTE(EN)=E=EDELQF(EN)*(PQTE(EN)**SIGME("QF",EN))
                                 *(PCEN(EN)**(-SIGME("QF",EN)));






 DEFDPXDPCT(J)..    DPXDPCT(J)=E=A(J)*DPLLPYT(J)*DPYDPPU(J)*DPUDPVT(J)*DPVDPRT(J)*DPRDPQT(J)
                                *DPQDPCT(J);

 DEFDPXDPFT(J)..    DPXDPFT(J)=E=A(J)*DPLLPYT(J)*DPYDPPU(J)*DPUDPVT(J)*DPVDPRT(J)*DPRDPQT(J)
                                *DPQDPFT(J);

 DEFDPXDPFeT(J)..   DPXDPFeT(J)=E=A(J)*DPLLPYT(J)*DPYDPPU(J)*DPUDPVT(J)*DPVDPRT(J)*DPRDPFeT(J);

 DEFDPXDPET(J)..    DPXDPET(J)=E=A(J)*DPLLPYT(J)*DPYDPPU(J)*DPUDPVT(J)*DPVDPET(J);

 DEFDPXDPKT(J)..    DPXDPKT(J)=E=A(J)*DPLLPYT(J)*DPYDPPU(J)*DPUDPKT(J);

 DEFDPXDPYT(J)..    DPXDPYT(J)=E=A(J)*DPLLPYT(J)*DPYDPLH(J);
 
 DEFDPXDPLL(J)..    DPXDPLL(J)=E=A(J)*DPLLPLL(J);



 DEFDPXDPCTE(EN)..    DPXDPCTE(EN)=E=AE(EN)*DPLLPYTE(EN)*DPYDPPUE(EN)*DPUDPVTE(EN)*DPVDPRTE(EN)
                                *DPRDPQTE(EN)*DPQDPCTE(EN);

 DEFDPXDPFTE(EN)..    DPXDPFTE(EN)=E=AE(EN)*DPLLPYTE(EN)*DPYDPPUE(EN)*DPUDPVTE(EN)*DPVDPRTE(EN)
                                *DPRDPQTE(EN)*DPQDPFTE(EN);

 DEFDPXDPFeTE(EN)..   DPXDPFeTE(EN)=E=AE(EN)*DPLLPYTE(EN)*DPYDPPUE(EN)*DPUDPVTE(EN)*DPVDPRTE(EN)
                                *DPRDPFeTE(EN);

 DEFDPXDPETE(EN)..    DPXDPETE(EN)=E=AE(EN)*DPLLPYTE(EN)*DPYDPPUE(EN)*DPUDPVTE(EN)*DPVDPETE(EN);

 DEFDPXDPKTE(EN)..    DPXDPKTE(EN)=E=AE(EN)*DPLLPYTE(EN)*DPYDPPUE(EN)*DPUDPKTE(EN);

 DEFDPXDPYTE(EN)..    DPXDPYTE(EN)=E=AE(EN)*DPLLPYTE(EN)*DPYDPLHE(EN);
 
 DEFDPXDPLLE(EN)..    DPXDPLLE(EN)=E=AE(EN)*DPLLPLLE(EN);




* Sectoral demand for input factors

 CDEMAND(J)..      CT(J) =E=DPXDPCT(J)*X(J)/EXP(ENEFF(J,"COALe")*TIME);

 EDEMAND(J)..      ET(J) =E=DPXDPET(J)*X(J)/EXP(ENEFF(J,"ELECe")*TIME);

 FDEMAND(J)..      FT(J) =E=DPXDPFT(J)*X(J)/EXP(ENEFF(J,"FUEL")*TIME);

 FeDEMAND(J)..     FeT(J) =E=DPXDPFeT(J)*X(J)/EXP(ENEFF(J,"FUELe")*TIME);

 KDEMAND(J)..      KT(J) =E=DPXDPKT(J)*X(J)/EXP(LAMBDAK(J)*TIME);

 LHDEMAND(J)..     LHT(J)=E=DPXDPYT(J)*X(J)/EXP(LAMBDALH(J)*TIME);

 LLDEMAND(J)..     LLT(J)=E=DPXDPLL(J)*X(J)/EXP(LAMBDALL(J)*TIME);


 CDEMANDE(EN)..      CTE(EN) =E=DPXDPCTE(EN)*XEE(EN)/EXP(ENEFFE(EN,"COALe")*TIME);

 EDEMANDE(EN)..      ETE(EN) =E=DPXDPETE(EN)*XEE(EN)/EXP(ENEFFE(EN,"ELECe")*TIME);

 FDEMANDE(EN)..      FTE(EN) =E=DPXDPFTE(EN)*XEE(EN)/EXP(ENEFFE(EN,"FUEL")*TIME);

 FeDEMANDE(EN)..     FeTE(EN) =E=DPXDPFeTE(EN)*XEE(EN)/EXP(ENEFFE(EN,"FUELe")*TIME);

 KDEMANDE(EN)..      KTE(EN) =E=DPXDPKTE(EN)*XEE(EN)/EXP(LAMBDAKE(EN)*TIME);

 LHDEMANDE(EN)..     LHTE(EN)=E=DPXDPYTE(EN)*XEE(EN)/EXP(LAMBDALHE(EN)*TIME);

 LLDEMANDE(EN)..     LLTE(EN)=E=DPXDPLLE(EN)*XEE(EN)/EXP(LAMBDALLE(EN)*TIME);



* Sector and consumer prices of coal, gas and fuel depend on the emission-
* coefficients and the price of an emission permit.

 DEFPF..            PF=E=PFX*(1+FTAX)+NFTAX;

 DEFPFC..           PFC=E=PF;

 DEFPFT(J)..        PFT(J)=E=PF;

 DEFPCC..           PCC=E=PC*(1+CTAX)+NCTAX+(EMISTAX+PEM)*CEMH;

 DEFPCT(J)..        PCT(J)=E=PC*(1+CTAX)+NCTAX+(EMISTAX+PEM)*SUM(EN,CEMS(EN,J));

 DEFPFeC..          PFeC=E=PFe*(1+FeTAX)+NFeTAX;

 DEFPFeT(J)..       PFeT(J)=E=PFe*(1+FeTAX)+NFeTAX;

 DEFPEC..           PEC=E=PE*(1+ETAX)+NETAX;

 DEFPET(J)..        PET(J)=E=PE*(1+ETAX)+NETAX;



 DEFPFE(EN)..        PFEN(EN)=E=PF;

 DEFPCE(EN)..        PCEN(EN)=E=PC*(1+CTAX)+NCTAX+(EMISTAX+PEM)*SUM(EN1,(CEME(EN1,EN)));

 DEFPFeE(EN)..       PFeEN(EN)=E=PFe*(1+FeTAX)+NFeTAX;

 DEFPEE(EN)..        PEEN(EN)=E=PE*(1+ETAX)+NETAX;


 
 DEFPNC(I)..          PNC(I)=E=PN(I)*(1+NTAX(I))+NNTAX(I)+(EMISTAX+PEM)*NEMH(I);

 DEFPNT(I,J)..        PNT(I,J)=E=PN(I)*(1+NTAX(I))+NNTAX(I)+(EMISTAX+PEM)*NEMS(I,J);

 DEFPNE(I,EN)..        PNEN(I,EN)=E=PN(I)*(1+NTAX(I))+NNTAX(I)+(EMISTAX+PEM)*NEME(I,EN);

 

* Market equilibrium conditions

 NMARKET(I)$N(I)..  X(I)=E=SUM(J,AT(I,J)*X(J))+ED(I)+DGV(I)+DHR(I)+DHP(I)+ZN(I)+EXDEM(I);


EMARKET..            E=E=SUM(J,ET(J))+SUM(EN,ETE(EN))
                               +DEGV+DEHR+DEHP+EXDEME+DUME
                              +ZNE;

CMARKET..            C=E=SUM(J,CT(J))+SUM(EN,CTE(EN))
                               +DCGV+DCHR+DCHP+EXDEMC
                              +ZNC; 

FMARKET..            F=E=SUM(J,FT(J))+SUM(EN,FTE(EN))
                               +DFGV+DFHR+DFHP+EXDEMF
                              +ZNF;

FeMARKET..           Fe=E=SUM(J,FeT(J))+SUM(EN,FeTE(EN))
                               +DFeGV+DFeHR+DFeHP+EXDEMFe
                              +ZNFe;                     


 EEN..              XEE("ELECe")=E=E;
 FEN..              XEE("FUEL")=E=F;
 FeEN..             XEE("FUELe")=E=Fe;
 CEN..              XEE("COALe")=E=C;


* Demand in energy sectors 
* XEE+EED it is total demand on energy products. 


 ENDEMSEC(I,EN)..   EDS(I,EN)=E=EIO(I,EN)*XEE(EN);

 ENDEMTOT(I)..      ED(I)=E=SUM(EN,EDS(I,EN));


 DUMELEC..          DUM1=E=DUME;


*-----------------------------------------------------------------------

* Domestic governmment final demand, given by a linear expenditure system

 DEFINGV..        INGV=E=SUM(J,TTAX(J)*X(J))+SUM(EN,ENTAX(EN)*XEE(EN))+VAT
                              +KPFTAX+LABTAXTOT; 


 DEFEXGOV..       EXGOV=L=INGV-SBAR-SUBTOT;


 DEFSUBTOT..    SUBTOT=E=LABSUB+SUM(I,SUB(I)*X(I)*PN(I))+
                         SUBK*K*PK/(1+SUBK)+SUBC*C*PC+SUBF*F*PF+
                         SUBFe*Fe*PFe+SUBE*E*PE
                          +INGV*(SUBHR+SUBHP);


 DEFIN..       IN=E=K*PK+LL*PLLOWX+LH*PLHIGHX+SUM(J,PROF(J)*X(J))
                           +SUM(EN,PROFEN(EN)*XEE(EN));

 DEFINH..      IN=E=INHR+INHP;
 DEFEXH..      EXH=E=EXHR+EXHP;
 DEFEX..       EXH=L=IN-LABTAXTOT-LABSOCTOT+INGV*(SUBHR+SUBHP);

 DEFINHR..       INHR+INGV*SUBHR=G=EXHR;
 DEFINHP..       INHP+INGV*SUBHP=G=EXHP;

* DEFHPTOT..       INHP=E=(IN-LABTAXTOT)/2+INGV*SUBHP;


 DEFGOVTOT..     EXGOV=E=SUM(I$N(I),DGV(I)*MCT(I))+DCGV*PC*(1+SUBC)+DFGV*PF*(1+SUBF)+DFeGV*PFe*(1+SUBFe)+DEGV*PE*(1+SUBE)+DMGV*PM;
 DEFEXHR..       EXHR=E=SUM(I$N(I),DHR(I)*PCON(I))+DCHR*PCONC+DFHR*PCONF+DFeHR*PCONFe+DEHR*PCONE;
 DEFEXHP..       EXHP=E=SUM(I$N(I),DHP(I)*PCON(I))+DCHP*PCONC+DFHP*PCONF+DFeHP*PCONFe+DEHP*PCONE;



* Domestic household final demand, given by a linear expenditure system

 DEFDTGVBAR..      DTGVBAR=E=
                            SUM(J$N(J),MCT(J)*DGVBAR(J))
                            +PE*(1+SUBE)*DEGVBAR
                            +PFX*(1+SUBF)*DFGVBAR
                            +PFe*(1+SUBFe)*DFeGVBAR
                            +PM*DMGVBAR
                            +PC*(1+SUBC)*DCGVBAR ;

 DEFDTHRBAR..      DTHRBAR=E=
                            SUM(J$N(J),PNC(J)*RR*ALFAVAT(J)*DHRBAR(J))
                            +PEC*RR*ALFAVATE("ELECe")*DEHRBAR
                            +PFC*RR*ALFAVATE("FUEL")*DFHRBAR
                            +PFeC*RR*ALFAVATE("FUELe")*DFeHRBAR
                            +PM*RR*ALFAVATM*DMHRBAR
                            +PCC*RR*ALFAVATE("COALe")*DCHRBAR ;

 DEFDTHPBAR..      DTHPBAR=E=
                            SUM(J$N(J),PNC(J)*RR*ALFAVAT(J)*DHPBAR(J))
                            +PEC*RR*ALFAVATE("ELECe")*DEHPBAR
                            +PFC*RR*ALFAVATE("FUEL")*DFHPBAR
                            +PFeC*RR*ALFAVATE("FUELe")*DFeHPBAR
                            +PM*RR*ALFAVATM*DMHPBAR
                            +PCC*RR*ALFAVATE("COALe")*DCHPBAR ;


 NCONSGV(I)$N(I).. DGV(I)=E=DGVBAR(I)+(BETGV(I)/MCT(I))*(EXGOV-DTGVBAR);


 NCONSHR(I)$N(I).. DHR(I)=E=DHRBAR(I)+(BETHR(I)/(PNC(I)*RR*ALFAVAT(I))) *(EXHR-DTHRBAR);
 NCONSHP(I)$N(I).. DHP(I)=E=DHPBAR(I)+(BETHP(I)/(PNC(I)*RR*ALFAVAT(I))) *(EXHP-DTHPBAR);


 CCONSGV..         DCGV=E=DCGVBAR+(BETCGV/(PC*(1+SUBC)))*(EXGOV-DTGVBAR);


 CCONSHR..         DCHR=E=DCHRBAR+(BETCHR/(PCC*RR*ALFAVATE("COALe")))*(EXHR-DTHRBAR);
 CCONSHP..         DCHP=E=DCHPBAR+(BETCHP/(PCC*RR*ALFAVATE("COALe")))*(EXHP-DTHPBAR);


 ECONSGV..         DEGV=E=DEGVBAR+(BETEGV/(PE*(1+SUBE)))*(EXGOV-DTGVBAR);


 ECONSHR..         DEHR=E=DEHRBAR+(BETEHR/(PEC*RR*ALFAVATE("ELECe")))*(EXHR-DTHRBAR);
 ECONSHP..         DEHP=E=DEHPBAR+(BETEHP/(PEC*RR*ALFAVATE("ELECe")))*(EXHP-DTHPBAR);


 FCONSGV..         DFGV=E=DFGVBAR+(BETFGV/(PFX*(1+SUBF)))*(EXGOV-DTGVBAR);


 FCONSHR..         DFHR=E=DFHRBAR+(BETFHR/(PFC*RR*ALFAVATE("FUEL")))*(EXHR-DTHRBAR);
 FCONSHP..         DFHP=E=DFHPBAR+(BETFHP/(PFC*RR*ALFAVATE("FUEL")))*(EXHP-DTHPBAR);


 FeCONSGV..        DFeGV=E=DFeGVBAR+(BETFeGV/(PFe*(1+SUBFe)))*(EXGOV-DTGVBAR);


 FeCONSHR..        DFeHR=E=DFeHRBAR+(BETFeHR/(PFeC*RR*ALFAVATE("FUELe")))*(EXHR-DTHRBAR);
 FeCONSHP..        DFeHP=E=DFeHPBAR+(BETFeHP/(PFeC*RR*ALFAVATE("FUELe")))*(EXHP-DTHPBAR);


 MCONSGV..         DMGV=E=DMGVBAR+(BETMGV/PM)*(EXGOV-DTGVBAR);


 MCONSHR..         DMHR=E=DMHRBAR+(BETMHR/(PM*RR*ALFAVATM))*(EXHR-DTHRBAR);
 MCONSHP..         DMHP=E=DMHPBAR+(BETMHP/(PM*RR*ALFAVATM))*(EXHP-DTHPBAR);


* VAT income for the government

 DEFVATR..           VATR=E=SUM(I$N(I),VATNR(I)*PNC(I))+VATCR*PCC+VATFR*PFC
                        +VATFeR*PFeC+VATER*PEC+VATMR*PM;

 DEFVATP..           VATP=E=SUM(I$N(I),VATNP(I)*PNC(I))+VATCP*PCC+VATFP*PFC
                        +VATFeP*PFeC+VATEP*PEC+VATMP*PM;

 DEFVAT..            VAT=E=VATR+VATP;


 DEFVATNR(I)$N(I)..  VATNR(I)$N(I) =E=(RR*ALFAVAT(I)-1)*DHR(I);

 DEFVATCR..          VATCR =E=(RR*ALFAVATE("COALe")-1)*DCHR;

 DEFVATFR..          VATFR =E=(RR*ALFAVATE("FUEL")-1)*DFHR;

 DEFVATFeR..         VATFeR =E=(RR*ALFAVATE("FUELe")-1)*DFeHR; 

 DEFVATER..          VATER =E=(RR*ALFAVATE("ELECe")-1)*DEHR;

 DEFVATMR..          VATMR =E=(RR*ALFAVATM-1)*DMHR;


 DEFVATNP(I)$N(I)..  VATNP(I)$N(I) =E=(RR*ALFAVAT(I)-1)*DHP(I);

 DEFVATCP..          VATCP =E=(RR*ALFAVATE("COALe")-1)*DCHP;

 DEFVATFP..          VATFP =E=(RR*ALFAVATE("FUEL")-1)*DFHP;

 DEFVATFeP..         VATFeP =E=(RR*ALFAVATE("FUELe")-1)*DFeHP; 

 DEFVATEP..          VATEP =E=(RR*ALFAVATE("ELECe")-1)*DEHP;

 DEFVATMP..          VATMP =E=(RR*ALFAVATM-1)*DMHP;


 DEFVATN(I)$N(I)..  VATN(I)$N(I) =E=(RR*ALFAVAT(I)-1)*(DHR(I)+DHP(I));

 DEFVATC..          VATC =E=(RR*ALFAVATE("COALe")-1)*(DCHR+DCHP);

 DEFVATF..          VATF =E=(RR*ALFAVATE("FUEL")-1)*(DFHR+DFHP);

 DEFVATFe..         VATFe =E=(RR*ALFAVATE("FUELe")-1)*(DFeHR+DFeHP); 

 DEFVATE..          VATE =E=(RR*ALFAVATE("ELECe")-1)*(DEHR+DEHP);

 DEFVATM..          VATM =E=(RR*ALFAVATM-1)*(DMHR+DMHP);


* Price of consumption by households including VAT

 DEFPCONN(I)$N(I).. PCON(I)$N(I)  = PNC(I)*(DHR(I)+DHP(I)+VATN(I))/(DHR(I)+DHP(I));

 DEFPCONC..         PCONC = PCC*(DCHR+DCHP+VATC)/(DCHR+DCHP);

 DEFPCONF..         PCONF = PFC*(DFHR+DFHP+VATF)/(DFHR+DFHP);

 DEFPCONFe..        PCONFe = PFeC*(DFeHR+DFeHP+VATFe)/(DFeHR+DFeHP);

 DEFPCONE..         PCONE = PEC*(DEHR+DEHP+VATE)/(DEHR+DEHP);




* SECTORAL WAGES [thousand EKK/person]


DEFWAGEL(I)..        WAGEL(I)=E=WAGEL0(I)*PLLOWX;

DEFWAGELE(EN)..      WAGELE(EN)=E=WAGELE0(EN)*PLLOWX;

DEFWAGEH(J)..        WAGEH(J)=E=WAGEH0(J)*PLHIGHX;

DEFWAGEHE(EN)..      WAGEHE(EN)=E=WAGEHE0(EN)*PLHIGHX;

DEFWAGELL..          WAGELL=E=WAGELL0*PLLOWX;

DEFWAGELH..          WAGELH=E=WAGELH0*PLHIGHX;

*Market clearing for labour
*Supply=Demand [measured in EKK]

*growth rate of supply (((1+LABLGRO)**TIME)*IO("LABOURL","TOTDEL")) is equivalent to 
*(((1+LABLGRO)**TIME)*IO("LABOURL","TOTDEL")/PLLOW0), beacuse PLLOW0=1, We can drop it

 LLMARKET..         LL*LSWITCH+(1-LSWITCH)*(((1+LABLGRO)**TIME)*IO("LABOURL","TOTDEL"))
                    =E=SUM(J,LLT(J))+SUM(EN,LLTE(EN)) ;

 LHMARKET..         LH*HSWITCH+(1-HSWITCH)*(((1+LABHGRO)**TIME)*IO("LABOURH","TOTDEL"))
                    =E=SUM(J,LHT(J))+SUM(EN,LHTE(EN));


 DEFLABHTAX(I)..    LABHTAX(I)=E=LHT(I)*PLHIGHX*GAMMAHTAX(I);

 DEFLABHSOC(I)..    LABHSOC(I)=E=LHT(I)*PLHIGHX*GAMMAHSOC(I);

 DEFLABLTAX(I)..    LABLTAX(I)=E=LLT(I)*PLLOWX*GAMMALTAX(I);

 DEFLABLSOC(I)..    LABLSOC(I)=E=LLT(I)*PLLOWX*GAMMALSOC(I);

 DEFLABHTAXE(EN)..  LABHTAXE(EN)=E=LHTE(EN)*PLHIGHX*GAMMAEHTAX(EN);

 DEFLABHSOCE(EN)..  LABHSOCE(EN)=E=LHTE(EN)*PLHIGHX*GAMMAEHSOC(EN);

 DEFLABLTAXE(EN)..  LABLTAXE(EN)=E=LLTE(EN)*PLLOWX*GAMMAELTAX(EN);

 DEFLABLSOCE(EN)..  LABLSOCE(EN)=E=LLTE(EN)*PLLOWX*GAMMAELSOC(EN);
 

 DEFLABTAXTOT..     LABTAXTOT=E=SUM(I,LABHTAX(I))+SUM(EN,LABHTAXE(EN))+ 
                                SUM(I,LABLTAX(I))+SUM(EN,LABLTAXE(EN));

 DEFLABSOCTOT..     LABSOCTOT=E=SUM(I,LABHSOC(I))+SUM(EN,LABHSOCE(EN))+ 
                                SUM(I,LABLSOC(I))+SUM(EN,LABLSOCE(EN));



 DEFLL..            LL=E=LL*LSWITCH+(1-LSWITCH)*(((1+LABLGRO)**TIME)
                        *IO("LABOURL","TOTDEL"));


 DEFLH..            LH=E=LH*HSWITCH+(1-HSWITCH)*(((1+LABHGRO)**TIME)
                        *IO("LABOURH","TOTDEL"));


 DEFPLLOWX..        PLLOWX=E=((1+PLLGRO)**TIME)*LSWITCH+(1-LSWITCH)*PLLOWX;

 DEFPLHIGHX..       PLHIGHX=E=((1+PLHGRO)**TIME)*HSWITCH+(1-HSWITCH)*PLHIGHX;

 DEFPLLOW..     PLLOW=E=PLLOWX/(1+SUBLL);

 DEFPLHIGH..    PLHIGH=E=PLHIGHX/(1+SUBLH);


 DEFLABSUB..    LABSUB=E=SUBLL*EMPLTOT*PLLOW+SUBLH*EMPHTOT*PLHIGH;

                       
 DEFUNL..           UNL=E=EMPLBAR-EMPLTOT;

 DEFUNH..           UNH=E=EMPHBAR-EMPHTOT;

* EMPL(I)=LLT(I)*PLLOW/WAGEL(I), where WAGEL(I)=WAGEL0(I)*PLLOW
* [thousand of people]
 DEFEMPL(I)..       EMPL(I)=E=LLT(I)/WAGEL0(I);

 DEFEMPLE(EN)..     EMPLE(EN)=E=LLTE(EN)/WAGELE0(EN);

 DEFEMPH(I)..       EMPH(I)=E=LHT(I)/WAGEH0(I);

 DEFEMPHE(EN)..     EMPHE(EN)=E=LHTE(EN)/WAGEHE0(EN);

 DEFEMPLT..         EMPLTOT=E=SUM(I,EMPL(I))+SUM(EN,EMPLE(EN));

 DEFEMPHT..         EMPHTOT=E=SUM(I,EMPH(I))+SUM(EN,EMPHE(EN));

 DEFEMPT..          EMPTOT=E=EMPLTOT+EMPHTOT;

 DEFEMPLB..         EMPLBAR=E=(EMPLTOT0+UNL0)*((1+LLFGRO)**TIME);

 DEFEMPHB..         EMPHBAR=E=(EMPHTOT0+UNH0)*((1+LHFGRO)**TIME);

 DEFEMPB..          EMPBAR=E=EMPLBAR+EMPHBAR;

 DEFEMPHG(I)..      EMPHGRO(I)=E=EMPH(I)/EMPH0(I);

 DEFEMPHEG(EN)..    EMPHEGRO(EN)=E=EMPHE(EN)/EMPHE0(EN);

 DEFEMPLG(I)..      EMPLGRO(I)=E=EMPL(I)/EMPL0(I);

 DEFEMPLEG(EN)..    EMPLEGRO(EN)=E=EMPLE(EN)/EMPLE0(EN);


* Emissions
* To calculate unit abatement costs an abatement cost curve can be used.
* First the share of emission reduction has to be calculated, after that
* the corresponding abatement cost can be found.


 DEFEMRED..     EMRED=E=100*(CLEAN/REF);


*[bln EKK / mln t]
 EMCOSTCO2..        ABAT=E=0.1;


* Calculation of total emission [mln t] - measured in base year prices

 DEFEMIS..          SUM(J,PEMS(J)*X(J))+ SUM(EN,PEME(EN)*XEE(EN))

                    +SUM((I,J),NEMS(I,J)*AT(I,J)*X(J))

                    +SUM(I,(NEMH(I)*(DGV(I)+DHR(I)+DHP(I))))+SUM((I,EN),NEME(I,EN)*EDS(I,EN))

                    +SUM((EN,J),CEMS(EN,J)*CT(J))

                    +CEMH*(DCHR+DCHP+DCGV)+SUM((EN,EN1),CEME(EN,EN1)*CTE(EN))

                    -TOTEM=E=0;

 EMISCALC(J)..   EMIS(J)=E=PEMS(J)*X(J) + SUM(I,NEMS(I,J)*AT(I,J)*X(J)) + SUM(EN,CEMS(EN,J))*CT(J);

 DEFENEMIS(EN)..  ENEMIS(EN)=E=SUM(I,NEME(I,EN)*EDS(I,EN)) + SUM(EN1,CEME(EN,EN1))*CTE(EN)+PEME(EN)*XEE(EN);

 EMISHOUS..       HOUSEMIS=E=SUM(I,(NEMH(I)*(DGV(I)+DHR(I)+DHP(I))))+CEMH*(DCHR+DCHP+DCGV);

*EMISTAX not for houeholds
*[bln EKK]
 FUNDCALC..       FUND=E=(0+EMISTAX)*(SUM(J,EMIS(J))+SUM(EN,ENEMIS(EN)));


*  The market for emission permits
* PK*ABAT= MC of cleaning
*PEM[thousand EKK/t]

 PERMARKET..    TOTEM-CLEAN-EMLIM=L=0;

 COMPLEM..      PEM*(TOTEM-CLEAN-EMLIM)=E=0;

 PROFCLEAN..    PEM-PK/(1+SUBK)*ABAT=L=0;

 DEFCLEANSEC..  CLEAN=E=SUM(J,CLEANSEC(J))+SUM(EN,CLEANSECEN(EN));

 COMPCLEANJ(J)..  CLEANSEC(J)*(PEM-PK/(1+SUBK)*ABAT)=E=0;

 COMPCLEANE(EN)..  CLEANSECEN(EN)*(PEM-PK/(1+SUBK)*ABAT)=E=0;

 WARUNEK(J)..    CLEANSEC(J)=L=EMIS(J) ;
 
 WARUNEKE(EN)..  CLEANSECEN(EN)=L=ENEMIS(EN);



* Definition of export and import


 DEFEXPN(I)$N(I)..   ZN(I)=E=ZN0(I)*((PN(I)/PWN(I))**(EPSI(I)*EPSIX))*EXP(EXPGRO(I)*TIME);

 DEFEXPNE..   ZNE=E=EZN0("ELECe")*((PE/PWEN("ELECe"))**(EPSIEN("ELECe")*EPSIX))*EXP(EXPGROEN("ELECe")*TIME);
 DEFEXPNC..   ZNC=E=EZN0("COALe")*((PC/PWEN("COALe"))**(EPSIEN("COALe")*EPSIX))*EXP(EXPGROEN("COALe")*TIME);
 DEFEXPNF..   ZNF=E=EZN0("FUEL")*((PFX/PWEN("FUEL"))**(EPSIEN("FUEL")*EPSIX))*EXP(EXPGROEN("FUEL")*TIME);
 DEFEXPNFe..  ZNFe=E=EZN0("FUELe")*((PFe/PWEN("FUELe"))**(EPSIEN("FUELe")*EPSIX))*EXP(EXPGROEN("FUELe")*TIME);


 DEFIMPORT..        M=E=SUM(J,AM(J)*X(J))+IMPENS("COALe")+IMPENS("ELECe")+IMPENS("FUELe")
                        +DMGV+DMHR+DMHP+EXDEMM ;

 DEFIMPTOT..         MTOT=E=M+IMPENS("FUEL");

 IMPORSEC(EN)..     IMPENS(EN)=E=ME(EN)*XEE(EN);

 IMPOREN..          IMPEN=E=SUM(EN,IMPENS(EN));


 DEFIMPORTS..  IMPORTS=E=PM*M+PFX*IMPENS("FUEL");     

 DEFEXPORTS..  EXPORTS=E=SUM(I$N(I),PN(I)*ZN(I))
                     +PC*ZNC+PFX*ZNF
                     +PFe*ZNFe+PE*ZNE;

 ACCOUNT..     EXPORTS-IMPORTS=E= SBAR;



*Market clearing for capital
*We have to add the value of shift scale into equation KPFTAX  

 KMARKET..          K=E=SUM(J,KT(J))+SUM(EN,KTE(EN))+CLEAN*ABAT; 

 DEFKFTAX..         KFTAX=E=F*FTAX+C*CTAX+Fe*FeTAX+E*ETAX+F*NFTAX+C*NCTAX+Fe*NFeTAX+E*NETAX
                        +SUM(I,X(I)*(NTAX(I)+NNTAX(I)));
 DEFKPFTAX..        KPFTAX=E=MCEN("FUEL")*F*FTAX+MCEN("COALe")*C*CTAX+MCEN("FUELe")*Fe*FeTAX
                                +MCEN("ELECe")*E*ETAX+F*NFTAX+C*NCTAX+Fe*NFeTAX+E*NETAX
                                +SUM(I,X(I)*NNTAX(I))+SUM(I,MCT(I)*X(I)*NTAX(I))+EMISTAX*TOTEM;

* SECTORAL RATE [capital units]

DEFRATE(I)..        RATE(I)=E=RATE0(I)*PK;

DEFRATEE(EN)..      RATEE(EN)=E=RATEE0(EN)*PK;

DEFRATTOT..         RATTOT=E=RATTOT0*PK;

* SECTORAL durable resources [bln EKK]
 DEFTRWAL(I)..       TRWAL(I)=E=RATE(I)*(KT(I)+CLEANSEC(I)*ABAT);

 DEFTRWALE(EN)..     TRWALE(EN)=E=RATEE(EN)*(KTE(EN)+CLEANSECEN(EN)*ABAT);

 DEFTRWALTOT..       TRWALTOT=E=RATTOT*(K+EMISTAX*TOTEM);

* Calculation of trade balance and GDP


 DEFGDP..          GDP=E=SUM(I$N(I),PCON(I)*(DHR(I)+DHP(I)))+SUM(I,MCT(I)*DGV(I))
                        +PEM*CLEAN
                        +PCONC*(DCHR+DCHP)+PCONE*(DEHR+DEHP)
                        +PCONF*(DFHR+DFHP)+PCONFe*(DFeHR+DFeHP)
                        +MCEN("COALe")*DCGV+MCEN("ELECe")*DCGV+MCEN("FUEL")*DFGV+MCEN("FUELe")*DFeGV
                        +IO("GROSSPROD","EXTDEM")+IO("GROSSPROD","STOCK")+SBAR;



 DEFGDP0..         GDP0=E=IO("GROSSPROD","CONHR")+IO("GROSSPROD","CONHP")
                         +IO("GROSSPROD","CONGOV")+IO("GROSSPROD","EXTDEM")
                         +IO("GROSSPROD","EXPORT")-IO("IMPORT","TOTDEL")
                         +IO("GROSSPROD","STOCK");


 DEFPKB..        PKB=E=PK/(1+SUBK)*K +PLLOW*LL + PLHIGH*LH 
                       +SUM(J,TTAX(J)*X(J))+SUM(EN,ENTAX(EN)*XEE(EN))+VAT
                       +EMISTAX*TOTEM+KPFTAX-SUBTOT  
                       +PEM*CLEAN
                       +SUM(J,PROF(J)*X(J))+SUM(EN,PROFEN(EN)*XEE(EN));
                
              
 DEFPKB0..       PKB0=E=IO("CAPITAL","TOTDEL")+IO("LABOURL","TOTDEL")+IO("LABOURH","TOTDEL")
                       +IO("TAXES","TOTDEL")+IO("PROFIT","TOTDEL");

 


 DEFPKKB..       PKKB=E=SUM(I,X(I)*PN(I))+C*PC+F*PF+Fe*PFe+E*PE
                      -SUM((I,J),AT(I,J)*X(J)*PNC(J))-SUM((I,EN1),EDS(I,EN1)*PNEN(I,EN1))
                      -SUM(EN,CTE(EN)*PCEN(EN)+FTE(EN)*PFEN(EN)+FeTE(EN)*PFeEN(EN)+ETE(EN)*PEEN(EN))
                      -SUM(J, CT(J)*PCT(J)+FT(J)*PFT(J)+FeT(J)*PFeT(J)+ET(J)*PET(J))
                      -IMPORTS+VAT;
                     
      
 DEFPKKB0..       PKKB0=E=SUM(J,IO("GROSSPROD",J))+SUM(EN1,IO("GROSSPROD",EN1))
                       -SUM((I,J),IO(I,J))-SUM((I,EN1),IO(I,EN1))
                       -SUM((EN,J),IO(EN,J))-SUM((EN,EN1),IO(EN,EN1))
                       -IO("IMPORT","TOTDEL")+IO("TAXES","CONHP")+IO("TAXES","CONHR");

                  
 DEFXPGRO(J)..      XPGRO(J)=E=X(J)*MCT(J)/X0(J);

 DEFXPEEGRO(EN)..   XPEEGRO(EN)=E=XEE(EN)*MCEN(EN)/XE0(EN);

 DEFXGRO(J)..       XGRO(J)=E=X(J)/X0(J);

 DEFXEEGRO(EN)..    XEEGRO(EN)=E=XEE(EN)/XE0(EN);


**Lower bounds for prices

        PN.LO("METALe")=0.10;

        PN.LO("MINERALe")=0.10;

        PN.LO("WOODPAP")=0.10;

        PN.LO("ELENGEE")=0.10;

        PN.LO("LIGHT")=0.10;

        PN.LO("FOOD")=0.10;

        PN.LO("CHEMICe")=0.10;

        PN.LO("CONSTR")=0.10;

        PN.LO("AGRICUL")=0.10;

        PN.LO("TRANSPe")=0.10;

        PN.LO("SERVICES")=0.10;

        PN.LO("PUBLSERV")=0.10;

        PN.LO("PUBLSERVe")=0.10;

        PLL.LO(I)=0.10;

        PYT.LO(I)=0.10;

        PPU.LO(I)=0.10;

        PVT.LO(I)=0.10;

        PRT.LO(I)=0.10;

        PQT.LO(I)=0.10;


        PLLE.LO(EN)=0.10;

        PYTE.LO(EN)=0.10;

        PPUE.LO(EN)=0.10;

        PVTE.LO(EN)=0.10;

        PRTE.LO(EN)=0.10;

        PQTE.LO(EN)=0.10;


        PK.LO=0.1000;

        PLLOW.LO=0.1000;

       PLHIGH.LO=0.1000;

        PLLOWX.LO=0.1000;



       PLHIGHX.LO=0.1000;


        PE.LO=0.1000;
        
        PC.LO=0.1000;

        PF.LO=0.100;

        PFe.LO=0.1000;

** Initial values
* Prices


     PN.L("METALe")=1.000;

     PN.L("MINERALe")=1.000;

     PN.L("WOODPAP")=1.000;

     PN.L("ELENGEE")=1.000;

     PN.L("LIGHT")=1.000;

     PN.L("FOOD")=1.000;

     PN.L("CHEMICe")=1.000;

     PN.L("CONSTR")=1.000;

     PN.L("AGRICUL")=1.000;

     PN.L("TRANSPe")=1.000;

     PN.L("SERVICES")=1.000;

     PN.L("PUBLSERV")=1.000;

     PN.L("PUBLSERVe")=1.000;

     PLL.L(I)=1.000;

     PYT.L(I)=1.000;

     PPU.L(I)=1.000;

     PVT.L(I)=1.000;
     
     PRT.L(I)=1.000;

     PQT.L(I)=1.000;


     PLLE.L(EN)=1.000;

     PYTE.L(EN)=1.000;

     PPUE.L(EN)=1.000;

     PVTE.L(EN)=1.000;
     
     PRTE.L(EN)=1.000;

     PQTE.L(EN)=1.000;


     PK.L=1.000;

     PLLOW.L=1.000;

     PLHIGH.L=1.000;

     PLLOWX.L=1.000;

     PLHIGHX.L=1.000;

     PE.L=1.000;

     PC.L=1.000;

     PF.L=1.000;

     PFe.L=1.000;

     PFC.L=1.000;

     PFT.L(I)=1.000;
     PFEN.L(EN)=1.000;

     PCC.L=1.000;

     PCT.L(I)=1.000;
     PCEN.L(EN)=1.000;

     PFeC.L=1.000;

     PFeT.L(I)=1.000;
     PFeEN.L(EN)=1.000;

     PEC.L=1.000;

     PET.L(I)=1.000;
     PEEN.L(EN)=1.000;

     MCT.L(I)=1.000;

     MCEN(EN).L=1.000;
 
PNC.L(I)=1.000000;
PNT.L(I,J)=1.000000;
PNEN.L(I,EN)=1.000000;


     PCON.L("METALe")=1.18;

     PCON.L("MINERALe")=1.18;

     PCON.L("WOODPAP")=1.18;

     PCON.L("ELENGEE")=1.18;

     PCON.L("LIGHT")=1.18;

     PCON.L("FOOD")=1.18;

     PCON.L("CHEMICe")=1.18;

     PCON.L("CONSTR")=1.18;

     PCON.L("AGRICUL")=1.18;

     PCON.L("TRANSPe")=1.18;

     PCON.L("SERVICES")=1.18;

     PCON.L("PUBLSERV")=1.0;

     PCON.L("PUBLSERVe")=1.18;

     PCONE.L=1.18;

     PCONC.L=1.0;

     PCONF.L=1.18;

     PCONFe.L=1.18;



     WAGEL.L(I)=WAGEL0(I);

     WAGELE.L(EN)=WAGELE0(EN);

     WAGEH.L(I)=WAGEH0(I);

     WAGEHE.L(EN)=WAGEHE0(EN);

     WAGELL.L=WAGELL0;

     WAGELH.L=WAGELH0;


* Quantities

     X.L(I)=X0(I);
XEE.L(EN)=XE0(EN);

     XGRO.L(I)=1.000;

     XPGRO.L(I)=1.000;

     XEEGRO.L(EN)=1.000;

     XPEEGRO.L(EN)=1.000;



     ZN.L(I)$N(I)=ZN0(I);
     ZNC.L=EZN0("COALe");
ZNF.L=EZN0("FUEL");
ZNFe.L=EZN0("FUELe");
ZNE.L=EZN0("ELECe");
EXPORTS.L=IO("EXPORT","GROSSPROD");

LABHTAX.L(I)=GAMMAHTAX(I)*LHIGH0(I);
LABHSOC.L(I)=GAMMAHSOC(I)*LHIGH0(I);
LABLTAX.L(I)=GAMMALTAX(I)*LLOW0(I);
LABLSOC.L(I)=GAMMALSOC(I)*LLOW0(I);
LABHTAXE.L(EN)=GAMMAEHTAX(EN)*IO("LABOURH",EN);
LABHSOCE.L(EN)=GAMMAEHSOC(EN)*IO("LABOURH",EN);
LABLTAXE.L(EN)=GAMMAELTAX(EN)*IO("LABOURL",EN);
LABLSOCE.L(EN)=GAMMAELSOC(EN)*IO("LABOURL",EN);
LABTAXTOT.L=SUM(I,LABHTAX.L(I))+SUM(EN,LABHTAXE.L(EN))+SUM(I,LABLTAX.L(I))+SUM(EN,LABLTAXE.L(EN));
LABSOCTOT.L=SUM(I,LABHSOC.L(I))+SUM(EN,LABHSOCE.L(EN))+SUM(I,LABLSOC.L(I))+SUM(EN,LABLSOCE.L(EN));


VATNR.L(I)=(PCON.L(I)-1)*D0HR(I);
VATCR.L=0.0*DC0HR;
VATFR.L=0.18*DF0HR;
VATFeR.L=0.18*DFe0HR; 
VATER.L=0.18*DE0HR;
VATMR.L=(ALFAVATM-1)*DM0HR;

VATNP.L(I)=(PCON.L(I)-1)*D0HP(I);
VATCP.L=0.0*DC0HP;
VATFP.L=0.18*DF0HP;
VATFeP.L=0.18*DFe0HP; 
VATEP.L=0.18*DE0HP;
VATMP.L=(ALFAVATM-1)*DM0HP;

VATR.L=IO("TAXES","CONHR");
VATP.L=IO("TAXES","CONHP");
VAT.L=VATR.L+VATP.L;


INGV.L=IO("TAXES","TOTDEL")+LABTAXTOT.L;
EXGOV.L=IO("GROSSPROD","CONGOV");

IN.L=IO("CAPITAL","TOTDEL")+IO("LABOURL","TOTDEL")+IO("LABOURH","TOTDEL")+IO("PROFIT","TOTDEL");
*INHR.L=
*INHP.L=EXHP.L+
EXH.L=IO("GROSSPROD","CONHR")+IO("GROSSPROD","CONHP");
EXHR.L=IO("GROSSPROD","CONHR");
EXHP.L=IO("GROSSPROD","CONHP");

DTGVBAR.L=DCGVBAR+DFGVBAR+DFeGVBAR+DEGVBAR+SUM(I,DGVBAR(I))+DMGVBAR;
DTHRBAR.L=SUM(J$N(J),PCON.L(J)*DHRBAR(J))
                            +PCONE.L*DEHRBAR
                            +PCONF.L*DFHRBAR
                            +PCONFe.L*DFeHRBAR
                            +PM*ALFAVATM*DMHRBAR
                            +PCONC.L*DCHRBAR ;
DTHPBAR.L=SUM(J$N(J),PCON.L(J)*DHPBAR(J))
                            +PCONE.L*DEHPBAR
                            +PCONF.L*DFHPBAR
                            +PCONFe.L*DFeHPBAR
                            +PM*ALFAVATM*DMHPBAR
                            +PCONC.L*DCHPBAR ;

     DGV.L(I)=D0GV(I);

     DHR.L(I)=D0HR(I);
     DHP.L(I)=D0HP(I);

     DCGV.L=DC0GV;

     DCHR.L=DC0HR;
     DCHP.L=DC0HP;

     DEGV.L=DE0GV;

     DEHR.L=DE0HR;
     DEHP.L=DE0HP;

     DFGV.L=DF0GV;

     DFeGV.L=DFe0GV;

     DFHR.L=DF0HR;
     DFHP.L=DF0HP;

     DFeHR.L=DFe0HR;
     DFeHP.L=DFe0HP;

     DMGV.L=DM0GV;

     DMHR.L=DM0HR;
     DMHP.L=DM0HP;



     KT.L(I)=K0(I);
     KTE.L(EN)=CAPENS0(EN);

RATE.L(I)=RATE0(I);
RATEE.L(EN)=RATEE0(EN);
RATTOT=RATTOT0;
TRWAL.L(I)=TRWAL0(I);
TRWALE.L(EN)=TRWALE0(EN);
TRWALTOT.L=SUM(I,TRWAL0(I))+SUM(EN,TRWALE0(EN));



     LLT.L(I)=LLOW0(I);

     LHT.L(I)=LHIGH0(I);

     LLTE.L(EN)=IO("LABOURL",EN);

     LHTE.L(EN)=IO("LABOURH",EN);


     LL.L=IO("LABOURL","TOTDEL");

     LH.L=IO("LABOURH","TOTDEL");

  

     ET.L(I)=E0(I);

     FT.L(I)=F0(I);

     FeT.L(I)=Fe0(I);

     CT.L(I)=C0(I);

     ETE.L(EN)=EN0(EN);

     FTE.L(EN)=FN0(EN);

     FeTE.L(EN)=FeN0(EN);

     CTE.L(EN)=CN0(EN);


     

* Derivatives

     DPXDPKT.L(I)=KX0(I);

     DPXDPYT.L(I)=LHIGHX0(I);

     DPXDPLL.L(I)=LLOWX0(I);

     DPXDPET.L(I)=EX0(I);

     DPXDPFT.L(I)=FX0(I);

     DPXDPCT.L(I)=CX0(I);

     DPXDPFeT.L(I)=FeX0(I);


     DPXDPKTE.L(EN)=KXE0(EN);

     DPXDPYTE.L(EN)=LHIGHXE0(EN);

     DPXDPLLE.L(EN)=LLOWXE0(EN);

     DPXDPETE.L(EN)=EXE0(EN);

     DPXDPFTE.L(EN)=FXE0(EN);

     DPXDPCTE.L(EN)=CXE0(EN);

     DPXDPFeTE.L(EN)=FeXE0(EN);



     DPLLPLL.L(J)=DELLL(J);

     DPLLPYT.L(J)=1-DELLL(J); 

     DPYDPLH.L(J)=DELYL(J);


     DPYDPPU.L(J)=1-DELYL(J);

     DPUDPVT.L(J)=1-DELUK(J);

     DPUDPKT.L(J)=DELUK(J);

     DPVDPET.L(J)=DELVE(J);

     DPVDPRT.L(J)=1-DELVE(J);

     DPRDPFeT.L(J)=DELRFe(J);

     DPRDPQT.L(J)=1-DELRFe(J);

     DPQDPFT.L(J)=1-DELQF(J);

     DPQDPCT.L(J)=DELQF(J);


     DPLLPLLE.L(EN)=EDELLL(EN);

     DPLLPYTE.L(EN)=1-EDELLL(EN); 

     DPYDPLHE.L(EN)=EDELYL(EN);

     DPYDPPUE.L(EN)=1-EDELYL(EN);

     DPUDPVTE.L(EN)=1-EDELUK(EN);

     DPUDPKTE.L(EN)=EDELUK(EN);

     DPVDPETE.L(EN)=EDELVE(EN);

     DPVDPRTE.L(EN)=1-EDELVE(EN);

     DPRDPFeTE.L(EN)=EDELRFe(EN);

     DPRDPQTE.L(EN)=1-EDELRFe(EN);

     DPQDPFTE.L(EN)=1-EDELQF(EN);

     DPQDPCTE.L(EN)=EDELQF(EN);




* Other variables

F.L=IO("GROSSPROD","FUEL");
E.L=IO("GROSSPROD","ELECe");
C.L=IO("GROSSPROD","COALe");
Fe.L=IO("GROSSPROD","FUELe");


    IMPEN.L=IMPEN0;

     IMPENS.L(EN)=IMPENS0(EN);

     ED.L(I)=ED0(I);

     EDS.L(I,EN)=EDS0(I,EN);


M.L=SUM(J,AM(J)*X.L(J))+DMGV.L+DMHR.L+DMHP.L+IMPEN0-IMPENS.L("FUEL")
    +IO("IMPORT","EXPORT")+IO("IMPORT","EXTDEM");

MTOT.L=M.L+IMPENS("FUEL");

IMPORTS.L=IO("IMPORT","TOTDEL");     


TOTEM.L= SUM(J,PEMS(J)*X.L(J))+ SUM(EN,PEME(EN)*XEE.L(EN))

                    +SUM((I,J),NEMS(I,J)*AT(I,J)*X.L(J))

                    +SUM(I,(NEMH(I)*(DGV.L(I)+DHR.L(I)+DHP.L(I))))+SUM((I,EN),NEME(I,EN)*EDS.L(I,EN))

                    +SUM((EN,J),CEMS(EN,J)*CT.L(J))

                    +CEMH*(DCHR.L+DCHP.L+DCGV.L)+SUM((EN,EN1),CEME(EN,EN1)*CTE.L(EN));

                  
CLEAN.L=TOTEM.L-EMLIM;

EMIS.L(J)=EMIS(J)=E=PEMS(J)*X.L(J) + SUM(I,NEMS(I,J)*AT(I,J)*X.L(J)) + SUM(EN,CEMS(EN,J)*CT.L(J));

ENEMIS.L(EN)=SUM(I,NEME(I,EN)*EDS.L(I,EN))  + SUM(EN1,CEME(EN,EN1)*CTE.L(EN))+PEME(EN)*XEE.L(EN);

HOUSEMIS.L=SUM(I,(NEMH(I)*(DGV.L(I)+DHR.L(I)+DHP.L(I))))+CEMH*(DCHR.L+DCHP.L+DCGV.L);



*====== Employment =====

EMPL.L(I)=EMPL0(I);

EMPLE.L(EN)=EMPLE0(EN);

EMPH.L(I)=EMPH0(I);

EMPHE.L(EN)=EMPHE0(EN);

EMPLTOT.L=EMPLTOT0;

EMPHTOT.L=EMPHTOT0;

EMPTOT.L=EMPTOT0;

UNL.L=UNL0;

UNH.L=UNH0;

EMPLBAR.L=EMPLTOT.L+UNL.L;

EMPHBAR.L=EMPHTOT.L+UNH.L;

EMPBAR.L=EMPLBAR.L+EMPHBAR.L;

EMPHGRO.L(I)=1.000 ;

EMPHEGRO.L(EN)=1.000 ;

EMPLGRO.L(I)=1.000 ;

EMPLEGRO.L(EN)=1.000 ;


* The model contains all equations mentioned before.

 MODEL ESTONIA /ALL/;

 ESTONIA.OPTFILE=1;


* Solve the model until a feasible solution is found

 SET ITER /1*10/;

 ESTONIA.MODELSTAT =0;

 LOOP(ITER $ (ESTONIA.MODELSTAT NE 2),

      SOLVE ESTONIA USING NLP MINIMIZING DUM1;)



* Output generation

 PARAMETER PIO (S, S1);

 PIO(I,J)=AT(I,J)*X.L(J);

 PIO(I,EN)=EDS.L(I,EN);

 PIO(I,"CONGOV")=DGV.L(I);

 PIO(I,"CONHR")=DHR.L(I);
 PIO(I,"CONHP")=DHP.L(I);

 PIO(I,"STOCK")=IO(I,"STOCK");

 PIO(I,"EXPORT")$N(I)=ZN.L(I);

 PIO(I,"EXTDEM")=IO(I,"EXTDEM");



 PIO("COALe",I)=CT.L(I);

 PIO("ELECe",I)=ET.L(I);

 PIO("FUEL",I)=FT.L(I);

 PIO("FUELe",I)=FeT.L(I);


 PIO("COALe",EN)=CTE.L(EN);

 PIO("ELECe",EN)=ETE.L(EN);

 PIO("FUEL",EN)=FTE.L(EN);

 PIO("FUELe",EN)=FeTE.L(EN);


 PIO("COALe","CONGOV")=DCGV.L;

 PIO("COALe","CONHR")=DCHR.L;
 PIO("COALe","CONHP")=DCHP.L;

 PIO("FUEL","CONGOV")=DFGV.L;

 PIO("FUEL","CONHR")=DFHR.L;
 PIO("FUEL","CONHP")=DFHP.L;

 PIO("FUELe","CONGOV")=DFeGV.L;

 PIO("FUELe","CONHR")=DFeHR.L;
 PIO("FUELe","CONHP")=DFeHP.L;

 PIO("ELECe","CONGOV")=DEGV.L;

 PIO("ELECe","CONHR")=DEHR.L;
 PIO("ELECe","CONHP")=DEHP.L;

 PIO(EN,"STOCK")=IO(EN,"STOCK");

 PIO("COALe","EXPORT")=ZNC.L;

 PIO("ELECe","EXPORT")=ZNE.L;

 PIO("FUELe","EXPORT")=ZNFe.L;

 PIO("FUELe","EXTDEM")=IO("FUELe","EXTDEM");

 PIO("FUEL","EXPORT")=ZNF.L;

 PIO("FUEL","EXTDEM")=IO("FUEL","EXTDEM");

 PIO("ELECe","EXTDEM")=IO("ELECe","EXTDEM");

 PIO("COALe","EXTDEM")=IO("COALe","EXTDEM");


 PIO(I,"TOTDEL")=X.L(I);

 PIO("FUEL","TOTDEL")=XEE.L("FUEL");

 PIO("ELECe","TOTDEL")=XEE.L("ELECe")  ;

 PIO("COALe","TOTDEL")=XEE.L("COALe");

 PIO("FUELe","TOTDEL")=XEE.L("FUELe");


 PIO("IMPORT",J)=AM(J)*X.L(J);

 PIO("IMPORT",EN)=IMPENS.L(EN);

 PIO("IMPORT","CONGOV")=DMGV.L;

 PIO("IMPORT","CONHR")=DMHR.L;
 PIO("IMPORT","CONHP")=DMHP.L;

 PIO("IMPORT","EXPORT")=IO("IMPORT","EXPORT");

 PIO("IMPORT","STOCK")=IO("IMPORT","STOCK");

 PIO("IMPORT","EXTDEM")=IO("IMPORT","EXTDEM");

 PIO("IMPORT","TOTDEL")=MTOT.L;

 PIO("COALe","STOCK")=IO("COALe","STOCK");

 PIO("FUELe","STOCK")=IO("FUELe","STOCK");

 PIO("CAPITAL",J)=KT.L(J);

 PIO("CAPITAL",EN)=KTE.L(EN);

 PIO("CAPITAL","CONHR")=IO("CAPITAL","CONHR");
 PIO("CAPITAL","CONHP")=IO("CAPITAL","CONHP");

 PIO("CAPITAL","CONGOV")=IO("CAPITAL","CONGOV");

 PIO("CAPITAL","TOTDEL")=SUM(J,KT.L(J))+SUM(EN,KTE.L(EN))+CLEAN.L*ABAT; 

 PIO("LABOURL",J)=LLT.L(J);

 PIO("LABOURL",EN)=LLTE.L(EN);

 PIO("LABOURL","CONHR")=IO("LABOURL","CONHR");
 PIO("LABOURL","CONHP")=IO("LABOURL","CONHP");

 PIO("LABOURL","CONGOV")=IO("LABOURL","CONGOV");

 PIO("LABOURL","TOTDEL")=LL.L;

 PIO("LABOURH",J)=LHT.L(J);

 PIO("LABOURH",EN)=LHTE.L(EN);

 PIO("LABOURH","CONHR")=IO("LABOURH","CONHR");
 PIO("LABOURH","CONHP")=IO("LABOURH","CONHP");

 PIO("LABOURH","CONGOV")=IO("LABOURH","CONGOV");

 PIO("LABOURH","TOTDEL")=LH.L;


 PIO("TAXES",EN)=ENTAX(EN)*XEE.L(EN);

 PIO("TAXES","CONHR")=VATR.L;
 PIO("TAXES","CONHP")=VATP.L;

 PIO("TAXES","TOTDEL")=SUM(I,PIO("TAXES",I))+PIO("TAXES","ELECe")+PIO("TAXES","FUEL")+

                     PIO("TAXES","FUELe")+PIO("TAXES","COALe")+

                        PIO("TAXES","CONHR")+PIO("TAXES","CONHP");


                     
 PIO("PROFIT",J)=PROF(J)*X.L(J);

 PIO("PROFIT",EN)=PROFEN(EN)*XEE.L(EN);

 PIO("PROFIT","CONHR")=IO("PROFIT","CONHR");
 PIO("PROFIT","CONHP")=IO("PROFIT","CONHP");

 PIO("PROFIT","TOTDEL")=SUM(I,PIO("PROFIT",I))+SUM(EN,PIO("PROFIT",EN));



 PIO("GROSSPROD",J)=SUM(I,PIO(I,J))+SUM(EN,PIO(EN,J))+PIO("IMPORT",J)+PIO("LABOURL",J)
                        +PIO("LABOURH",J)+PIO("CAPITAL",J)+PIO("PROFIT",J)+PIO("TAXES",J);


 PIO("GROSSPROD",EN)=SUM(I,PIO(I,EN))+SUM(EN1,PIO(EN1,EN))+PIO("IMPORT",EN)+PIO("LABOURL",EN)
                        +PIO("LABOURH",EN)+PIO("CAPITAL",EN)+PIO("PROFIT",EN)+PIO("TAXES",EN);


 PIO("GROSSPROD","CONHR")=SUM(I,DHR.L(I))+DCHR.L+DFHR.L+DFeHR.L+DEHR.L+

                          DMHR.L+VATR.L;
 
 PIO("GROSSPROD","CONHP")=SUM(I,DHP.L(I))+DCHP.L+DFHP.L+DFeHP.L+DEHP.L+

                          DMHP.L+VATP.L;

 PIO("GROSSPROD","STOCK")=SUM(I,PIO(I,"STOCK"))+PIO("ELECe","STOCK")+PIO("FUEL","STOCK")+

                         PIO("FUELe","STOCK")+PIO("COALe","STOCK")+

                         PIO("IMPORT","STOCK")+PIO("TAXES","STOCK") ;

 PIO("GROSSPROD","CONGOV")=SUM(I,DGV.L(I))+DFGV.L+DFeGV.L+DCGV.L+DEGV.L+DMGV.L;

 PIO("GROSSPROD","EXPORT")=SUM(I$N(I),ZN.L(I))+

                          PIO("COALe","EXPORT")+PIO("FUEL","EXPORT")+PIO("FUELe","EXPORT")+

                          PIO("ELECe","EXPORT");

PIO("GROSSPROD","EXTDEM")=IO("GROSSPROD","EXTDEM");

PIO("GROSSPROD","TOTDEL")=SUM(I,PIO("GROSSPROD",I))+SUM(EN,PIO("GROSSPROD",EN))+

                          PIO("GROSSPROD","CONGOV")+

                          PIO("GROSSPROD","CONHR")+

                          PIO("GROSSPROD","CONHP")+

                          PIO("GROSSPROD","STOCK")+

                          PIO("GROSSPROD","EXTDEM")+

                          PIO("GROSSPROD","EXPORT");




 PARAMETER OUT(SV, SH);

 OUT(SV,SH) = PIO(SV,SH);


 PARAMETER OUTDIF(SV, SH);

 OUTDIF(SV,SH)=OUT(SV,SH)-IO(SV,SH);

 PARAMETER SBARX;
           SBARX=PIO("GROSSPROD","EXPORT")-PIO("IMPORT","TOTDEL");


SCALAR NEWYEAR;

NEWYEAR=1997+TIME;

 PARAMETER  INCOMDIF;
            INCOMDIF=PLLOW/PLHIGH;



FILE FF /POLCES2.OUT/;

PUT FF;



*Printing of the I-O-table

PARAMETER TELELR;

PARAMETER ITER1;

PARAMETER ITER2;

PARAMETER MULT;

PARAMETER NUMBER;

ITER1=1;

ITER2=1;

MULT=1;

SET TAB/1*5/;

PUT "I-O-TABEL FOR" NEWYEAR::0 ;

PUT//;

LOOP(TAB,

    NUMBER=5*MULT;

    PUT ' ':20;

    LOOP(SH, PUT $ ((ITER1 <= NUMBER) AND

    (ORD(SH) EQ ITER1)) SH.TL;

       ITER1=ITER1+1 $ ((ITER1 <= NUMBER) AND

    (ORD(SH) EQ ITER1)));

    PUT/;

    LOOP(SV, PUT SV.TL;

      ITER2=1+5*(MULT-1);

      LOOP(SH, PUT $ ((CARD(SH) GE ITER2)

        AND (ITER2 <= NUMBER)

        AND (ORD(SH) EQ ITER2))  (1*PIO(SV,SH))::5;

          ITER2=ITER2+1 $((CARD(SH) GE ITER2)

        AND (ITER2 <= NUMBER)

        AND (ORD(SH) EQ ITER2)));

        PUT/);

        MULT=MULT+1;

PUT//;)

PUT ///;

PUT"***********************************************************************";

PUT ///;


*Printing of the I-O-table differences

PARAMETER ITER1a;

PARAMETER ITER2a;

PARAMETER MULTa;

PARAMETER NUMBERa;

ITER1a=1;

ITER2a=1;

MULTa=1;

SET TAB/1*5/;

PUT "DIFFERECSES IN I-O-TABEL FOR" NEWYEAR::0 ;

PUT//;

LOOP(TAB,

    NUMBER=5*MULTa;

    PUT ' ':20;

    LOOP(SH, PUT $ ((ITER1a <= NUMBER) AND

    (ORD(SH) EQ ITER1a)) SH.TL;

       ITER1a=ITER1a+1 $ ((ITER1a <= NUMBER) AND

    (ORD(SH) EQ ITER1a)));

    PUT/;

    LOOP(SV, PUT SV.TL;

      ITER2a=1+5*(MULTa-1);

      LOOP(SH, PUT $ ((CARD(SH) GE ITER2a)

        AND (ITER2a <= NUMBER)

        AND (ORD(SH) EQ ITER2a))  (1*(PIO(SV,SH)-IO(SV,SH)))::5;

          ITER2a=ITER2a+1 $((CARD(SH) GE ITER2a)

        AND (ITER2a <= NUMBER)

        AND (ORD(SH) EQ ITER2a)));

        PUT/);

        MULTa=MULTa+1;

PUT//;)

PUT ///;

PUT"***********************************************************************";

PUT ///;




*Printing of growth in production

PARAMETER GR(I);

GR(I)=X.L(I)/IO(I,"TOTDEL");

PUT "GROWTH OF SECTORS IN 2010  ";PUT /;

PUT"                          (1992=1.00) ";PUT /;

LOOP (I,PUT I.TL; PUT GR(I)::2; PUT/;)

PUT//;

PUT"***********************************************************************";

PUT //;

*Printing of prices

PUT "PRICES IN " NEWYEAR:4:0;

PUT //;

PUT "PRICES OF GOODS"; PUT /;

PUT"                   (1992=1.00) WORLDPRICE";PUT /;


LOOP (I$N(I), PUT I.TL; PUT MCT.L(I)::2;

PUT PWN(I)::2; PUT/);

PUT /;

PUT" PRICES OF PRODUCTION FACTORS (1992=1.00) "; PUT /;

PUT" PRICE OF CAPITAL      "; PUT PK.L::2; PUT /;

PUT" PRICE OF LABOURL      "; PUT PLLOW.L::2; PUT /;

PUT" PRICE OF LABOURH      "; PUT PLHIGH.L::2; PUT /;

PUT" PRICE OF COALe         "; PUT PC::2; PUT /;

PUT" PRICE OF FUEL         "; PUT PF::2; PUT /;

PUT" PRICE OF FUELe         "; PUT PFe::2; PUT /;

PUT" PRICE OF ELECTRICITY  "; PUT PE.L::2; PUT /;

PUT" PRICE OF IMPORTS      "; PUT PM::2; PUT /;

PUT //;

PUT"***********************************************************************";

PUT //;


*Printing of employment variables

PUT "                  EMPL         EMPL0         EMPH         EMPH0   "; PUT / ;

LOOP (I$N(I), PUT I.TL; PUT EMPL.L(I)::2; PUT EMPL0(I)::2;
                        PUT EMPH.L(I)::2; PUT EMPH0(I)::2; PUT/);

PUT /;

PUT "                EMPLE        EMPLE0         EMPHE         EMPHE0   "; PUT / ;

LOOP (EN, PUT EN.TL; PUT EMPLE.L(EN)::2; PUT EMPLE0(EN)::2;
                     PUT EMPHE.L(EN)::2; PUT EMPHE0(EN)::2; PUT/);

PUT /;


PUT  "EMPLBAR  "  EMPLBAR.L ; PUT /;

PUT  "EMPHBAR  "  EMPHBAR.L ; PUT /;

PUT  "UNL  " UNL.L ; PUT /;

PUT  "UNH  " UNH.L ; PUT /;

PUT  "VATR  " VATR.L ; PUT /;

PUT  "VATP  " VATP.L ; PUT /;

PUT  "VATC  " VATC.L ; PUT /;

PUT  "VATE  " VATE.L ; PUT /;

PUT  "VATF  " VATF.L ; PUT /;

PUT  "VATFe  " VATFe.L ; PUT /;

PUT  "VATM  " VATM.L ; PUT //;

PUT  "              VATN   "; PUT / ;

LOOP (I, PUT I.TL; PUT VATN.L(I)::2; PUT/);

PUT /;

PUT  "LABSUB  " LABSUB.L ; PUT /;

PUT  "PLHIGHX  " PLHIGHX.L ; PUT /;

PUT  "PLLOWX  " PLLOWX.L ; PUT /;


 
PUT "                 LABHTAX       LABHSOC      LABLTAX        LABLSOC   "; PUT / ;

LOOP (I, PUT I.TL; PUT LABHTAX.L(I)::2; PUT LABHSOC(I)::2;
                     PUT LABLTAX.L(I)::2; PUT LABLSOC(I)::2; PUT/);

PUT /;

PUT "                 LABHTAXE      LABHSOCE     LABLTAXE       LABLSOCE  "; PUT / ;

LOOP (EN, PUT EN.TL; PUT LABHTAXE.L(EN)::2; PUT LABHSOCE(EN)::2;
                     PUT LABLTAXE.L(EN)::2; PUT LABLSOCE(EN)::2; PUT/);

PUT /;

PUT  "LABTAXTOT   " LABTAXTOT.L  ; PUT /;

PUT  "LABSOCTOT   " LABSOCTOT.L  ; PUT / ;



*Printing of other variables

*SCALAR GDP90;

*GDP90 = GDP0;

*SCALAR GDP1;

*YY1=YY.L*1;

*SCALAR GDPGRO;

*GNPGRO=YY1/GNP90;

*PUT "GDP IS " YY1::0; PUT/;

*PUT "                   (1990=1.00)  "; PUT/;

*PUT "GROWTH IN GDP IS "; PUT GNPGRO::2; PUT / ;


*PUT "GDP_P92  " GDP_P92*1::0; PUT / ;


*PUT //;
                


PUT "****************************************************************";

PUT //;

*SCALAR FUELGRO;

*FUELGRO =(((F.L/IO("FUELe","TOTDEL"))**(1/TIME))-1);



*SCALAR COALGRO;

*COALGRO =(((C.L/IO("COALe","TOTDEL"))**(1/TIME))-1);



*Printing assumptions

PUT "ASSUMPTIONS MADE IN THE MODEL  "; PUT //;

PUT "GROWTH OF EXPORT DEMAND IS         " EXPGRO("METALe")::6; PUT /;

PUT "LABOURL PRODUCTIVITY GROWTH IS      " LAMBDALL("METALe")::6; PUT /;

PUT "LABOURH PRODUCTIVITY GROWTH IS      " LAMBDALH("METALe")::6; PUT /;

PUT "CAPITAL PRODUCTIVITY GROWTH IS     " LAMBDAK("METALe")::6; PUT /;

PUT "ENERGY EFFICIENCY GROWTH (COALe) IS " ENEFF("METALe","COALe")::6; PUT /;

PUT "ENERGY EFFICIENCY GROWTH (FUELe) IS " ENEFF("METALe","FUELe")::6; PUT /;

PUT "ENERGY EFFICIENCY GROWTH (FUEL) IS " ENEFF("METALe","FUEL")::6; PUT /;

PUT "ENERGY EFFICIENCY GROWTH (ELECe) IS " ENEFF("METALe","ELECe")::6; PUT /;

PUT "CAPITAL ENDOWMENT IS            " (1*K)::0;

PUT "    (GROWTH " CAPGRO " % ANNUALLY )"; PUT /;

PUT "LABOUR ENDOWMENT IS             " (1*LL)::0;

PUT "    (GROWTL " LABLGRO " % ANNUALLY )"; PUT /;

PUT "    (GROWTH " LABHGRO " % ANNUALLY )"; PUT /;

PUT /;

PUT "SIGMA: BASEIND BULKIND AGRI FOOD EXPTRAD CONSTR TRADE SERVICES TRANSPORT PUBLIC"; PUT /;


PUT "OF     "

LOOP(I, PUT SIGMA("QF",I))

PUT "RFe     "

LOOP(I, PUT SIGMA("RFe",I))

PUT "VE     "

LOOP(I, PUT SIGMA("VE",I))

PUT "UK     "

LOOP(I, PUT SIGMA("UK",I))

PUT "YL     "

LOOP(I, PUT SIGMA("YL",I))

PUT ///;

PUT "EXPORT DEMAND ELASTICITY " PUT /;

LOOP (I$N(I), PUT I.TL; PUT (EPSI(I)*EPSIX); PUT /)


PUT //;

PUT ///;