double precision sc(mxpart,mxpart) common/scprods/sc