From e36ee2a19a91460dfb1e4be703cb33907bfd20f4 Mon Sep 17 00:00:00 2001
From: "patrick.lambert" <patrick.lambert@irstea.fr>
Date: Wed, 2 Oct 2019 11:45:50 +0200
Subject: [PATCH] correction of temperature at sea

---
 exploration/scriptR/BruchPrg.R | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/exploration/scriptR/BruchPrg.R b/exploration/scriptR/BruchPrg.R
index 7b00104..9507fc8 100644
--- a/exploration/scriptR/BruchPrg.R
+++ b/exploration/scriptR/BruchPrg.R
@@ -126,11 +126,11 @@ vonBertalaffyIncrement = function(nStep, L0, Linf, K, deltaT, sigma, withTempEff
 # parametres
 # ----------------------------------------------------
 L0 = 2
-Linf = 80
+Linf = 70
 #kopt = 0.3900707
 #koptMale = 0.2
 #koptFemale = 0.3
-timestep = 0.25
+timestep = 0.25 # time step of the simumation
 sigma = 0.2
 
 tempEffect = mean(temperatureEffect( TrefAtSea, 3, 17, 26))
@@ -203,8 +203,9 @@ extractTemp<-as.data.frame(tempData %>%
 )
 
 plot(extractTemp$Lat, extractTemp$meanSpring, pch=20)
+abline(v=45.25) # Garonne
 abline(v=49.25) # vire
-abline(v=53.75) # elbe
+abline(v=52.25) # Rhine
 
 
 abline(h=11)
@@ -212,7 +213,7 @@ abline(h=12.5)
 # ======================================================================
 # exploration of growth for male and female
 # ======================================================================
-correction=mean(temperatureEffect(Tref, 3, 17, 26))
+correction=mean(temperatureEffect(TrefAtSea, 3, 17, 26))
 
 age=seq(0,10,.25)
 present = vonBertalaffyGrowth(age, 2, 60, 0.3900707 * correction)
-- 
GitLab