Field sigma_c should be removed if we use only the internal value computed as sigma_c = 2*Dn.n/cross_section.