anchor.diff:      !-----------------------------------------------------------------------
   
       integer                                               :: i
  -
   
       i = max(1, index_fcp)
   
  
positive.code:       write(fid_tot,'(3(A16))') 'Surface_density', 'Temperature', 'Optical_depth'
  
        do j = 1, nb_it
          temp_real(j)  =  temp_real(j) 
          sigma_real(j) =  sigma_real(j) 
  
          call variables(temp_real(j), sigma_real(j), omega, f, optical_depth, tau_eff)
  
  
negative.code:     real(kind = x_precision),intent(in) :: sigma_c_thick
      real(kind = x_precision),intent(in) :: temp_c_thick
      integer                 ,intent(in) :: k
  
      !------------------------------------------------------------------------
  
      write(*,*)'**** Critical Point',k,'********'
      write(*,*)'****************************************'
      write(*,"(' Optically thin (T,sigma) :',1p,E12.4,4x,1p,E12.4)")temp_c_thin,sigma_c_thin
  