UMAT Non linear Viscoelasticity - Algorithm verification?

Hello everybody,

           My task is to write a UMAT for Non-linear Viscoelasicity based on Schapery's model in ABAQUS. I have coded the routine, but I am still not able to predict correct response. I tried to debug it as much as I can. I understand that it would be difficult for anybody to understand my code. So, I have attached a flowchart (pdf) of the whole algorithm that I follow. It would be nice if you could check that and correct me if I am wrong.

Thanks in advance for your time and consideration

Regards

Shree

PS. I have also attached a copy of the paper that I am referring for the model.