Dear Stefaan, Thank you for your time and clear answer, i want to know whether the results depend on case.clmcorup/dn, case.clmvalup/dn and case.clmscup/cn files. I think they shouldn't, but the results seem to depend on them. I'll do calculations with different case.dmatup/dn files and compare the results. yours sincerely, hongxia