parameters kn instead of k_n Modal ChemArithmetic

models.pl

@@ -276,7 +276,7 @@ inherits(Ancestors) :-
 
 parametrize :-
   biocham_command,
-  doc('Replace numeric constants in a model by parameters initialized with their value.'),
+  doc('Replace numeric constants in a model by parameters k1, k2,... initialized with their value.'),
   parametrize_initial_concentration,
   parametrize_reaction.
 

util.pl

@@ -658,7 +658,7 @@ stutter(N, In, Out) :-
 
 %! new_parameter(-Name)
 %
-% Generate the name of a new (and non-existent) parameter of the form k_n
+% Generate the name of a new (and non-existent) parameter of the form kn
 
 new_parameter(Name) :-
   (
@@ -670,7 +670,7 @@ new_parameter(Name) :-
   ),
   Np is N+1,
   nb_setval(parameter_index, Np),
-  atom_concat('k_',N,Try),
+  atom_concat('k',N,Try),
   (
     single_model(ModelId),
     identifier_kind(ModelId, Try, parameter)