$INCLUDE GTAP_input_header

SET
SIM  simulations
;

$SET ds foodpris51
*$include GTAP_input_header
*$GDXIN circe_13x16_6_forDirk.gdx
*$GDXIN india3x3.gdx
$GDXIN %ds%_resPD.gdx
*$INCLUDE circe_12x16_6.map


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

$LOAD
$LOAD r
$LOAD f
$LOAD y
*$LOAD ghga
$LOAD SIM


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

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

*         reu(r) /EU /

*        Sector groupings: energy
         icrop(y)   the sector of agriculture   /PDR,WHT,GRO/
         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

* Consumers
         jcon(j) consumption sectors /PRI,PUB/
         jpub(j) public consumers  /PUB/
         jpri(j) /PRI/
         jinv(j) /INV/
         jr(j,r)

*         Kghga(ghga)  /C2F6,C3F8,C4F10,C5F12,C6F14,C7F16,CF4,CH4,CO2NOC
*                         HFC125,HFC134a,HFC143a,HFC152a,HFC227ea,HFC236fa,HFC23,HFC245fa,HFC32,HFC365mfc,HFC4310mee,N2O
*                         SF6,cC4F8/

;

$LOAD jr

i(y) = yes$(not ycpg(y));

ine(i) = yes$(not ie(i));

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

SETS
*         a     all technologies or activities to produce goods   /ALL, NUKE,RENE/
         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/

*         g technologies to generate elesticity  /NUKE,    RENE,    TGAS,    TCOL,    TREF,    HYDR,    BIOM,    SOLA,    WIND/
*         rweo  regions in LDE project model  / US,EU28,JPN,RUS,CHINA,INDIA,BRAZIL,ROW/
*         ga(g,a)  mapping technology g to a  /NUKE.NUKE,RENE.RENE,TGAS.RENE,TCOL.RENE,TREF.RENE,HYDR.RENE,BIOM.RENE,SOLA.RENE,WIND.RENE/
*         ga(g,a)  mapping technology g to a  /NUKE.NUKE,RENE.RENE,TGAS.GASE,TCOL.COLE,TREF.REFE,HYDR.RENE,BIOM.RENE,SOLA.RENE,WIND.RENE/
         fuel(y)   fuel inputs sectors for electricity /COL,GAS,REF,PRO,AGR/
         ay(a,y)  mapping technology a to good y
         aynoC(a,y)  mapping non-coal technology a to good y
         aynoF(a,y)  mapping non-fossil technology a to good y
*         afuel(a,y)  mapping technology to fuels used by the technology  /RENE.COL, RENE.GAS, RENE.REF,NUKE.NMM, RENE.AGR/
         afuel(a,y)  mapping technology to fuels used by the technology  /COLE.COL, GASE.GAS, REFE.REF,NUKE.PRO, 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)) ) ;
aynoC(a,i) = yes$(ay(a,i) and (not icol(i)) and (not sameas(a,'COLE')) ) ;
aynoF(a,i) = yes$(ay(a,i) and (not ief(i)) and (not aef(a)) ) ;

$include time_declare.gms

$ontext

$LOAD GDPr_T GDP_I
$LOAD PINDEX_T
$LOAD PINDEXW_T
$LOAD PD_T
$LOAD XDr_T Xr_T
$LOAD P_T
$LOAD IO_ABSr_T
$LOAD EXRr_T
$LOAD TRANRr_T
$LOAD CTOP_T
$LOAD FD_ABSr_T INC_T

parameter
 GDPr_m(r,SIM,t)
 GDP_I_m(f,y,a,r,SIM,t)
 PINDEX_m(r,SIM,t)
 PINDEXW_m(SIM,t)
 PD_m(y,r,SIM,t)
 XDr_m(y,r,SIM,t)
 Xr_m(y,r,SIM,t)
 P_m(y,r,SIM,t)
 IO_ABSr_m(yy,a,y,r,SIM,t)
 EXRr_m(y,r,rr,SIM,t)
 TRANRr_m(y,r,SIM,t)
 CTOP_m(y,r,SIM,t)
 FD_ABSr_m(y,j,r,SIM,t)
 INC_m(j,r,SIM,t)
  ;
$GDXIN %ds%_resPD.gdx
$LOAD GDPr_m    =     GDPr_T
$LOAD  GDP_I_m  =      GDP_I
$LOAD PINDEX_m   =    PINDEX_T
$LOAD PINDEXW_m  =    PINDEXW_T
$LOAD PD_m      =     PD_T
$LOAD XDr_m     =     XDr_T
$LOAD Xr_m      =     Xr_T
$LOAD P_m       =     P_T
$LOAD IO_ABSr_m =     IO_ABSr_T
$LOAD EXRr_m    =     EXRr_T
$LOAD TRANRr_m  =     TRANRr_T
$LOAD CTOP_m    =     CTOP_T
$LOAD FD_ABSr_m =     FD_ABSr_T
$LOAD INC_m =     INC_T

 GDPr_T(r,SIM,t)$GDPr_m(r,SIM,t)=  GDPr_m(r,SIM,t)  ;
 GDP_I(f,y,a,r,SIM,t)$sameas(y,'wht') =  GDP_I_m(f,y,a,r,SIM,t) ;
 PINDEX_T(r,SIM,t)$PINDEX_m(r,SIM,t)=    PINDEX_m(r,SIM,t)  ;
 PINDEXW_T(SIM,t)$PINDEXW_m(SIM,t)=    PINDEXW_m(SIM,t)  ;
 PD_T(y,r,SIM,t)$PD_m(y,r,SIM,t)=   PD_m(y,r,SIM,t)   ;
 XDr_T(y,r,SIM,t)$XDr_m(y,r,SIM,t)=  XDr_m(y,r,SIM,t)  ;
 Xr_T(y,r,SIM,t)$Xr_m(y,r,SIM,t) =   Xr_m(y,r,SIM,t)  ;
 P_T(y,r,SIM,t)$P_m(y,r,SIM,t)=  P_m(y,r,SIM,t) ;
 IO_ABSr_T(yy,a,y,r,SIM,t)$IO_ABSr_m(yy,a,y,r,SIM,t)= IO_ABSr_m(yy,a,y,r,SIM,t) ;
 EXRr_T(y,r,rr,SIM,t)$EXRr_m(y,r,rr,SIM,t)=   EXRr_m(y,r,rr,SIM,t) ;
 TRANRr_T(y,r,SIM,t)$TRANRr_m(y,r,SIM,t) =  TRANRr_m(y,r,SIM,t) ;
 CTOP_T(y,r,SIM,t)$CTOP_m(y,r,SIM,t)=           CTOP_m(y,r,SIM,t) ;
 FD_ABSr_T(y,j,r,SIM,t)$FD_ABSr_m(y,j,r,SIM,t)= FD_ABSr_m(y,j,r,SIM,t) ;
 INC_T(j,r,SIM,t)$INC_m(j,r,SIM,t) =   INC_m(j,r,SIM,t)  ;
$GDXIN

$GDXIN ..\FoodPris20191010A\%ds%_resPD.gdx
$LOADm GDPr_m    =     GDPr_T
$LOADm  GDP_I_m  =      GDP_I
$LOADm PINDEX_m   =    PINDEX_T
$LOADm PINDEXW_m  =    PINDEXW_T
$LOADm PD_m      =     PD_T
$LOADm XDr_m     =     XDr_T
$LOADm Xr_m      =     Xr_T
$LOADm P_m       =     P_T
$LOADm IO_ABSr_m =     IO_ABSr_T
$LOADm EXRr_m    =     EXRr_T
$LOADm TRANRr_m  =     TRANRr_T
$LOADm CTOP_m    =     CTOP_T
$LOADm FD_ABSr_m =     FD_ABSr_T
$LOADm INC_m =     INC_T

 GDPr_T(r,SIM,t)$GDPr_m(r,SIM,t)=  GDPr_m(r,SIM,t)  ;
 GDP_I(f,y,a,r,SIM,t)$sameas(y,'wht') =  GDP_I_m(f,y,a,r,SIM,t) ;
 PINDEX_T(r,SIM,t)$PINDEX_m(r,SIM,t)=    PINDEX_m(r,SIM,t)  ;
 PINDEXW_T(SIM,t)$PINDEXW_m(SIM,t)=    PINDEXW_m(SIM,t)  ;
 PD_T(y,r,SIM,t)$PD_m(y,r,SIM,t)=   PD_m(y,r,SIM,t)   ;
 XDr_T(y,r,SIM,t)$XDr_m(y,r,SIM,t)=  XDr_m(y,r,SIM,t)  ;
 Xr_T(y,r,SIM,t)$Xr_m(y,r,SIM,t) =   Xr_m(y,r,SIM,t)  ;
 P_T(y,r,SIM,t)$P_m(y,r,SIM,t)=  P_m(y,r,SIM,t) ;
 IO_ABSr_T(yy,a,y,r,SIM,t)$IO_ABSr_m(yy,a,y,r,SIM,t)= IO_ABSr_m(yy,a,y,r,SIM,t) ;
 EXRr_T(y,r,rr,SIM,t)$EXRr_m(y,r,rr,SIM,t)=   EXRr_m(y,r,rr,SIM,t) ;
 TRANRr_T(y,r,SIM,t)$TRANRr_m(y,r,SIM,t) =  TRANRr_m(y,r,SIM,t) ;
 CTOP_T(y,r,SIM,t)$CTOP_m(y,r,SIM,t)=           CTOP_m(y,r,SIM,t) ;
 FD_ABSr_T(y,j,r,SIM,t)$FD_ABSr_m(y,j,r,SIM,t)= FD_ABSr_m(y,j,r,SIM,t) ;
 INC_T(j,r,SIM,t)$INC_m(j,r,SIM,t) =   INC_m(j,r,SIM,t)  ;
$GDXIN

$GDXIN ..\FoodPris20191010B\%ds%_resPD.gdx
$LOADm GDPr_m    =     GDPr_T
$LOADm  GDP_I_m  =      GDP_I
$LOADm PINDEX_m   =    PINDEX_T
$LOADm PINDEXW_m  =    PINDEXW_T
$LOADm PD_m      =     PD_T
$LOADm XDr_m     =     XDr_T
$LOADm Xr_m      =     Xr_T
$LOADm P_m       =     P_T
$LOADm IO_ABSr_m =     IO_ABSr_T
$LOADm EXRr_m    =     EXRr_T
$LOADm TRANRr_m  =     TRANRr_T
$LOADm CTOP_m    =     CTOP_T
$LOADm FD_ABSr_m =     FD_ABSr_T
$LOADm INC_m =     INC_T

 GDPr_T(r,SIM,t)$GDPr_m(r,SIM,t)=  GDPr_m(r,SIM,t)  ;
 GDP_I(f,y,a,r,SIM,t)$sameas(y,'wht') =  GDP_I_m(f,y,a,r,SIM,t) ;
 PINDEX_T(r,SIM,t)$PINDEX_m(r,SIM,t)=    PINDEX_m(r,SIM,t)  ;
 PINDEXW_T(SIM,t)$PINDEXW_m(SIM,t)=    PINDEXW_m(SIM,t)  ;
 PD_T(y,r,SIM,t)$PD_m(y,r,SIM,t)=   PD_m(y,r,SIM,t)   ;
 XDr_T(y,r,SIM,t)$XDr_m(y,r,SIM,t)=  XDr_m(y,r,SIM,t)  ;
 Xr_T(y,r,SIM,t)$Xr_m(y,r,SIM,t) =   Xr_m(y,r,SIM,t)  ;
 P_T(y,r,SIM,t)$P_m(y,r,SIM,t)=  P_m(y,r,SIM,t) ;
 IO_ABSr_T(yy,a,y,r,SIM,t)$IO_ABSr_m(yy,a,y,r,SIM,t)= IO_ABSr_m(yy,a,y,r,SIM,t) ;
 EXRr_T(y,r,rr,SIM,t)$EXRr_m(y,r,rr,SIM,t)=   EXRr_m(y,r,rr,SIM,t) ;
 TRANRr_T(y,r,SIM,t)$TRANRr_m(y,r,SIM,t) =  TRANRr_m(y,r,SIM,t) ;
 CTOP_T(y,r,SIM,t)$CTOP_m(y,r,SIM,t)=           CTOP_m(y,r,SIM,t) ;
 FD_ABSr_T(y,j,r,SIM,t)$FD_ABSr_m(y,j,r,SIM,t)= FD_ABSr_m(y,j,r,SIM,t) ;
 INC_T(j,r,SIM,t)$INC_m(j,r,SIM,t) =   INC_m(j,r,SIM,t)  ;
$GDXIN

$GDXIN ..\FoodPris20191010C\%ds%_resPD.gdx
$LOADm GDPr_m    =     GDPr_T
$LOADm  GDP_I_m  =      GDP_I
$LOADm PINDEX_m   =    PINDEX_T
$LOADm PINDEXW_m  =    PINDEXW_T
$LOADm PD_m      =     PD_T
$LOADm XDr_m     =     XDr_T
$LOADm Xr_m      =     Xr_T
$LOADm P_m       =     P_T
$LOADm IO_ABSr_m =     IO_ABSr_T
$LOADm EXRr_m    =     EXRr_T
$LOADm TRANRr_m  =     TRANRr_T
$LOADm CTOP_m    =     CTOP_T
$LOADm FD_ABSr_m =     FD_ABSr_T
$LOADm INC_m =     INC_T

 GDPr_T(r,SIM,t)$GDPr_m(r,SIM,t)=  GDPr_m(r,SIM,t)  ;
 GDP_I(f,y,a,r,SIM,t)$sameas(y,'wht') =  GDP_I_m(f,y,a,r,SIM,t) ;
 PINDEX_T(r,SIM,t)$PINDEX_m(r,SIM,t)=    PINDEX_m(r,SIM,t)  ;
 PINDEXW_T(SIM,t)$PINDEXW_m(SIM,t)=    PINDEXW_m(SIM,t)  ;
 PD_T(y,r,SIM,t)$PD_m(y,r,SIM,t)=   PD_m(y,r,SIM,t)   ;
 XDr_T(y,r,SIM,t)$XDr_m(y,r,SIM,t)=  XDr_m(y,r,SIM,t)  ;
 Xr_T(y,r,SIM,t)$Xr_m(y,r,SIM,t) =   Xr_m(y,r,SIM,t)  ;
 P_T(y,r,SIM,t)$P_m(y,r,SIM,t)=  P_m(y,r,SIM,t) ;
 IO_ABSr_T(yy,a,y,r,SIM,t)$IO_ABSr_m(yy,a,y,r,SIM,t)= IO_ABSr_m(yy,a,y,r,SIM,t) ;
 EXRr_T(y,r,rr,SIM,t)$EXRr_m(y,r,rr,SIM,t)=   EXRr_m(y,r,rr,SIM,t) ;
 TRANRr_T(y,r,SIM,t)$TRANRr_m(y,r,SIM,t) =  TRANRr_m(y,r,SIM,t) ;
 CTOP_T(y,r,SIM,t)$CTOP_m(y,r,SIM,t)=           CTOP_m(y,r,SIM,t) ;
 FD_ABSr_T(y,j,r,SIM,t)$FD_ABSr_m(y,j,r,SIM,t)= FD_ABSr_m(y,j,r,SIM,t) ;
 INC_T(j,r,SIM,t)$INC_m(j,r,SIM,t) =   INC_m(j,r,SIM,t)  ;
$GDXIN

$GDXIN ..\FoodPris20191010D\%ds%_resPD.gdx
$LOADm GDPr_m    =     GDPr_T
$LOADm  GDP_I_m  =      GDP_I
$LOADm PINDEX_m   =    PINDEX_T
$LOADm PINDEXW_m  =    PINDEXW_T
$LOADm PD_m      =     PD_T
$LOADm XDr_m     =     XDr_T
$LOADm Xr_m      =     Xr_T
$LOADm P_m       =     P_T
$LOADm IO_ABSr_m =     IO_ABSr_T
$LOADm EXRr_m    =     EXRr_T
$LOADm TRANRr_m  =     TRANRr_T
$LOADm CTOP_m    =     CTOP_T
$LOADm FD_ABSr_m =     FD_ABSr_T
$LOADm INC_m =     INC_T

 GDPr_T(r,SIM,t)$GDPr_m(r,SIM,t)=  GDPr_m(r,SIM,t)  ;
 GDP_I(f,y,a,r,SIM,t)$sameas(y,'wht') =  GDP_I_m(f,y,a,r,SIM,t) ;
 PINDEX_T(r,SIM,t)$PINDEX_m(r,SIM,t)=    PINDEX_m(r,SIM,t)  ;
 PINDEXW_T(SIM,t)$PINDEXW_m(SIM,t)=    PINDEXW_m(SIM,t)  ;
 PD_T(y,r,SIM,t)$PD_m(y,r,SIM,t)=   PD_m(y,r,SIM,t)   ;
 XDr_T(y,r,SIM,t)$XDr_m(y,r,SIM,t)=  XDr_m(y,r,SIM,t)  ;
 Xr_T(y,r,SIM,t)$Xr_m(y,r,SIM,t) =   Xr_m(y,r,SIM,t)  ;
 P_T(y,r,SIM,t)$P_m(y,r,SIM,t)=  P_m(y,r,SIM,t) ;
 IO_ABSr_T(yy,a,y,r,SIM,t)$IO_ABSr_m(yy,a,y,r,SIM,t)= IO_ABSr_m(yy,a,y,r,SIM,t) ;
 EXRr_T(y,r,rr,SIM,t)$EXRr_m(y,r,rr,SIM,t)=   EXRr_m(y,r,rr,SIM,t) ;
 TRANRr_T(y,r,SIM,t)$TRANRr_m(y,r,SIM,t) =  TRANRr_m(y,r,SIM,t) ;
 CTOP_T(y,r,SIM,t)$CTOP_m(y,r,SIM,t)=           CTOP_m(y,r,SIM,t) ;
 FD_ABSr_T(y,j,r,SIM,t)$FD_ABSr_m(y,j,r,SIM,t)= FD_ABSr_m(y,j,r,SIM,t) ;
 INC_T(j,r,SIM,t)$INC_m(j,r,SIM,t) =   INC_m(j,r,SIM,t)  ;
$GDXIN

execute_unload '%ds%_resallPD.gdx ', y, a, r, SIM,t,GDPr_T, GDP_I,PINDEX_T,PINDEXW_T,PD_T,XDr_T,Xr_T,P_T,IO_ABSr_T,EXRr_T,TRANRr_T,CTOP_T,FD_ABSr_T,INC_T;
$exit
$offtext

$GDXIN %ds%_resallPD.gdx
*$GDXIN %ds%_resPD.gdx
$LOAD GDPr_T GDP_I
$LOAD PINDEX_T
$LOAD PINDEXW_T
$LOAD PD_T
$LOAD XDr_T Xr_T
$LOAD P_T
$LOAD IO_ABSr_T FD_ABSr_T
$LOAD EXRr_T
$LOAD TRANRr_T
$LOAD INC_T
$GDXIN


PARAMETER
         VA_T(y,*,SIM,t)
         GDP_D_T(r,SIM,t)  GDP deflated by local price index
         EXRrV_T(y,r,rr,SIM,t)  total export from r to rr
         EXPORT_T(y,r,SIM,t)  net export to other regions
         ShWhtFD(r,SIM,t) wheat consumption percentage share in total HH consumption
         ShFoodFD(r,SIM,t) wheat consumption percentage share in total HH consumption
         FDwWht(r,SIM,t)  consumer goods that needs input of the most wheat in their production
         ;

         GDP_D_T(r,SIM,t)$GDPr_T(r,SIM,t) =  GDPr_T(r,SIM,t) / PINDEX_T(r,SIM,t) * PINDEXW_T(SIM,t) ;
         VA_T(i,r,SIM,t)$PINDEX_T(r,SIM,t) =  ( PD_T(i,r,SIM,t)*XDr_T(i,r,SIM,t) - sum((ii,a)$ay(a,i),P_T(ii,r,SIM,t)*IO_ABSr_T(ii,a,i,r,SIM,t)) ) ;
         VA_T(i,'World',SIM,t) = sum(r, VA_T(i,r,SIM,t) ) ;
         EXRrV_T(i,r,rr,SIM,t)$sameas(i,'WHT') = PD_T(i,r,SIM,t)*EXRr_T(i,r,rr,SIM,t)  ;
         EXPORT_T(i,r,SIM,t) = sum(rr,PD_T(i,r,SIM,t)*EXRr_T(i,r,rr,SIM,t)) +  TRANRr_T(i,r,SIM,t)*PD_T(i,r,SIM,t)
                               - sum(rr,PD_T(i,rr,SIM,t)*EXRr_T(i,rr,r,SIM,t))  ;
         ShWhtFD(r,SIM,t)$INC_T('PRI',r,SIM,t) = P_T('WHT',r,SIM,t)*FD_ABSr_T('WHT','PRI',r,SIM,t)/INC_T('PRI',r,SIM,t)* 100 ;
*         ShFoodFD(r,SIM,t)$INC_T('PRI',r,SIM,t) = P_T('FOD',r,SIM,t)*FD_ABSr_T('FOD','PRI',r,SIM,t)/INC_T('PRI',r,SIM,t)* 100 ;
         FDwWht(r,SIM,t) = P_T('WHT',r,SIM,t)*FD_ABSr_T('WHT','PRI',r,SIM,t)+P_T('AGR',r,SIM,t)*FD_ABSr_T('AGR','PRI',r,SIM,t)+P_T('FOD',r,SIM,t)*FD_ABSr_T('FOD','PRI',r,SIM,t) ;
parameters
         CTOP(y,r,t)    detailed crop production output data from FAO (tonnes)
         CTOPX(y,r,t)    detailed crop production output data from FAO (tonnes)
         CTOPC_T(y,r,SIM,t)    changes in detailed crop production output data from FAO (tonnes)
         CTOPP_T(y,r,SIM,t)    per capita detailed crop production output data from FAO (kg per capita)
         CTOPCP_T(y,r,SIM,t)    percentage changes in detailed crop production output data from FAO (kg per capita)

         CTOPX_T(y,r,SIM,t)    crop consumption (tonnes)
         CTOPXP_T(y,r,SIM,t)   percentage changes in crop consumption (kg)
         SelfS(y,r,SIM,t)  changes in grain self sufficiency rate
         SelfSP(y,r,SIM,t)  percentage changes in grain self sufficiency rate
         ;
$GDXIN 'FAOcrop51.gdx'
$LOAD CTOP=CTOPA CTOPX=CTOPXA
$GDXIN
         CTOP_T(y,r,SIM,t)$(sameas(y,'WHT') and XDr_T(y,r,'BAU',t)) = CTOP(y,r,t)*XDr_T(y,r,SIM,t)/XDr_T(y,r,'BAU',t) ;
*         CTOPC_T(y,r,SIM,t)$CTOP_T(y,r,SIM,t)= CTOP_T(y,r,SIM,t) - CTOP_T(y,r,'BAU',t) ;
*         CTOPC_T('agr',r,SIM,t)= sum(i$icrop(i), CTOPC_T(i,r,SIM,t) ) ;
*         CTOP_T('agr',r,SIM,t)= sum(i$icrop(i), CTOP_T(i,r,SIM,t) ) ;

         CTOPX_T(y,r,SIM,t)$(sameas(y,'WHT') and Xr_T(y,r,'BAU',t)) = CTOPX(y,r,t)*Xr_T(y,r,SIM,t)/Xr_T(y,r,'BAU',t) ;
*         CTOPX_T('agr',r,SIM,t)= sum(i$icrop(i), CTOPX_T(i,r,SIM,t) ) ;

*Next calculate per capita variables
*set age, sex;
parameter
          pop(r) population from UN population projection 2010 revision (Thousand)

         ;
$GDXIN foodpris51.gdx
$LOAD pop
$GDXIN

$ontext
         CTOPP_T(y,r,SIM,t) = CTOP_T(y,r,SIM,t)*1E3/(pop(r)*1E9) ;
         CTOPCP_T(y,r,SIM,t)$CTOP_T(y,r,'BAU',t) = CTOP_T(y,r,SIM,t)/ CTOP_T(y,r,'BAU',t)*100-100 ;
         CTOPXP_T(y,r,SIM,t)$CTOPX_T(y,r,'BAU',t) = CTOPX_T(y,r,SIM,t)/CTOPX_T(y,r,'BAU',t)*100-100 ;
*         CTOPXP_T(y,r,SIM,t)$pop(r) = CTOPX_T(y,r,SIM,t)*1E3/(pop(r)*1E9) ;
         SelfS(y,r,SIM,t)$CTOPX_T(y,r,SIM,t) = CTOP_T(y,r,SIM,t)/ CTOPX_T(y,r,SIM,t) ;
         SelfSP(y,r,SIM,t)$SelfS(y,r,'BAU',t) = SelfS(y,r,SIM,t)/SelfS(y,r,'BAU',t)*100-100 ;
$offtext

SETS
sr  scenarios  /BAU,s01r00*s01r24,s02r00*s02r24,s03r00*s03r24,s04r00*s04r24,s05r00*s05r24,s06r00*s06r24,s07r00*s07r24,s08r00*s08r24,s09r00*s09r24,s10r00*s10r24,s11r00*s11r24,s12r00*s12r24,s13r00*s13r24,s14r00*s14r24,s15r00*s15r24,s16r00*s16r24/
*tC(t) years /2011,2063*2066/
tC(t) years /2036*2039, 2063*2066/
;
set mapSIM(SIM,sr,t) ;

parameter
YCh(sr,t,r) percentage changes from average yields of region rr
PW_T(y,*,sr,t)     Deviations from BAU: Armington price globally weighted by production of regional crops

;
*$GDXIN 'wheatPD2C.gdx'
*$LOAD YCh
*$GDXIN
$GDXIN %ds%_resPD.gdx
$LOAD mapSIM
$GDXIN

mapSIM('BAU','BAU','2011') = Yes ;

*PW_T(y,'World',sr,t)$(sameas(y,'WHT') and sum((r,SIM)$mapSIM(SIM,sr,t),P_T(y,r,SIM,'2011')))= sum((r,SIM)$mapSIM(SIM,sr,t),P_T(y,r,SIM,'2011')*CTOP(y,r,'2011'))/sum(r,CTOP(y,r,'2011'))*100-100 ;
PW_T(y,r,sr,t)$(sameas(y,'WHT') and sum(SIM$mapSIM(SIM,sr,t),P_T(y,r,SIM,'2011'))) = sum(SIM$mapSIM(SIM,sr,t),P_T(y,r,SIM,'2011'))/P_T(y,r,'BAU','2011')*100-100 ;

set VR  variable names to report in results  /
         CTOP  production output in tonnes
         CTOPPC  production output in tonnes per capita
         CTOPX   consumption in tonnes
         CTOPXPC  consumption in tonnes per capita
         PW  percentage changes of prices compared to Baseline level
         X  percentage changes of consumption comparecd to baseline level or
         XD  percentage changes of production output comparecd to baseline level
         GDP
         VA  primary income in 1e9 USD2011
         VA_LAB labor income
         VA_CAP capital income
         VA_LND land income
         VA_TAX taxes to the government
         Export net export in tonnes
         ExportV net export in 1e9 USD2011
         SSR percentage changes of self-sufficiency rates compared to baseline
         SW         Share of wheat consumption in total (%)
         SWT        share of consumption of wht agr and fod in total consumption (%)
         TC  total consumption value in 1e9 USD2011 ('INC')
         PW_EX world price weighted by export
         pW_pop world price weighted by population
         /
    VRC(VR)  per capita data /CTOP,CTOPX/
         ;
Parameter
ResVR(y,*,sr,t,VR)  results related to changes in grain self-sufficiency rates
ResEXR(r,rr,sr,t)  total export volume of wheat from r to rr in billion USD 2011
ResEXRV(r,rr,sr,t)  total export values of wheat from r to rr in billion USD 2011
*ResOth(r,sr,t,VR)  other results for analysis

         ;
ResVR(y,r,sr,t,'CTOP')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t),CTOP_T(y,r,SIM,'2011')) ;
*ResVR(y,r,sr,t,'CTOPPC')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t),CTOPP_T(y,r,SIM,'2011')) ;
ResVR(y,r,sr,t,'CTOPX')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t),CTOPX_T(y,r,SIM,'2011')) ;
*ResVR(y,r,sr,t,'CTOPXPC')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t),CTOPXP_T(y,r,SIM,'2011')) ;
ResVR(y,r,sr,t,'PW') = PW_T(y,r,sr,t) ;

*ResVR(y,r,sr,t,'X')$(sameas(y,'WHT') and tC(t) and (not sameas(sr,'BAU'))) = sum(SIM$mapSIM(SIM,sr,t),Xr_T(y,r,SIM,'2011'))/Xr_T(y,r,'BAU','2011')*100-100 ;       ResVR(y,r,sr,t,'X')$(ResVR(y,r,sr,t,'X')=-100) = 0 ;
*ResVR(y,r,sr,t,'XD')$(sameas(y,'WHT') and tC(t) and (not sameas(sr,'BAU'))) = sum(SIM$mapSIM(SIM,sr,t),XDr_T(y,r,SIM,'2011'))/XDr_T(y,r,'BAU','2011')*100-100 ;    ResVR(y,r,sr,t,'XD')$(ResVR(y,r,sr,t,'XD')=-100) = 0 ;
*ResVR('CPG',r,sr,t,'GDP') = sum(SIM$mapSIM(SIM,sr,t),GDP_D_T(r,SIM,'2011')) ;
ResVR(y,r,sr,t,'VA')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t),VA_T(y,r,SIM,'2011')) ;  ResVR(y,'World',sr,t,'VA')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t),VA_T(y,'World',SIM,'2011')) ;
ResVR(y,r,sr,t,'VA_LAB')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t), GDP_I('LAB',y,'ALL',r,SIM,'2011')) ;
ResVR(y,r,sr,t,'VA_CAP')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t), GDP_I('CAP',y,'ALL',r,SIM,'2011')+GDP_I('RES',y,'ALL',r,SIM,'2011')) ;
*ResVR(y,r,sr,t,'VA_LND')$sameas(y,'WHT') = sum(SIM$mapSIM(SIM,sr,t), GDP_I('RES',y,'ALL',r,SIM,'2011')) ;
ResVR(y,r,sr,t,'VA_TAX')$sameas(y,'WHT') = ResVR(y,r,sr,t,'VA') - ResVR(y,r,sr,t,'VA_LAB')-ResVR(y,r,sr,t,'VA_CAP')-ResVR(y,r,sr,t,'VA_LND') ;


ResVR(y,'World',sr,t,VR)$VRC(VR)=sum(r$(not sameas(r,'World')), ResVR(y,r,sr,t,VR));
ResVR(y,'World',sr,t,'PW')$ResVR(y,'World',sr,t,'CTOP')= sum(r$(not sameas(r,'World')),ResVR(y,r,sr,t,'PW')*ResVR(y,r,sr,t,'CTOP'))/sum(r$(not sameas(r,'World')),ResVR(y,r,sr,t,'CTOP')) ;
*ResVR(y,r,sr,t,'CTOPX')$ResVR(y,'World',sr,t,'CTOP') = ResVR(y,r,sr,t,'CTOPX') * ResVR(y,'World',sr,t,'CTOP')/ResVR(y,'World',sr,t,'CTOPX') ;

ResVR(y,r,sr,t,'Export') = ResVR(y,r,sr,t,'CTOP') - ResVR(y,r,sr,t,'CTOPX') ;
*ResVR(y,r,sr,t,'SSR')$sameas(y,'WHT')  =  sum(SIM$mapSIM(SIM,sr,t), SelfSP(y,r,SIM,'2011')) ;   ResVR(y,r,sr,t,'SSR')$(ResVR(y,r,sr,t,'SSR')=-100) = 0 ;

EXRr_T(i,r,rr,SIM,t)$(not sameas(i,'WHT'))  = 0 ;
ResEXR(r,rr,sr,t)$tC(t) = sum(SIM$mapSIM(SIM,sr,t),EXRr_T('wht',r,rr,SIM,'2011')) ;
ResEXRV(r,rr,sr,t)$tC(t) = sum(SIM$mapSIM(SIM,sr,t),EXRrV_T('wht',r,rr,SIM,'2011')) ;

ResVR('CPG',r,sr,t,'ExportV') =  sum((rr,SIM)$mapSIM(SIM,sr,t),PD_T('wht',r,SIM,'2011')*EXRr_T('wht',r,rr,SIM,'2011')) ;
ResVR('CPG',r,sr,t,'GDP') =  sum(SIM$mapSIM(SIM,sr,t),GDPr_T(r,SIM,'2011')) ;
ResVR('CPG',r,sr,t,'SW') =  sum(SIM$mapSIM(SIM,sr,t),ShWhtFD(r,SIM,'2011')) ;
ResVR('CPG',r,sr,t,'SWT') =  sum(SIM$mapSIM(SIM,sr,t),FDwWht(r,SIM,'2011')) ;
ResVR('CPG',r,sr,t,'TC') =  sum(SIM$mapSIM(SIM,sr,t),INC_T('PRI',r,SIM,'2011')) ;

ResVR(y,r,sr,'2011','X')$(not sameas(sr,'BAU')) = 0 ;
ResVR(y,r,sr,'2011','XD')$(not sameas(sr,'BAU')) = 0 ;
ResVR(y,'World',sr,'2011','X')$(not sameas(sr,'BAU')) = 0 ;
ResVR(y,'World',sr,'2011','XD')$(not sameas(sr,'BAU')) = 0 ;

ResVR(y,'World',sr,t,'PW_EX')$ResVR(y,'World',sr,t,'CTOP') = sum(r$(not sameas(r,'World') and (ResVR('CPG',r,sr,t,'ExportV')>0)),ResVR(y,r,sr,t,'PW')*ResVR('CPG',r,sr,t,'ExportV'))
                                                 /sum(r$(not sameas(r,'World') and (ResVR('CPG',r,sr,t,'ExportV')>0)),ResVR('CPG',r,sr,t,'ExportV')) ;
ResVR(y,'World',sr,t,'PW_pop')$ResVR(y,'World',sr,t,'CTOP')= sum(r$(not sameas(r,'World')),ResVR(y,r,sr,t,'PW')*pop(r))/sum(r$(not sameas(r,'World')),pop(r)) ;

ResVR('CPG',r,sr,t,'ExportV') =  sum(SIM$mapSIM(SIM,sr,t),EXPORT_T('wht',r,SIM,'2011')) ;

*execute_unload 'report51PD.gdx', y, r,f, j, SIM, mapSIM,  GDP_I;
execute_unload 'report51PD.gdx', y, r,f, j, SIM, VR, mapSIM, ResVR;

*execute 'gdxxrw.exe report44.gdx o=report44.xls par=ResVR' ;