write MA(k) kinetics

gpac.pl

@@ -1639,19 +1639,19 @@ monomial_to_reaction([R,Exp], Product, Name_list) :-
    (
       R > 0
    ->
-      exponant_to_solution(Exp, Name_list, Sol, Kin),
+      exponant_to_solution(Exp, Name_list, Sol, _Kin),
       Reac = '_',
       Rate = R,
       Prod = Product
     ;
-      exponant_to_solution(Exp, Name_list, _Sol, Kin),
+      exponant_to_solution(Exp, Name_list, _Sol, _Kin),
       nth1(N, Name_list, Product),
       length(Exp, Len),
       Left is N-1,
       Right is Len-N,
       displace_exponent([-1],Left,Right,EmY),
       add_list(Exp, EmY, Exp_cor),
-      exponant_to_solution(Exp_cor, Name_list, Sol, _Kin),
+      exponant_to_solution(Exp_cor, Name_list, Sol, _Kinetics),
       Reac = Product,
       Rate is -R,
       Prod = '_'
@@ -1659,9 +1659,9 @@ monomial_to_reaction([R,Exp], Product, Name_list) :-
    (
       Sol = '_'
    ->
-      add_reaction(Rate*Kin for Reac=>Prod)
+      add_reaction('MA'(Rate) for Reac=>Prod)
    ;
-      add_reaction(Rate*Kin for Reac=[Sol]=>Prod)
+      add_reaction('MA'(Rate) for Reac=[Sol]=>Prod)
    ).
 
 
@@ -1701,7 +1701,7 @@ negate_name([], _N, [], []) :- !.
 negate_name([Name|TailName], N, [N|TailNeg], [Name_p,Name_m|TailNewName]) :-
   !,
   name_p_m(Name,Name_p,Name_m),
-  add_reaction(fast*Name_p*Name_m for Name_p+Name_m=>_),
+  add_reaction('MA'(fast) for Name_p+Name_m=>_),
   NN is N+1,
   negate_name(TailName, NN, TailNeg, TailNewName).