Jump to content

User:Splee/Sandbox

From Wikipedia, the free encyclopedia
< User:Splee

This is an old revision of this page, as edited by Splee (talk | contribs) at 18:33, 13 March 2015. The present address (URL) is a permanent link to this revision, which may differ significantly from the current revision.

  • ===============================================================================
  • import data

$if not declared Debt $include data3


Parameters XDE(i) eXport Demand Elasticity, 

          XXD(i,t) eXtra eXport Demand
          WPM(i,t)
          WPX(i,t)
          WPC(i,t)

PKMULT(j,t);

  • use 100 as a value to signify homogeneity; absolute values.

XDE(i) = 2;


XXD(i,t) = 0; WPM(i,t) = 1; WPX(i,t) = 1; WPC(i,t) = 1; PKMULT(j,t) =1;


set s /1*100/;

$ontext $model:belgium $sectors: 

        q(i,t)$VOM(i)
        costs(i,t)$(VOA(i)*mk0(i))
        qq(i,t)$(VO(i)-VX(i))
        qa(i,t)$(VO(i)+VM(i)-VX(i))
        qdebt(i,t)$VD(i)
        qequity(i,t)$VE(i)
        invest(j,t)
        invprod(t)
        capital(j,t)
        choose(t)$(sum(i,L0T(i))*LM)
        train(i,t)$(LTT(i)*LM)
        labour(i,t)$(L0T(i)*LM)
        m(i,t)$VM(i)
        Export(i,t)$(VX(i) and XDE(i) eq 100)
        u(t)
        qmr(t)$(sum(i,MR(i)))
        v(i,t)$GVA0(i)
        timealloc(t)
        margins(t)

$commodities: 

        pp(i,t)$(VO(i)+VX(i)) ! producer price
        uc(i,t)$(VOA(i)*mk0(i))
        dp(i,t)$VO(i) ! domestic price
        ap(i,t)$(VO(i)-VX(i)) ! armington price
        xp(i,t)$VX(i) ! export price
        mp(i,t)$VM(i) ! import price
        pmr(t)$(sum(i,MR(i)))  ! price of minimum requirements

        pul(t)$LM                !price of unskilled labour
        wl(i,t)$(VL(i)*LM)    !wage of sector specific labour
        pl(i,t)$(VL(i)*LM)    !price of labour
        plt(i)$(VL(i)*LM)
        ip(t)
        w(t)$(LM=0)
        rk(i,t)
        re(i,t)$VE(i)
        rd(i,t)$VD(i)
        pk(j,t)
        pkt(j)
        fe(t)
        gtran(t)
        gtran2(t)$(Closure = HARB  or Closure = PSBR)
        p(t)
        vp(i,t)$GVA0(i)

        PKA(j,t)$PHI
        time(t)
        tmp(t)

$consumers: 

        GOVT(t)
        PRIV(s)
        FC(i,t)$n0(i)
        Exports(i,t)$(VX(i) and XDE(i) < 100)

$auxiliary: 

        nmult(i,t)$n0(i)
        dmshare(i,t)$VO(i)
        dxshare(i,t)$mk0(i)
        mk(i,t)$mk0(i)
        SNP(i,t)$(VOA(i)*mk0(i) and FEE(i)=0)
        QLAB(i,t)$(VL(i) and LM)
        termcon
         QX(i,t)$(XDE(i)<100 and VX(i)>0)
        gdebt(t)$(Closure=PSBR)
        drisk(t)$(Closure=PSBR)
        govprice(t)



$prod:q(i,t)$VOM(i) s:0 t:1 

     O:pp(j,t)    Q:(Eff(i,t)*Make(j,i))  P:((1-mk0(i))/Eff(i,t)) A:FC(i,t)$N0(i) N:mk(i,t)$mk0(i) M:1$mk0(i)
     I:dp(j,t)    Q:(Use(j,i)/InputEff(j,i,t)*(1-mk0(i)))  P:(InputEff(j,i,t)*(1+ut0(j,i)))  A:GOVT(t)  T:ut(j,i,t)
     I:vp(i,t)    Q:(GVA0(i)*(1-mk0(i)))

$prod:costs(i,t)$(VOA(i)*mk0(i)) s:ElasG(i) 

     O:uc(i,t)#(s)    Q:(VOM(i)*mk0(i)/card(s)) A:PRIV(s) N:SNP(i,t)$(FEE(i)=0)
     I:dp(j,t)    Q:(Use(j,i)    /InputEff(j,i,t)*mk0(i)) P:(InputEff(j,i,t)*(1+ut0(j,i)))  A:GOVT(t)  T:ut(j,i,t)
     I:vp(i,t)    Q:(GVA0(i)*mk0(i))

$report: 

     V:IntUse(j,i,t)        I:dp(j,t)      PROD:q(i,t)

$prod:v(i,t)$GVA0(i) s:ElasSub(i) cap:eDebtEquity(i) 

     O:vp(i,t)          Q:GVA0(i)     P:(1-PT0(i))  A:GOVT(t) T:PT(i,t)
     I:wl(i,t)$LM       Q:VL(i)  A:GOVT(t) P:(1+ET0(i))  T:ET(i,t)
     I:w(t)$(LM=0)      Q:VL(i)   A:GOVT(t) P:(1+ET0(i))  T:ET(i,t)
  • I:rk(i,t) Q:VK(i) A:GOVT(t) P:(1+CT0(i)) T:CT(i,t)
  • Debt:
     I:rd(i,t)   Q:VD(i)    cap:
  • Equity:
     I:re(i,t)   Q:VE(i)   P:(1+ct0(i))  A:GOVT(t) T:ct(i,t) cap:


$demand:FC(i,t)$n0(i) 

       D:uc(i,t)

$prod:qdebt(i,t)$VD(i) 

        O:rd(i,t)    Q:(VD(i)*debtprod(i,t))
        I:rk(i,t)    Q:VD(i)

$prod:qequity(i,t)$VE(i) 

        O:re(i,t)    Q:(VE(i)*equityprod(i,t))
        I:rk(i,t)    Q:VE(i)

$prod:qq(i,t)$(VO(i)-VX(i)) t:1 

       O:ap(i,t) Q:(VO(i)-VX(i))
       O:xp(i,t) Q:VX(i)
       I:pp(i,t) Q:VO(i)

$prod:qa(i,t)$(VO(i)+VM(i)-VX(i)) s:Elasm(i) 

     O:dp(i,t) Q:(VO(i)+VM(i)-VX(i)+max(0,supply(i,"margin")) )
     I:ap(i,t) Q:(VO(i)-VX(i))
     I:mp(i,t) Q:VM(i)
     I:tmp(t)   Q:(max(0,supply(i,"margin")))


$prod:invest(j,t) s:0 ADJ:1 gds(ADJ):0 

     O:pk(j,t)    Q:I0(j)
     I:ip(t)    Q:(I0(j)/J0(j)*kshare(j)*TI)
     I:ip(t)    Q:(I0(j)/J0(j)*kshare(j)*TI*(PHI/2)*(d0+g0)/(1+g0)) gds:
     I:pka(j,t)   Q:(I0(j)*(PHI/2)*(d0+g0)/(1+g0)) ADJ:


$prod:invprod(t) 

     O:ip(t)     Q:TI
     I:dp(i,t)   Q:VI(i)  A:GOVT(t) P:(1+ut0(i,"inv")) T:ut(i,"inv",t)

$report: 

     V:inv_q(i,t)  I:dp(i,t)  PROD:invprod(t)

$prod:capital(j,t) s:0 

     O:pk(j,t+1)        Q:((K0(j)+I0(j))*(1-d0))
     O:pkt(j)$tlast(t)  Q:((K0(j)+I0(j))*(1-d0))
     O:rk(j,t)          Q:((K0(j)+I0(j))*RK0*CapEff(j,t))
     O:pka(j,t)         Q:((K0(j)+I0(j))*PKA0)
     I:pk(j,t)          Q:(( K0(j)+I0(j))*PKMULT(j,t))

$prod:choose(t)$(sum(i,L0T(i))*LM) t:3 

     O:pl(i,t)   Q:(L0T(i)*(LU0(i)+g0))
     I:pul(t)    Q:(sum(i,L0T(i)*(LU0(i)+g0)))

$report: 

     V:choosesec(i,t)$LM O:pl(i,t) PROD:choose(t)

$prod:train(i,t)$(LTT(i)*LM) 

     O:pl(i,t+1)       Q:LTT(i)
     O:plt(i)$tlast(t) Q:LTT(i)
     I:wl(i,t)         Q:(VL(i)*.2)

$prod:labour(i,t)$(L0T(i)*LM) 

     O:pl(i,t+1)       Q:(L0T(i)*(1-LD0(i)))
     O:plt(i)$tlast(t) Q:(L0T(i)*(1-LD0(i)))
     O:wl(i,t)         Q:(VL(i)*1.2)
     I:pl(i,t)         Q:L0T(i)

$prod:m(i,t)$VM(i) 

     O:mp(i,t)     Q:VM(i)
     I:FE(t)       Q:(VM(i)*WPM(i,t))

$prod:Export(i,t)$(VX(i) and XDE(i)=100) s:0 

     O:FE(t)       Q:(VX(i)*WPX(i,t))
     I:xp(i,t)     Q:VX(i)

$demand:Exports(i,t)$(XDE(i)<100 and VX(i)>0) 

     E:xp(i,t)     Q:(-VX(i)*QREF(t))    R:QX(i,t)
     E:xp(i,t)     Q:(-XXD(i,t))
     D:FE(t)       Q:(VX(i)*QREF(t)+XXD(i,t))
     E:FE(t)       Q:(2*(VX(i)*QREF(t)+XXD(i,t)))

$constraint:QX(i,t)$(XDE(i)<100 and VX(i)>0) 

        QX(i,t) =e= (xp(i,t)/fe(t))**(-XDE(i));

$prod:Margins(t) s:0 

        O:tmp(t)  Q:(sum(i,max(0,-supply(i,"margin"))))
        I:dp(i,t) Q:(max(0,-supply(i,"margin")))

$demand:GOVT(t) s:1 

     E:dp(i,t)$(Closure <> SPEND) Q:(-VG(i)*QREF(t))
     D:dp(i,t)$(Closure =  SPEND) Q:( VG(i)*QREF(t))  P:PREF(t)

     D:GTRAN2(t)$(Closure = HARB or Closure = PSBR) Q:(1e4/scale*QREF(t))
     E:GTRAN(t)$(Closure = HARB or Closure = PSBR)  Q:(1e4/scale*QREF(t))
     E:GTRAN(t)                   Q:((GovtDef-GovtSurp)*QREF(t))

$report: 

     V:govt_d(i,t)  D:dp(i,t) DEMAND:GOVT(t)

$constraint:gdebt(t)$(Closure=PSBR) 

     gdebt(t) =e= (govdebt*1000*govprice(t))$(ord(t)=1)
                  +gdebt(t-1)$(ord(t)>1) *(1+govrate)
                  - GOVT(t)+GTRAN(t)*(10+GovtDef-GovtSurp)*QREF(t)
                  + DebtShock(t)*1000*govprice(t)/scale;


$constraint:drisk(t)$(Closure=PSBR)

  • this is calibrated so that every £100bn increases risk premia by 1% times elasdebtrisk.
        drisk(t) =e= [max[0,gdebt(t)/govprice(t)-BaseDebt(t)]/scale/100000]*elasdebtrisk;

$constraint:govprice(t) 

        govprice(t) =e= sum[i,dp(i,t)*VI(i)]/sum[i,VI(i)] /PREF(t);

$prod:u(t) s:LabLeis gds:elasG("hh") 

     O:p(t)   Q:(sum(i,(VH(i)-MR(i))*(1+ut0(i,"hse")))+Leisure)
     I:dp(i,t)   Q:(VH(i)-MR(i)) P:(1+ut0(i,"hse")*WPC(i,t))  A:GOVT(t)  T:ut(i,"hse",t)
  • leisure demand
     I:time(t)  Q:Leisure

$report: 

    V:con_d(i,t)   I:dp(i,t) PROD:u(t)
    V:con_m(i,t)   I:mp(i,t) PROD:u(t)
    V:qleisure(t)  I:time(t) PROD:u(t)
    V:qlabour(t)   I:time(t) PROD:timealloc(t)

$prod:timealloc(t) 

     O:pul(t)$LM     Q:TL
     O:w(t)$(LM=0)   Q:TL
     I:time(t)       Q:TL


$prod:qmr(t)$(sum(i,abs(MR(i)))) s:0 

    O:pmr(t)  Q:(sum(i,  MR(i)*ut0(i,"hse") ))
    I:dp(i,t)   Q:MR(i)  P:(1+ut0(i,"hse"))  A:GOVT(t)  T:ut(i,"hse",t)


$demand:PRIV(s) s:inter

  • consumption
     D:p(t)     Q:((sum(i,(VH(i)-MR(i))*(1+ut0(i,"hse")))+Leisure)*QREF(t)/card(s))    P:PREF(t)
  • time endowment
     E:time(t)     Q:((Leisure+TL)*QREF(t)/card(s))
  • minimum requirements
     E:pmr(t)    Q:(-sum(i,  MR(i)*ut0(i,"hse") )    *QREF(t)/card(s))

     E:pul(t)$LM    Q:(-TL*QREF(t)/card(s))
  • Labour earnings:
  • E:pl(i,tfirst)$LM Q:((L0(i)*(1-LU0(i)-g0))/card(s))
  • E:pul(tfirst)$LM Q:(sum(i,L0(i)*LU0(i)*(1-lost))/card(s))
  • E:plt(i)$LM Q:(-(L0(i)*(1-LD0(i)) + LT0(i)) /(1+g0)/card(s)) R:TERMCON
  • E:pl(i,t)$(tnotfirst(t) and LM) Q:(L0(i)*(1-LU0(i)-g0)*(1+g0) - (LT0(i) + L0(i)*(1-LD0(i)))) R:Qlab(i,t)
  • E:pul(t)#(i)$(tnotfirst(t) and LM) Q:(L0(i)*LU0(i)*(1+g0)*(1-lost)) R:Qlab(i,t)
  • E:pul(t)$LM Q:(sum(i,L0(i)*(g0+LU0(i)*lost))*QREF(t))


  • Factor earnings:Capital
     E:pk(j,tfirst)     Q:(K0(j)/card(s))
     E:pkt(j)           Q:(-K0(j)/card(s))  R:TERMCON
  • Fixed demands/endowments:
     E:dp(i,t)    Q:(-VS(i)*QREF(t)/card(s))
  • Trade balance:
  • D:FE(t) Q:(TradeBal*QREF(t)/card(s))
     E:FE(t)      Q:(-TradeBal*QREF(t)/card(s))
  • Varietal scaling:
     E:pp(j,t)#(i)$(VSE(i)*n0(i))  Q:(Eff(i,t)*Make(j,i)/card(s))  R:nmult(i,t)
     E:pp(j,t)#(i)$(VSE(i)*n0(i))  Q:(-Eff(i,t)*Make(j,i)/card(s))
  • Government deficit
     E:GTRAN(t)                  Q:((GovtSurp-GovtDef)*QREF(t)/card(s))
     E:GTRAN(t)$(Closure = HARB or Closure = PSBR) Q:(-1e4/scale*QREF(t)/card(s))
     E:GTRAN2(t)$(Closure = HARB or Closure = PSBR)   Q:(1e4/scale*QREF(t)/card(s))



$report: 

    V:WELF(s)  W:PRIV(s)


$constraint:SNP(i,t)$(VOA(i)*mk0(i) and FEE(i)=0) 

        QREF(t) =E= costs(i,t);

$constraint:nmult(i,t)$n0(i) 

        nmult(i,t) =E= costs(i,t) / QREF(t);

$constraint:dxshare(i,t)$(VO(i)-VX(i)) 

       dxshare(i,t) =E= (VO(i)-VX(i))/VO(i)* PP(i,t)**(ElasM(i)  -1) * XP(i,t)**(1-ElasM(i)  );

$constraint:dmshare(i,t)$VDB(i) 

       dmshare(i,t) =E= (VO(i)-VX(i))/VDB(i)
       * [DP(i,t)/MP(i,t)]**(ElasM(i)-1);

$constraint:mk(i,t)$mk0(i) 

       mk(i,t) =E= (1+CV(i))*[(1/ElasF(i)) + (1/(n0(i)*nmult(i,t)*PercN(i,t)))
               * sum[j, Make(j,i)/VOM(i) * dxshare(j,t)
               *[(1/ElasM(j)) +dmshare(j,t) * (-1/ElasP(j)) - (1/ElasM(j))]]
               - (1/ElasF(i))];

$constraint:QLAB(i,t)$(VL(i) and LM) 

       QLAB(i,t) =E= labour(i,t-1);



  • ===============================================================================
  • TERMINAL CONDITION
  • TC=1 terminal growth path
  • TC=2 fixed terminal stock
  • TC=3 fixed terminal growth path

$constraint:termcon$(TC=1) 

sum[(i,t)$tlast(t), k0(i)* invest(i,t)] /sum[(i,t)$tlast(t), k0(i)* invest(i,t-1)]
 =E= sum[t$tlast(t), u(t)] /sum[t$tlast(t), u(t-1)];

$constraint:termcon$(TC=2) 

      termcon =E= qreft*TCM;

$constraint:termcon$(TC=3) 

         [sum[(i,t)$tlast(t), k0(i)*capital(i,t)]
         /sum[(i,t)$tlast(t), k0(i)*capital(i,t-1)] ]=E= 1+g0*TCM ;


$offtext $sysinclude mpsgeset belgium $offwarning

  • numeraire

fe.fx("%firstyear%") = 1;


  • initial levels of auxilliary variables:

nmult.l(i,t) = 1; dxshare.l(i,t)$VO(i) = (VO(i)-VX(i))/VO(i); dmshare.l(i,t)$VDB(i) = (VO(i)-VX(i))/VDB(i); mk.l(i,t) = mk0(i); termcon.l = qreft; qLab.l(i,t) = QREF(t-1); qlab.fx(i,tfirst) = 0;

qx.l(i,t) = 1; dp.lo(j,t) = .01; gdebt.l(t) = BaseDebt(t); govprice.l(t) = 1;

  • initial levels of quantity variables:

q.l(i,t) = QREF(t); u.l(t) = QREF(t); costs.l(i,t) = QREF(t); qa.l(i,t) = QREF(t); invest.l(j,t) = QREF(t); capital.l(j,t) = QREF(t); choose.l(t) = QREF(t); train.l(i,t) = QREF(t); labour.l(i,t) = QREF(t); m.l(i,t) = QREF(t); Export.l(i,t) = QREF(t); v.l(i,t) = QREF(t); timealloc.l(t) = QREF(t); qmr.l(t) = QREF(t); invprod.l(t) = QREF(t); qq.l(i,t) = QREF(t); qdebt.l(i,t) = QREF(t); qequity.l(i,t) = QREF(t); margins.l(t) = QREF(t);

  • initial levels of price variables:

ip.l(t) = pref(t); pp.l(i,t) = pref(t); p.l(t) = pref(t); uc.l(i,t) = pref(t); dp.l(i,t) = pref(t); xp.l(i,t) = pref(t); mp.l(i,t) = pref(t); pul.l(t) = pref(t); wl.l(i,t) = pref(t); pl.l(i,t) = pref(t); rk.l(j,t) = pref(t); pmr.l(t) = pref(t); fe.l(t) = pref(t); w.l(t) = pref(t); vp.l(i,t) = pref(t); gtran.l(t) = pref(t); gtran2.l(t) = pref(t); pka.l(i,t)= pref(t); time.l(t) = pref(t); ap.l(i,t) = pref(t); re.l(i,t) = pref(t); rd.l(i,t) = pref(t); tmp.l(t) = pref(t);

plt.l(i) = preft; plt.l(i) = preft; pk.l(j,t) = pref(t)*pk0; pkt.l(j) = preft*pk0;

$libinclude store these q v fe rk re rd qdebt qequity u w dp p invest capital labour m Export xp mp qlabour qleisure

$libinclude store calc Welfare = sum[s,Welf.l(s)]/card(s)*(TL+Leisure)*scale;

  • for these calculations, R_ = real; N_ = nominal; _PI = price index
  • _C = consumption; _G = government; _I = investment; _X = exports; _M = imports
  • _Y = GDP
  • in some cases domestic (_D) and imported (_D) parameters are calculated.

$libinclude store calc RC (t) = scale*sum[i,[con_d.l(i,t)+MR(i)*QREF(t)]]; $libinclude store calc S_RC (i,t) = scale*[con_d.l(i,t)+MR(i)*QREF(t)]; $libinclude store calc NC (t) = scale*sum[i,[con_d.l(i,t)+MR(i)*QREF(t)]*dp.l(i,t)]; $libinclude store calc CPI (t)= NC(t)/RC(t);

$libinclude store calc RG (t) = scale*sum[i,VG(i)$(Closure<>Spend)*QREF(t)+govt_d.l(i,t)$(Closure=Spend)]; $libinclude store calc S_RG (i,t) = scale*[VG(i)$(Closure<>Spend)*QREF(t)+govt_d.l(i,t)$(Closure=Spend)]; $libinclude store calc NG (t) = scale*sum[i,dp.l(i,t)*VG(i)*QREF(t)$(Closure<>Spend)+govt_d.l(i,t)$(Closure=Spend)]; $libinclude store calc GPI (t)= NG(t)/RG(t);

$libinclude store calc RI (t) = scale*sum[i,inv_q.l(i,t)]; $libinclude store calc NI (t) = scale*sum[i,inv_q.l(i,t)*dp.l(i,t)*(1+ut(i,"inv",t))]; $libinclude store calc S_RI (i,t) = scale*[inv_q.l(i,t)]; $libinclude store calc IPI (t) = NI(t)/RI(t);


$libinclude store calc RM (t)= scale*sum[i,m.l(i,t)*VM(i)]; $libinclude store calc NM (t)= scale*sum[i,m.l(i,t)*VM(i)*mp.l(i,t)]; $libinclude store calc S_RM (i,t)= scale*[m.l(i,t)*VM(i)]; $libinclude store calc MPI (t)=NM(t)/RM(t);

$libinclude store calc QtyX (i,t)= scale*[(qx.l(i,t)*VX(i)*QREF(t)+XXD(i,t))$(XDE(i)<100)+QREF(t)*VX(i)*export.l(i,t)$(XDE(i)=100)]; $libinclude store calc RX (t)= sum[i,QtyX(i,t)]; $libinclude store calc S_RX (i,t)= QtyX(i,t); $libinclude store calc NX (t)= sum[i,QtyX(i,t)*xp.l(i,t)]; $libinclude store calc XPI (t)=NX(t)/RX(t); $libinclude store calc TOT (t)=XPI(t)/MPI(t);


$libinclude store calc RY (t) = RC(t)+RG(t)+RI(t)+RX(t)-RM(t); $libinclude store calc S_RY (i,t) = S_RC(i,t)+S_RG(i,t)+S_RI(i,t)+S_RX(i,t)-S_RM(i,t); $libinclude store calc NY (t) = NC(t)+NG(t)+NI(t)+NX(t)-NM(t); $libinclude store calc YPI (t) = NY(t)/RY(t); $libinclude store calc RER (t)=YPI(t)/FE.L(t);

$libinclude store calc disc_GDP = sum[t, RY(t)];

$libinclude store calc Output (i,t) = v.l(i,t)*GVA0(i,r)*scale; $libinclude store calc LabGVA (i,t) = w.l(t)*VL(i)*SCALE;


belgium.workspace = 100; belgium.WORKSPACE$(CARD(t) gt 5) =250; belgium.WORKSPACE$(CARD(t) gt 10) = 250; belgium.workspace$(CARD(t)*CARD(i) > 200) = 250;

$libinclude store setmodel belgium

belgium.iterlim = 0; $include belgium.gen solve belgium using mcp;

$setglobal solve1 $include belgium.gen $setglobal solve2 solve belgium using mcp;


abort$(belgium.objval gt 30/scale) "Benchmark replication failed!!!" ,belgium.objval, q.m,costs.m,qa.m,capital.m,invest.m,qequity.m,qdebt.m,choose.m,train.m,labour.m,m.m, export.m,u.m,v.m,invprod.m,ip.m,pp.m,uc.m,dp.m,xp.m,mp.m,pul.m,margins.m,tmp.m, wl.m,pl.m,plt.m,w.m,rk.m,re.m,rd.m,pk.m,pkt.m,gtran.m,p.m,time.m,timealloc.m, priv.m,nmult.m,dmshare.m,dxshare.m,mk.m,vp.m,pka.m,mk.m,termcon.m, snp.m,qlab.m,w.m,gtran2.m,gdebt.m,GOVT.m,pk.m,pmr.m,qmr.m,belgium.objval;


belgium.iterlim = 10000; $libinclude store base $libinclude store setbase base

  • insert baseline scenario here:


$libinclude store baseline $libinclude store setbase baseline

$exit

  • if this file is included in another, exit here

$if not "%system.incparent%" == "" $exit


  • standard simulations that will run when this file is run:

fe.fx("%firstyear%") = 1.1;


  • $ontext
  • initial levels of price variables:

ip.l(t) = pref(t)*1.1; pp.l(i,t) = pref(t)*1.1; p.l(t) = pref(t)*1.1; uc.l(i,t) = pref(t)*1.1; dp.l(i,t) = pref(t)*1.1; xp.l(i,t) = pref(t)*1.1; mp.l(i,t) = pref(t)*1.1; pul.l(t) = pref(t)*1.1; wl.l(i,t) = pref(t)*1.1; pl.l(i,t) = pref(t)*1.1; rk.l(j,t) = pref(t)*1.1; pmr.l(t) = pref(t)*1.1; fe.l(t) = pref(t)*1.1; w.l(t) = pref(t)*1.1; vp.l(i,t) = pref(t)*1.1; gtran.l(t) = pref(t)*1.1; gtran2.l(t) = pref(t)*1.1; pka.l(i,t)= pref(t)*1.1; time.l(t) = pref(t)*1.1; ap.l(i,t) = pref(t)*1.1; re.l(i,t) = pref(t)*1.1; rd.l(i,t) = pref(t)*1.1; tmp.l(t) = pref(t)*1.1; GOVPRICE.l(t) = 1.1; plt.l(i) = preft*1.1; plt.l(i) = preft*1.1; pk.l(j,t) = pref(t)*pk0*1.1; pkt.l(j) = preft*pk0*1.1;

  • $offtext

belgium.iterlim = 0; $libinclude store numeraire0

belgium.iterlim = 1000; $libinclude store numeraire1


fe.fx("%firstyear%") = 1; $libinclude store basecheck

set fss (i) /FINS, INSU, FAUX/;


XXD(fss,t) = -.1 * VX(fss)*QREF(t);

$libinclude store lessdemand XXD(i,t) = 0; $libinclude store backtobase

XXD(fss,t) = .1 * VX(fss)*QREF(t); $libinclude store extrademand


$libinclude store display percent changes

Retrieved from "https://en.wikipedia.org/w/index.php?title=User:Splee/Sandbox&oldid=651225194"