diff --git a/src/SubModel.cc b/src/SubModel.cc
index 6f504424a006ea0e00828ad40593e41f3fcdd30f..bacd904a1cb6b56f5dabb675207f2e796754293a 100644
--- a/src/SubModel.cc
+++ b/src/SubModel.cc
@@ -1716,7 +1716,7 @@ PacModelTable::writeOutput(ostream& output) const
         string pac_model_name = "M_.pac." + name + ".mce.";
         output << pac_model_name << "z = NaN(" << get<2>(val).size() << ",1);" << endl
                << pac_model_name << "w = NaN(" << get<2>(val).size() << ",1);" << endl
-               << pac_model_name << "b = NaN(" << get<2>(val).size() << ",1);" << endl;
+               << pac_model_name << "b = true(" << get<2>(val).size() << ",1);" << endl;
         for (int component_idx {1};
              auto& [component, growth_component, auxname, kind, coeff, growth_neutrality_param,
                     h_indices, original_growth_component, growth_component_info] : get<2>(val))
@@ -1733,7 +1733,7 @@ PacModelTable::writeOutput(ostream& output) const
               {
                 ostringstream mce_w;
                 coeff->writeOutput(mce_w);
-                output << pac_model_name << "b(" << component_idx << ") = 0;" << endl
+                output << pac_model_name << "b(" << component_idx << ") = false;" << endl
                        << pac_model_name << "w(" << component_idx << ") = " << mce_w.str() << ";"
                        << endl;
               }
@@ -1741,7 +1741,7 @@ PacModelTable::writeOutput(ostream& output) const
                the parameter index
             */
             if ((params.size() == 1) && (typeid(*coeff) == typeid(VariableNode)))
-              output << pac_model_name << "b(" << component_idx << ") = 1;" << endl
+              output << pac_model_name << "b(" << component_idx << ") = true;" << endl
                      << pac_model_name << "w(" << component_idx
                      << ") = " << symbol_table.getTypeSpecificID(*params.begin()) + 1 << ";"
                      << endl;