An error occurred while loading the file. Please try again.
-
Gaetano Raffaele authored5285094a
/**
*
* @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.);