Commit d43fcf3e authored by Poulet Camille's avatar Poulet Camille
Browse files

Correction_GR3D_function

parent baffec1a
...@@ -101,7 +101,7 @@ spawnerSurvivalPreReproduction <- function(Triver,Tmin, Topt, Tmax){ ...@@ -101,7 +101,7 @@ spawnerSurvivalPreReproduction <- function(Triver,Tmin, Topt, Tmax){
#River survival #River survival
survProbOptRiver = survivePar$survProbOptRiver survProbOptRiver = survivePar$survProbOptRiver
spRiver = survProbOptRiver * TemperatureEffect(Triver, Tmin, Topt, Tmax) spRiver = survProbOptRiver * temperatureEffect(Triver, Tmin, Topt, Tmax)
#StockAfterSurv = S * SpRiver #StockAfterSurv = S * SpRiver
...@@ -159,7 +159,7 @@ stockRecruitementRelationship <- function(Triver,Tmin, Topt, Tmax, survivalStock ...@@ -159,7 +159,7 @@ stockRecruitementRelationship <- function(Triver,Tmin, Topt, Tmax, survivalStock
alleeEffect = 1/(1 + exp(-log(19)*(survivalStock - n/ratioS95_S50*surfaceWatershed)/(n*surfaceWatershed - n/ratioS95_S50*surfaceWatershed))) alleeEffect = 1/(1 + exp(-log(19)*(survivalStock - n/ratioS95_S50*surfaceWatershed)/(n*surfaceWatershed - n/ratioS95_S50*surfaceWatershed)))
Rj = (alphaj * SurvivalStock * alleeEffect)/(betaj + survivalStock * alleeEffect) Rj = (alphaj * survivalStock * alleeEffect)/(betaj + survivalStock * alleeEffect)
#Rj = ((aj * S) * p)/(Bj +S * p) #Rj = ((aj * S) * p)/(Bj +S * p)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment