Saturation: need to have N_cm and tau variable in the parameter file