Bilag 4. Gekko kode

Top  Previous 

 

1ymul.gcm:

CLOSE *; CLEAR;

MODE sim;

MODEL Jun19 ;

 

DATE per0 = 2018;  // Last historical period  ;

DATE per1 = 2019;  // First Simulation period ;

DATE per2 = 2100;  // Last simulation period  ;

 

READ lang100;

RUN iweu ;

RUN fVo1;

RUN qo1 ;

RUN lo ;

RUN fIbo ;

RUN fImo ;

RUN Sya ;

RUN tsysp ;

RUN Syv ;

RUN Syc ;

RUN Tirk ;

RUN Sywp ;

RUN Tyrrs ;

RUN Tyrrr ;

RUN Spg ;

RUN Spp_ce ;

RUN Spp_cf ;

RUN Spp_cg ;

RUN sppu ;

RUN Spzej ;

RUN Spzuqr ;

RUN Spr ;

 

EXIT;

 

iweu.gcm:

READ lang100 ;

SERIES <%per1 %per2> iweu   + -0.01 ;

SERIES <%per1 %per2> iwbus  + -0.01 ;

SERIES <%per1 %per2> iuwsd  + -0.01 ;

SERIES <%per1 %per2> iuwse  + -0.01 ;

TIME %per1 %per2 ;

SIM ;

 

fVo1.gcm:

READ lang100;

TIME %per1 %per1  ;

ENDO jrfvmo1; EXO fvmo1;

SERIES fvmo1 = @fvmo1 + 0.001*@Y/@pvmo1 ;

SIM <fix>;

UNFIX;

TIME %per1 %per2  ;

SIM;

 

Qo1.gcm:

READ lang100;

TIME %per1 %per2;

SERIES zQo1 = zQo1 + 0.001*Y[%per1]/(lnakk[%per1]*hgo1[%per1]*.001*klo1[%per1]) ;

SIM ;

 

lo.gcm:

READ lang100;

TIME %per1 %per2;

SERIES klo1 = klo1 + 0.001*Y[%per1]/(lnakk[%per1]*hqo1[%per1]) ;

SIM ;

 

fIbo.gcm:

READ lang100 ;

TIME %per1 %per2;

SERIES zfibo1 = @zfibo1 + 0.001*@y[%per1]/pibo1[%per1] ;

SIM ;

 

fImo.gcm:

READ lang100;

TIME %per1 %per2;

SERIES zfImro1 = zfImro1 + 0.001*@Y[%per1]/pimro1[%per1] ;

SIM ;

 

Sya.gcm:

READ lang100 ;

TIME %per1 %per2 ;

SERIES tsyae = @tsyae - 0.001*Y[%per1]/(ysya[%per1]*ksya[%per1]) ;

SIM ;

 

tsysp.gcm:

READ lang100 ;

TIME %per1 %per2 ;

SERIES tsysp1 = @tsysp1 - 0.001*@Y[%per1]/@Ysp1[%per1] ;

SIM ;

 

Syv.gcm:

READ lang100 ;

TIME %per1 %per2;

SERIES tsyv = @tsyv - 0.001*Y[%per1]/fkncb[%per1] ;

SIM ;

 

Syc.gcm:

READ lang100 ;

TIME %per1 %per2 ;

SERIES tsyc = @tsyc * (@Syc[%per1]-0.001*@Y[%per1])/@Syc[%per1] ;

SIM ;

 

Tirk.gcm:

READ lang100 ;

TIME %per1 %per2 ;

SERIES ttirk = @ttirk - 0.001*Y[%per1]/(ktirk[%per1]*((1-tsyc[%per1])*Yre[%per1]-Tire_o[%per1]-ivsk[%per1])-(1-d4703)*Tiro[%per1])  ;

SIM ;

 

Sywp.gcm:

READ lang100 ;

TIME %per1 %per2 ;

ENDO; EXO;

SERIES Sywp = @Sywp - 0.001*@Y[%per1]  ;

SERIES dtsywp = 1;

EXO Sywp ;

ENDO ztsywp ;

SIM <fix> ;

UNFIX;

 

Tyrrs.gcm:

READ lang100;

TIME %per1 %per2;

SERIES Tyrrrd = Tyrrrd + 0.001*@Y[%per1]/pttyl[%per1] ;

SIM ;

 

Tyrrr.gcm:

READ lang100 ;

TIME %per1 %per2 ;

SERIES Tyrrsd = Tyrrsd + 0.001*@Y[%per1]/pttyl[%per1] ;

SIM ;

 

Spg.gcm:

READ lang100;

TIME %per1 %per2 ;

ENDO;EXO;

SERIES spg = spg - 0.001*y[%per1]  ;

EXO spg ;

ENDO tg ;

SIM <fix> ;

UNFIX ;

 

Spp_ce.gcm:

READ lang100;

TIME %per1 %per2;

SERIES dtp1 = 1;

SERIES dtp2 = 1;

SERIES ztpcex = @ztpcex - 0.001*@Y[%per1]/@fCe[%per1]  ;

SIM ;

 

Spp_cf.gcm:

READ lang100;

TIME %per1 %per2 ;

SERIES dtp1 = 1;

SERIES dtp2 = 1;

SERIES ztpcf = ztpcf - 0.001*Y[%per1]/fCf[%per1] ;

SIM ;

 

Spp_cg.gcm:

READ lang100;

TIME %per1 %per2;

SERIES dtp1 = 1;

SERIES dtp2 = 1;

SERIES ztpcg = @ztpcg - 0.001*@Y[%per1]/@fCg[%per1]  ;

SIM ;

 

Sppu.gcm:

READ lang100 ;

TIME %per1 %per2 ;

SERIES dfcg = 1 ;

SERIES Sppukto = @Sppukto + 0.001*@Y[%per1]  ;

TIME %per1 %per1 ;

SERIES JDtpcs = JDtpcs - 0.001*@Y[%per1]/fcs[%per1]  ;

TIME %per1 %per2 ;

SIM ;

 

Spzej.gcm:

READ lang100 ;

TIME %per1 %per2 ;

SERIES tspzej  = @tspzej*(1 - 0.001*Y[%per1]/Spzej[%per1]) ;

SIM ;

 

Spzuqr.gcm:

READ lang100 ;

TIME %per1 %per2 ;

SERIES Spzuqr = @Spzuqr+0.001*@Y[%per1] ;

SIM ;

 

Spr.gcm:

READ lang100;

TIME %per1 %per2;

SERIES trcb = @trcb - 0.001*@Y[%per1]/(fcb[%per1]*pcb[%per1]) ;

SIM ;