anchor.diff:           index_fcp, index_scp, sigma_c_thin, temp_c_thin, nb_it)
       implicit none
   
  +    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
       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
   
  
positive.code:   subroutine build_s_curve(sigma_real_thick, sigma_real_thin, temp_real_thick, nb_it, index_scp, sigma_real, temp_real)
      implicit none
  
      real(kind = x_precision),dimension(nb_it),intent(in) :: sigma_real_thick
      real(kind = x_precision),dimension(nb_it),intent(in) :: sigma_real_thin
      real(kind = x_precision),dimension(nb_it),intent(in) :: temp_real_thick
      integer                                  ,intent(in) :: nb_it
      integer                                  ,intent(in) :: index_scp
      !real(kind = x_precision),dimension(:),allocatable,intent(out) :: sigma_real
      !real(kind = x_precision),dimension(:),allocatable,intent(out) :: temp_real
  
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)
  