anchor.diff:  
           call variables(temp_real(j), sigma_real(j), omega, f, optical_depth, tau_eff)
   
  -        call save_data(fid_tot, fname_tot, nb_it, file_length, sigma_real(j), temp_real(j), tau_eff)
  +        ! call save_data(fid_tot, fname_tot, nb_it, file_length, sigma_real(j), temp_real(j), tau_eff)
  +
  +        write(fid_tot,fmt = '(3(e16.6e2))') sigma_real(j), temp_real(j), tau_eff
  +
           file_length = file_length + 1
         enddo
  -
  +      
  +      close(fid_tot)
   
         ! ADIMENSIONED critical T and S [OUTPUT]
   
  
positive.code:  !____________________________________________________________________________________
  
      do k = 1, n_cell
         file_length = 0
  
         r     = r_state%r(k)
         omega = x_state%Omega(k)
  
negative.code:     character(len = 8)                                     :: number_of_cell
      character(len = 64)                                    :: fname_tot
      integer                                                :: fid_tot
      integer                                                :: ios
   ! For the critical points
      integer                                                :: index_fcp
      real(kind = x_precision)                               :: sigma_c_thick
  