Merge branch 'development'

Author: zingale

CHANGES

@@ -1,3 +1,8 @@
+# changes since last release
+
+  -- A new network, subch2, was added that combines aprox13 and the
+     subch networks. (#184)
+
 # 19.05
 
   -- The aprox21 network was missing the analytic Jacobian term for

integration/BS/bs_integrator.F90

@@ -59,7 +59,7 @@ subroutine bs_integrator(state_in, state_out, dt, time, status)
     type (bs_t) :: bs
 
     real(rt) :: ener_offset
-    real(rt) :: edot, t_enuc, t_sound, limit_factor
+    real(rt) :: t_enuc, t_sound, limit_factor
 
     logical :: success
 

integration/BS/bs_jac.F90

@@ -16,14 +16,12 @@ subroutine jac(bs)
     use extern_probin_module, only: jacobian, integrate_temperature, integrate_energy, react_boost
     use burn_type_module, only: burn_t, net_ienuc, net_itemp
     use bs_type_module, only: bs_t, bs_to_burn, burn_to_bs
-    use rpar_indices, only: irp_y_init, irp_t_sound
+    use rpar_indices, only: irp_y_init
 
     implicit none
 
     type (bs_t) :: bs
 
-    real(rt) :: limit_factor, t_sound, t_enuc
-
     integer :: n
 
     ! Initialize the Jacobian to zero.

integration/BS/bs_rhs.F90

@@ -18,14 +18,12 @@ subroutine f_rhs(bs)
     use extern_probin_module, only: integrate_temperature, integrate_energy, react_boost
     use bs_type_module, only: bs_t, clean_state, renormalize_species, update_thermodynamics, &
                               burn_to_bs, bs_to_burn
-    use rpar_indices, only: irp_y_init, irp_t_sound
+    use rpar_indices, only: irp_y_init
 
     implicit none
 
     type (bs_t) :: bs
 
-    real(rt) :: limit_factor, t_sound, t_enuc
-
     ! We are integrating a system of
     !
     ! y(1:nspec_evolve) = dX/dt

integration/jacobian_sparsity.F90

@@ -192,8 +192,6 @@ subroutine get_jac_entry(state, row, col, val)
     integer, intent(in) :: row, col
     real(rt), intent(out) :: val
 
-    integer :: csr_loc
-
     !$gpu
 
 #ifdef REACT_SPARSE_JACOBIAN

networks/aprox13/actual_rhs.F90

@@ -1133,7 +1133,7 @@ subroutine screen_aprox13(btemp, bden, y, rr)
     double precision, intent(in   ) :: y(nspec)
     type (rate_t),    intent(inout) :: rr
 
-    integer          :: i, jscr
+    integer          :: jscr
     double precision :: sc1a,sc1adt,sc2a,sc2adt,sc3a,sc3adt
     double precision :: sc1add,sc2add
 !    double precision :: sc3add