anchor.diff:  
        dens_ambient = model_S_state(npts_model_S,idens_model)
   
  -     ! check that the cutoff temperature match
  +     ! Check that the cutoff temperature match
        if (model_S_state(npts_model_S,itemp_model) /= &
            model_P_state(npts_model_P,itemp_model)) then
           call bl_error("ERROR: Cutoff temperatures for primary and secondary models do not agree.")
  
positive.code:      use model_parser_module
       use bl_constants_module
       use fundamental_constants_module
       use eos_module
       use com, only: mass_p, mass_s, com_loc_p, com_loc_s
  
  
negative.code: 
       call close_model_file()
  
       ! Compute the mass of the primary
       mass_p_initial = ZERO
       do i = 1, npts_model_P-1
          if (i == 1) then
  