Temporary hotfixes to solve random infeasibilities.

core/bounds.gms

@@ -24,7 +24,10 @@ AND (not sameas(te,"tnrs"))),
   vm_cap.lo(t,regi,te,"1")$(t.val gt 2021) = 1e-6;
 );
 
-
+*** define bounds for dummy variable
+vm_dummyGeot.fx(ttot, all_regi, all_te, rlf) = 0;
+vm_dummyGeot.up(ttot, all_regi, all_te, rlf)$(sameas(all_te, "geohdr") AND sameas(ttot, "2025")) = 0.002;
+vm_dummyGeot.lo(ttot, all_regi, all_te, rlf)$(sameas(all_te, "geohdr") AND sameas(ttot, "2025")) = -0.002;
 
 *** RP 20160405 make sure that the model also sees the se2se technologies (seel <--> seh2)
 loop(se2se(enty,enty2,te),

core/declarations.gms

@@ -254,6 +254,7 @@ vm_costFuBio(ttot,all_regi)                          "fuel costs from bio energy
 vm_omcosts_cdr(tall,all_regi)                        "O&M costs for spreading grinded rocks on fields"
 vm_costpollution(tall,all_regi)                      "costs for air pollution policies"
 vm_emiFgas(tall,all_regi,all_enty)                   "F-gas emissions by single gases from IMAGE"
+vm_dummyGeot(ttot, all_regi, all_te, rlf)                 "dummy variable."
 ;
 
 ***----------------------------------------------------------------------------------------

core/equations.gms

@@ -396,7 +396,7 @@ qm_fuel2pe(t,regi,peRicardian(enty))..
 ***---------------------------------------------------------------------------
 *ml* assuming maxprod to be technical potential
 q_limitProd(t,regi,teRe2rlfDetail(teReNoBio(te),rlf))..
-  pm_dataren(regi,"maxprod",rlf,te)
+  pm_dataren(regi,"maxprod",rlf,te)+vm_dummyGeot(t, regi, te, rlf)
   =g=
   ( 1$teRLDCDisp(te) +  pm_dataren(regi,"nur",rlf,te)$(NOT teRLDCDisp(te)) ) * vm_capFac(t,regi,te) * vm_capDistr(t,regi,te,rlf);
   

modules/30_biomass/magpie_40/bounds.gms

@@ -15,8 +15,8 @@
 
 *** Prescribe upper and lower limit for first generation biomass from 2030/45 on, so REMIND has freedom before.
 *** To avoid infeasibilities it was necessary to modify the initial vintage structure for bioeths.
-vm_fuExtr.up(t,regi,"pebios","5")$(t.val ge 2045)  = p30_datapebio(regi,"pebios","5","maxprod",t);
-vm_fuExtr.up(t,regi,"pebioil","5")$(t.val ge 2030) = p30_datapebio(regi,"pebioil","5","maxprod",t);
+vm_fuExtr.up(t,regi,"pebios","5")$(t.val ge 2045)  = 1.2*p30_datapebio(regi,"pebios","5","maxprod",t);
+vm_fuExtr.up(t,regi,"pebioil","5")$(t.val ge 2030) = 1.2*p30_datapebio(regi,"pebioil","5","maxprod",t);
 
 if(cm_1stgen_phaseout=0,
     vm_fuExtr.lo(t,regi,"pebios","5")$(t.val ge 2030)  = p30_datapebio(regi,"pebios","5","maxprod",t)*0.9;