Computational domain in the (r,z) plane showing active and ghost grid zones with solid and dashed lines, respectively. The zone centers are marked with filled/open circles for active/ghost zones. The arrows indicate the backward finite-difference scheme used to discretize spatial derivatives. The innermost grid cell highlighted with a backslash palette refers to the seed value of the gas surface density n_0,0.