Here is a summary of the provided article:

## Infiltration

The surface moisture flux that remains after surface runoff is removed is divided into inputs to surface water (q_in,h2osfc) and the soil (q_in,soil). If q_in,soil exceeds the maximum soil infiltration capacity (q_infl,max), then infiltration excess (Hortonian) runoff is generated (q_infl,excess) and transferred from q_in,soil to q_in,h2osfc. 

The final moisture fluxes are:
- q_in,h2osfc = f_h2osfc * q_in,surface + q_infl,excess - q_evap,h2osfc
- q_in,soil = (1-f_h2osfc) * q_in,surface - q_infl,excess - (1 - f_sno - f_h2osfc) * q_evap,soil

The balance of surface water is then calculated as the change in surface water store (ΔW_sfc) over the time step.

Bottom drainage from the surface water store (q_drain,h2osfc) is added to q_in,soil to give the total infiltration into the surface soil layer (q_infl).

Infiltration (q_infl) and explicit surface runoff (q_over) are not allowed for glaciers.