Missing rotation intensity fix.

GAMS/IntExtOpt.gms

@@ -381,6 +381,7 @@ $gdxin
  importAmount.L(all_types) = previousImportAmount(all_types);
  exportAmount.L(all_types) = previousExportAmount(all_types);
  rotationIntensity.L(location) = previousRotationIntensity(location);
+ rotationIntensity.L(location) $ (previousRotationIntensity(location) = 0) = 0.02;
  woodYieldRota.L(location) = woodYieldMax(location) * (1 - exp(-woodYieldParam(location) / rotationIntensity.L(location))) * rotationIntensity.L(location);
  woodSupply.L(wood_type) = demand(wood_type);
 
@@ -448,8 +449,8 @@ $gdxin
  parameter rotationPeriod(location);
  parameter irrigWaterWithdrawn(crop, location);
  parameter fertiliserUsed(crop, location);
-* totalLandCoverChange(land_cover, land_cover_after) = sum(location, landCoverChange.L(land_cover, land_cover_after, location));
-* landCoverDiff(land_cover) = sum(location, landCoverArea.L(land_cover, location)) - sum(location, previousLandCoverArea(land_cover, location));
+ totalLandCoverChange(land_cover, land_cover_after) = sum((protection, location), landCoverChange.L(land_cover, land_cover_after, protection, location));
+ landCoverDiff(land_cover) = sum((protection, location), landCoverArea.L(land_cover, protection, location)) - sum((protection, location), previousLandCoverArea(land_cover, protection, location));
  rotationPeriod(location) = 1 / rotationIntensity.L(location);
  irrigWaterWithdrawn(crop, location) = sum(farming_type, irrigMaxRate(crop, location) * irrigI.L(crop, farming_type, location) * cropArea.L(crop, farming_type, location) * 0.01);
  fertiliserUsed(crop, location) = sum(farming_type, fertI.L(crop, farming_type, location) * cropArea.L(crop, farming_type, location));