Fixes in rc after test runs for forestry and default settings

CHANGELOG.md

@@ -13,12 +13,13 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 - **38_factor_costs** Added the new "sticky" realization to the factor costs module. The realization "sticky_feb18" favors expansion in cells with preexisting farmland and capital based on capital investment decisions.
 
  - **modules** added endogenous implementation of local biophysical (bph) impacts of afforestation to existing realizations in modules 32_forestry (dynamic_oct19) and 56_ghg_policy (price_jan20). default = off
- - **73_timber** Added timber module which brings the ability of producing woody biomass for timber plantations and natural vegetation. Default = off.
- - **32_forestry** New realization for timber production from timber plantations. This builds up on previous forestry realization for afforestation.
- - **35_natveg** New realization for timber production from natural vegetation.
+ - **73_timber** Added timber module which brings the ability of producing woody biomass for timber plantations and natural vegetation. Default = off. New switch: `s73_foresight`. New scalars : `s73_timber_prod_cost`, `s73_timber_harvest_cost`,`s73_cost_multiplier`,`s73_free_prod_cost`
+ - **32_forestry** New realization `dynamic_may20` for forestry land use dynamics. This builds up on previous forestry realization for afforestation. New switches: `s32_fix_plant`, `c32_interest_rate`. New scalars : `s32_plant_share`, `s32_forestry_int_rate`, `s32_investment_cost`, `s52_plantation_threshold`.
+ - **35_natveg** New realization `dynamic_may20` for natural vegetation land use dynamics. New forest protection scenario.
+ - **52_carbon** Added interface which is used for calculating additional investment needed in plantations when carbon stocks are lower than a specified threshold. New scalar: `s52_plantation_threshold`.
  - **57_maccs** Added MACCs from Harmsen PBL 2019
  - **15_food** Added the option to include calories from alcohol consumption in healthy and sustainable EAT-Lancet diets.
- - **scripts** added start script for making timber production runs (forestry.R)
+ - **scripts** added start script for making timber production runs (forestry.R).
 
 ### changed
  - **scripts** updated selection routine for start and output scripts
@@ -31,6 +32,9 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
  - **scripts** removed GP_final start script due to obsolescence
 
 ### fixed
+ - **32_forestry** Rotation length calculation based on correct marginals of growth function in timber plantations. Clearer calculations for harvested area for timber production.
+ - **35_natveg** Clearer calculations for harvested area for timber production.
+ - **52_carbon** Fix to the Carbon densities received from LPJmL for timber plantations.
 
 ## [4.2.1] - 2020-05-15
 

config/default.cfg

@@ -24,10 +24,10 @@ cfg$model <- "main.gms"   #def = "main.gms"
 # which input data sets should be used?
 
 cfg$input <- c("isimip_rcp-IPSL_CM5A_LR-rcp2p6-co2_rev48_c200_690d3718e151be1b450b394c1064b1c5.tgz",
-         "rev4.51_h12_magpie.tgz",
-         "rev4.51_h12_validation.tgz",
- 	 "calibration_H12_c200_26Feb20.tgz",
-         "additional_data_rev3.85.tgz")
+         "rev4.52_h12_magpie.tgz",
+         "rev4.52_h12_validation.tgz",
+         "calibration_H12_c200_26Feb20.tgz",
+         "additional_data_rev3.86.tgz")
 
 #a list of repositories (please pay attention to the list format!) in which the
 #files should be searched for. Files will be searched in all repositories until
@@ -314,8 +314,8 @@ cfg$gms$c15_EAT_scen <- "FLX"                # def = FLX
 # * Scenario target for the inclusion of alcohol in the EAT-Lancet diet
 # * only activated if s15_exo_diet is set to 1
 # * (0):      no alcohol consumption, as in the original version of the EAT-Lancet diet
-# * (0.014):  maximum target for alcohol consumption is 1.4% of total calorie consumption 
-# *           (see Lassen et al., 2020) 
+# * (0.014):  maximum target for alcohol consumption is 1.4% of total calorie consumption
+# *           (see Lassen et al., 2020)
 cfg$gms$s15_alc_scen <- 0.014                # def = 0.014
 
 # * Switch and specification of countries for which exogenous food scenarios
@@ -456,6 +456,33 @@ cfg$gms$c32_tcre_ctrl <- "ann_TCREmean"
 # * 0=off
 cfg$gms$s32_initial_distribution         <- 0     # def = 0
 
+# Switch to fix timber plantations area after sm_fix_SSP2 value
+# This makes sure that cellular timber plantations are not changed after 2020.
+# This switch does not impact the model behavior before sm_fix_SSP2 as plantations are
+# not fixed for that period in the dynamic realization.
+# * 1= Plantation area is fixed after 2020
+# * 0= Plantation area is not fixed after 2020
+cfg$gms$s32_fix_plant <- 0                          # def = 0
+
+# Percentage of overall timber production which guides plantations establishment
+# decision. The model sees a proportion of wood demand from future and establishes
+# new plantations accordingly in the current time step (accepted values 0 to 1).
+cfg$gms$s32_plant_share <- 0.25                     # def = 0.25
+
+# Switch fore regional or global interest rate for rotation length calculations.
+# Using the global setting would mean that the timber plantation decisions are
+# decoupled from other decisions in the model which are based on regional interest
+# rates.
+# * ("regional")  = Regionally differentiated interest rates
+# * ("global")    = One global interest rate
+cfg$gms$c32_interest_rate <- "regional"             # def = "regional"
+
+# Global interest rate for plantations in case c32_interest_rate switch is global
+# Accepted values between 0 and 1.
+cfg$gms$s32_forestry_int_rate   <- 0.05              # def = 0.05
+
+# Additional investment cost in plantations for using unproductive cells
+cfg$gms$s32_investment_cost     <- 200              # def = 200
 
 # ***---------------------    34_urban   ---------------------------------------
 # * (static): static urban
@@ -644,6 +671,9 @@ cfg$gms$carbon <- "normal_dec17"                       # def = normal_dec17
 # *             nocc (no climate change)
 cfg$gms$c52_carbon_scenario  <- "nocc"   # def = "nocc"
 
+# * Minimum threshold of carbon density (tC/ha) in timber plantations
+cfg$gms$s52_plantation_threshold <- 8   # def = 8
+
 # ***---------------------    53_methane   -------------------------------------
 # * (ipcc2006_flexreg_apr16): IPCC 1996 methodology
 # * (off):                    methane calculations deactivated
@@ -925,6 +955,32 @@ cfg$gms$timber <- "biomass_mar20"                  # def = biomass_mar20
 # * 0=off
 cfg$gms$s73_demand_switch         <- 0     # def = 0
 
+# Setting to define if the model should be forward looking or not in terms of
+# seeing the future demand for current timestep establishment of new plantations
+# in the forward setting, the model sees only current demand for establishment
+# decisions in the historical time period but then sees future demand for non
+# historic time periods. In myopic setting, the model sees always the current
+# timber demand for establishment decisions.
+# * 1 = forward looking. Model sees future demand for establishment in current step
+# * 0 = myopic. Model sees current demand for establishment in current step
+cfg$gms$s73_foresight <- 0                 # def = 0
+
+# harvesting cost per ton of dry matter produced (USD/tDM)
+s73_timber_prod_cost <- 2000            # def = 2000
+
+# harvesting cost per ha of forests (USD/ha)
+s73_timber_harvest_cost <- 2000         # def = 2000
+
+# Cost multiplier for harvesting costs to make natural vegetation harvest expensive
+# than timber plantation harvst. This provides a signal to the model to harvest
+# timber plantations first.
+s73_cost_multiplier <- 1.5              # def = 1.5
+
+# Cost of production without using any land in case the model is running into infeasibilities.
+# This is a last ditch effort for the model and the variable assocaited with this cost
+# should not be used in a normally feasible model run (USD/tDM)
+s73_free_prod_cost <- 1000000           # def = 1000000
+
 
 # ***-------------------    80_optimization ------------------------------------
 # * (nlp_apr17):    solve procedure solving the whole, nonlinear problem
@@ -1001,7 +1057,7 @@ cfg$files2export$start <- c("input/info.txt",
                             ".Rprofile",
                             "input/*.spam",
                             "input/lpj_envflow_total_c200.mz",
-                            "input/lpj_watavail_total_c200.mz",                            
+                            "input/lpj_watavail_total_c200.mz",
                             "input/validation.mif",
                             "calib_*.cs3",
                             "input/spamplot_*.pdf")

core/sets.gms

@@ -253,7 +253,8 @@ sets
 ***Forestry**
    ac Age classes  / ac0,ac5,ac10,ac15,ac20,ac25,ac30,ac35,ac40,ac45,ac50,
                     ac55,ac60,ac65,ac70,ac75,ac80,ac85,ac90,ac95,ac100,
-                    ac105,ac110,ac115,ac120,ac125,ac130,ac135,ac140,ac145,acx /
+                    ac105,ac110,ac115,ac120,ac125,ac130,ac135,ac140,ac145,
+                    ac150,ac155,acx /
 
   ac_est(ac) Dynamic subset of age classes for establishment
 

main.gms

@@ -145,60 +145,60 @@ $title magpie
 *'  * Always try to access model outputs through the corresponding magpie package instead of accessing them directly with readGDX. It cannot be guaranteed that your script will work in the future if you do otherwise (as only the corresponding magpie package will be continuously adapted to changes in the GAMS code).
 
 *##################### R SECTION START (VERSION INFO) ##########################
-* 
-* Used data set: isimip_rcp-IPSL_CM5A_LR-rcp2p6-co2_rev48_c200_690d3718e151be1b450b394c1064b1c5.tgz
-* md5sum: b07787909a9339225fd097c813277b19
-* Repository: /p/projects/landuse/data/input/archive
-* 
-* Used data set: rev4.51_h12_magpie.tgz
-* md5sum: 7d380cccd273baf44d169492b08cd98e
-* Repository: /p/projects/rd3mod/inputdata/output
-* 
-* Used data set: rev4.51_h12_validation.tgz
-* md5sum: 99e4aaca789892fc44896d4fcd0e0459
-* Repository: /p/projects/rd3mod/inputdata/output
-* 
+*
+* Used data set: isimip_rcp-IPSL_CM5A_LR-rcp2p6-co2_rev46_c200_690d3718e151be1b450b394c1064b1c5.tgz
+* md5sum: NA
+* Repository: https://rse.pik-potsdam.de/data/magpie/public
+*
+* Used data set: rev4.52_h12_magpie.tgz
+* md5sum: NA
+* Repository: scp://cluster.pik-potsdam.de/p/projects/rd3mod/inputdata/output
+*
+* Used data set: rev4.52_h12_validation.tgz
+* md5sum: NA
+* Repository: scp://cluster.pik-potsdam.de/p/projects/rd3mod/inputdata/output
+*
 * Used data set: calibration_H12_c200_26Feb20.tgz
-* md5sum: 0fc75955deb5916d4ac078e8d817feda
-* Repository: /p/projects/rd3mod/mirror/rse.pik-potsdam.de/data/magpie/intern
-* 
-* Used data set: additional_data_rev3.85.tgz
-* md5sum: 23bc2b7ff07da84b798d741c3d2ad3d2
-* Repository: /p/projects/rd3mod/mirror/rse.pik-potsdam.de/data/magpie/intern
-* 
+* md5sum: NA
+* Repository: https://rse.pik-potsdam.de/data/magpie/public
+*
+* Used data set: additional_data_rev3.86.tgz
+* md5sum: NA
+* Repository: scp://cluster.pik-potsdam.de/p/projects/landuse/data/input/archive
+*
 * Low resolution: c200
 * High resolution: 0.5
-* 
+*
 * Total number of cells: 200
-* 
+*
 * Number of cells per region:
 *   CAZ  CHA  EUR  IND  JPN  LAM  MEA  NEU  OAS  REF  SSA  USA
 *    28   24   10    7    3   53   17    8   22    7   11   10
-* 
+*
 * Regionscode: 690d3718e151be1b450b394c1064b1c5
-* 
-* Regions data revision: 4.51
-* 
+*
+* Regions data revision: 4.52
+*
 * lpj2magpie settings:
 * * LPJmL data folder: /p/projects/landuse/data/input/lpj_input/isimip_rcp/IPSL_CM5A_LR/rcp2p6/co2
-* * Additional input folder: /p/projects/landuse/data/input/other/rev48
-* * Revision: 48
+* * Additional input folder: /p/projects/landuse/data/input/other/rev46
+* * Revision: 46
 * * Call: lpj2magpie(input_folder = path(cfg$lpj_input_folder, gsub("-",     "/", cfg$input)), input2_folder = path(cfg$additional_input_folder,     paste("rev", floor(cfg$revision), sep = "")), output_file = lpj2magpie_file,     rev = cfg$revision)
-* 
+*
 * aggregation settings:
 * * Input resolution: 0.5
 * * Output resolution: c200
-* * Input file: /p/projects/landuse/data/input/archive/isimip_rcp-IPSL_CM5A_LR-rcp2p6-co2_rev48_0.5.tgz
-* * Output file: /p/projects/landuse/data/input/archive/isimip_rcp-IPSL_CM5A_LR-rcp2p6-co2_rev48_c200_690d3718e151be1b450b394c1064b1c5.tgz
+* * Input file: /p/projects/landuse/data/input/archive/isimip_rcp-IPSL_CM5A_LR-rcp2p6-co2_rev46_0.5.tgz
+* * Output file: /p/projects/landuse/data/input/archive/isimip_rcp-IPSL_CM5A_LR-rcp2p6-co2_rev46_c200_690d3718e151be1b450b394c1064b1c5.tgz
 * * Regionscode: 690d3718e151be1b450b394c1064b1c5
 * * (clustering) n-repeat: 5
 * * (clustering) n-redistribute: 0
 * * Call: aggregation(input_file = lpj2magpie_file, regionmapping = paste0("../",     cfg$regionmapping), output_file = aggregation_file, rev = cfg$revision,     res_high = cfg$high_res, res_low = cfg$low_res, hcells = cfg$highres_cells,     weight = cfg$cluster_weight, nrepeat = cfg$nrepeat, nredistribute = cfg$nredistribute,     sum_spam_file = cfg$spamfile, debug = FALSE, seed = cfg$seed)
-* 
-* 
-* 
-* Last modification (input data): Thu Sep  3 10:35:15 2020
-* 
+*
+*
+*
+* Last modification (input data): Thu Sep 10 09:59:16 2020
+*
 *###################### R SECTION END (VERSION INFO) ###########################
 
 $offupper

modules/14_yields/dynamic_aug18/presolve.gms

@@ -43,7 +43,7 @@ p14_growing_stock(t,j,ac,land_natveg,"natveg") =
 p14_growing_stock(t,j,ac,land_natveg,"natveg") = p14_growing_stock(t,j,ac,land_natveg,"natveg")$(p14_growing_stock(t,j,ac,land_natveg,"natveg")>0)+0.0001$(p14_growing_stock(t,j,ac,land_natveg,"natveg")=0);
 p14_growing_stock(t,j,ac,"forestry","plantations") = p14_growing_stock(t,j,ac,"forestry","plantations")$(p14_growing_stock(t,j,ac,"forestry","plantations")>0)+0.0001$(p14_growing_stock(t,j,ac,"forestry","plantations")=0);
 
-** Used in equations -- Annual value hence division by timestep
+** Used in equations
 ***************************************************************
 ** If the plantation yield switch is on, forestry yields are treated are plantation yields
 pm_timber_yield(t,j,ac,"forestry")$(s14_timber_plantation_yield = 1) = p14_growing_stock(t,j,ac,"forestry","plantations") ;

modules/20_processing/off/sets.gms

@@ -16,5 +16,5 @@ kpr(kall) Products that can be processed
 knpr(kall) Plant-based products that cannot be processed
    /alcohol,distillers_grain,ethanol,fibres,fish,livst_chick,livst_egg,livst_milk,
    livst_pig,livst_rum,oils,oilcakes,pasture,puls_pro,res_cereals,res_fibrous,
-   res_nonfibrous,scp/
+   res_nonfibrous,scp,wood,woodfuel/
 ;

modules/20_processing/substitution_dec18/sets.gms

@@ -16,7 +16,7 @@ kpr(kall) Products that can be processed
 knpr(kall) Products that cannot be processed
    /alcohol,distillers_grain,ethanol,fibres,fish,livst_chick,livst_egg,livst_milk,
    livst_pig,livst_rum,oilcakes,pasture,puls_pro,res_cereals,res_fibrous,
-   res_nonfibrous,scp/
+   res_nonfibrous,scp,wood,woodfuel/
 
 processing_subst20 Processing activities
 / milling, refining, extracting, distilling, fermentation, breeding, ginning, substitutes /

modules/21_trade/selfsuff_reduced/input.gms

@@ -5,7 +5,7 @@
 *** |  MAgPIE License Exception, version 1.0 (see LICENSE file).
 *** |  Contact: magpie@pik-potsdam.de
 
-$setglobal c21_trade_liberalization  l908080r807070
+$setglobal c21_trade_liberalization  l909090r808080
 *   options are "regionalized" and "globalized" and "fragmented"
 
 scalars

modules/21_trade/selfsuff_reduced/preloop.gms

@@ -19,3 +19,5 @@ elseif (s21_trade_tariff=0),
 
 pm_selfsuff_ext(t_ext,i,kforestry) = f21_self_suff("y2150",i,kforestry);
 pm_selfsuff_ext(t_all,i,kforestry) = f21_self_suff(t_all,i,kforestry);
+*Implausible jumps in selfsuff for wood in some regions.
+pm_selfsuff_ext(tstart21,i,kforestry) = f21_self_suff("y2010",i,kforestry);