anchor.diff:      real(kind = x_precision)                               :: f     = 0.0d0
       integer                                                :: optical_depth =0
   
  -    real(kind = x_precision),dimension(nb_it)              :: temp_t_thick  
  -    real(kind = x_precision),dimension(nb_it)              :: sigma_t_thick 
  -    real(kind = x_precision),dimension(nb_it)              :: temp_t_thin   
  -    real(kind = x_precision),dimension(nb_it)              :: sigma_t_thin  
  -    real(kind = x_precision),dimension(nb_it)              :: temp_real     
  -    real(kind = x_precision),dimension(nb_it)              :: sigma_real    
  +    real(kind = x_precision),dimension(nb_it)              :: temp_t_thick
  +    real(kind = x_precision),dimension(nb_it)              :: sigma_t_thick
  +    real(kind = x_precision),dimension(nb_it)              :: temp_t_thin
  +    real(kind = x_precision),dimension(nb_it)              :: sigma_t_thin
  +    real(kind = x_precision),dimension(nb_it)              :: temp_real
  +    real(kind = x_precision),dimension(nb_it)              :: sigma_real
   
       character(len = 8)                                     :: number_of_cell
       !character(len = 64)                                    :: fname_thick
  
positive.code:   subroutine write_critical_points(n,radius, sigma_c_thin, temp_c_thin,sigma_c_thick, temp_c_thick)
  
      implicit none
       
      integer                                  ,intent(in):: n 
      real(kind = x_precision),dimension(n)    ,intent(in):: radius
  
      real(kind = x_precision),dimension(n)    ,intent(in):: sigma_c_thin
  
negative.code:       temp_thick(k)  = log10( temp_c_thick * state_0%T_0 )
        sigma_thick(k) = log10( sigma_c_thick * state_0%S_0 )
        temp_thin(k)   = log10( temp_c_thin * state_0%T_0 )
        sigma_thin(k)  = log10( sigma_c_thin * state_0%S_0 )
  
  
      enddo
  