ok Rob, 

after doing some variations in the settings system ran without error but it is not converging what could be the possible reasons behind that also see the attached picture to understand the behavior of scaled residuals currently epsilon is increasing but in such a way that it is looking constant. after 1000 iteration it was on 1e-1 i don't know what should i change now should i adapt a region... the system in the picture is without patching without adaption.. results are not accurate ... also if you have any information on how to compare the DSL number graph with the simulation number density graph then please guide me ...thank you 

best regards,