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: 			  
  			  !save a model and photo
  			  call star_write_model(id, s% job% save_model_filename, ierr)
  	          call get_name_for_restart_file(s% model_number, s% photo_digits, num_str)
  	          filename = trim(s% photo_directory) // '/' // trim(num_str) 
  	          call output_to_file(filename, id, ierr) 
  	          if (ierr /= 0) return   
  			  
  			  !turn off burning
  			  do r=1,g% num_reactions
  
negative.code:       real(dp), parameter :: Zsol = 0.0142
        type (star_info), pointer :: s
        type (Net_General_Info), pointer :: g
        character (len=strlen) :: photoname
        
        ierr = 0	 
        call star_ptr(id, s, ierr)
  