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: 
      implicit none
  
      integer                                  ,intent(in):: n
      real(kind = x_precision),dimension(n)    ,intent(in):: radius
      real(kind = x_precision),dimension(n)    ,intent(in):: sigma_c_thin
      real(kind = x_precision),dimension(n)    ,intent(in):: temp_c_thin
      real(kind = x_precision),dimension(n)    ,intent(in):: sigma_c_thick
      real(kind = x_precision),dimension(n)    ,intent(in):: temp_c_thick
      real(kind = x_precision),dimension(n)    ,intent(in):: tau_thin
      real(kind = x_precision),dimension(n)    ,intent(in):: tau_thick
      !-----------------------------------------------------------------------
      integer                                             :: i
      integer                                             :: fid_4 = 11
  
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)
  