FLR

Unit: these are like separate FLBiols so all we need do is apply the current method over the 3rd dim

season: calculations depend on whether it is the last season in a year, if so you have to do the plusgroup calculations and calculate F as the ratio of log(N(a+1,t+1,1)/ N(a,t,n)}, otherwise it is log(N(a,t,i+1)/ N(a,t,1)}

area: this is more difficult as it’s not simply partial F. If you have an infinite F in one area and zero F in another F doesn’t sum to infinity, plus next season/year the area with infinite F will be repopulated.

I think we can implement harvest for unit and season but should turn it off for area for the moment.


Laurie