Bugfix in grasslands_apr22 pasture realization.

CHANGELOG.md

@@ -12,6 +12,8 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/).
 - **14_yields_and_config** The new default is to not use yield calibration factors from a calibration run. The switch s14_use_yield_calib can optionally reenable the use of yield calibration factors.
 - **scripts** LUH2_disaggregation output script was modified. Specifically, flooded area was made compatible with the LUH definition, cropland and grazing land were added to the states.nc file, and specific naming/details (datatype,  zname, xname, and yname) were added when creating the .nc files.
 - **scripts** For the emulator scripts select a different bioenergy demand variable that excludes bioenergy sources other than second generation bioenergy crops. Set the minimal bioenergy demand to zero. Both avoid artificial clustering of data points and allow for better fits.
+- **default.cfg** Removed description of cfg$gms$c31_past_suit_scen since no longer needed due to changes in 31_past described below. Its function is now done by cfg$gms$c31_grassl_yld_scenario.
+- **31_past** in grasslands_apr22 realization: changed structure of f31_pastr_suitability to align with ssp-rcp specific input data formulation. Changed input filename from cs3 to cs2. Added `cc`, `nocc` and `nocc_hist` options for `i31_manpast_suit` and `i31_manpast_suit` and changed input gams code from table to parameter. Climate scenario assignment moved from preloop.gms to input.gms. Removed pastSuit set in sets.gms as not needed anymore. Adjusted not_used.txt in both grasslands_apr22 and static realizations.
 - **36_employment** regression between hourly labor regression and GDP pc changed from linear to log-log
 
 ### added

CITATION.cff

@@ -112,7 +112,13 @@ authors:
     affiliation: "Potsdam Institute for Climate Impact Research"
     email: [email protected]
 
-  - family-names: Lotze-Campen
+  - family-names: Köberle
+    given-names: Alexandre
+    orcid: https://orcid.org/0000-0003-0328-4750
+    affiliation: "Potsdam Institute for Climate Impact Research"
+    email: [email protected]
+
+- family-names: Lotze-Campen
     given-names: Hermann
     orcid: https://orcid.org/0000-0002-0003-5508
     affiliation: "Potsdam Institute for Climate Impact Research"

config/default.cfg

@@ -718,19 +718,12 @@ cfg$gms$s31_fac_req_past  <- 1             # def = 1
 # * reflecting that not all pasture area is grazed.
 cfg$gms$s31_unequal  <- 1             # def = 1
 
-# * Scenario for grassland and pasture yields (only take effect for realization grasslands_apr22)
+# * Scenario for grassland, pasture suitability areas and pasture yields (only take effect for realization grasslands_apr22)
 # *   options:   cc        (climate change)
 # *             nocc       (no climate change)
 # *             nocc_hist  (no climate change after year defined by sm_fix_cc)
 cfg$gms$c31_grassl_yld_scenario  <- "cc"   # def = "cc"
 
-
-# * SSP scenario switch for pasture suitability areas (only take effect for realization grasslands_apr22)
-# * options: ssp126, ssp245, ssp370, ssp460, ssp585
-# *          nocc  (fixed to values from 1995),
-# *          nocc_hist (ssp245 until year defined in sm_fix_cc, fixed thereafter)
-cfg$gms$c31_past_suit_scen  <- "ssp370"   # def = "ssp370"
-
 # * Switch that allows selecting how yield calibration factors will be calculated.
 # * If 0, grassland yield calibration will be calculated as multiplicative relative values.
 # * If 1, calibration values will be limited to additive absolute values when modeled yields strongly

modules/31_past/endo_jun13/not_used.txt

@@ -2,5 +2,4 @@ name,type,reason
 pcm_land, input, not needed
 vm_tau,input,questionnaire
 fm_pastr_tau_hist,input,questionnaire
-sm_fix_SSP2, input, not needed
 sm_fix_cc, input, not needed

modules/31_past/grasslands_apr22/input.gms

@@ -11,26 +11,26 @@ scalars
   s31_cost_grass_prod  Grasslands factor costs (USD05MER per tDM)                 / 1 /
 ;
 
-$setglobal c31_past_suit_scen  ssp370
+$setglobal c31_grassl_yld_scenario  cc
+*   options:  cc        (climate change)
+*             nocc      (no climate change)
+*             nocc_hist (no climate change after year defined by sm_fix_cc)
 
-table f31_pastr_suitability(t_all,j,ssp_past) Areas suitable for pasture management (mio. ha)
+parameter f31_pastr_suitability(t_all,j) Areas suitable for pasture management (mio. ha)
+/
 $ondelim
-$include "./modules/31_past/input/f31_pastr_suitability.cs3"
+$include "./modules/31_past/input/f31_pastr_suitability.cs2"
 $offdelim
-;
+/;
+
+i31_manpast_suit(t_all,j) = f31_pastr_suitability(t_all,j)
+* set values to 1995 if nocc scenario is used, or to sm_fix_cc after sm_fix_cc if nocc_hist is used
+$if "%c31_grassl_yld_scenario%" == "nocc" i31_manpast_suit(t_all,j) = f31_pastr_suitability("y1995",j);
+$if "%c31_grassl_yld_scenario%" == "nocc_hist" i31_manpast_suit(t_all,j)$(m_year(t_all) > sm_fix_cc) = i31_manpast_suit(t_all,j)$(m_year(t_all) = sm_fix_cc);
 
-table f31_LUH2v2(t_all,j,f31_luh) Different land type areas (mio. ha)
-$ondelim
-$include "./modules/31_past/input/f31_LUH2v2.cs3"
-$offdelim
-;
 
 scalar s31_limit_calib   Relative managament calibration switch (1=limited 0=pure relative) / 1 /;
 
-$setglobal c31_grassl_yld_scenario  cc
-*   options:  cc        (climate change)
-*             nocc      (no climate change)
-*             nocc_hist (no climate change after year defined by sm_fix_cc)
 
 table f31_grassl_yld(t_all,j,grassland,w) LPJmL potential yields per cell (rainfed only) (tDM per ha)
 $ondelim

modules/31_past/grasslands_apr22/preloop.gms

@@ -5,15 +5,6 @@
 *** |  MAgPIE License Exception, version 1.0 (see LICENSE file).
 *** |  Contact: [email protected]
 
-$ifthen "%c31_past_suit_scen%" == "nocc" 
-  i31_manpast_suit(t_all,j) = f31_pastr_suitability("y1995",j,"ssp245");
-$elseif "%c31_past_suit_scen%" == "nocc_hist" 
-  i31_manpast_suit(t_all,j) = f31_pastr_suitability(t_all,j,"ssp245");
-  i31_manpast_suit(t_all,j)$(m_year(t_all) > sm_fix_cc) = f31_pastr_suitability(t_all,j,"ssp245")$(m_year(t_all) = sm_fix_cc);
-$else 
-  i31_manpast_suit(t_all,j) = f31_pastr_suitability(t_all,j,"%c31_past_suit_scen%");
-  i31_manpast_suit(t_all,j)$(m_year(t_all) <= sm_fix_SSP2) = f31_pastr_suitability(t_all,j,"ssp245");
-$endif
 
 pc31_grass(j,grassland) = f31_LUH2v2("y1995",j,grassland);
 

modules/31_past/grasslands_apr22/sets.gms

@@ -19,7 +19,4 @@ sets
  grass_from31(grassland) pasture management options
    / pastr,range /
 
- ssp_past SSP scenarios for pasture suitability areas
-   / ssp126, ssp245, ssp370, ssp460, ssp585 /
-
 ;

modules/31_past/input/files

@@ -1,5 +1,5 @@
 * list of files that are required here
-f31_pastr_suitability.cs3
+f31_pastr_suitability.cs2
 f31_LUH2v2.cs3
 f31_grassl_yld.cs3
 f31_grass_bio_hist.cs3

modules/31_past/static/not_used.txt

@@ -4,5 +4,4 @@ vm_yld, input, not needed
 pm_land_conservation,input,questionnaire
 vm_tau,input,questionnaire
 fm_pastr_tau_hist,input,questionnaire
-sm_fix_SSP2, input, not needed
 sm_fix_cc, input, not needed