Merge pull request #6 from FelicitasBeier/regional_scen

Regional scen
magpiemodel · Mar 28, 2020 · 755e14f · 755e14f
2 parents 7a0a28c + 140bf0e
commit 755e14f
Show file tree

Hide file tree

Showing 33 changed files with 927 additions and 208 deletions.
diff --git a/config/default.cfg b/config/default.cfg
@@ -96,6 +96,34 @@ cfg$gms$s_use_gdx <- 2  # def = 2
 # ***                           MODULES
 # ***---------------------------------------------------------------------------
 
+#### Useful shortcuts ####
+# Vector of all iso countries (used for regional policy implementations)
+all_iso_countries <- "ABW,AFG,AGO,AIA,ALA,ALB,AND,ARE,ARG,ARM,
+                          ASM,ATA,ATF,ATG,AUS,AUT,AZE,BDI,BEL,BEN,
+                          BES,BFA,BGD,BGR,BHR,BHS,BIH,BLM,BLR,BLZ,
+                          BMU,BOL,BRA,BRB,BRN,BTN,BVT,BWA,CAF,CAN,
+                          CCK,CHN,CHE,CHL,CIV,CMR,COD,COG,COK,COL,
+                          COM,CPV,CRI,CUB,CUW,CXR,CYM,CYP,CZE,DEU,
+                          DJI,DMA,DNK,DOM,DZA,ECU,EGY,ERI,ESH,ESP,
+                          EST,ETH,FIN,FJI,FLK,FRA,FRO,FSM,GAB,GBR,
+                          GEO,GGY,GHA,GIB,GIN,GLP,GMB,GNB,GNQ,GRC,
+                          GRD,GRL,GTM,GUF,GUM,GUY,HKG,HMD,HND,HRV,
+                          HTI,HUN,IDN,IMN,IND,IOT,IRL,IRN,IRQ,ISL,
+                          ISR,ITA,JAM,JEY,JOR,JPN,KAZ,KEN,KGZ,KHM,
+                          KIR,KNA,KOR,KWT,LAO,LBN,LBR,LBY,LCA,LIE,
+                          LKA,LSO,LTU,LUX,LVA,MAC,MAF,MAR,MCO,MDA,
+                          MDG,MDV,MEX,MHL,MKD,MLI,MLT,MMR,MNE,MNG,
+                          MNP,MOZ,MRT,MSR,MTQ,MUS,MWI,MYS,MYT,NAM,
+                          NCL,NER,NFK,NGA,NIC,NIU,NLD,NOR,NPL,NRU,
+                          NZL,OMN,PAK,PAN,PCN,PER,PHL,PLW,PNG,POL,
+                          PRI,PRK,PRT,PRY,PSE,PYF,QAT,REU,ROU,RUS,
+                          RWA,SAU,SDN,SEN,SGP,SGS,SHN,SJM,SLB,SLE,
+                          SLV,SMR,SOM,SPM,SRB,SSD,STP,SUR,SVK,SVN,
+                          SWE,SWZ,SXM,SYC,SYR,TCA,TCD,TGO,THA,TJK,
+                          TKL,TKM,TLS,TON,TTO,TUN,TUR,TUV,TWN,TZA,
+                          UGA,UKR,UMI,URY,USA,UZB,VAT,VCT,VEN,VGB,
+                          VIR,VNM,VUT,WLF,WSM,YEM,ZAF,ZMB,ZWE"
+
 # ***---------------------    09_drivers   ----------------------------------------
 # * (aug17): default drivers
 cfg$gms$drivers <- "aug17"          # def = aug17
@@ -117,13 +145,29 @@ cfg$gms$costs <- "default"          # def = default
 # ***---------------------    12_interest_rate ---------------------------------
 # * (glo_jan16): global interest rates - annuity calculated with macro m_annuity
 # * (reg_feb18): regional interest rates, dependent on development state
-
 cfg$gms$interest_rate <- "reg_feb18"          # def = reg_feb18
 
-# * interest rate scenario: low, medium, high or mixed
+# * interest rate scenario: coupling, low, medium, high or mixed
 # * This only has effects on glo_jan16 realization, but for the coupling in both realizations
 cfg$gms$c12_interest_rate <- "medium"          # def = "medium"
 
+# * Regional interest rate policy scenario switch
+# * Specification of countries for which GDP-dependent interest rates apply.
+# * Options: list of iso-codes of countries where GDP-dependent interest rates
+# * should be in effect.
+# * In all other countries, the interest rate defined in s12_default_interest
+# * is in effect.
+# *
+# * Note: must be written in the format: "IND, BRA, DEU"
+# * Default: all iso countries
+# * Note: Does not take effect when c12_interest_rate="coupling"
+cfg$gms$gdp_countries12  <- all_iso_countries      # def = all_iso_countries
+
+# * Alternative interest rate that applies to the non-selected (!) countries
+# * Note: Does not take effect when c12_interest_rate="coupling"
+cfg$gms$s12_interest_noselect <- 0.02           # def = 0.02
+
+
 # ***---------------------    13_tc    -----------------------------------------
 # * (endo_jun18): endogenous technological change with full cost accounting and
 # *               stepwise updated crop area information
@@ -169,7 +213,8 @@ cfg$gms$s15_convergence <- 0.005            # def = 0.005
 # *  options: 0 (exogenous food demand) and 1 (endogenous food demand)
 cfg$gms$s15_elastic_demand <- 1            # def =1
 
-# * food scenario
+# * food scenario for selected (and respectively not selected) countries
+# * in scen_countries15
 # *   options:   SSP: "SSP1", "SSP2", "SSP3", "SSP4", "SSP5"
 # *             SRES: "A1", "A2", "B1", "B2"
 # *            OTHER: "PB" (planetary boundaries)
@@ -236,6 +281,14 @@ cfg$gms$c15_kcal_scen <- "healthy_BMI"       # def = healthy_BMI
 # *   options:   BMK, FLX, PSC, VEG, VGN, FLX_hmilk, FLX_hredmeat
 cfg$gms$c15_EAT_scen <- "FLX"                # def = FLX
 
+# * Switch and specification of countries for which exogenous food scenarios
+# * (calorie, diet, food waste scenarios) shall be applied
+# * Options: list of iso-codes of countries where exogenous scenario should be
+# * in effect
+# * Note: must be written in the format: "IND, BRA, DEU"
+# * Default: all iso countries
+cfg$gms$scen_countries15  <- all_iso_countries
+
 
 # ***---------------------    16_demand    -------------------------------------
 # * (sector_may15): default for flexible regions
@@ -423,6 +476,13 @@ cfg$gms$s42_irrigation_efficiency <- 0.66       # def=0.66
 # * (mixed): EFP policy only in hic regions
 cfg$gms$c42_env_flow_policy <- "off"             # def = "off"
 
+# * Switch and specification of countries for which environmental flow policy
+# * shall apply.
+# * Options: list of iso-codes of countries where EFP should be in effect.
+# * Note: must be written in the format: "IND, BRA, DEU"
+# * Default: all iso countries
+cfg$gms$EFP_countries  <- all_iso_countries # def = all_iso_countries
+
 # * Environmental flow protection scenario
 # * (0): do not consider environmental flows.
 # *      s42_env_flow_base_fraction and
@@ -465,14 +525,29 @@ cfg$gms$climate <- "static"    # def = static
 # * (exoeff_aug16): exogenous nr efficiency
 cfg$gms$nr_soil_budget <- "exoeff_aug16"    # def = exoeff_aug16
 
-# * scenario for nr efficiency on croplands or pastures
-# * options: constant,
+# * Scenario for nr efficiency on croplands or pastures for selected (and
+# * respectively non-selected) countries in cropneff_countries and pastneff_countries
+# * Options: constant,
 # *          neff55_55_starty1990,neff60_60_starty1990,neff65_70_starty1990,
 # *          neff65_70_starty2010,neff60_60_starty2010,neff55_60_starty2010,
 # *          neff70_75_starty2010,neff75_80_starty2010,neff80_85_starty2010
 # *          neff85_85_starty2010
-cfg$gms$c50_scen_neff <- "neff60_60_starty2010"   # def = neff60_60_starty2010
-cfg$gms$c50_scen_neff_pasture <- "constant"   	  # def = constant
+cfg$gms$c50_scen_neff_select <- "neff60_60_starty2010"   # def = neff60_60_starty2010
+cfg$gms$c50_scen_neff_noselect <- "neff60_60_starty2010"   # def = neff60_60_starty2010
+cfg$gms$c50_scen_neff_pasture_select <- "constant"   	  # def = constant
+cfg$gms$c50_scen_neff_pasture_noselect <- "constant"   	  # def = constant
+
+# * Switch and specification of countries for which above defined (c50_scen_neff_select,
+# * c50_scen_neff_pasture_select) apply, respectively.
+# * For all other countries c50_scen_neff_noselect / c50_scen_neff_pasture_noselect
+# * apply.
+# * Options: list of iso-codes of countries where selected neff scenario should
+# * be in effect.
+# * Note: must be written in the format: "IND, BRA, DEU"
+# * Default: above defined neff scenario applies to all iso countries
+cfg$gms$cropneff_countries  <- all_iso_countries
+cfg$gms$pastneff_countries  <- all_iso_countries
+
 
 # ***---------------------    51_nitrogen   ------------------------------------
 # * (ipcc2006_sep16): IPCC based nitrogen implementation
@@ -556,7 +631,22 @@ cfg$gms$s56_cprice_red_factor <- 1 		# def = 1
 # *   	SSPDB-SSP5-45-REMIND-MAGPIE,SSPDB-SSP5-60-REMIND-MAGPIE,
 # *   	SSPDB-SSP5-Ref-REMIND-MAGPIE,
 # *     coupling
-cfg$gms$c56_pollutant_prices <- "R2M41-SSP2-NPi"         # def = "R2M41-SSP2-NPi"
+# * Note: c56_pollutant_prices_select applies to countries selected in policy_countries56
+# * c56_pollutant_prices_noselect applies to all other countries. In default, all countries
+# * are selected and only c56_pollutant_prices_select has to be set.
+# * Available scenarios for c56_pollutant_prices_noselect are identical to c56_pollutant_prices_select
+# * (see above) except for emulator and coupling (which can only be chosen for c56_pollutant_prices_select)
+cfg$gms$c56_pollutant_prices_select <- "R2M41-SSP2-NPi"     # def = R2M41-SSP2-NPi
+cfg$gms$c56_pollutant_prices_noselect <- "R2M41-SSP2-NPi"     # def = R2M41-SSP2-NPi
+
+# * Switch and specification of countries for which pollutant pricing in
+# * c56_pollutant_prices_select applies.
+# * Options: list of iso-codes of countries where ghg policy should be applied
+# * Note: must be written in the format: "IND, BRA, DEU"
+# * Default: all iso countries
+# * Note: Only for exogenous scenarios. Does not take any effect when "coupling" /
+# * "emulator" is selected above.
+cfg$gms$policy_countries56  <- all_iso_countries
 
 # * Switch for C price driven afforestation (1=on 0=off)
 cfg$gms$s56_c_price_induced_aff <- 1
@@ -662,7 +752,20 @@ cfg$gms$c60_1stgen_biodem <- "const2020"               # def = const2020
 # *   	SSPDB-SSP5-45-REMIND-MAGPIE,SSPDB-SSP5-60-REMIND-MAGPIE,
 # *   	SSPDB-SSP5-Ref-REMIND-MAGPIE,
 # *     coupling
-cfg$gms$c60_2ndgen_biodem <- "R2M41-SSP2-NPi"     # def = R2M41-SSP2-NPi
+# * Note: c60_2ndgen_biodem_select applies to countries selected in scen_countries60
+# * c60_2ndgen_biodem_noselect applies to all other countries. In default, all countries
+# * are selected and only c60_2ndgen_biodem_select has to be set.
+cfg$gms$c60_2ndgen_biodem_select <- "R2M41-SSP2-NPi"     # def = R2M41-SSP2-NPi
+cfg$gms$c60_2ndgen_biodem_noselect <- "R2M41-SSP2-NPi"     # def = R2M41-SSP2-NPi
+
+# * Switch and specification of countries for which 2nd gen bioenergy demand
+# * scenario applies.
+# * Options: list of iso-codes of countries where scenario should be applied
+# * Note: must be written in the format: "IND, BRA, DEU"
+# * Default: all iso countries
+# * Note: Only for exogenous scenarios. Does not take any effect when "coupling" /
+# * "emulator" selected.
+cfg$gms$scen_countries60  <- all_iso_countries
 
 # * residue demand for 2nd generation bioenergy scenarios
 # * options:    ssp1, ssp2, ssp3, ssp4, ssp5, off

diff --git a/main.gms b/main.gms
@@ -148,23 +148,23 @@ $title magpie
 * 
 * Used data set: isimip_rcp-IPSL_CM5A_LR-rcp2p6-co2_rev42_c200_690d3718e151be1b450b394c1064b1c5.tgz
 * md5sum: 94c214b0a7b46f403dd0aab57b7e476c
-* Repository: /p/projects/landuse/data/input/archive
+* Repository: http://rse.pik-potsdam.de/data/magpie/intern
 * 
-* Used data set: rev4.37_690d3718e151be1b450b394c1064b1c5_magpie.tgz
-* md5sum: 3e22e59029e3cd043e247b435cc773ce
-* Repository: /p/projects/rd3mod/inputdata/output
+* Used data set: rev4.42_690d3718e151be1b450b394c1064b1c5_magpie.tgz
+* md5sum: 531f893581b6e3ea3ca871d76e0dd8a3
+* Repository: http://rse.pik-potsdam.de/data/magpie/intern
 * 
-* Used data set: rev4.37_690d3718e151be1b450b394c1064b1c5_validation.tgz
-* md5sum: 7573940b148f440a5f2d2ceedf400afc
-* Repository: /p/projects/rd3mod/inputdata/output
+* Used data set: rev4.42_690d3718e151be1b450b394c1064b1c5_validation.tgz
+* md5sum: ad8da81da4ddbb3cf4b341b40ceec3fb
+* Repository: http://rse.pik-potsdam.de/data/magpie/intern
 * 
-* Used data set: calibration_H12_c200_12Sep18.tgz
-* md5sum: 0a7d88e902918eb6a5263faaf066cc5d
-* Repository: /p/projects/landuse/data/input/calibration
+* Used data set: calibration_H12_c200_26Feb20.tgz
+* md5sum: 0fc75955deb5916d4ac078e8d817feda
+* Repository: http://rse.pik-potsdam.de/data/magpie/intern
 * 
 * Used data set: additional_data_rev3.77.tgz
 * md5sum: 28184d92028972c171145399ab21fdca
-* Repository: /p/projects/landuse/data/input/archive
+* Repository: http://rse.pik-potsdam.de/data/magpie/intern
 * 
 * Low resolution: c200
 * High resolution: 0.5
@@ -177,7 +177,7 @@ $title magpie
 * 
 * Regionscode: 690d3718e151be1b450b394c1064b1c5
 * 
-* Regions data revision: 4.37
+* Regions data revision: 4.42
 * 
 * lpj2magpie settings:
 * * LPJmL data folder: /p/projects/landuse/data/input/lpj_input/isimip_rcp/IPSL_CM5A_LR/rcp2p6/co2
@@ -197,7 +197,7 @@ $title magpie
 * 
 * 
 * 
-* Last modification (input data): Mon Feb 17 14:01:16 2020
+* Last modification (input data): Wed Mar 04 13:53:55 2020
 * 
 *###################### R SECTION END (VERSION INFO) ###########################
 
@@ -222,7 +222,7 @@ $offlisting
 *                    Key parameters during model runs
 
 $setglobal c_timesteps  coup2100
-$setglobal c_title  SSP2_rev42_rev4p37
+$setglobal c_title  default
 
 scalars
   s_use_gdx   use of gdx files                                       / 2 /

diff --git a/modules/09_drivers/aug17/declarations.gms b/modules/09_drivers/aug17/declarations.gms
@@ -7,7 +7,7 @@
 parameters
 
 * Population
-  im_pop_iso(t,iso)                              Population (mio. per yr)
+  im_pop_iso(t_all,iso)                              Population (mio. per yr)
   im_pop(t,i)                                    Population (mio. per yr)
 
 * GDP in MER

diff --git a/modules/09_drivers/aug17/preloop.gms b/modules/09_drivers/aug17/preloop.gms
@@ -12,7 +12,7 @@ im_demography(t,iso,sex,age) = f09_demography(t,iso,"%c09_gdp_scenario%",sex,age
 
 i09_gdp_ppp_iso(t,iso) = f09_gdp_ppp_iso(t,iso,"%c09_gdp_scenario%");
 i09_gdp_mer_iso(t,iso) = f09_gdp_mer_iso(t,iso,"%c09_gdp_scenario%");
-im_pop_iso(t,iso) = f09_pop_iso(t,iso,"%c09_pop_scenario%");
+im_pop_iso(t_all,iso) = f09_pop_iso(t_all,iso,"%c09_pop_scenario%");
 * get the scenario GDP & Populaiton data for MAgPIE regions
 i09_gdp_mer(t,i) = sum(i_to_iso(i,iso),i09_gdp_mer_iso(t,iso));
 i09_gdp_ppp(t,i) = sum(i_to_iso(i,iso),i09_gdp_ppp_iso(t,iso));

diff --git a/modules/12_interest_rate/glo_jan16/declarations.gms b/modules/12_interest_rate/glo_jan16/declarations.gms
@@ -6,6 +6,9 @@
 *** |  Contact: magpie@pik-potsdam.de
 
 parameters
-  pm_interest(i)        Current interest rate in each region (% per yr)
-  p12_interest(t,i)	    Interest rate (% per yr)
+ pm_interest(i)        Current interest rate in each region (% per yr)
+ p12_interest(t,i)	    Interest rate (% per yr)
+* country-specific scenario switch
+ p12_country_dummy(iso)               Dummy parameter indicating whether country is affected by interest rate scenario (1)
+ p12_reg_shr(t_all,i) Weighted share of region with regards to interest rate scenario of countries (1)
 ;
diff --git a/modules/12_interest_rate/glo_jan16/input.gms b/modules/12_interest_rate/glo_jan16/input.gms
@@ -7,6 +7,23 @@
 
 $setglobal c12_interest_rate  medium
 
+* Set-switch for countries affected by country-specific interest rate scenario
+* Default: all iso countries selected
+sets
+  gdp_countries12(iso) countries to be affected by chosen interest rate scenario / ABW,AFG,AGO,AIA,ALA,ALB,AND,ARE,ARG,ARM,ASM,ATA,ATF,ATG,AUS,AZE,
+BDI,BEN,BES,BFA,BGD,BHR,BHS,BIH,BLM,BLR,BLZ,BMU,BOL,BRA,BRB,BRN,BTN,BVT,BWA,CAF,CAN,
+CCK,CHN,CHE,CHL,CIV,CMR,COD,COG,COK,COL,COM,CPV,CRI,CUB,CUW,CXR,CYM,DJI,DMA,DOM,DZA,
+ECU,EGY,ERI,ESH,ETH,FJI,FLK,FRO,FSM,GAB,GEO,GGY,GHA,GIB,GIN,GLP,GMB,GNB,GNQ,GRD,GRL,
+GTM,GUF,GUM,GUY,HKG,HMD,HND,HTI,IDN,IMN,IND,IOT,IRN,IRQ,ISL,ISR,JAM,JEY,JOR,JPN,KAZ,
+KEN,KGZ,KHM,KIR,KNA,KOR,KWT,LAO,LBN,LBR,LBY,LCA,LIE,LKA,LSO,MAC,MAF,MAR,MCO,MDA,MDG,
+MDV,MEX,MHL,MKD,MLI,MMR,MNE,MNG,MNP,MOZ,MRT,MSR,MTQ,MUS,MWI,MYS,MYT,NAM,NCL,NER,NFK,
+NGA,NIC,NIU,NOR,NPL,NRU,NZL,OMN,PAK,PAN,PCN,PER,PHL,PLW,PNG,PRI,PRK,PRY,PSE,PYF,QAT,
+REU,RUS,RWA,SAU,SDN,SEN,SGP,SGS,SHN,SJM,SLB,SLE,SLV,SMR,SOM,SPM,SRB,SSD,STP,SUR,SWZ,
+SXM,SYC,SYR,TCA,TCD,TGO,THA,TJK,TKL,TKM,TLS,TON,TTO,TUN,TUR,TUV,TWN,TZA,UGA,UKR,UMI,
+URY,USA,UZB,VAT,VCT,VEN,VGB,VIR,VNM,VUT,WLF,WSM,YEM,ZAF,ZMB,ZWE /
+;
+
+
 table f12_interest(t_all,scen12)  Interest rate scenarios (% per yr)
 $ondelim
 $include "./modules/12_interest_rate/input/f12_interest_rate.cs3"

diff --git a/modules/12_interest_rate/glo_jan16/preloop.gms b/modules/12_interest_rate/glo_jan16/preloop.gms
@@ -12,7 +12,27 @@ t_to_i_to_dev(t,i,"lic") = yes$(im_gdp_pc_ppp(t,i) <= 1045);
 t_to_i_to_dev(t,i,"mic") = yes$(im_gdp_pc_ppp(t,i) > 1045 AND im_gdp_pc_ppp(t,i) < 12746);
 t_to_i_to_dev(t,i,"hic") = yes$(im_gdp_pc_ppp(t,i) >= 12746);
 
-$ifthen "%c12_interest_rate%" == "coupling" p12_interest(t,i) = f12_interest_coupling(t);
-$elseif "%c12_interest_rate%" == "mixed" p12_interest(t,i) = sum(t_to_i_to_dev("y1995",i,dev), sum(scen12_to_dev(scen12,dev), f12_interest(t,scen12)));
-$else p12_interest(t,i) = f12_interest(t,"%c12_interest_rate%");
+* Country switch to determine countries for which chosen interest rate scenario
+* shall be applied.
+* In the default case, the interest rate scenario affects all countries when
+* activated.
+p12_country_dummy(iso) = 0;
+p12_country_dummy(gdp_countries12) = 1;
+* Because MAgPIE is not run at country-level, but at region level, a region
+* share is calculated that translates the countries' influence to regional level.
+* Countries are weighted by their population size.
+p12_reg_shr(t_all,i) = sum(i_to_iso(i,iso), p12_country_dummy(iso) * im_pop_iso(t_all,iso)) / sum(i_to_iso(i,iso), im_pop_iso(t_all,iso));
+
+
+$ifthen "%c12_interest_rate%" == "coupling"
+ p12_interest(t,i) = f12_interest_coupling(t)
+ ;
+$elseif "%c12_interest_rate%" == "mixed"
+ p12_interest(t,i) = sum(t_to_i_to_dev("y1995",i,dev), sum(scen12_to_dev(scen12,dev), f12_interest(t,scen12))) * p12_reg_shr(t,i)
+                      + f12_interest(t,"medium") * (1-p12_reg_shr(t,i))
+ ;
+$else
+ p12_interest(t,i) =  f12_interest(t,"%c12_interest_rate%") * p12_reg_shr(t,i)
+                      + f12_interest(t,"medium") * (1-p12_reg_shr(t,i))
+ ;
 $endif
diff --git a/modules/12_interest_rate/reg_feb18/declarations.gms b/modules/12_interest_rate/reg_feb18/declarations.gms
@@ -6,13 +6,17 @@
 *** |  Contact: magpie@pik-potsdam.de
 
 scalars
-  s12_min_dev  Minimum development_state of all regions in 1995 (1)
-  s12_max_dev  Maximum development_state of all regions in 1995 (1)
-  s12_slope_a        Slope of the linear relationship between development state and interest rate (1)
-  s12_intercept_b    Intercept of the linear relationship between development state and interest rate (1)
+  s12_min_dev         Minimum development_state of all regions in 1995 (1)
+  s12_max_dev         Maximum development_state of all regions in 1995 (1)
+  s12_slope_a         Slope of the linear relationship between development state and interest rate (1)
+  s12_intercept_b     Intercept of the linear relationship between development state and interest rate (1)
 ;
 
 parameters
-  pm_interest(i)        Current interest rate in each region (% per yr)
-  p12_interest(t,i)     Interest rate in each region and timestep (% per yr)
+ pm_interest(i)                      Current interest rate in each region (% per yr)
+ p12_interest(t,i)                   Interest rate in each region and timestep (% per yr)
+* country-specific scenario switch
+ p12_interest_reg_policy(t_all)       Regional interest rate policy fader
+ p12_country_dummy(iso)               Dummy parameter indicating whether country is affected by interest rate scenario (1)
+ p12_reg_shr(t_all,i)                 Weighted share of region with regards to interest rate scenario of countries (1)
 ;
diff --git a/modules/12_interest_rate/reg_feb18/input.gms b/modules/12_interest_rate/reg_feb18/input.gms
@@ -7,12 +7,37 @@
 
 $setglobal c12_interest_rate  medium
 
+* Set-switch for countries affected by country-specific interest rate scenario
+* Default: all iso countries selected
+sets
+  gdp_countries12(iso) countries to be affected by chosen interest rate scenario / ABW,AFG,AGO,AIA,ALA,ALB,AND,ARE,ARG,ARM,ASM,ATA,ATF,ATG,AUS,AZE,
+BDI,BEN,BES,BFA,BGD,BHR,BHS,BIH,BLM,BLR,BLZ,BMU,BOL,BRA,BRB,BRN,BTN,BVT,BWA,CAF,CAN,
+CCK,CHN,CHE,CHL,CIV,CMR,COD,COG,COK,COL,COM,CPV,CRI,CUB,CUW,CXR,CYM,DJI,DMA,DOM,DZA,
+ECU,EGY,ERI,ESH,ETH,FJI,FLK,FRO,FSM,GAB,GEO,GGY,GHA,GIB,GIN,GLP,GMB,GNB,GNQ,GRD,GRL,
+GTM,GUF,GUM,GUY,HKG,HMD,HND,HTI,IDN,IMN,IND,IOT,IRN,IRQ,ISL,ISR,JAM,JEY,JOR,JPN,KAZ,
+KEN,KGZ,KHM,KIR,KNA,KOR,KWT,LAO,LBN,LBR,LBY,LCA,LIE,LKA,LSO,MAC,MAF,MAR,MCO,MDA,MDG,
+MDV,MEX,MHL,MKD,MLI,MMR,MNE,MNG,MNP,MOZ,MRT,MSR,MTQ,MUS,MWI,MYS,MYT,NAM,NCL,NER,NFK,
+NGA,NIC,NIU,NOR,NPL,NRU,NZL,OMN,PAK,PAN,PCN,PER,PHL,PLW,PNG,PRI,PRK,PRY,PSE,PYF,QAT,
+REU,RUS,RWA,SAU,SDN,SEN,SGP,SGS,SHN,SJM,SLB,SLE,SLV,SMR,SOM,SPM,SRB,SSD,STP,SUR,SWZ,
+SXM,SYC,SYR,TCA,TCD,TGO,THA,TJK,TKL,TKM,TLS,TON,TTO,TUN,TUR,TUV,TWN,TZA,UGA,UKR,UMI,
+URY,USA,UZB,VAT,VCT,VEN,VGB,VIR,VNM,VUT,WLF,WSM,YEM,ZAF,ZMB,ZWE /
+;
+
 table f12_interest_bound(t,bound12) Lower and higher bounds of interest rates (% per yr)
 $ondelim
 $include "./modules/12_interest_rate/input/f12_interest_rate_bound.cs3"
 $offdelim
 ;
 
+table f12_interest(t_all,scen12)  Interest rate scenarios (% per yr)
+$ondelim
+$include "./modules/12_interest_rate/input/f12_interest_rate.cs3"
+$offdelim
+;
+
+scalars
+  s12_interest_noselect   Interest rate scenario chosen for regional interest rate scenario switch (1) / 0.02 /
+;
 
 $if "%c12_interest_rate%" == "coupling" parameter f12_interest_coupling(t_all) Interest rate (% per yr)
 $if "%c12_interest_rate%" == "coupling" /