Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/SedimentLoadCalculation.java @ 6970:7be97faf5848

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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.
author Sascha L. Teichmann <teichmann@intevation.de>
date Thu, 05 Sep 2013 17:15:04 +0200
parents 733c35239fe7
children 2b022ca95b3b
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;
6498
686bbfeacfe5 issue1235: sort km when running over them for total sediment load.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6411
diff changeset
14 import java.util.TreeSet;
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
15 import java.util.List;
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
16
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
17 import org.apache.log4j.Logger;
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
18
5831
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5788
diff changeset
19 import org.dive4elements.river.artifacts.access.SedimentLoadAccess;
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5788
diff changeset
20 import org.dive4elements.river.artifacts.model.Calculation;
bd047b71ab37 Repaired internal references
Sascha L. Teichmann <teichmann@intevation.de>
parents: 5788
diff changeset
21 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
22 import org.dive4elements.river.artifacts.model.Range;
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
23
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
24
4483
cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4373
diff changeset
25 /** Calculate sediment load. */
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
26 public class SedimentLoadCalculation
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
27 extends Calculation
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
28 {
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
29
4483
cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4373
diff changeset
30 /** Private logger. */
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
31 private static final Logger logger = Logger
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
32 .getLogger(SedimentLoadCalculation.class);
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
33
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
34 protected String river;
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
35 protected String yearEpoch;
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
36 protected double kmUp;
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
37 protected double kmLow;
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
38 protected int[] period;
6740
68c7b5811c4f Cosmetics, docs, in range of issue1393.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6678
diff changeset
39 /** Years of chosen epochs. */
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
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
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
50 String river = access.getRiver();
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
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
6948
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
131 /**
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
132 * @param[out] resLoad resulting SedimentLoad
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
133 */
6928
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
134 private void calculateEpochKm(
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
135 List<SedimentLoad> epochLoads,
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
136 SedimentLoad resLoad,
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
137 double km
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
138 ) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
139 int cSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
140 int fmSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
141 int sSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
142 int ssSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
143 int ssbSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
144 int sseSum = 0;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
145 for (SedimentLoad load : epochLoads) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
146 SedimentLoadFraction f = load.getFraction(km);
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
147 if (f.getCoarse() > 0d) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
148 double c = resLoad.getFraction(km).getCoarse();
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
149 resLoad.setCoarse(km, c + f.getCoarse(), f.getCoarseRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
150 cSum++;
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 if (f.getFineMiddle() > 0d) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
153 double fm = resLoad.getFraction(km).getFineMiddle();
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
154 resLoad.setFineMiddle(km, fm + f.getFineMiddle(), f.getFineMiddleRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
155 fmSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
156 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
157 if (f.getSand() > 0d) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
158 double s = resLoad.getFraction(km).getSand();
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
159 resLoad.setSand(km, s + f.getSand(), f.getSandRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
160 sSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
161 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
162 if (f.getSuspSand() > 0d) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
163 double s = resLoad.getFraction(km).getSuspSand();
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
164 resLoad.setSuspSand(km, s + f.getSuspSand(), f.getSuspSandRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
165 ssSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
166 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
167 if (f.getSuspSandBed() > 0d) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
168 double s = resLoad.getFraction(km).getSuspSandBed();
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
169 resLoad.setSuspSandBed(km, s + f.getSuspSandBed(), f.getSuspSandBedRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
170 ssbSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
171 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
172 if (f.getSuspSediment() > 0d) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
173 double s = resLoad.getFraction(km).getSuspSediment();
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
174 resLoad.setSuspSediment(km, s + f.getSuspSediment(), f.getSuspSedimentRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
175 sseSum++;
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
176 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
177 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
178
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
179 SedimentLoadFraction fr = resLoad.getFraction(km);
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
180 // Prevent divisions by zero, the fraction defaults to 0d.
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
181 if (cSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
182 resLoad.setCoarse(km, fr.getCoarse()/cSum, fr.getCoarseRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
183 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
184 if (fmSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
185 resLoad.setFineMiddle(km, fr.getFineMiddle()/fmSum,
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
186 fr.getFineMiddleRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
187 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
188 if (sSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
189 resLoad.setSand(km, fr.getSand()/sSum, fr.getSandRange());
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 if (ssSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
192 resLoad.setSuspSand(km, fr.getSuspSand()/ssSum,
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
193 fr.getSuspSandRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
194 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
195 if (ssbSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
196 resLoad.setSuspSandBed(km, fr.getSuspSandBed()/ssbSum,
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
197 fr.getSuspSandBedRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
198 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
199 if (sseSum != 0) {
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
200 resLoad.setSuspSediment(km, fr.getSuspSediment()/sseSum, fr.getSuspSedimentRange());
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
201 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
202 }
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
203
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
204 private SedimentLoadResult calculateEpoch(int i) {
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
205 List<SedimentLoad> epochLoads = new ArrayList<SedimentLoad>();
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
206 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
207 epochLoads.add(SedimentLoadFactory.getLoadWithData(
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
208 this.river,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
209 this.yearEpoch,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
210 this.kmLow,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
211 this.kmUp,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
212 j,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
213 j));
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
214 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
215
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
216 SedimentLoad resLoad = new SedimentLoad();
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
217 TDoubleArrayList kms = new TDoubleArrayList();
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
218
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
219 for (SedimentLoad load : epochLoads) {
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
220 for (double km : load.getKms()) {
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
221 if (!kms.contains(km)) {
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
222 kms.add(km);
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
223 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
224 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
225 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
226
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
227 for (int j = 0; j < kms.size(); j++) {
6928
2119f13f95ab SedimentLoadCalculation: Refactoring, extracted method.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6740
diff changeset
228 calculateEpochKm(epochLoads, resLoad, kms.get(j));
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
229 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
230 resLoad.setDescription("");
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
231 resLoad.setEpoch(true);
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
232
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
233 SedimentLoadResult result;
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
234 SedimentLoad sl = calculateTotalLoad(resLoad, this.epoch[i][0]);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
235 if (this.unit.equals("m3_per_a")) {
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
236 SedimentLoad slu = calculateUnit(sl, this.epoch[i][0]);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
237 result = new SedimentLoadResult(
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
238 this.epoch[i][0],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
239 this.epoch[i][1],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
240 slu);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
241 }
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
242 else {
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
243 result = new SedimentLoadResult(
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
244 this.epoch[i][0],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
245 this.epoch[i][1],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
246 sl);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
247 }
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
248
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
249 return result;
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
250 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
251
6948
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
252 /**
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
253 * Calculate/Fetch values at off. epochs.
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
254 * @param i index in epochs.
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
255 */
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
256 private SedimentLoadResult calculateOffEpoch(int i) {
5661
b6297a67823e SedimentLoadFactory: Renamed getLoadwithData to getLoadWithData, update callers.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4543
diff changeset
257 SedimentLoad load = SedimentLoadFactory.getLoadWithData(
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
258 this.river,
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
259 this.yearEpoch,
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
260 kmLow,
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
261 kmUp,
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
262 this.epoch[i][0],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
263 this.epoch[i][1]);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
264 SedimentLoadResult result;
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
265 SedimentLoad sl = calculateTotalLoad(load, this.epoch[i][0]);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
266 if (unit.equals("m3_per_a")) {
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
267 SedimentLoad slu = calculateUnit(sl, epoch[i][0]);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
268 result = new SedimentLoadResult(
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
269 this.epoch[i][0],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
270 this.epoch[i][1],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
271 slu);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
272 }
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
273 else {
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
274 result = new SedimentLoadResult(
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
275 this.epoch[i][0],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
276 this.epoch[i][1],
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
277 sl);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
278 }
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
279
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
280 return result;
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
281 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
282
6678
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
283 /**
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
284 * Fetch loads for a single year, calculate total and
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
285 * return the result containing both.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
286 * @param y year, e.g. 1980
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
287 */
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
288 private SedimentLoadResult calculateYear(int y) {
5661
b6297a67823e SedimentLoadFactory: Renamed getLoadwithData to getLoadWithData, update callers.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4543
diff changeset
289 SedimentLoad load = SedimentLoadFactory.getLoadWithData(
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
290 this.river,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
291 this.yearEpoch,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
292 this.kmLow,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
293 this.kmUp,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
294 y,
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
295 y);
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
296
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
297 SedimentLoadResult result;
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
298 SedimentLoad sl = calculateTotalLoad(load, y);
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
299 if (unit.equals("m3_per_a")) {
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
300 SedimentLoad slu = calculateUnit(sl, y);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
301 result = new SedimentLoadResult(y, 0, slu);
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
302 }
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
303 else {
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
304 result = new SedimentLoadResult(y, 0, sl);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
305 }
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
306 return result;
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
307 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
308
6374
48e92ff57f23 SedimentLoad*: Set range to fraction/fractions values.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5994
diff changeset
309 /** Add up the loads of a year. */
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
310 private SedimentLoad calculateTotalLoad(SedimentLoad load, int year) {
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
311 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
312 boolean problemThisYear = false;
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
313 if (!load.hasCoarse()) {
5742
a4a894b15c35 Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5685
diff changeset
314 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
315 problemThisYear = true;
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
316 }
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
317 if (!load.hasFineMiddle()) {
5742
a4a894b15c35 Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5685
diff changeset
318 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
319 problemThisYear = true;
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
320 }
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
321 if (!load.hasSand()) {
5742
a4a894b15c35 Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5685
diff changeset
322 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
323 problemThisYear = true;
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
324 }
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
325 if (!load.hasSuspSand()) {
5742
a4a894b15c35 Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5685
diff changeset
326 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
327 problemThisYear = true;
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
328 }
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
329 if (!load.hasSuspSediment()) {
5742
a4a894b15c35 Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5685
diff changeset
330 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
331 problemThisYear = true;
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
332 }
5685
756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible).
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5661
diff changeset
333 if (problemThisYear) {
756df79274e1 SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible).
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5661
diff changeset
334 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
335 return load;
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
336 }
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
337 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
338 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
339
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
340
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
341
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
342 /**
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
343 * Set total values in load.
6678
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
344 *
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
345 * 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
346 * 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
347 * Merge these ranges and add (maybe new) respective fractions to
6678
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
348 * load. In the process, remember any 'unfished' ends from measurements
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
349 * where the km-ranges did not completely match.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
350 *
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
351 * @param load SedimentLoad to add total values (and ranges) to.
6678
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
352 * @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
353 */
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
354 private SedimentLoad partialTotal(SedimentLoad 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
355 SedimentLoad fairLoad = 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
356
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 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
358 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
359
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
360 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
361 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
362
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
363 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
364 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
365 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
366 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
367 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
368 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
369 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
370 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
371 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
372 // 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
373 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
374 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
375 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
376 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
377 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
378 }
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 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
380 // 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
381 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
382 < 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
383 // 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
384 // 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
385 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
386 lastSuspRange.setStart(fraction.getCoarseRange().getEnd());
6411
a21f7a3e4be5 Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6405
diff changeset
387 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
388 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
389 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
390 }
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 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
392 // 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
393 // Adjust and remember other 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
394 lastOtherRange = (Range) fraction.getSuspSedimentRange().clone();
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 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
396 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
397 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
398 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
399 }
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 }
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
401 }
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
402 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
403 // 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
404 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
405 // 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
406 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
407 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
408 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
409 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
410 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
411
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
412 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
413 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
414 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
415 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
416 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
417 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
418 lastSuspRange.setStart(lastOtherRange.getEnd());
6411
a21f7a3e4be5 Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6405
diff changeset
419 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
420 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
421 }
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 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
423 // 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
424 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
425 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
426 }
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
427 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
428 && 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
429 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
430 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
431 }
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 fairLoad.setTotal(km, total + fraction.getSuspSediment(), 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
433 }
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 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
435 lastSuspRange = suspSedimentRange;
6411
a21f7a3e4be5 Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6405
diff changeset
436 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
437 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
438 }
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
439 }
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
440 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
441 // 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
442 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
443 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
444 fraction.getSand() +
6411
a21f7a3e4be5 Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6405
diff changeset
445 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
446 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
447 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
448 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
449 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
450 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
451 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
452 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
453 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
454 // 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
455 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
456 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
457 }
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 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
459 // 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
460 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
461 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
462 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
463 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
464 }
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
465 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
466 && 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
467 && 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
468 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
469 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
470 }
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 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
472 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
473 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
474 // 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
475 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
476 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
477 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
478 }
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 }
df1a3826c691 SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6377
diff changeset
480 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
481 // 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
482 // 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
483 if (fraction.hasButSuspValues()) {
6500
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
484 double total = fraction.getCoarse() +
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
485 fraction.getFineMiddle() +
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
486 fraction.getSand() +
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
487 fraction.getSuspSand();
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
488 lastOtherRange = fraction.getCoarseRange();
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
489 lastOtherValue = total;
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
490 lastSuspRange = null;
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
491 }
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
492 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
493 lastSuspRange = fraction.getSuspSedimentRange();
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
494 lastSuspValue = fraction.getSuspSediment();
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
495 lastOtherRange = null;
f3b3742fbdad SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6499
diff changeset
496 }
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
497 }
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
498 }
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
499 return fairLoad;
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
500 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
501
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
502
6678
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
503 /**
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
504 * Transform values in load.
6948
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
505 * Background is to transform values measured in
6951
f5f60b3d9180 Doc fix.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6950
diff changeset
506 * 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
507 *
6678
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
508 * @param load The load of which values should be transformed.
6948
d4108d6c4000 Cosmetics, docs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6928
diff changeset
509 * @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
510 *
6678
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
511 * @return parameter load with transformed values.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6543
diff changeset
512 */
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
513 private SedimentLoad calculateUnit(SedimentLoad load, int year) {
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
514 SedimentDensity density =
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
515 SedimentDensityFactory.getSedimentDensity(river, kmLow, kmUp, year);
6951
f5f60b3d9180 Doc fix.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6950
diff changeset
516
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
517 for (double km: load.getKms()) {
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
518 double dens = density.getDensity(km, year);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
519 SedimentLoadFraction fraction = load.getFraction(km);
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
520 double coarse = fraction.getCoarse();
6376
8d6f23243d1d SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6375
diff changeset
521 double fineMiddle = fraction.getFineMiddle();
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
522 double sand = fraction.getSand();
6376
8d6f23243d1d SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6375
diff changeset
523 double suspSand = fraction.getSuspSand();
8d6f23243d1d SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6375
diff changeset
524 double bedSand = fraction.getSuspSandBed();
8d6f23243d1d SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 6375
diff changeset
525 double sediment = fraction.getSuspSediment();
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
526 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
527 double loadTotal = fraction.getLoadTotal();
6374
48e92ff57f23 SedimentLoad*: Set range to fraction/fractions values.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 5994
diff changeset
528 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
529 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
530 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
531 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
532 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
533 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
534 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
535 load.setLoadTotal(km, (loadTotal * dens), fraction.getLoadTotalRange());
4522
a9753f717b3d Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents: 4483
diff changeset
536 }
4373
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
537 return load;
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
538 }
1fb224bb2c6b Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff changeset
539 }
4483
cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents: 4373
diff changeset
540 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf-8 :
http://dive4elements.wald.intevation.org