diff --git a/src/binary_c_derivatives.h b/src/binary_c_derivatives.h
index e44e9b395fc98031925b733e0e5a381acc329f0c..95837943accd0dddae4e7399ebc263d62a736e29 100644
--- a/src/binary_c_derivatives.h
+++ b/src/binary_c_derivatives.h
@@ -148,6 +148,8 @@
         "Central helium",                       \
         "Central carbon",                       \
         "He core mass (no 3DUP)",               \
+        "Number of thermal pulses",             \
+        "Number of thermal pulses since mcmin", \
         "nova H layer mass",                    \
         "Roche radius"                          \
     }
diff --git a/src/binary_star_functions/limit_accretion_rates.c b/src/binary_star_functions/limit_accretion_rates.c
index ee6b3fa39619478d638d7b71f7de9c7913ca4d8d..ead3808741db789ee737c11c89a3aed9607db013 100644
--- a/src/binary_star_functions/limit_accretion_rates.c
+++ b/src/binary_star_functions/limit_accretion_rates.c
@@ -409,7 +409,7 @@ void limit_accretion_rates(struct stardata_t * RESTRICT const stardata)
                                      * Remove all the accreted layer : anything 
                                      * remaining becomes part of the white dwarf.
                                      */
-                                    accretor->derivative[DERIVATIVE_STELLAR_H_LAYER_MASS] -=
+                                    accretor->derivative[DERIVATIVE_STELLAR_H_LAYER_MASS] +=
                                         dm_loss_rate;
                                     
                                     /*