adding two more components, which are the microbial species into the chemical components list.
The x-z coordinate setting will cause trouble. Switch back to x-y direction.
changing the mesh to x-y direction.
removing the source term for GROUNDWATER_FLOW, using fixed head to control the flow rate.
fullfill the initial conditions and adding variables to the output.