fixed bug in cSIR.c