calc_smolt function

Smolt production