The variables 'biomass' and 'ind' was removed from stomate.f90 and stomate_lpj.f90. Now biomass is only calculated in slowproc.f90 just after stomate_main subroutine call.