**                                               **
*        Load set contents from the gdx file      *
**                                               **

** The source of input GDX file:
* GTAP data is so detail that we generally need to aggregate to less sectors, regions, and factos. This can be done by
* a modified version of GTAP6inGAMS package provided by Thomas F. Rutherford at http://www.mpsge.org/gtap6/.
** the end of the source of input GDX file.

*$include GTAP_input_header
$GDXIN %ds%.gdx

**                                               **
*        Load set contents from the gdx file      *
**                                               **
$LOAD
$LOAD r=R
$LOAD f=F
$LOAD y=I
$LOAD ty         =        ty
$LOAD ti         =        ti
$LOAD tf         =        tf
$LOAD tx         =        tx
$LOAD tm         =        tm
$LOAD tg         =        tg
$LOAD tp         =        tp
$LOAD vafm       =        vafm
$LOAD vfm        =        vfm
$LOAD vxmd       =        vxmd
* If an error occurs in the next line, try to replace the right "vtwr" with "aggvtwr"
$LOAD vtwr       =        aggvtwr
$LOAD vst        =        vst
$LOAD vdgm       =        vdgm
$LOAD vigm       =        vigm
$LOAD vdpm       =        vdpm
$LOAD vipm       =        vipm
$LOAD vdep       =        vdep
$LOAD vkb       =        vkb
$load  edf edg edp eif eig eip
$load mdf mdg mdp mif mig mip

$GDXIN

**                                               **
*        Define subsets of factors and sectors    *
**                                               **
* General valid settings
SETS
                 flab(f) labour /LAB/
                 fcap(f) capital /CAP/
                 fnat(f) nature /RES/

         ycpg(y) capital good /CPG/
         i(y)

*        Sector groupings: energy
         ie(y)     all energy   /REF,GAS,COL,CRU,ELC/
         ief(y)    Nonrenewable fossil energy   /REF,GAS,COL,CRU/
         ipep(y) primary energy production /CRU,GAS,COL/
         iref(y)  ref oil        /REF/
         ielc(y)  elc            /ELC/
         icol(y)  coal           /COL/
         igas(y)  gas            /GAS/
         icru(y)  crude          /CRU/
* non-energy sectors
         ine(y) all non-energy sectors

         iagr(y)   the sector of agriculture   /AGR/
         iser(y)   services sector  /SER/
         itran(y)  transport sectors  /tran/

* Consumers
         jcon(j) consumption sectors /RHH,UHH,PRI,PUB/
         jpub(j) public consumers  /PUB/
         jpri(j) /RHH,UHH,PRI/
         jrhh(j) rural households /RHH/
         jinv(j) /INV/
         jr(j,r)  map of final deliveries to regions

* Alway defin a region as rest of the world "ROW"
         rusa(r) /USA/
         reu(r) /EU/
         rjpn(r) /JPN/
         rrus(r) /RUS/
         rchn(r) /CHN/
         rind(r) /IND/
         rbra(r) /BRA/
         rrow(r) /ROW/
         rdvp(r) /CHN,IND,BRA,ROW/
         rrhh(r) regions with RHH and UHH divided
;

i(y)$(not ycpg(y)) = yes;
ine(i)$(not ie(i)) = yes;
jr(j,r)$(ord(j)>2) = yes; rrhh(r) = NO ;
display jr;

Alias (i,ii),(i,iii);

* the end of general valid settings

* Below project-specific settings can be added.
SETS
         a     all technologies or activities to produce goods   /ALL, COLE,GASE,REFE,HYDR,NUKE,BIOM,SOLA,WIND,RENE/
         aef(a)  technologies in electricity sector     /COLE,GASE,REFE/
         aren(a)  technologies in electricity sector     /HYDR,BIOM,SOLA,WIND,RENE/
         anuk(a)  technologies in electricity sector     /NUKE/
         aint(a) intermittent technologies in electricity sector /BIOM,SOLA,WIND,RENE/

         fuel(y)   fuel inputs sectors for electricity /COL,GAS,REF,NMM,AGR/
         ay(a,y)  mapping technology a to good y
         afuel(a,y)  mapping technology to fuels used by the technology  /COLE.COL, GASE.GAS, REFE.REF,NUKE.NMM, BIOM.AGR, RENE.SER/
         ;
alias(a,aa) ;
*ay(a,y) = yes$(ord(a)=1) ;
ay(a,i) = yes$( (ord(a)=1 and not ielc(i)) or (ord(a)>1 and ielc(i)) ) ;