Fix eps_nuc rate output for some reverse rates

Sometimes b_r had the negative built, sometimes it didn't

Author: rjfarmer

net/private/net_derivs.f90

@@ -716,7 +716,7 @@ subroutine get_general_2_to_2_derivs(i,ierr)
             n% eps_nuc_categories(icat_r) = n% eps_nuc_categories(icat_r) + b_r
             if (show_eps_nuc .and. abs(b) > 1d2) &
                write(*,1) trim(reaction_Name(ir)) // ' eps_nuc',  b, b_f, b_r
-                        
+
             n% eps_nuc_rate(i) = b_f * Qconv
             n% eps_nuc_rate(r_i) = b_r  * Qconv
             n% eps_neu_rate(i) = 0d0
@@ -863,7 +863,7 @@ subroutine get_basic_2_to_2_derivs(i,ierr)
                write(*,1) trim(reaction_Name(ir)) // ' eps_nuc',  b, b_f, b_r
 
             n% eps_nuc_rate(i) = b_f * Q  * Qconv
-            n% eps_nuc_rate(r_i) = b_r * Q  * Qconv
+            n% eps_nuc_rate(r_i) = -b_r * Q  * Qconv
             n% eps_neu_rate(i) = 0d0
             n% eps_neu_rate(r_i) = 0d0
 
@@ -1005,7 +1005,7 @@ subroutine get_basic_2_to_1_derivs(i,ierr)
                write(*,1) trim(reaction_Name(ir)) // ' eps_nuc',  b, b_f, b_r
 
             n% eps_nuc_rate(i) = b_f * Q  * Qconv
-            n% eps_nuc_rate(r_i) = b_r * Q  * Qconv
+            n% eps_nuc_rate(r_i) = -b_r * Q  * Qconv
             n% eps_neu_rate(i) = 0d0
             n% eps_neu_rate(r_i) = 0d0