* variables in the main model
* Numerairre to make price consistent over time
         PINDEXW_T(SIM,t)   =          PINDEXW.L   ;

* sectors
         FD_T(j,r,SIM,t)$jr(j,r)   =            FD.L(j,r)   ;
         FD_ABS_T(y,j,r,SIM,t)$(jr(j,r) and FDZ(y,j,r))   =      FD_ABS.L(y,j,r)     ;
         IFAF_T(f,a,y,r,SIM,t)$(IFAZ(f,a,y,r) and (FreeF_SH(f,r)<>0) and ay(a,y))   =         IFAF.L(f,a,y,r)   ;
*         IFAF_noC_T(f,a,y,r,SIM,t)$(i(y) and IFAZ(f,a,y,r) and (FreeF_SH_noC(f,r)<>0) and fcap(f) and ay(a,y) and (not icol(y))) = IFAF_noC.L(f,a,y,r) ;
*         IFAF_noF_T(f,a,y,r,SIM,t)$(i(y) and IFAZ(f,a,y,r) and (FreeF_SH_noF(f,r)<>0) and fcap(f) and ay(a,y) and (not ief(y))) = IFAF_noF.L(f,a,y,r) ;
         IM_T(y,r,SIM,t)$IMZ(y,r)   =            IM.L(y,r)   ;
         IO_ABS_T(y,a,yy,r,SIM,t)$(IOAZ(y,a,yy,r) and ay(a,yy))   =     IO_ABS.L(y,a,yy,r)   ;
         IF_ABS_T(f,a,y,r,SIM,t)$(IFAZ(f,a,y,r) and ay(a,y))   =  IF_ABS.L(f,a,y,r)  ;
         TRANS_T(SIM,t)$(SUM((i,r),TRANRZ(i,r)))   =             TRANS.L   ;
         X_T(y,r,SIM,t)$XZ(y,r)   =             X.L(y,r)   ;
         XD_T(y,r,SIM,t)$XDZ(y,r)   =            XD.L(y,r)   ;
         XDA_T(a,y,r,SIM,t)$(ay(a,y) and XDAZ(a,y,r))   =            XDA.L(a,y,r)   ;
* COMMODITIES
         P_C_T(y,j,r,SIM,t)$(jr(j,r) and FDZ(y,j,r))   =         P_C.L(y,j,r)  /PINDEXW.L   ;
         P_IO_T(y,a,yy,r,SIM,t)$(IOAZ(y,a,yy,r) and ay(a,yy))   =       P_IO.L(y,a,yy,r)  /PINDEXW.L   ;
         P_T(y,r,SIM,t)$XZ(y,r)   =             P.L(y,r)  /PINDEXW.L   ;
         PCAP_T(SIM,t)   =              PCAP.L  /PINDEXW.L   ;
         PD_T(y,r,SIM,t)$XDZ(y,r)   =            PD.L(y,r)  /PINDEXW.L   ;
         PDA_T(a,y,r,SIM,t)$(ay(a,y) and XDAZ(a,y,r))   =            PDA.L(a,y,r)  /PINDEXW.L  ;
         PF_T(f,r,SIM,t)$FTOT(f,r)   =            PF.L(f,r)  /PINDEXW.L   ;
         PFF_T(f,r,SIM,t)$(FreeF_SH(f,r)<>0)   =            PFF.L(f,r)  /PINDEXW.L  ;
         PFF_noC_T(f,r,SIM,t)$(FreeF_SH_noC(f,r)<>0 and fcap(f)) = PFF_noC.L(f,r) /PINDEXW.L  ;
         PFF_noF_T(f,r,SIM,t)$(FreeF_SH_noF(f,r)<>0 and fcap(f)) = PFF_noF.L(f,r)  /PINDEXW.L ;
         PFNF_T(f,r,SIM,t)$(1-FreeF_SH(f,r)<>0)   =            PFNF.L(f,r)   /PINDEXW.L ;
         PFA_T(f,a,y,r,SIM,t)$(IFAZ(f,a,y,r) and ay(a,y))   =         PFA.L(f,a,y,r)  /PINDEXW.L  ;
         PFD_T(j,r,SIM,t)$jr(j,r)   =           PFD.L(j,r)  /PINDEXW.L   ;
         PI_T(SIM,t)   =                PI.L  /PINDEXW.L   ;
         PIM_T(y,r,SIM,t)$IMZ(y,r)   =           PIM.L(y,r)  /PINDEXW.L   ;
         PS_T(r,SIM,t)   =              PS.L(r)  /PINDEXW.L   ;
         PT_T(j,r,SIM,t)$(jcon(j) and jr(j,r))   =            PT.L(j,r)  /PINDEXW.L   ;
         PTRANS_T(SIM,t)$(SUM((i,r),TRANRZ(i,r)))   =            PTRANS.L   /PINDEXW.L  ;
*         PCAPA_T(a,y,r,SIM,t)$(ielc(y) and ay(a,y) and XDAZ(a,y,r))   =          PCAPA.L(a,y,r) ;

         PP_T(r,SIM,t)$C_toggle(r)    =   PP.L(r) /PINDEXW.L  ;

* CONSUMERS:
         GLOBANK_T(SIM,t)   =           GLOBANK.L   /PINDEXW.L  ;
         GLOBTRST_T(SIM,t)   =          GLOBTRST.L   /PINDEXW.L  ;
         INC_T(j,r,SIM,t)   =           INC.L(j,r)  /PINDEXW.L   ;
         REGINC_T(r,SIM,t)   =          REGINC.L(r)  /PINDEXW.L   ;
* AUXILIARY
         INV_T(r,SIM,t)   =             INV.L(r)   ;
         RORE_T(r,SIM,t)   =            RORE.L(r)   ;
         ROREZ_T(r,SIM,t)   =            ROREZ(r)   ;
         VKE_T(r,SIM,t)   =             VKE.L(r)   ;

         GDP_T(r,SIM,t)   =             GDP.L(r)   ;
         EMSTOT_T(r,SIM,t)   =          EMSTOT.L(r)   ;
         PINDEX_T(r,SIM,t)   =          PINDEX.L(r)   ;
         FOSS_T(y,r,SIM,t)   =          FOSS.L(y,r)   ;
*         EFFTFP_T(f,r,SIM,t)$(ord(t)=1)  =         1;
         EFFTFP_S_T(f,i,r,SIM,t)$(ord(t)=1 and fnat(f) and IFZ(f,i,r))  = 1;
         EFFTFP_T(f,r,SIM,t)$(not fnat(f))  =         EFFTFP.L(f,r);
         EFFTFP_S_T(f,i,r,SIM,t)$(fnat(f) and IFZ(f,i,r))  =         EFFTFP_S.L(f,i,r);
         CAPA_T(a,y,r,SIM,t)$(ielc(y) and ay(a,y) and XDAZ(a,y,r)) =         CAPA.L(a,y,r);
         TS_Primen_T(y,r,SIM,t)$(ipep(y) or ielc(y)) =  TS_Primen.L(y,r);

         KHAT_T(r,SIM,t)  = KHAT0.L(r)$calib_KHAT(r) + KHAT(r)$(not calib_KHAT(r)) ;
         FreeF_SH_T(f,r,SIM,t) = FreeF_SH(f,r) ;

         PERMITr_T(r,SIM,t)   =  ( SUM((ghga,ii,a,i)$(ioz(ii,i,r) ),eghga_io(ghga,ii,i,r)*IO_ABS.L(ii,a,i,r))
                                    + sum((ghga,ie,j)$(jr(j,r)),eghga_hh(ghga,ie,j,r)*FD_ABS.L(ie,j,r)) ) ;

         PERMIT_T(r,SIM,t)$(not C_toggle(r) and ord(t)=1)   =  ( SUM((ghga,ii,a,i)$(ioz(ii,i,r)),eghga_io(ghga,ii,i,r)*IO_ABS.L(ii,a,i,r))
                                    + sum((ghga,ie,j)$(jr(j,r)),eghga_hh(ghga,ie,j,r)*FD_ABS.L(ie,j,r)) )
                                         / (SUM((ghga,ii,i)$ioz(ii,i,r),eghga_io(ghga,ii,i,r))+ sum((ghga,i,j)$jr(j,r),eghga_hh(ghga,i,j,r)) ) ;
         PERMIT_T(r,SIM,t)$(not C_toggle(r) and ord(t)>1)   =  ( SUM((ghga,ii,a,i)$ioz(ii,i,r),eghga_io(ghga,ii,i,r)*IO_ABS.L(ii,a,i,r))
                                    + sum((ghga,ie,j)$jr(j,r),eghga_hh(ghga,ie,j,r)*FD_ABS.L(ie,j,r)) )
                                         / (SUM((ghga,ii,i)$ioz(ii,i,r),eghga_io(ghga,ii,i,r))+ sum((ghga,i,j)$jr(j,r),eghga_hh(ghga,i,j,r)) )
                                         *eghga_abz_T(r,t)   ;
         PERMIT_T(r,SIM,t)$C_toggle(r)   =          PERMITZ(r)*eghga_abz_T(r,t)   ;

*         PERMITrN_T(r,SIM,t)$(ord(SIM)>1)   =  PERMITr_T(r,SIM,t) - PERMITZ(r)* PERMITr_T(r,'BAU','2004')  ;

         tflnd_T(f,y,r,SIM,t)$(C_toggle(r) and fnat(f) and iagr(y) and rrhh(r)) =  -(tflnd.L(f,y,r)+tf(f,y,r));


* Reported levels of variables
* sectors
         FDr_T(j,r,SIM,t)   =           FD.L(j,r)   *   SUM(i,FDZ(i,j,r)*PFDZ(i,j,r))   ;
         FD_ABSr_T(y,j,r,SIM,t)   =     FD_ABS.L(y,j,r) * FDZ(y,j,r) ;
         IFANFr_T(f,a,y,r,SIM,t)$(((1-FreeF_SH(f,r))<>0) and ay(a,y))   =         IFANFr.L(f,a,y,r)   ;
         IFAr_T(f,a,y,r,SIM,t)$IFAZ(f,a,y,r)   =        IFAZ(f,a,y,r)*IFAF.L(f,a,y,r) * FreeF_SH(f,r) + IFANFr.L(f,a,y,r)     ;
         IMr_T(y,r,SIM,t)   =           IM.L(y,r)   *    IMTOTZ(y,r) ;
         IO_ABSr_T(y,a,yy,r,SIM,t)$(ay(a,yy) and IOAZ(y,a,yy,r))   =    IO_ABS.L(y,a,yy,r) * IOAZ(y,a,yy,r)  ;
         IF_ABSr_T(f,a,y,r,SIM,t)$(IFAZ(f,a,y,r) and ay(a,y))   =  IF_ABS.L(f,a,y,r)*IFAZ(f,a,y,r)  ;
         TRANSr_T(SIM,t)   =            TRANS.L   *     (SUM((i,r),TRANRZ(i,r)))   ;
         Xr_T(y,r,SIM,t)   =            X.L(y,r)   *     XZ(y,r)   ;
         XDr_T(y,r,SIM,t)   =           XD.L(y,r)   *    XDZ(y,r)   ;
         XDAr_T(a,y,r,SIM,t)$(ay(a,y) and XDAZ(a,y,r))   =          XDAZ(a,y,r)*  XDA.L(a,y,r)   ;
* COMMODITIES
         P_IOr_T(y,a,yy,r,SIM,t)$(IOAZ(y,a,yy,r) and ay(a,yy))   =      P_IO.L(y,a,yy,r) *  P_IOZ(y,yy,r) /PINDEXW.L    ;
         PFAr_T(f,a,y,r,SIM,t)$ay(a,y)   =        PFA.L(f,a,y,r)   * PFAZ(f,y,r) /PINDEXW.L  ;
* AUXILIARY
         INVr_T(r,SIM,t)   =         INVZ(r)*INV.L(r)*PI.L/ PINDEXW.L ;
         VKEr_T(r,SIM,t)   =            VKE.L(r)   *     VKEZ(r)   ;
         GDPr_T(r,SIM,t)   =            GDP.L(r)   *     GDPZ(r)   ;
         EMSTOTr_T(r,SIM,t)   =         EMSTOT.L(r)   *  EMSTOTZ(r)   ;
         FOSSr_T(y,r,SIM,t)   =         FOSS.L(y,r)   *  FOSSZ(y,r)   ;
         FTOT_T(f,r,SIM,t)    =    FTOT(f,r)  ;

* Reported levels of variables
         EXRr_T(y,rr,r,SIM,t) =  EXRr.L(y,rr,r) ;
         VTWRr_T(y,r,SIM,t)  =   VTWRr.L(y,r) ;
         XDDr_T(y,r,SIM,t)  = XDDr.L(y,r) ;
         TRANRr_T(y,r,SIM,t)   = TRANRr.L(y,r)  ;

* Calculated parameters
         SAVr_T(r,SIM,t)  = SAVZ(r)*PS.L(r) /PINDEX.L(r)  ;
* Import is calculated at  c.i.f. prices
         IMc_T(y,r,SIM,t) = [ sum(rr,(PD.L(y,rr)*(1+tx(y,rr,r)))*EXRr.L(y,rr,r)) + PTRANS.L*VTWRr.L(y,r) ] /PINDEX.L(r);
* Export is calculated at f.o.b. prices (exports plus export taxes) plus internatinal transpartation margins provided by the region (vst/TRANRZ(i,r))
         EXc_T(y,r,SIM,t) = [ sum(rr,(PD.L(y,r)*(1+tx(y,r,rr)))*EXRr.L(y,r,rr)) + PD.L(y,r)*TRANRr.L(y,r) ] /PINDEX.L(r);
* net export
         EX_IM_BAL_T(r,SIM,t) = sum(i,EXc_T(i,r,SIM,t) - IMc_T(i,r,SIM,t));

* other parameters
         GLOBTRUSTSHARE_T(r,SIM,t) = GLOBTRUSTSHARE(r) ;
         eff_if_T(f,a,y,r,SIM,t)$ay(a,y) =  eff_if(f,a,y,r);
         eff_T(yy,a,y,r,SIM,t)$ay(a,y) =eff(yy,a,y,r);
         eff_f_T(y,j,r,SIM,t) =  eff_f(y,j,r);
         eff_xd_T(a,i,r,SIM,t)$(XDAZ(a,i,r) and ay(a,i)) =  eff_xd(a,i,r) ;
         NFreeF_SH_T(f,a,y,r,SIM,t)=         NFreeF_SH(f,a,y,r) ;

*         eghga_xd_T(ghga,i,r,SIM,t) =   eghga_xd(ghga,i,r)*XD.L(i,r) ;
         eghga_hh2_T(ghga,r,j,SIM,t) =    sum((ie)$(FDZ(ie,j,r)),eghga_hh(ghga,ie,j,r)*FD_ABS.L(ie,j,r)) ;
*         eghga_FD_T(ghga,r,j,SIM,t) =    eghga_FD(ghga,j,r)*FD.L(j,r) ;
         eghga_io_T(ghga,ii,a,i,r,SIM,t)$(IOAZ(ii,a,i,r) and ay(a,i)) =   eghga_io(ghga,ii,i,r)*IO_ABS.L(ii,a,i,r)*IOAZ(ii,a,i,r)/sum(aa$ay(aa,i),IOAZ(ii,aa,i,r)) ;
*         eghga_f_T(ghga,f,i,r,SIM,t) =   eghga_f(ghga,f,i,r)*IF_ABS.L(f,i,r) ;

*          c_LAB_T(r,SIM,t) = c_LAB(r);
*          c_CAP_T(r,SIM,t) = c_CAP(r);
          eKLE_T(y,r,SIM,t) = eKLE(y,r);
          eNEL_T(y,r,SIM,t) = eNEL(y,r);

*         CO2_ffA_T(r,SIM,t) =  PERMITr_T(r,SIM,t) / pop(r,t) ;
*         CO2_ffA_T('World',SIM,t) =  sum(r,PERMITr_T(r,SIM,t)) / sum(r,pop(r,t)) ;

* macroeconomic variables
         GDP_I(f,i,a,r,SIM,t)$ay(a,i) = P_IF.L(f,a,i,r)*IF_ABS.L(f,a,i,r)*IFAZ(f,a,i,r)/PINDEX.L(r);
         GDP_I('ti',i,a,r,SIM,t)$ay(a,i) = sum(ii$ti(ii,i,r),ti(ii,i,r)*P_IO.L(ii,a,i,r)*IO_ABS.L(ii,a,i,r)*IOAZ(ii,a,i,r))/PINDEX.L(r);
         GDP_I('tf',i,a,r,SIM,t)$ay(a,i) = sum(f$tf(f,i,r),tf(f,i,r)*P_IF.L(f,a,i,r)*IF_ABS.L(f,a,i,r)*IFAZ(f,a,i,r))/PINDEX.L(r);
         GDP_I('ty',i,a,r,SIM,t)$ay(a,i) = ty(i,r)*PDA.L(a,i,r)*XDA.L(a,i,r)*XDAZ(a,i,r)/PINDEX.L(r);
*         GDP_I('PP',i,a,r,SIM,t)$(C_toggle(r) and ay(a,i)) = sum(ii$IOAZ(ii,a,i,r),sum(ghga,eghga_IO(ghga,ii,i,r))*IO_ABS.L(ii,a,i,r)*IOAZ(ii,a,i,r)**2/sum(aa$ay(aa,i,r),IOAZ(ii,aa,i,r)) )*PP.L(r)/PINDEX.L(r);
         GDP_I('tflnd',i,a,r,SIM,t)$(C_toggle(r) and ay(a,i)) = sum(f$fnat(f),-tflnd.L(f,i,r)*P_IF.L(f,a,i,r)*IF_ABS.L(f,a,i,r)*IFAZ(f,a,i,r))/PINDEX.L(r);
         GDP_I('TS_Primen',i,a,r,SIM,t)$ay(a,i) =  sum(ii$ipep(ii),P.L(ii,r)*IO_ABS.L(ii,a,i,r)*IOAZ(ii,a,i,r)*TS_Primen.L(ii,r)) ;
         GDP_I('VA',i,a,r,SIM,t)$ay(a,i)= ( PDA.L(a,i,r)*XDA.L(a,i,r)*XDAZ(a,i,r) - sum(ii,P.L(ii,r)*IO_ABS.L(ii,a,i,r)*IOAZ(ii,a,i,r)) ) / PINDEX.L(r)   ;

* Consumption taxes should be excluded from GDP
         GDP_E(j,r,SIM,t) = SUM(i$FDZ(i,j,r),P.L(i,r)*FD_ABS.L(i,j,r)*FDZ(i,j,r))/ PINDEX.L(r);
*         GDP_E(j,a,r,SIM,t)$((ord(j)>2) and aj(a,j,r)) = SUM(i$FDAZ(i,a,j,r),P.L(i,r)*FD_ABS.L(i,a,j,r)*FDAZ(i,a,j,r))/ PINDEX.L(r);
*         GDP_E(j,a,r,SIM,t)$((ord(j)<=2) and aj(a,j,r) ) =  sum[(i,age)$(ageajr(i,age,a,j,r) and (ord(j)<=2) ),P.L(i,r)*FD_ABSCG.L(i,age,a,j,r)*FDCG(i,age,'tot',j,a)]/ PINDEX.L(r) ;
*         GDP_E(age,a,r,SIM,'2012')$rchn(r) = sum[(i,j)$(ageajr(i,age,a,j,r) and sameas(j,'RHH') ),P.L(i,r)*FD_ABSCG.L(i,age,a,j,r)*FDCG(i,age,'tot',j,a)]/ PINDEX.L(r) ;
*         GDP_E(age,a,r,SIM,'2013')$rchn(r) = sum[(i,j)$(ageajr(i,age,a,j,r) and sameas(j,'UHH') ),P.L(i,r)*FD_ABSCG.L(i,age,a,j,r)*FDCG(i,age,'tot',j,a)]/ PINDEX.L(r) ;

* all international tariff (export and import taxes to the country) should be excluded from GDP calculation.
         GDP_E('EX-tariff',r,SIM,t) = sum(i,EXc_T(i,r,SIM,t) - sum(rr$EXRZ(i,r,rr),(PD.L(i,r)*tx(i,r,rr))*EXRr.L(i,r,rr))/PINDEX.L(r));
         GDP_E('IM+tariff',r,SIM,t) = -sum(i,IMc_T(i,r,SIM,t) + sum(rr$EXRZ(i,rr,r),(PD.L(i,rr)*(1+tx(i,rr,r))*tm(i,rr,r))*EXRr.L(i,rr,r) + PTRANS.L*vtwr(i,rr,r)*EXRr.L(i,rr,r)/EXRZ(i,rr,r)*tm(i,rr,r) ) /PINDEX.L(r) );

         CTOP_T(y,r,SIM,t)$XD_T(y,r,'BAU',t) = CTOP(y,r,t)*XD_T(y,r,SIM,t)/XD_T(y,r,'BAU',t) ;
         CTOP_T(y,'World',SIM,t) = sum(r$CTOP(y,r,t),CTOP_T(y,r,SIM,t)) ;