In the init.c file the disc is made to terminate on the surface on the star. In pluto.ini ALPHA is changed to 0.01 from 0.7 and in eta_visc.c file disc is made to terminate on the surface of star and the there is a factor of 0.1 multipled to the analytic expression to supress the numerical oscillation.
corrected typos