In the files ets_visc.c and init.c I have made the disc terminate on the star, and figured out that the spurious numerical oscillations at the end of 37 % and beginning of 38 % of integration did not occur as in the case of disc not terminating on the surface of star.