anchor.diff:        use crlibm_lib
         use rates_def
         use net_def
  -	  use photo_out, only: output_star_photo
  -	  use star_utils, only: get_name_for_restart_file
         
         implicit none
         
  
positive.code:       real(dp), parameter :: Zsol = 0.0142
        type (star_info), pointer :: s
        type (Net_General_Info), pointer :: g
        character (len=strlen) :: filename, num_str, fstring
        integer :: num_digits
        
        ierr = 0	 
        call star_ptr(id, s, ierr)
  
negative.code:       
        ierr = 0
        call get_net_ptr(s% net_handle, g, ierr)
        if (ierr /= 0) stop 'bad handle'	  
        
  !     increase VARCONTROL and MDOT: increase varcontrol and Mdot when the model hits the TPAGB phase
        if ((s% initial_mass < 10) .and. (s% center_h1 < 1d-4) .and. (s% center_he4 < 1d-4)) then
  