FishNutrient.java 20.06 KiB
/**
 * 
 * @author Camille Poulet, Patrick Lambert
 * @copyright Copyright (c) 2019, Irstea
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
package species;
import environment.SeaBasin;
import environment.Time.Season;
import fr.cemagref.simaqualife.pilot.Pilot;
import species.DiadromousFish.Gender;
import species.DiadromousFish.Stage;
import species.ReproduceAndSurviveAfterReproduction;
import java.util.ArrayList;
import java.util.HashMap;
import java.util.Hashtable;
import java.util.List;
import java.util.Map;
import com.thoughtworks.xstream.XStream;
import com.thoughtworks.xstream.io.xml.DomDriver;
/**
 * @author camille.poulet
public class FishNutrient {
	private static enum SpawningPosition  {PRE,POST}; // on crer un static pour rserver une mme classe mmoire pour toutes les instances 
	private ArrayList<String> nutrientsOfInterest;
	private double residenceTime;
	private Map<String, Double> excretionRate; 
	/**
	 * Main feature for weight (g) computation before spawning i.e. gametes expelling according to gender, for a given length (cm) 
	 * //Voir pour un retour  la ligne lors du commentaire 
	 * key gender
	 * value
	 * 		key feature
	 * 		value value
	private Map <Gender,Map<String, Double>> fishFeaturesPreSpawning;
	private Map <Gender, Map<String, Double>> fishFeaturesPostSpawning;
	private Map <String, Double> juvenileFeatures;
	/**
	 *  Weight  of gametes spawned for both males and females 
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
* key gender * value g * usually computed as the difference between unspawned gonad (inbound) and spawned gonad (outbound; "spent gonad") */ private Map <Gender, Double> spawnedGametesWeight; /** * chemical composition of carcass before gametes expelling (before spawning) i.e. soma + gonads + gametes * <key> gender * <value> * <key> chemical element * <value> value ratio element / (total wet weight) g/g */ private Map<DiadromousFish.Gender,Map<String,Double>> compoCarcassPreSpawning; //package permettant la cration d'une table de hachage ie fonctionnant en cl -valeur. Cl unique, mais valeur peut tre associe plusieurs cls; //La class d'objet Map a pour point faible la taille des donnes stocker. Plus on a de valeur dans la table, plus c'est lourd et lent! Donc, trouver un compromis entre temps de calcul et espace. //N'accepte pas la valeur nulle et thread safe i.e. utilisable simultanment par plusieurs lments du programme. /** * chemical composition of carcass after spawning i.e. soma + spent gonad (without gametes) * <value> * <key> chemical element * <value> value */ private Map<DiadromousFish.Gender, Map<String, Double>> compoCarcassPostSpawning; /** * chemical composition of gametes estimated from the MIGADO dataset (BDalosesBruch) * <key> gender * <value> * <key> chemical element * <value> value */ private Map<DiadromousFish.Gender, Map<String,Double>> compoGametes; // For juveniles - Based on Taverny (1991) /** * chemical composition of juveniles * <key> stage * <value> * <key> chemical element * <value> value */ private Map<String,Double> compoJuvenile; /** * */ public FishNutrient() { // TODO Auto-generated constructor stub } /** * Constructor based on the 5 Map of fish composition * @param fishFeaturesPreSpawning * @param compoCarcassPreSpawning * @param compoCarcassPostSpawning * @param compoGametes * @param compoJuvenile */ public FishNutrient(ArrayList<String> nutrientsOfInterest, double residenceTime, Map <String, Double> excretionRate,
141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210
Map<Gender, Map<String, Double>> fishFeaturesPreSpawning, Map<Gender, Map<String, Double>> fishFeaturesPostSpawning, Map<Gender, Double> gameteSpawned, Map<Gender, Map<String, Double>> compoCarcassPreSpawning, Map<Gender, Map<String, Double>> compoCarcassPostSpawning, Map<Gender, Map<String, Double>> compoGametes, Map<String, Double> juvenileFeatures, Map<String, Double> compoJuvenile) { super(); this.nutrientsOfInterest = nutrientsOfInterest; this.excretionRate = excretionRate; this.residenceTime = residenceTime; this.fishFeaturesPreSpawning = fishFeaturesPreSpawning; this.fishFeaturesPostSpawning = fishFeaturesPostSpawning; this.spawnedGametesWeight = gameteSpawned; this.compoCarcassPreSpawning = compoCarcassPreSpawning; this.compoCarcassPostSpawning = compoCarcassPostSpawning; this.compoGametes = compoGametes; this.juvenileFeatures = juvenileFeatures; this.compoJuvenile = compoJuvenile; } /** * compute the nutrient fluxes for a single fish (in the super individual) * that dies before spawning * @param fish */ public Map<String,Double> computeNutrientsInputForDeathBeforeSpawning(DiadromousFish fish, ArrayList<String> nutrientsOfInterest) { Map<String,Double> nutrientsInput = new Hashtable<String, Double>(); // On crer ici une Map, classe mre des hashtable (Homme = classe mere ie Map//Jules = hashtable) for (String nutrient : nutrientsOfInterest) { if (fish.getStage()== Stage.MATURE) { double totalWeightPre = this.getWeight(fish,SpawningPosition.PRE); double carcass = totalWeightPre * compoCarcassPreSpawning.get(fish.getGender()).get(nutrient); double excretion = totalWeightPre * residenceTime * excretionRate.get(nutrient) ; double nutrientImport = carcass + excretion; nutrientsInput.put(nutrient, nutrientImport); } else { nutrientsInput.put(nutrient, 0.); } } //TODO Multiply by fish amount return nutrientsInput; } public Map<String,Double> computeNutrientsInputForDeathBeforeSpawning(DiadromousFish fish) { return computeNutrientsInputForDeathBeforeSpawning(fish,this.nutrientsOfInterest); } /** * compute the nutrient fluxes for a single fish (in the super individual) * that dies after spawning (gametes expelling) * @param fish * @return nutrientsInput */ public Map<String, Double> computeNutrientsInputForDeathAfterSpawning(DiadromousFish fish, ArrayList<String> nutrientsOfInterest) { Map<String,Double> nutrientsInput = new Hashtable<String,Double>(); for (String nutrient : nutrientsOfInterest) {
211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280
if (fish.getStage()== Stage.MATURE) { double totalWeightPost = this.getWeight(fish, SpawningPosition.POST); double carcass = totalWeightPost * compoCarcassPostSpawning.get(fish.getGender()).get(nutrient); //double gametes = (totalWeightPre - totalWeightPost) FAUX car perte de poids somatique due a la reproduction double gametes = spawnedGametesWeight.get(fish.getGender()) *compoGametes.get(fish.getGender()).get(nutrient); //TODO: FAUX ! Revoir comment calculer les gamtes double excretion = totalWeightPost * residenceTime * excretionRate.get(nutrient); double nutrientImport = carcass + gametes + excretion; nutrientsInput.put(nutrient,nutrientImport); } else { nutrientsInput.put(nutrient,0.); } } return nutrientsInput; } public Map<String, Double> computeNutrientsInputForDeathAfterSpawning(DiadromousFish fish){ return computeNutrientsInputForDeathAfterSpawning(fish, this.nutrientsOfInterest); } /** * compute the nutrient fluxes for a single fish (in the super individual) * that survives after spawning * Map: model output = element of interest ie string + double ie the quantification of this fluxes. * @return nutrientsInput */ public Map<String,Double>computeNutrientsInputForSurvivalAfterSpawning(DiadromousFish fish, ArrayList<String> nutrientsOfInterest) { Map<String,Double> nutrientsInput = new Hashtable<String,Double>(); for (String nutrient: nutrientsOfInterest) { if (fish.getStage()==Stage.MATURE) { //TODO Fix with new data double totalWeightPost = this.getWeight(fish, SpawningPosition.POST); //Gamete compositions depends on sex. double gametes = spawnedGametesWeight.get(fish.getGender()) * compoGametes.get(fish.getGender()).get(nutrient); double excretion = totalWeightPost * residenceTime * excretionRate.get(nutrient); double nutrientImport = gametes + excretion; nutrientsInput.put(nutrient, nutrientImport); } else { nutrientsInput.put(nutrient,0.); } } return nutrientsInput; } public Map<String,Double>computeNutrientsInputForSurvivalAfterSpawning(DiadromousFish fish) { return computeNutrientsInputForSurvivalAfterSpawning(fish, this.nutrientsOfInterest); } public Map<String,Double> computeNutrientsExportForJuveniles(DiadromousFish juvenileFish, ArrayList<String>nutrientsOfInterest) { Map<String,Double> nutrientsExport = new Hashtable<String,Double>(); for(String nutrient: nutrientsOfInterest) { if(juvenileFish.getStage()==Stage.IMMATURE) {
281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350
double JuvenileMass = this.getWeight(juvenileFish); nutrientsExport.put(nutrient, JuvenileMass * compoJuvenile.get(nutrient)); } } return nutrientsExport; } public Map<String,Double> computeNutrientsExportForJuveniles(DiadromousFish juvenileFish){ return computeNutrientsExportForJuveniles(juvenileFish, this.nutrientsOfInterest); } /** * Compute the weight for a fish with length (cm) * @param fish * @return weight (g) */ public double getWeight (DiadromousFish fish, SpawningPosition spawningPosition) { double weight = 0.; if (fish.getStage()==Stage.IMMATURE) weight = juvenileFeatures.get("aLW") * Math.pow(fish.getLength(),juvenileFeatures.get("bLW")); else //Stage.MATURE if (spawningPosition == SpawningPosition.PRE) weight = fishFeaturesPreSpawning.get(fish.getGender()).get("aLW") * Math.pow(fish.getLength(), fishFeaturesPreSpawning.get(fish.getGender()).get("bLW") ); else weight = fishFeaturesPostSpawning.get(fish.getGender()).get("aLW") * Math.pow(fish.getLength(), fishFeaturesPostSpawning.get(fish.getGender()).get("bLW")); return weight; } /** * Compute the weight for a fish with length (cm) * @param fish * @return weight (g) */ public double getWeight (DiadromousFish fish) { return getWeight (fish, SpawningPosition.PRE); } /** * @param args */ /** * @param args */ /** * @param args */ public static void main(String[] args) { double aResidenceTime =30; System.out.println("aResidenceTime: " + aResidenceTime); // Map <String, Double> anExcretionRate = new Hashtable <String, Double>(); anExcretionRate.put("N", 24.71E-6); //values from Barber et al, Alosa sapidissima in ug/g wet mass/hour : convertit en g anExcretionRate.put("P", 2.17E-6); //values from Barber et al, Alosa sapidissima in ug/g wet mass/hour: convertit en g System.out.println("anExcretionRate: " + anExcretionRate.toString()); // /* * A feature pre spawning */ Map<Gender, Map<String, Double>> aFeaturePreSpawning = new Hashtable<DiadromousFish.Gender, Map<String,Double>>();
351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420
/* * For females */ Map<String,Double> aFeature = new Hashtable<String,Double>(); aFeature.put("aLW", Math.exp(-4.9078)); //weight size relationship computed from BDalosesBruch aFeature.put("bLW", 3.147); //aFeature.put("bLW",3.3429);// parametre "b" de la relation taille/poids - Coefficient d'allometrie //aFeature.put("aLW",1.2102E-6 * Math.pow(10., aFeature.get("bLW"))); // parametre "a" de la relation taille/poids en kg/cm- Traduit la condition //aFeature.put("GSI",0.15); aFeaturePreSpawning.put(Gender.FEMALE, aFeature); /* * For males */ aFeature = new Hashtable<String,Double>(); aFeature.put("aLW", Math.exp(-1.304)); aFeature.put("bLW", 2.1774); //aFeature.put("aLW",2.4386E-6 * Math.pow(10, aFeature.get("bLW"))); // Conversion des g/mm en g.cm (from Taverny, 1991) //aFeature.put("GSI",.08); aFeaturePreSpawning.put(Gender.MALE,aFeature); System.out.println("aFeaturePreSpawning: " + aFeaturePreSpawning.toString()); // /* * a Feature post Spawning */ Map<Gender, Map<String, Double>> aFeaturePostSpawning = new Hashtable<DiadromousFish.Gender, Map<String,Double>>(); /* * For females */ aFeature = new Hashtable<String,Double>(); aFeature.put("aLW", Math.exp(-4.3276)); //weight size relationship computed from BDalosesBruch aFeature.put("bLW", 2.9418); //aFeature.put("GSI",0.10); //From BDalosesBruch //aFeature.put("aLW",aFeaturePreSpawning.get(Gender.FEMALE).get("aLW")/(1+aFeature.get("GSI"))); // parametre "a" de la relation taille/poids avec Lt en cm - Traduit la condition //aFeature.put("bLW",aFeaturePreSpawning.get(Gender.FEMALE).get("bLW"));// parametre "b" de la relation taille/poids - Coefficient d'allometrie aFeaturePostSpawning.put(Gender.FEMALE, aFeature); /* * For males */ aFeature = new Hashtable<String,Double>(); aFeature.put("aLW", Math.exp(-4.5675));// parametre "a" de la relation taille/poids - Coefficient d'allometrie aFeature.put("bLW", 2.9973); //aFeature.put("GSI",.05); From BDalosesBruch //aFeature.put("aLW",aFeaturePreSpawning.get(Gender.MALE).get("aLW")/(1+aFeature.get("GSI"))); //aFeature.put("bLW",aFeaturePreSpawning.get(Gender.MALE).get("bLW")); aFeaturePostSpawning.put(Gender.MALE,aFeature); System.out.println("aFeaturePostSpawning: " + aFeaturePostSpawning.toString()); Map<Gender, Double> aGameteSpawned = new Hashtable <DiadromousFish.Gender,Double>(); aGameteSpawned.put(Gender.FEMALE, 131.); // Compute from the difference between spawned and unspawned ovaries ie correspond to a mean weight of eggs spawned aGameteSpawned.put(Gender.MALE, 44.8); // Compute from the difference between spawned and unspawned testes ie correspond to a mean weight of sperm spawned System.out.println("aGameteSpawned: " + aGameteSpawned.toString()); // carcass composition for fish before spawning Map<Gender, Map<String, Double>> aCompoCarcassPreSpawning = new Hashtable<DiadromousFish.Gender,Map<String,Double>>(); Map<String,Double> aCompo = new Hashtable<String,Double>(); aCompo.put("N", 2.958 / 100.); //On remplit une collection avec un put. aCompo.put("P", 0.673 / 100.); aCompoCarcassPreSpawning.put(Gender.FEMALE,aCompo); aCompo = new Hashtable<String,Double>(); aCompo.put("N", 2.941 / 100.); aCompo.put("P", 0.666 / 100.);