Fix mass flux units in "-surface ...,forcing".
Add a reference.