anchor.diff:      real(kind = x_precision),dimension(n_cell),intent(out) :: s              ! Surface density for the critical point
       !------------------------------------------------------------------------
       integer                                                :: i,j,k
  -    integer                                                :: file_length    ! Number of lines of the data saving file
   
       ! For the temperature and s
       real(kind = x_precision)                               :: temp
       real(kind = x_precision)                               :: Smin
       real(kind = x_precision)                               :: Smax
  -    real(kind = x_precision)                               :: Sigma          ! S
  +    real(kind = x_precision)                               :: Sigma          ! Sigma
   
       real(kind = x_precision)                               :: omega          ! Angular velocity
       real(kind = x_precision)                               :: r              ! Radius
  
positive.code:   use mod_read_parameters
    use mod_variables
    implicit none
  
  ! _____________________________________________________________________________________________
  
    ! BE SURE TO CREATE TWO FOLDERS IN YOUR WORKING DIRECTORY : s_curves  & critical_points
  !_____________________________________________________________________________________________
  
  contains
    !-------------------------------------------------------------------------
  
negative.code:         sigma_real(j) =  sigma_real(j) * state_0%S_0
  
          write(fid_tot,fmt = '(3(e16.6e2))') sigma_real(j), temp_real(j), tau_eff
        enddo
        
        close(fid_tot)
  