anchor.diff:  
     subroutine set_small
   
  -    use meth_params_module, only: small_temp, small_pres, small_dens !, small_ener
  +    use meth_params_module, only: small_temp, small_pres, small_dens, small_ener
   
       implicit none
   
  
positive.code:   double precision :: ambient_density
  
    ! Bulk system motion
  
    double precision :: bulk_velx, bulk_vely, bulk_velz
  
  contains
  
    ! This routine calls all of the other subroutines at the beginning
    ! of a simulation to fill in the basic problem data.
  
    subroutine initialize(name, namlen, init)
  
      use bl_constants_module, only: ZERO
      use bl_error_module, only: bl_error
  
      implicit none
  
      integer :: namlen, i, init
      integer :: name(namlen)
  
      ! Build "probin" filename -- the name of the file containing the fortin namelist.
      allocate(character(len=namlen) :: probin)
  
negative.code:        call init_1d(model_S_r, model_S_state, npts_model, dx, radius_S_initial, mass_S, &
                      central_density_S, stellar_temp, stellar_comp, ambient_state)
  
         if (ioproc == 1 .and. init == 1) then
            print *, "Generated initial model for secondary WD of mass", mass_S, &
                     ", central density", central_density_S, ", and radius", radius_S_initial
         endif
  