anchor.diff:    subroutine curve(temperature, s)
       implicit none
   
  -    real(kind = x_precision),dimension(n_cell),intent(out)  :: temperature
  -    real(kind = x_precision),dimension(n_cell),intent(out)  :: s
  +    real(kind = x_precision),dimension(n_cell),intent(out) :: temperature
  +    real(kind = x_precision),dimension(n_cell),intent(out) :: s
   
       integer                                                :: i     = 0
       integer                                                :: j     = 0
  
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
  
      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) :: sigma_real_thin
      real(kind = x_precision),dimension(nb_it),intent(in) :: temp_real_thick
      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
      real(kind = x_precision),dimension(nb_it),intent(out) :: sigma_real
  
negative.code: 
      enddo
  
      !call write_critical_points(r_state%r, rs, sigma_thin, temp_thin, sigma_thick, temp_thick)
  
    end subroutine curve
  
  