.. _oyster.benthos:

OYSTER : Filtrage benthique
===========================

* Ce module complémetnaire permet de prendre en compte la présence d'huitres sur le fond qui filtrent le matériel phytoplanctonique.
* Dans cette option, les huitres ne sont pas simulées en tant que variable d'état, mais leur densité par maille agit sur le terme de "filtrage benthique" (:ref:`broutage.benthos`)
* Ce "filtrage benthique" correspond à la filtration des huitres qui dépend des MES et de la température.
* Les paramètres qui sont utilisés dans ce module sont donnnés dans la namelist **namoysterfiltbent** dans *parabiolo.txt* (:ref:`oyster.parameters`)

.. note::
   
   Ce module est activé en utilisant la clé **key_oyster_benthos** (+ **key_biolo** )
   

La distribution des huitres
---------------------------
* La distribution spatiale des huitres est determinée par le nombre d'huitres *nbhuitre(i,j)* 
* *nbhuitre* doit être initialisé, estimé ou lu dans un fichier. 
* Cette définition est réalisée par l'utilisateur dans *bioloinit.F90*
* la hauteur des tables y peuvent être également définies  (*hautable*) (en m)

La filtration 
-------------

* L'intensité de la filtration dépend de la température et de la quantité de matière en suspension (MES) dans l'eau.
* Un effet de colmatage est pris en compte :

	| :math:`Colmat_{huitre} = \min (0.,X_{colmat_thr} - MES)`
	| :math:`Effettemp_{huitre} = X_{Filt-temp1} ∗ [TEMPERATURE−X_{Tempopt-Filt}]^2`

* Si MES < *p_huitre_thr_mes*  

	| :math:`Filt= X_{filtmax} - Effettemp_{huitre}`
	
* sinon  

	| :math:`Filt= X_{filt1}*MES + X_{filt2} - Effettemp_{huitre}`
	| :math:`Filt= Filt * (X_{DW}^{X_{expfilt}}`
	| :math:`Filt= Filt * \exp ( X_{colmat} * Colmat_{huitre} )`
	| Si Filt < 0  ==> Filt=0.
	| :math:`T_{filtbenth} = Filt * nbhuitre(i,j) /surf(i,j)`
	
* avec :

	* :math:`X_{colmat_thr}` = *p_huitre_thr_colmat*
	* :math:`X_{Filt-temp1}` = *p_huitre_temp1*	
	* :math:`X_{Tempopt-Filt}` = *p_huitre_temp_opt*	
	* :math:`X_{filtmax}` = *p_huitre_filt_std*
	* :math:`X_{filt1}` = *p_huitre_loi_mes1*
	* :math:`X_{filt2}` = *p_huitre_loi_mes2*
	* :math:`X_{DW}` = *p_huitre_DW*
	* :math:`X_{expfilt}` = *p_huitre_exp_allom*
	* :math:`X_{colmat}` = *p_huitre_loi_colmat*