Add mod file to compute the boundaries (steady states).

models/boundaries.mod

+@#define steady = false
+
+@#include "init.mod"
+
+@#include "symdecls.mod"
+
+@#include "modeqs.mod"
+
+load_params_and_steady_state('eagle_steady.txt');
+
+options_.noprint=true;
+steady(solve_algo = 5);
+
+write_steady_state('initial-condition.inc', M_, oo_);
+
+endval;
+EAA_thetatbar = 4.125;
+RW_thetatbar  = 4.125;
+US_thetatbar  = 4.125;
+EAB_thetatbar = 4.125;
+CH_thetatbar  = 4.125;
+JP_thetatbar  = 4.125;
+LAC_thetatbar = 4.125;
+EAA_thetat = 4.125;
+RW_thetat  = 4.125;
+US_thetat  = 4.125;
+EAB_thetat = 4.125;
+CH_thetat  = 4.125;
+JP_thetat  = 4.125;
+LAC_thetat = 4.125;
+end;
+
+options_.noprint=true;
+steady(solve_algo = 5);
+
+write_steady_state('terminal-condition.inc', M_, oo_);