diff --git a/exploration/NEA_calibration_offline/6-Growth_Calibration.Rmd b/exploration/NEA_calibration_offline/6-Growth_Calibration.Rmd
index d73cc0bbc2507f3ec7ed127a0f5b28f171ac3274..18ef26e8ce9738bf4df118038b53bb001ac93532 100644
--- a/exploration/NEA_calibration_offline/6-Growth_Calibration.Rmd
+++ b/exploration/NEA_calibration_offline/6-Growth_Calibration.Rmd
@@ -26,6 +26,7 @@ library(officedown)
 library(flextable)
 library(stringr)
 library(tictoc)
+library(tibble)
 
 ```
 
@@ -344,9 +345,8 @@ write_csv(growParOptim, "growParOptim.csv")
 
 load("growth_optimal.RData")
 
-
 growthCurves = computeGrowOpt(growth_optimal$par, growthInOffshore)
-
+save(growthCurves, file = "growthCurves.RData")
 
 #growth curve according to metapopulations 
 growthCurves %>% 
@@ -357,7 +357,8 @@ growthCurves %>%
   #scale_x_continuous(limits = c(0, 10)) +
   geom_hline(yintercept = 45, color = "pink", linetype =1)+ #female
   geom_hline(yintercept = 40, color = "lightblue", linetype = 1)+
-  facet_wrap(.~ metapop, ncol =3)
+  geom_vline(aes(age), xintercept = 8, linetype = 1)+
+  facet_wrap(.~ metapop, ncol =1)
 
 #temperature effect 
 temp <- seq(0,45,0.1)
diff --git a/exploration/NEA_calibration_offline/growParOptim.csv b/exploration/NEA_calibration_offline/growParOptim.csv
index fbb7cbc1caa6ed4bbbae4e4fcbe6913d7badefad..722ac8edb25100e683386c355bc585f6ef68371e 100644
--- a/exploration/NEA_calibration_offline/growParOptim.csv
+++ b/exploration/NEA_calibration_offline/growParOptim.csv
@@ -1,7 +1,7 @@
 parameter,value
 Tmin,0
-Topt,4.4903872865652845
-Tmax,51.96494839598789
-L0,5.778843636796766
-Linf,47.707177806062305
-K,0.5061897626434646
+Topt,4.5125688473103995
+Tmax,76.23141515590089
+L0,6.5184840221819975
+Linf,47.53851690994118
+K,0.4964834345987972