Consider the Armington aggregation of domestic and imported commodities in the GTAP-EG CGE model
(Rutherford, T. and Paltsev, S. [2000]: GTAPinGAMS and GTAP-EG: Global Datasets for Economic Research and Illustrative Models. Available online: https://www.mpsge.org/papers/gtaptext.pdf):
$prod:a(i,r)$a0(i,r) s:4 m:8 s.tl(m):0
o:pa(i,r) q:a0(i,r)
i:py(i,r) q:d0(i,r)
i:py(i,s) q:vxmd(i,s,r) p:pmx0(i,s,r) s.tl:
+ a:ra(s) t:tx(i,s,r) a:ra(r) t:(tm(i,s,r)*(1+tx(i,s,r)))
i:pt#(s) q:vtwr(i,s,r) p:pmt0(i,s,r) s.tl: a:ra(r) t:tm(i,s,r)
Is there a way to compute a reporting variable for i:pt#(s) with indices (i,s,r) directly in MPSGE?
As far as I understand, one can only get a Reporting Variable for this input with indices (i,r):
V:VA_PT(i,r) i:PT prod:A(i,r)
In order to get the values of the reporting variable with indices (i,s,r), one would have to refer to the reporting variable for i:py(i,s):
V:VA_PYxm(i,s,r) i:PY(i,s) prod:A(i,r)
And then multiply it with vtwr(i,s,r)/vxmd(i,s,r).