-
|
Dear MAgPIE community, I'm trying to extract pasture yield and plot it on a global map with Panoply. Could you please help me check if my scripts below are right? Also, could you please tell me which variable I should write for Many thanks! past_pro<- collapseNames(readGDX(gdx="fulldata.gdx","oq14_yield_past")[,,"level"]) grid_pastproduc <- gdxAggregate(gdx="fulldata.gdx",x = past_pro, to="grid",weight=pasture, absolute=TRUE,dir=".") write.magpie(grid_pastproduc,"grid_pastproduc.nc") |
Beta Was this translation helpful? Give feedback.
Replies: 1 comment 1 reply
-
|
Reading Here, the weight argument has to be a 0.5 degree spatial scale magpie object most likely. Other MAgPIE members might have a good suggestion for that. Note that pasture yields are relative measures, so the choice of weight is critical in interpreting the yield results on a 0.5 degree scale. |
Beta Was this translation helpful? Give feedback.
-
|
Many thanks Abhimishr for your suggestion! I will have another try and keep that in mind. |
Beta Was this translation helpful? Give feedback.
Reading
oq14_yield_pastis not the right way to read pasture yields as you are reading the values of the equation.What you need in this case is
vm_yldvariable forpasture.This should be read in as
ov_yldand then you can filter forpasture.Here, the weight argument has to be a 0.5 degree spatial scale magpie object most likely. Other MAgPIE members might have a good suggestion for that. Note that pasture yields are relative measures, so the choice of weight is critical in interpreting the yield results on a 0.5 degree scale.