sprat-xsa-experimental

1
sr.flsr.bevholt <- function (ssb, a, b) {
2
  r <- (ssb * a) / (ssb + b)
3
  return(r)  
4
}
5

6
sr.eqsim.bevholt <- function (ssb, a, b) {
7
  r <- ((ssb*a)/(1 + ssb*b))
8
  return(r)
9
}
10

11
sr.segreg <- function(ssb, a, b) {
12
  if (ssb > b)
13
    return (a*b)
14
  
15
  return (a*ssb)
16
}
17

18
# test values
19
a <- 754.26
20
b <- 1.0042e+8
21
ssb <- 89255193
22

23

24
rec <- sr.segreg(ssb, a, b)
25
print(rec/1e+9)
26

27
# segreg: 170t.t -> 75.7 m.rec
28
# bevholt: 170t.t -> 84.1 m.rec
29
# bevholt: 89t.t -> 59.5 m.rec
30
# segreg: 89t.t -> 67.3