wfcreate u 300
smpl 1 300

' Model parameters and starting values
genr alfa1 = 0.0004

genr drm = 0
genr EX = 0
genr ran1 = 0
genr ZED = 0
genr u = 0.6925
genr a1 = 0.75
genr A1 = 0.0004
genr a2 = 0.064
genr a4 = 0.05
genr add3 = 0.02
genr B1 = 0.08
genr B2 = 0.05
genr delta = 0.10667
genr eta = 0.6
genr gik = 0.03001
genr go = -0.0018
genr lam = 0.27988
genr lam0 = 0.13259
genr lam01 = 0.67132
genr lam02 = -0.04341
genr lam1 = 6.6
genr lam11 = 2.2
genr lam12 = 2.2
genr lam2 = 2.2
genr lam21 = 2.2
genr lam22 = 2.2
genr lam3 = 0.1
genr lam31 = 0.1
genr lam32 = 0.1
genr lam4 = 2.2
genr lam41 = 6.6
genr lam42 = 2.2
genr lam5 = 2.2
genr lam51 = 2.2
genr lam52 = 6.6
genr lamsta = 0.27988
genr mu1 = 0.0153
genr mu2 = 0.05239
genr N = 80.92
genr Nsta = 80.92
genr Nt = 80
genr omega = 0.15
genr pb = 18.182
genr phi1 = 0.15255
genr phi3 = 0.92
genr pi = 0.00609
genr prodg = 0.03
genr psi = 0.05
genr rb = 0.035
genr rbl = 0.055
genr rl = 0.06279
genr rlc = 0.07
genr rm = 0.0104
genr ro = 0.05
genr ro0 = -0.23322
genr ro1 = 1
genr ro2 = 2
genr ro3 = 0.45621
genr rrle = 0.05636
genr theta = 0.22844
genr theta1 = 0.05
genr upsilon = 0.8
genr xi = 0.2
genr zeta = 0.16667
genr zeta1 = 0.16666
genr rk = 0
genr LR = 0.05998
genr p = 12.542
genr NUC = 9.6966
genr UC = 9.6966
genr NHUC = 9.799
genr ed = 1.1988
genr es = 1.1988
genr pe = 7036
genr wsta = 6021.6
genr vk = 1154.639
genr pr = 7.4313
genr consk = 365.081
genr ypde = 391.115
genr gk = 114.580
genr ypdk = 391.114
genr ist = 119.670
genr yk = 601.341
genr tk = 115.799
genr sk = 598.349
genr sek = 598.349
genr dkk = 118.689
genr ik = 102.717
genr ista = 102.717
genr k = 894.419
genr Bs = 2201.551
genr Ld = 2233.399
genr Ls = 2233.399
genr Bpd = 1654.156
genr Bps = 1654.156
genr Bld = 90.979
genr Bls = 90.979
genr W = 72.059
genr Hbd = 125.470
genr Hbs = 125.470
genr Vn = 14252.440
genr V = 14481.380
genr Hpd = 228.941
genr Hps = 228.941
genr Md = 2509.390
genr Ms = 2509.390
genr Bcbd = 354.410
genr Bcbs = 354.410
genr H = 354.410
genr Cons = 4578.814
genr I = 996.012
genr WB = 5830.985
genr Yp = 6357.664
genr T = 1452.345
genr Ypd = 4905.319
genr Y = 7533.453
genr Se = 7504.444
genr G = 1437.048
genr S = 7504.442
genr Fbsta = 115.235
genr Fb = 115.233
genr Fd = 242.596
genr FTe = 1639.125
genr FT = 1647.966
genr Fue = 1321.548
genr Fu = 1330.389
genr DK = 1488.581
genr Bbd = 150.522
genr Bbs = 150.521
smpl 2 300
genr LR = LR(-1)*(1+0.000166638893517757)
genr p = p(-1)*(1+0.00610710930159719)
genr NUC = NUC(-1)*(1+0.0060963260015392)
genr UC = UC(-1)*(1+0.00608613758114007)
genr NHUC = NHUC(-1)*(1+0.00611283695761089)
genr ed = ed(-1)*(1+0.0137688120397055)
genr es = es(-1)*(1+0.0137688120397055)
genr pe = pe(-1)*(1+0.0221453287197231)
genr wsta = wsta(-1)*(1+0.0300004559339808)
genr vk = vk(-1)*(1+0.0299941110312354)
genr pr = pr(-1)*(1+0.0299985222402839)
genr consk = consk(-1)*(1+0.0299977439651067)
genr ypde = ypde(-1)*(1+0.0299967944143256)
genr gk = gk(-1)*(1+0.0299990415641673)
genr ypdk = ypdk(-1)*(1+0.0299992044513289)
genr ist = ist(-1)*(1+0.0300001529449552)
genr yk = yk(-1)*(1+0.0299999543449181)
genr tk = tk(-1)*(1+0.0299993677688564)
genr sk = sk(-1)*(1+0.0300000611777949)
genr sek = sek(-1)*(1+0.0300000611777949)
genr dkk = dkk(-1)*(1+0.0300084042036439)
genr ik = ik(-1)*(1+0.0299976835765579)
genr ista = ista(-1)*(1+0.0299976835765579)
genr k = k(-1)*(1+0.0300096072487153)
genr Bs = Bs(-1)*(1+0.0356447013903545)
genr Ld = Ld(-1)*(1+0.036258857792866)
genr Ls = Ls(-1)*(1+0.036258857792866)
genr Bpd = Bpd(-1)*(1+0.0362608193005673)
genr Bps = Bps(-1)*(1+0.0362608193005673)
genr Bld = Bld(-1)*(1+0.0362607295607424)
genr Bls = Bls(-1)*(1+0.0362607295607424)
genr W = W(-1)*(1+0.0362647777722354)
genr Hbd = Hbd(-1)*(1+0.0362690960658345)
genr Hbs = Hbs(-1)*(1+0.0362690960658345)
genr Vn = Vn(-1)*(1+0.0362685316484268)
genr V = V(-1)*(1+0.0362685233449995)
genr Hpd = Hpd(-1)*(1+0.0362684335098344)
genr Hps = Hps(-1)*(1+0.0362684335098344)
genr Md = Md(-1)*(1+0.0362702742965679)
genr Ms = Ms(-1)*(1+0.0362702742965679)
genr Bcbd = Bcbd(-1)*(1+0.0362712955990152)
genr Bcbs = Bcbs(-1)*(1+0.0362712955990152)
genr H = H(-1)*(1+0.0362712955990152)
genr Cons = Cons(-1)*(1+0.0362714343245936)
genr I = I(-1)*(1+0.03627124852599)
genr WB = WB(-1)*(1+0.0362723374246043)
genr Yp = Yp(-1)*(1+0.0362724788646234)
genr T = T(-1)*(1+0.0362725013194669)
genr Ypd = Ypd(-1)*(1+0.0362724722162939)
genr Y = Y(-1)*(1+0.0362740778081103)
genr Se = Se(-1)*(1+0.0362741975270204)
genr G = G(-1)*(1+0.0362744423669779)
genr S = S(-1)*(1+0.0362740821253926)
genr Fbsta = Fbsta(-1)*(1+0.0362771474246502)
genr Fb = Fb(-1)*(1+0.0362777132251448)
genr Fd = Fd(-1)*(1+0.0362781049746268)
genr FTe = FTe(-1)*(1+0.0362813432320328)
genr FT = FT(-1)*(1+0.0362802678638128)
genr Fue = Fue(-1)*(1+0.0362827086646726)
genr Fu = Fu(-1)*(1+0.0362813675269096)
genr DK = DK(-1)*(1+0.0362818567173622)
genr Bbd = Bbd(-1)*(1+0.0364465236174427)
genr Bbs = Bbs(-1)*(1+0.0364345895887535)
smpl 1 300


model wgm
wgm.append Bpd=Vn(-1)*(lam0+lam1*rb(-1)-lam2*rm(-1)-lam4*rk(-1)-lam5*rbl(-1)-lam3*Ypd/V)
wgm.append pe=Vn(-1)*(lam01-lam11*rb(-1)-lam21*rm(-1)+lam41*rk(-1)-lam51*rbl(-1)-lam31*Ypd/V)/ed
wgm.append Bld=Vn(-1)*(lam02-lam12*rb(-1)-lam22*rm(-1)-lam42*rk(-1)+lam52*rbl(-1)-lam32*Ypd/V)/pb
wgm.append @IDENTITY rk=Fd/(pe(-1)*ed(-1))
wgm.append @IDENTITY pb=1/rbl
wgm.append @IDENTITY ed=es
wgm.append @IDENTITY yk=sek+ista-ik(-1)
wgm.append @IDENTITY Nsta=yk/pr
wgm.append @IDENTITY N=N(-1)+eta*(Nsta-N(-1))
wgm.append @IDENTITY pr=pr(-1)*(1+prodg)
wgm.append @IDENTITY WB=N*W
wgm.append sek=sk/2+sk(-1)*(1+(prodg+ran1))/2
wgm.append @IDENTITY ist=xi*sek
wgm.append ista=ik(-1)+omega*(ist-ik(-1))
wgm.append @IDENTITY ik=ista-(sk-sek)
wgm.append @IDENTITY sk=consk+gk+dkk
wgm.append @IDENTITY k=k(-1)*(1+gik)
wgm.append gik=go+psi*u(-1)-theta1*rrle
wgm.append @IDENTITY u=yk/k(-1)
wgm.append @IDENTITY rrle=(1+rl)/(1+pi)-1
wgm.append @IDENTITY dkk= k-k(-1)+delta*k(-1)
wgm.append @IDENTITY Se=sek*p
wgm.append @IDENTITY FT=S-WB+I-I(-1)-rl(-1)*I(-1)
wgm.append @IDENTITY S=sk*p
wgm.append @IDENTITY I=ik*UC
wgm.append @IDENTITY FTe=Fue+Fd+rl(-1)*(Ld(-1)-I(-1))
wgm.append lamsta=(Fd+Fue+rl(-1)*(Ld(-1)-I(-1)))/((1-zeta)*sek*UC+(1+rl(-1))*zeta*sek*UC(-1))
wgm.append @IDENTITY zeta=ik(-1)/sek
wgm.append Fue=phi3*DK(-1)
wgm.append Fd=phi1*(FT(-1))
wgm.append @IDENTITY UC=WB/yk
wgm.append lam=lam(-1)+upsilon*(lamsta(-1)-lam(-1))
wgm.append p=(1+lam)*NHUC
wgm.append @IDENTITY pi=(p-p(-1))/p(-1)
wgm.append NHUC=(1-zeta1)*NUC+zeta1*(1+rlc(-1))*NUC(-1)
wgm.append @IDENTITY NUC=W/pr
wgm.append @IDENTITY Ld=Ld(-1)+DK+I-I(-1)-Fu-(es-es(-1))*pe
wgm.append @IDENTITY DK=dkk*p 
wgm.append @IDENTITY Fu=FT-Fd-rl(-1)*(Ld(-1)-I(-1))
wgm.append es=es(-1)+((1-phi3)*DK(-1))/pe
wgm.append wsta=exp(ro0+ro1*log(pr)+ro2*log(N(-1)/Nt(-1)))
wgm.append W=W(-1)+ro3*(wsta*p(-1)-W(-1))
wgm.append @IDENTITY Yp=WB+Fd+Fb+rm(-1)*Md(-1)+rb(-1)*Bpd(-1)+Bls(-1)
wgm.append @IDENTITY Ypd=Yp*(1-theta)
wgm.append @IDENTITY T=theta*Yp 
wgm.append @IDENTITY ypdk=Ypd/p
wgm.append @IDENTITY vk=V/p
wgm.append @IDENTITY tk=T/p
wgm.append @IDENTITY consk=a1*ypde+a2*vk(-1) 
wgm.append @IDENTITY ypde=ypdk/2+(ypdk(-1)*(1+prodg))/2
wgm.append @IDENTITY V=V(-1)+Yp-T-CONS+(pe-pe(-1))*ed(-1)+(pb-pb(-1))*Bls(-1)
wgm.append @IDENTITY CONS=consk*p
wgm.append @IDENTITY Vn=V-Hpd
wgm.append @IDENTITY Hpd=a4*CONS
wgm.append @IDENTITY Fb=rl(-1)*Ls(-1)+rb(-1)*Bbd(-1)-rm(-1)*Ms(-1)
wgm.append @IDENTITY Md=Vn-Bpd-pe*ed-pb*Bld
wgm.append @IDENTITY Bs=Bs(-1)+G-T-(Bls-Bls(-1))*pb+rb(-1)*(Bps(-1)+Bbs(-1))+Bls(-1)
wgm.append @IDENTITY G=gk*p
wgm.append @IDENTITY Bbd=Ms-Ls-Hbd
wgm.append @IDENTITY Bbs=Bbs(-1)+(Bs-Bs(-1))-(Bps-Bps(-1))-(Bcbs-Bcbs(-1))
wgm.append @IDENTITY Hbd=ro*Ms
wgm.append @IDENTITY LR=Bbd/Ms
wgm.append @IDENTITY EX=LR(-1)>B1
wgm.append @IDENTITY ZED=LR(-1)<=B2
wgm.append @IDENTITY drm=ZED*alfa1-EX*alfa1
wgm.append @IDENTITY rm=rm(-1)+drm
wgm.append @IDENTITY Fbsta=mu1*Y
wgm.append @IDENTITY rl=rm+mu2
wgm.append @IDENTITY mu2=(1/Ls(-1))*(Fbsta-rb(-1)*Bbd(-1)+rm*(Ms(-1)-Ls(-1)))
wgm.append @IDENTITY rbl=rb+add3
wgm.append @IDENTITY H=Hbs+Hps
wgm.append @IDENTITY Bcbd=H
wgm.append @IDENTITY Bcbs=Bcbd
wgm.append @IDENTITY Ls=Ld
wgm.append @IDENTITY Ms=Md
wgm.append @IDENTITY Hps=Hpd
wgm.append @IDENTITY Bps=Bpd
wgm.append @IDENTITY Bls=Bld
wgm.append @IDENTITY Hbs=Hbd
wgm.append @IDENTITY Y=sk*p+(ik-ik(-1))*UC