Not exactly, this line should be: `snow->density = ((double)SNOWICE_THRESH + (double)NEW_SNOW_DENSITY)/2.0;`

This is resetting the snow pack density after removing the denser glacial ice.