As part of an effort to validate the edge plasma model in the SOLPS-ITER code suite under ITER-relevant divertor plasma and neutral conditions, we report on progress in the modeling of the Alcator C-Mod divertor plasma with the new code. We perform simulations with a complete drifts model and kinetic neutrals, including effects of neutral viscosity, ion-molecule collisions and Lyα-opaque conditions, but assuming a pure deuterium plasma. Through a series of simulations with varying divertor geometries, we show the importance of including neutal leakage paths through the divertor substructure on the divertor plasma solution. Moreover, the impact of drifts on inner-outer target asymmetries is assessed. Including both effects, we achieve excellent agreement between simulations and upstream and outer target Langmuir Probe data. In absence of strong volumetric losses due to e.g. impurity radiation in our simulations, the strong inner target detachment observed experimentally remains elusive in our modeling at present.