Commit 9538a4823a08026ab744111064b047bad820c4bf

Authored by dumoda01
1 parent 70b751d6
Exists in master and in 1 other branch snow

Fonction de flottabilite reparee mais non activee par defaut

Showing 1 changed file with 6 additions and 2 deletions   Show diff stats
src/extras/bio/bio_fasham.F90
... ... @@ -207,9 +207,12 @@
207 207 mu4 = mu4 /secs_pr_day
208 208 gmax = gmax /secs_pr_day
209 209 w_p = w_p /secs_pr_day
  210 + w_pmin = w_pmin /secs_pr_day !DD
  211 + w_pmax = w_pmax /secs_pr_day !DD
  212 + theta = theta /secs_pr_day !DD
210 213 w_d = w_d /secs_pr_day
211 214 alpha= alpha/secs_pr_day
212   - inib = inib /secs_pr_day !CHG1
  215 + inib = inib /secs_pr_day !CHG1
213 216  
214 217 out_unit=unit
215 218  
... ... @@ -574,7 +577,8 @@
574 577 dd(l,b,ci)=vb*cc(l,ci)/(k4+min67+cc(l,ci))*(cc(b,ci)+b0)
575 578  
576 579 ! Taux de chute du phytoplancton en fonction de la croissance
577   -! ws(p,ci)=w_pmin+(w_pmax-w_pmin)*exp(-theta*sigma)
  580 +! ws(p,ci)=w_pmin+(w_pmax-w_pmin)*exp(-ff/theta)
  581 +! ws(p,ci)=w_pmax-(w_pmax-w_pmin)*(ff/theta/(1.+ff/theta))
578 582 ws(p,ci)=w_p
579 583  
580 584 do i=1,numc
... ...