About the IFPRI Standard CGE model

GAMS
1

I have calibrated the IFPRI dynamic CGE model and the model runs successfully. I want to include the timber dynamics into the model but timber harvest is not coupling into the model. The SAM of the model does not include any harvest entry. I already modified CES function like that:

QA(A) =E= alphaa(A) * (
deltaa1(A)*HARVEST(A)(-rhoa(A))
+ deltaa2(A)*QVA(A)(-rhoa(A))
+ (1 - deltaa1(A) - deltaa2(A))*QINTA(A)**(-rhoa(A))
) ** (-1 / rhoa(A));

How can I include harvest dynamic from the dynsim.dat file into the model so that harvest dynamic impacts dynamically into the model. Harvest dynamic included in the dynsim.dat file which is the supporting file of the dynsim101.gms file.
Harvest dynamics in the dynsim.dat file included like bellow:

  • TIMBER GROWING STOCK MODEL WITH MORTALITY SHOCKS — GENERALIZED FOR SIM

PARAMETER
GSTOCK(SIM,YR) growing stock of timber
GROWTH(SIM,YR) timber growth
HARVDYN(SIM,YR) timber harvest (70% of growth)
MORTALITY(SIM,YR) timber loss due to natural mortality (shocked)
DRAINAGE(SIM,YR) total timber drainage (harvest + mortality)
GSTOCKDUM(SIM) dummy for recursion
GSTOCK0(SIM) initial stock value
HARVESTSIM(A,SIM,YR) Harvest simulation parameter
;

  • — Initial Timber Stock —
    GSTOCK0(SIM) = 3000;
    GSTOCK(SIM,YR1) = GSTOCK0(SIM);
    GSTOCKDUM(SIM) = GSTOCK0(SIM);

  • — Mortality shocks from 2020 (YR1) to 2050 (YR31) —
    MORTALITY(SIM, YR1) = 0.005;
    MORTALITY(SIM, YR2) = 0.005;
    MORTALITY(SIM, YR3) = 0.005;
    MORTALITY(SIM, YR4) = 0.005;
    MORTALITY(SIM, YR5) = 0.005;
    MORTALITY(SIM, YR6) = 0.005;
    MORTALITY(SIM, YR7) = 0.005;
    MORTALITY(SIM, YR8) = 0.005;
    MORTALITY(SIM, YR9) = 0.005;
    MORTALITY(SIM, YR10) = 0.005;
    MORTALITY(SIM, YR11) = 0.0075;
    MORTALITY(SIM, YR12) = 0.005;
    MORTALITY(SIM, YR13) = 0.005;
    MORTALITY(SIM, YR14) = 0.005;
    MORTALITY(SIM, YR15) = 0.005;
    MORTALITY(SIM, YR16) = 0.005;
    MORTALITY(SIM, YR17) = 0.005;
    MORTALITY(SIM, YR18) = 0.005;
    MORTALITY(SIM, YR19) = 0.005;
    MORTALITY(SIM, YR20) = 0.005;
    MORTALITY(SIM, YR21) = 0.01;
    MORTALITY(SIM, YR22) = 0.005;
    MORTALITY(SIM, YR23) = 0.005;
    MORTALITY(SIM, YR24) = 0.005;
    MORTALITY(SIM, YR25) = 0.005;
    MORTALITY(SIM, YR26) = 0.005;
    MORTALITY(SIM, YR27) = 0.005;
    MORTALITY(SIM, YR28) = 0.005;
    MORTALITY(SIM, YR29) = 0.005;
    MORTALITY(SIM, YR30) = 0.005;
    MORTALITY(SIM, YR31) = 0.005;

  • — Recursive Timber Stock Loop from 2020–2050 —
    LOOP((SIM,YR),
    GROWTH(SIM, YR)(NOT YR1(YR)) = 0.04 * GSTOCKDUM(SIM); HARVDYN(SIM,YR)(NOT YR1(YR)) = 0.70 * GROWTH(SIM, YR);
    DRAINAGE(SIM, YR)$(NOT YR1(YR)) = HARVDYN(SIM,YR) + MORTALITY(SIM, YR) * GSTOCKDUM(SIM);

    GSTOCK(SIM, YR)YR1(YR) = GSTOCK0(SIM); GSTOCK(SIM, YR)(NOT YR1(YR)) = GSTOCKDUM(SIM) + GROWTH(SIM, YR) - DRAINAGE(SIM, YR);
    GSTOCKDUM(SIM) = GSTOCK(SIM, YR);
    HARVESTSIM(‘AFOR’,SIM,YR) = HARVDYN(SIM,YR);

);

DISPLAY GSTOCK, GROWTH, HARVDYN, MORTALITY, DRAINAGE;

Give me a simple way so that i can coupled the harvest dynamic into the main model dynamically. It will really be helpful for me because I have been trying since 3 weeks but failed.
Advance thanks for your suggestions.