dive4elements/river: artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/SedimentLoadCalculation.java annotate

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/SedimentLoadCalculation.java @ 8098:09725b65955a

Add new and simplyfied SedimentLoadFacet The SedimentLoadFacet is intended to work with the Measurement stations. It uses the same mechanismn to access the Mesurement station values as the calculation does. SedimentLoadLS values need a different facet that will come soon.

author	Andre Heinecke <andre.heinecke@intevation.de>
date	Fri, 15 Aug 2014 18:27:19 +0200
parents	fdb26fe898dc
children

rev	line source
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	1 /* Copyright (C) 2011, 2012, 2013 by Bundesanstalt für Gewässerkunde
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	2 * Software engineering by Intevation GmbH
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	3 *
5994 af13ceeba52a Removed trailing whitespace. Sascha L. Teichmann <teichmann@intevation.de> parents: 5863 diff changeset	4 * This file is Free Software under the GNU AGPL (>=v3)
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	5 * and comes with ABSOLUTELY NO WARRANTY! Check out the
5994 af13ceeba52a Removed trailing whitespace. Sascha L. Teichmann <teichmann@intevation.de> parents: 5863 diff changeset	6 * documentation coming with Dive4Elements River for details.
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	7 */
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	8
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 5788 diff changeset	9 package org.dive4elements.river.artifacts.model.minfo;
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	10
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	11 import gnu.trove.TDoubleArrayList;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	12
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	13 import java.util.ArrayList;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	14 import java.util.List;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	15
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	16 import org.apache.log4j.Logger;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	17
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 5788 diff changeset	18 import org.dive4elements.river.artifacts.access.SedimentLoadAccess;
bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 5788 diff changeset	19 import org.dive4elements.river.artifacts.model.Calculation;
bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 5788 diff changeset	20 import org.dive4elements.river.artifacts.model.CalculationResult;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	21 import org.dive4elements.river.artifacts.model.Range;
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	22
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	23
4483 cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 4373 diff changeset	24 /** Calculate sediment load. */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	25 public class SedimentLoadCalculation
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	26 extends Calculation
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	27 {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	28
4483 cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 4373 diff changeset	29 /** Private logger. */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	30 private static final Logger logger = Logger
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	31 .getLogger(SedimentLoadCalculation.class);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	32
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	33 protected String river;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	34 protected String yearEpoch;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	35 protected double kmUp;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	36 protected double kmLow;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	37 protected int[] period;
7532 aeeabb580b12 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7462 diff changeset	38 /** Years of chosen epochs.
aeeabb580b12 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7462 diff changeset	39 * epoch[0][0] is the start year of first entered epoch, e[0][1] the end-year. */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	40 protected int[][] epoch;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	41 protected String unit;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	42
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	43 public SedimentLoadCalculation() {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	44 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	45
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	46 /** Returns CalculationResult with array of SedimentLoadResults. */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	47 public CalculationResult calculate(SedimentLoadAccess access) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	48 logger.info("SedimentLoadCalculation.calculate");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	49
7261 a56fe3bc6700 Refactoring: Let RiverAccess.getRiver return an River. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7254 diff changeset	50 String river = access.getRiverName();
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	51 String yearEpoch = access.getYearEpoch();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	52 String unit = access.getUnit();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	53 int[] period = null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	54 int[][] epoch = null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	55 double kmUp = access.getUpperKM();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	56 double kmLow = access.getLowerKM();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	57 if (yearEpoch.equals("year")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	58 period = access.getPeriod();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	59 epoch = null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	60 }
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	61 else if (yearEpoch.equals("epoch") \|\| yearEpoch.equals("off_epoch")) {
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	62 epoch = access.getEpochs();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	63 period = null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	64 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	65 else {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	66 addProblem("minfo.missing.year_epoch");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	67 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	68
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	69 if (river == null) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	70 // TODO: i18n
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	71 addProblem("minfo.missing.river");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	72 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	73
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	74 if (period == null && epoch == null) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	75 addProblem("minfo.missing.time");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	76 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	77
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	78 if (!hasProblems()) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	79 this.river = river;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	80 this.yearEpoch = yearEpoch;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	81 this.unit = unit;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	82 this.period = period;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	83 this.epoch = epoch;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	84 this.kmUp = kmUp;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	85 this.kmLow = kmLow;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	86 return internalCalculate();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	87 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	88
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	89 return new CalculationResult();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	90 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	91
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	92 /** Returns CalculationResult with array of SedimentLoadResults. */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	93 private CalculationResult internalCalculate() {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	94 logger.debug("internalCalulate; mode:" + yearEpoch);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	95 if (yearEpoch.equals("year")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	96 List<SedimentLoadResult> results =
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	97 new ArrayList<SedimentLoadResult>();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	98 for (int i = period[0]; i <= period[1]; i++) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	99 SedimentLoadResult res = calculateYear(i);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	100 results.add(res);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	101 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	102 return new CalculationResult(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	103 results.toArray(new SedimentLoadResult[results.size()]), this);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	104 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	105 else if (yearEpoch.equals("epoch")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	106 List<SedimentLoadResult> results =
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	107 new ArrayList<SedimentLoadResult>();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	108 for (int i = 0; i < epoch.length; i++) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	109 SedimentLoadResult res = calculateEpoch(i);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	110 results.add(res);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	111 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	112 return new CalculationResult(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	113 results.toArray(new SedimentLoadResult[results.size()]), this);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	114 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	115 else if (yearEpoch.equals("off_epoch")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	116 List<SedimentLoadResult> results =
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	117 new ArrayList<SedimentLoadResult>();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	118 for (int i = 0; i < epoch.length; i++) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	119 SedimentLoadResult res = calculateOffEpoch(i);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	120 results.add(res);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	121 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	122 return new CalculationResult(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	123 results.toArray(new SedimentLoadResult[results.size()]), this);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	124 }
6950 40da4a6e2cfb SedimentLoadCalculation: Log an error on unknown mode. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6948 diff changeset	125 else {
40da4a6e2cfb SedimentLoadCalculation: Log an error on unknown mode. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6948 diff changeset	126 logger.error("Unknown mode " + yearEpoch);
40da4a6e2cfb SedimentLoadCalculation: Log an error on unknown mode. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6948 diff changeset	127 }
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	128 return null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	129 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	130
7054 f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	131
f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	132 /** Returns val if not NaN, 0d otherwise. */
f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	133 private static double makeNaN0(double val) {
f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	134 return Double.isNaN(val)
f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	135 ? 0d
f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	136 : val;
f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	137 }
f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	138
6948 d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	139 /**
7031 7ee382b3f839 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7023 diff changeset	140 * Take Loads and build average of all fractions at given km.
7ee382b3f839 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7023 diff changeset	141 * The average fractions value is set in resLoad.
7ee382b3f839 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7023 diff changeset	142 *
7ee382b3f839 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7023 diff changeset	143 * @param km km at which to build fractions average.
7ee382b3f839 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7023 diff changeset	144 * @param epochLoads the loads to build average over.
7ee382b3f839 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7023 diff changeset	145 * @param[out] resLoad resulting SedimentLoad.
6948 d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	146 */
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	147 private void calculateEpochKm(
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	148 List<SedimentLoadLSData> epochLoads,
963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	149 SedimentLoadLSData resLoad,
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	150 double km
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	151 ) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	152 int cSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	153 int fmSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	154 int sSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	155 int ssSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	156 int ssbSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	157 int sseSum = 0;
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	158 for (SedimentLoadLSData load : epochLoads) {
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	159 SedimentLoadFraction f = load.getFraction(km);
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	160 if (f.getCoarse() > 0d) {
7054 f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	161 double c = makeNaN0(resLoad.getFraction(km).getCoarse());
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	162 resLoad.setCoarse(km, c + f.getCoarse(), f.getCoarseRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	163 cSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	164 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	165 if (f.getFineMiddle() > 0d) {
7054 f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	166 double fm = makeNaN0(resLoad.getFraction(km).getFineMiddle());
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	167 resLoad.setFineMiddle(km, fm + f.getFineMiddle(), f.getFineMiddleRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	168 fmSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	169 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	170 if (f.getSand() > 0d) {
7054 f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	171 double s = makeNaN0(resLoad.getFraction(km).getSand());
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	172 resLoad.setSand(km, s + f.getSand(), f.getSandRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	173 sSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	174 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	175 if (f.getSuspSand() > 0d) {
7054 f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	176 double s = makeNaN0(resLoad.getFraction(km).getSuspSand());
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	177 resLoad.setSuspSand(km, s + f.getSuspSand(), f.getSuspSandRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	178 ssSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	179 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	180 if (f.getSuspSandBed() > 0d) {
7054 f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	181 double s = makeNaN0(resLoad.getFraction(km).getSuspSandBed());
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	182 resLoad.setSuspSandBed(km, s + f.getSuspSandBed(), f.getSuspSandBedRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	183 ssbSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	184 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	185 if (f.getSuspSediment() > 0d) {
7054 f4d6ef50665a issue1436: Repair sedimentload epoch calculation, Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7031 diff changeset	186 double s = makeNaN0(resLoad.getFraction(km).getSuspSediment());
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	187 resLoad.setSuspSediment(km, s + f.getSuspSediment(), f.getSuspSedimentRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	188 sseSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	189 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	190 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	191
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	192 SedimentLoadFraction fr = resLoad.getFraction(km);
7031 7ee382b3f839 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7023 diff changeset	193 // Prevent divisions by zero.
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	194 if (cSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	195 resLoad.setCoarse(km, fr.getCoarse()/cSum, fr.getCoarseRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	196 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	197 if (fmSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	198 resLoad.setFineMiddle(km, fr.getFineMiddle()/fmSum,
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	199 fr.getFineMiddleRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	200 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	201 if (sSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	202 resLoad.setSand(km, fr.getSand()/sSum, fr.getSandRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	203 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	204 if (ssSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	205 resLoad.setSuspSand(km, fr.getSuspSand()/ssSum,
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	206 fr.getSuspSandRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	207 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	208 if (ssbSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	209 resLoad.setSuspSandBed(km, fr.getSuspSandBed()/ssbSum,
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	210 fr.getSuspSandBedRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	211 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	212 if (sseSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	213 resLoad.setSuspSediment(km, fr.getSuspSediment()/sseSum, fr.getSuspSedimentRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	214 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	215 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	216
7532 aeeabb580b12 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7462 diff changeset	217 /**
aeeabb580b12 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7462 diff changeset	218 * Calculate result for the ith given epoch.
aeeabb580b12 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7462 diff changeset	219 * @param i index of epoch (if multiple given).
aeeabb580b12 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7462 diff changeset	220 */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	221 private SedimentLoadResult calculateEpoch(int i) {
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	222 List<SedimentLoadLSData> epochLoads = new ArrayList<SedimentLoadLSData>();
7931 b4117b8adb79 Include last year of epoch in calculation. "Tom Gottfried <tom@intevation.de>" parents: 7532 diff changeset	223 for (int j = epoch[i][0]; j <= epoch[i][1]; j++) {
5661 b6297a67823e SedimentLoadFactory: Renamed getLoadwithData to getLoadWithData, update callers. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 4543 diff changeset	224 epochLoads.add(SedimentLoadFactory.getLoadWithData(
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	225 this.river,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	226 this.yearEpoch,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	227 this.kmLow,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	228 this.kmUp,
7064 fabb5e7a3b66 Doc, comment. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7054 diff changeset	229 j, //syear
fabb5e7a3b66 Doc, comment. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7054 diff changeset	230 j)); //eyear
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	231 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	232
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	233 SedimentLoadLSData resLoad = new SedimentLoadLSData();
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	234 TDoubleArrayList kms = new TDoubleArrayList();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	235
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	236 for (SedimentLoadLSData load : epochLoads) {
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	237 for (double km : load.getKms()) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	238 if (!kms.contains(km)) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	239 kms.add(km);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	240 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	241 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	242 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	243
7254 fe32a7f9655e Reduce usage of careless 'for (int i = 0; i < x.size(); i++)' pattern. Sascha L. Teichmann <teichmann@intevation.de> parents: 7065 diff changeset	244 for (int j = 0, J = kms.size(); j < J; j++) {
6928 2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6740 diff changeset	245 calculateEpochKm(epochLoads, resLoad, kms.get(j));
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	246 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	247 resLoad.setDescription("");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	248 resLoad.setEpoch(true);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	249
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	250 SedimentLoadLSData sl = calculateTotalLoad(resLoad, this.epoch[i][0]);
7065 b17634649e61 Simplify code a bit. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7064 diff changeset	251
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	252 if (this.unit.equals("m3_per_a")) {
7065 b17634649e61 Simplify code a bit. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7064 diff changeset	253 sl = calculateUnit(sl, this.epoch[i][0]);
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	254 }
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	255
7065 b17634649e61 Simplify code a bit. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7064 diff changeset	256 return new SedimentLoadResult(
b17634649e61 Simplify code a bit. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7064 diff changeset	257 this.epoch[i][0],
b17634649e61 Simplify code a bit. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7064 diff changeset	258 this.epoch[i][1],
b17634649e61 Simplify code a bit. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7064 diff changeset	259 sl);
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	260 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	261
6948 d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	262 /**
d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	263 * Calculate/Fetch values at off. epochs.
d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	264 * @param i index in epochs.
d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	265 */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	266 private SedimentLoadResult calculateOffEpoch(int i) {
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	267 SedimentLoadLSData load = SedimentLoadFactory.getLoadWithData(
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	268 this.river,
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	269 this.yearEpoch,
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	270 kmLow,
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	271 kmUp,
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	272 this.epoch[i][0],
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	273 this.epoch[i][1]);
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	274 SedimentLoadResult result;
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	275 SedimentLoadLSData sl = calculateTotalLoad(load, this.epoch[i][0]);
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	276 if (unit.equals("m3_per_a")) {
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	277 SedimentLoadLSData slu = calculateUnit(sl, epoch[i][0]);
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	278 result = new SedimentLoadResult(
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	279 this.epoch[i][0],
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	280 this.epoch[i][1],
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	281 slu);
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	282 }
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	283 else {
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	284 result = new SedimentLoadResult(
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	285 this.epoch[i][0],
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	286 this.epoch[i][1],
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	287 sl);
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	288 }
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	289
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	290 return result;
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	291 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	292
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	293 /**
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	294 * Fetch loads for a single year, calculate total and
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	295 * return the result containing both.
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	296 * @param y year, e.g. 1980
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	297 */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	298 private SedimentLoadResult calculateYear(int y) {
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	299 SedimentLoadLSData load = SedimentLoadFactory.getLoadWithData(
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	300 this.river,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	301 this.yearEpoch,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	302 this.kmLow,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	303 this.kmUp,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	304 y,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	305 y);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	306
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	307 SedimentLoadResult result;
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	308 SedimentLoadLSData sl = calculateTotalLoad(load, y);
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	309 if (unit.equals("m3_per_a")) {
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	310 SedimentLoadLSData slu = calculateUnit(sl, y);
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	311 result = new SedimentLoadResult(y, 0, slu);
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	312 }
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	313 else {
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	314 result = new SedimentLoadResult(y, 0, sl);
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	315 }
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	316 return result;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	317 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	318
6374 48e92ff57f23 SedimentLoad: Set range to fraction/fractions values. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 5994 diff changeset	319 /** Add up the loads of a year. */
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	320 private SedimentLoadLSData calculateTotalLoad(SedimentLoadLSData load, int year) {
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	321 logger.debug("calculateTotalLoad");
5685 756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible). Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5661 diff changeset	322 boolean problemThisYear = false;
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	323 if (!load.hasCoarse()) {
5742 a4a894b15c35 Minor part of issue1235: dont show decimals of years. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5685 diff changeset	324 addProblem("missing.fraction.coarse", Integer.toString(year));
5685 756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible). Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5661 diff changeset	325 problemThisYear = true;
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	326 }
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	327 if (!load.hasFineMiddle()) {
5742 a4a894b15c35 Minor part of issue1235: dont show decimals of years. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5685 diff changeset	328 addProblem("missing.fraction.fine_middle", Integer.toString(year));
5685 756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible). Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5661 diff changeset	329 problemThisYear = true;
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	330 }
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	331 if (!load.hasSand()) {
5742 a4a894b15c35 Minor part of issue1235: dont show decimals of years. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5685 diff changeset	332 addProblem("missing.fraction.sand", Integer.toString(year));
5685 756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible). Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5661 diff changeset	333 problemThisYear = true;
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	334 }
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	335 if (!load.hasSuspSand()) {
5742 a4a894b15c35 Minor part of issue1235: dont show decimals of years. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5685 diff changeset	336 addProblem("missing.fraction.susp_sand", Integer.toString(year));
5685 756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible). Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5661 diff changeset	337 problemThisYear = true;
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	338 }
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	339 if (!load.hasSuspSediment()) {
5742 a4a894b15c35 Minor part of issue1235: dont show decimals of years. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5685 diff changeset	340 addProblem("missing.fraction.susp_sediment", Integer.toString(year));
5685 756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible). Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5661 diff changeset	341 problemThisYear = true;
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	342 }
5685 756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible). Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5661 diff changeset	343 if (problemThisYear) {
756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible). Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 5661 diff changeset	344 logger.warn("Some problem, not calculating total load.");
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	345 return load;
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	346 }
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	347 return partialTotal(load);
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	348 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	349
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	350
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	351
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	352 /**
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	353 * Set total values in load.
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	354 *
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	355 * Therefore, run over the sorted kms and find ranges where either all
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	356 * or all Geschiebe or just the Schwebstoff fractions are set.
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	357 * Merge these ranges and add (maybe new) respective fractions to
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	358 * load. In the process, remember any 'unfished' ends from measurements
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	359 * where the km-ranges did not completely match.
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	360 *
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	361 * @param load SedimentLoad to add total values (and ranges) to.
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	362 * @return input param load, with total values set.
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	363 */
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	364 private SedimentLoadLSData partialTotal(SedimentLoadLSData load) {
7022 0f99aa385d99 Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6989 diff changeset	365 // The load with balanced ranges, will be returned.
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	366 SedimentLoadLSData fairLoad = load;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	367
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	368 Range lastOtherRange = null;
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	369 double lastOtherValue = 0d;
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	370
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	371 Range lastSuspRange = null;
e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	372 double lastSuspValue = 0d;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	373
6970 7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	374 for (double km: load.getKms()) { // kms are already sorted!
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	375 logger.debug ("Trying to add at km " + km);
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	376 SedimentLoadFraction fraction = load.getFraction(km);
6970 7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	377 if (fraction.isComplete()) {
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	378 double total = fraction.getCoarse() +
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	379 fraction.getFineMiddle() +
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	380 fraction.getSand() +
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	381 fraction.getSuspSand() +
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	382 fraction.getSuspSediment();
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	383 // Easiest case. Add values up and set'em.
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	384 if (fraction.getCoarseRange().equals(
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	385 fraction.getSuspSedimentRange())) {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	386 lastOtherRange = null;
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	387 lastSuspRange = null;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	388 fairLoad.setTotal(km, total, fraction.getCoarseRange());
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	389 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	390 else {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	391 // Need to split a range.
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	392 if (fraction.getCoarseRange().getEnd()
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	393 < fraction.getSuspSedimentRange().getEnd()) {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	394 // Schwebstoff is longer.
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	395 // Adjust and remember schwebstoffs range and value.
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	396 lastSuspRange = (Range) fraction.getSuspSedimentRange().clone();
e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	397 lastSuspRange.setStart(fraction.getCoarseRange().getEnd());
6411 a21f7a3e4be5 Revert to rev 6404 for issue1235 Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6405 diff changeset	398 lastSuspValue = fraction.getSuspSediment();
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	399 lastOtherRange = null;
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	400 fairLoad.setTotal(km, total, fraction.getCoarseRange());
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	401 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	402 else {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	403 // Geschiebe is longer.
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	404 // Adjust and remember other values.
7023 a18a43764793 issue1235: Do not add 'other' values twice to suspended sed. values, fix range. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7022 diff changeset	405 lastOtherRange = (Range) fraction.getCoarseRange().clone();
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	406 lastOtherRange.setStart(fraction.getSuspSedimentRange().getEnd());
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	407 lastOtherValue = (total - fraction.getSuspSediment());
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	408 lastSuspRange = null;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	409 fairLoad.setTotal(km, total, fraction.getSuspSedimentRange());
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	410 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	411 }
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	412 }
6970 7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	413 else if (fraction.hasOnlySuspValues() && lastOtherRange != null) {
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	414 // Split stuff.
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	415 Range suspSedimentRange = fraction.getSuspSedimentRange();
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	416 // if intersects with last other range, cool! merge and add!
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	417 if (lastOtherRange.contains(km)) {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	418 double maxStart = 0d;
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	419 double minEnd = 0d;
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	420 maxStart = Math.max(suspSedimentRange.getStart(),
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	421 lastOtherRange.getStart());
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	422
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	423 minEnd = Math.min(suspSedimentRange.getEnd(),
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	424 lastOtherRange.getEnd());
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	425 double total = lastOtherValue + fraction.getSuspSediment();
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	426 Range totalRange = new Range(maxStart, minEnd);
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	427 if (suspSedimentRange.getEnd() > lastOtherRange.getEnd()) {
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	428 lastSuspRange = (Range) suspSedimentRange.clone();
e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	429 lastSuspRange.setStart(lastOtherRange.getEnd());
6411 a21f7a3e4be5 Revert to rev 6404 for issue1235 Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6405 diff changeset	430 lastSuspValue = fraction.getSuspSediment();
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	431 lastOtherRange = null;
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	432 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	433 else {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	434 // Other is "longer".
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	435 lastOtherRange.setStart(suspSedimentRange.getEnd());
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	436 lastSuspRange = null;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	437 }
6970 7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	438 if (lastOtherRange != null
7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	439 && Math.abs(suspSedimentRange.getEnd() - lastOtherRange.getEnd()) < 0.1d) {
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	440 lastOtherRange = null;
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	441 lastSuspRange = null;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	442 }
7023 a18a43764793 issue1235: Do not add 'other' values twice to suspended sed. values, fix range. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 7022 diff changeset	443 fairLoad.setTotal(km, total, totalRange);
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	444 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	445 else {
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	446 lastSuspRange = suspSedimentRange;
6411 a21f7a3e4be5 Revert to rev 6404 for issue1235 Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6405 diff changeset	447 lastSuspValue = fraction.getSuspSediment();
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	448 lastOtherRange = null;
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	449 }
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	450 }
6970 7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	451 else if (fraction.hasButSuspValues() && lastSuspRange != null) {
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	452 // If intersects with last suspsed range, merge and add
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	453 double total = fraction.getCoarse() +
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	454 fraction.getFineMiddle() +
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	455 fraction.getSand() +
6411 a21f7a3e4be5 Revert to rev 6404 for issue1235 Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6405 diff changeset	456 fraction.getSuspSand() +
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	457 lastSuspValue;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	458 double maxStart = Math.max(fraction.getCoarseRange().getStart(),
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	459 lastSuspRange.getStart());
e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	460 if (lastSuspRange.contains(km)) {
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	461 double minEnd = Math.min(fraction.getCoarseRange().getEnd(),
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	462 lastSuspRange.getEnd());
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	463 Range totalRange = new Range(maxStart, minEnd);
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	464 if (lastSuspRange.getEnd() > fraction.getCoarseRange().getEnd()) {
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	465 // SuspSed longer.
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	466 lastSuspRange.setStart(fraction.getCoarseRange().getEnd());
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	467 lastOtherRange = null;
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	468 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	469 else {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	470 // Other longer
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	471 lastOtherRange = (Range) fraction.getCoarseRange().clone();
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	472 lastOtherRange.setStart(lastSuspRange.getEnd());
e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	473 lastSuspRange = null;
e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	474 lastOtherValue = total - lastSuspValue;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	475 }
6970 7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	476 if (lastSuspRange != null
7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	477 && lastOtherRange != null
7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	478 && Math.abs(lastSuspRange.getEnd() - lastOtherRange.getEnd()) < 0.1d) {
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	479 lastOtherRange = null;
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	480 lastSuspRange = null;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	481 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	482 fairLoad.setTotal(km, total, totalRange);
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	483 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	484 else {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	485 // Ranges are disjoint.
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	486 lastOtherRange = fraction.getCoarseRange();
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	487 lastOtherValue = total - fraction.getSuspSediment();
6404 e6ceab9e3091 prep issue1235: Rename functions and variable to indicate use of other susp values. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6381 diff changeset	488 lastSuspRange = null;
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	489 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	490 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	491 else {
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	492 // Some values are missing or no intersection with former values.
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	493 // Stay as we are.
6970 7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	494 if (fraction.hasButSuspValues()) {
6500 f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	495 double total = fraction.getCoarse() +
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	496 fraction.getFineMiddle() +
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	497 fraction.getSand() +
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	498 fraction.getSuspSand();
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	499 lastOtherRange = fraction.getCoarseRange();
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	500 lastOtherValue = total;
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	501 lastSuspRange = null;
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	502 }
6970 7be97faf5848 flys/issue1235: Same kicks against a few inconsistencies and bugs in the calculation of sediment loads. I believe it _do_not_ delivers the right results. Sascha L. Teichmann <teichmann@intevation.de> parents: 6953 diff changeset	503 else if (fraction.hasOnlySuspValues()) {
6500 f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	504 lastSuspRange = fraction.getSuspSedimentRange();
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	505 lastSuspValue = fraction.getSuspSediment();
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	506 lastOtherRange = null;
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6499 diff changeset	507 }
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	508 }
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	509 }
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	510 return fairLoad;
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	511 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	512
6381 df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6377 diff changeset	513
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	514 /**
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	515 * Transform values in load.
6948 d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	516 * Background is to transform values measured in
6951 f5f60b3d9180 Doc fix. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6950 diff changeset	517 * t/a to m^3/a using the specific measured densities.
6948 d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	518 *
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	519 * @param load The load of which values should be transformed.
6948 d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	520 * @param year The year at which to look at density (e.g. 2003).
d4108d6c4000 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6928 diff changeset	521 *
6678 0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	522 * @return parameter load with transformed values.
0a70a320bfca Doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6543 diff changeset	523 */
8024 963ede7b32bb Renamed SedimentLoad to SedimentLoadLSData to make place for SedimentLoad in backend. Sascha L. Teichmann <teichmann@intevation.de> parents: 7931 diff changeset	524 private SedimentLoadLSData calculateUnit(SedimentLoadLSData load, int year) {
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	525 SedimentDensity density =
8065 fdb26fe898dc Sediment load: Cached sediment values dont depend on year. The cached object is an index over year so there is not need to ask the factory for this. Sascha L. Teichmann <teichmann@intevation.de> parents: 8024 diff changeset	526 SedimentDensityFactory.getSedimentDensity(river, kmLow, kmUp);
6951 f5f60b3d9180 Doc fix. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6950 diff changeset	527
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	528 for (double km: load.getKms()) {
7462 906ed0b1f3f1 Fix for flys/issue1654: Divide through denisity, do not multiply. Sascha L. Teichmann <teichmann@intevation.de> parents: 7261 diff changeset	529 double dens = 1d/density.getDensity(km, year);
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	530 SedimentLoadFraction fraction = load.getFraction(km);
a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	531 double coarse = fraction.getCoarse();
6376 8d6f23243d1d SedimentLoad: Switch from snake_case to camelCase. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 6375 diff changeset	532 double fineMiddle = fraction.getFineMiddle();
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	533 double sand = fraction.getSand();
6376 8d6f23243d1d SedimentLoad: Switch from snake_case to camelCase. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 6375 diff changeset	534 double suspSand = fraction.getSuspSand();
8d6f23243d1d SedimentLoad: Switch from snake_case to camelCase. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 6375 diff changeset	535 double bedSand = fraction.getSuspSandBed();
8d6f23243d1d SedimentLoad: Switch from snake_case to camelCase. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 6375 diff changeset	536 double sediment = fraction.getSuspSediment();
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	537 double total = fraction.getTotal();
6953 733c35239fe7 issue1462/3: When changing units (t/a -> m^3/a), also change total load. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6951 diff changeset	538 double loadTotal = fraction.getLoadTotal();
6374 48e92ff57f23 SedimentLoad: Set range to fraction/fractions values. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 5994 diff changeset	539 load.setCoarse(km, (coarse * dens), fraction.getCoarseRange());
48e92ff57f23 SedimentLoad: Set range to fraction/fractions values. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 5994 diff changeset	540 load.setFineMiddle(km, (fineMiddle * dens), fraction.getFineMiddleRange());
48e92ff57f23 SedimentLoad: Set range to fraction/fractions values. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 5994 diff changeset	541 load.setSand(km, (sand * dens), fraction.getSandRange());
48e92ff57f23 SedimentLoad: Set range to fraction/fractions values. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 5994 diff changeset	542 load.setSuspSand(km, (suspSand * dens), fraction.getSuspSandRange());
48e92ff57f23 SedimentLoad: Set range to fraction/fractions values. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 5994 diff changeset	543 load.setSuspSandBed(km, (bedSand * dens), fraction.getSuspSandBedRange());
48e92ff57f23 SedimentLoad: Set range to fraction/fractions values. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 5994 diff changeset	544 load.setSuspSediment(km, (sediment * dens), fraction.getSuspSedimentRange());
48e92ff57f23 SedimentLoad: Set range to fraction/fractions values. Felix Wolfsteller <felix.wolfsteller@intevation.de>* parents: 5994 diff changeset	545 load.setTotal(km, (total * dens), fraction.getTotalRange());
6953 733c35239fe7 issue1462/3: When changing units (t/a -> m^3/a), also change total load. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 6951 diff changeset	546 load.setLoadTotal(km, (loadTotal * dens), fraction.getLoadTotalRange());
4522 a9753f717b3d Completed and fixed sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: 4483 diff changeset	547 }
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	548 return load;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	549 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	550 }
4483 cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 4373 diff changeset	551 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf-8 :

Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/SedimentLoadCalculation.java @ 8098:09725b65955a