anchor.diff:    !-------------------------------------------------------------------------
     subroutine first_critical_point(sigma_real_thick, temp_real_thick, index_fcp, sigma_c_thick, temp_c_thick, nb_it)
       implicit none
  -    integer,intent(in)                                   :: nb_it
  -
  +
  +    integer                                  ,intent(in) :: nb_it
       real(kind = x_precision),dimension(nb_it),intent(in) :: sigma_real_thick
       real(kind = x_precision),dimension(nb_it),intent(in) :: temp_real_thick
   
  
positive.code:          index_fcp, index_scp, sigma_c_thin, temp_c_thin, nb_it)
      implicit none
  
      real(kind = x_precision),dimension(nb_it),intent(in) :: sigma_real_thick
      real(kind = x_precision),dimension(nb_it),intent(in) :: temp_real_thick
      real(kind = x_precision),dimension(nb_it),intent(in) :: sigma_real_thin
  
      integer                                  ,intent(in) :: nb_it
      integer                                  ,intent(in) :: index_fcp
      integer                                  ,intent(out):: index_scp
  
  
negative.code: 
        do j = 1, nb_it
  
          temp_real(j)  = log10( temp_real(j) * state_0%T_0 )
          sigma_real(j) = log10( sigma_real(j) * state_0%S_0 )
  
          write(fid_tot,'(1p,E12.6,4x,1p,E12.6,4x,1p,E12.6)') sigma_real(j), temp_real(j)
        enddo
  
        !close(fid_thick)
  