anchor.diff:    !-------------------------------------------------------------------------
     ! Subroutine to save in a file the critical points
     !-------------------------------------------------------------------------
  -  subroutine write_critical_points(sigma_c_thin, temp_c_thin,sigma_c_thick, temp_c_thick, tau_thin, tau_thick)
  +  subroutine write_critical_points(sigma_c_thin, temp_c_thin, tau_thin, sigma_c_thick, temp_c_thick, tau_thick)
   
       implicit none
   
       real(kind = x_precision),dimension(n_cell),intent(in) :: sigma_c_thin
       real(kind = x_precision),dimension(n_cell),intent(in) :: temp_c_thin
  +    real(kind = x_precision),dimension(n_cell),intent(in) :: tau_thin
       real(kind = x_precision),dimension(n_cell),intent(in) :: sigma_c_thick
       real(kind = x_precision),dimension(n_cell),intent(in) :: temp_c_thick
  -    real(kind = x_precision),dimension(n_cell),intent(in) :: tau_thin
       real(kind = x_precision),dimension(n_cell),intent(in) :: tau_thick
       !-----------------------------------------------------------------------
       integer             :: i
  
positive.code: 
      enddo
  
      call write_critical_points(sigma_thin, temp_thin, sigma_thick, temp_thick, tau_thin, tau_thick)
  
    end subroutine curve
  
  
negative.code:                                      'Temp_thick','Sigma_thick','Tau_thin','Tau_thick'
  
      do i=1, n_cell
        write(fid_4,'(7(e16.6e2))') r_state%r(i), temp_c_thin(i), sigma_c_thin(i), tau_thin(i),&
                                                  temp_c_thick(i), sigma_c_thick(i), tau_thick(i)
      end do
  
      close(fid_4)
  