* 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)   ;
         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  ;
         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  ;

         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_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)   ;

         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)  = SAVr.L(r)*PS.L(r) /PINDEXW.L  ;

* 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_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)/FD_ABSZ(ie,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)/IO_ABSZ(ii,a,i,r)*IOAZ(ii,a,i,r)/sum(aa$ay(aa,i),IOAZ(ii,aa,i,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) ;