Mercurial > dive4elements > river
annotate artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/SedimentLoadCalculation.java @ 6664:4ae1fbaaacdd
SedimentLoadUnknownFacet: Removed another debugging obstacle (wrong class for logger).
author | Felix Wolfsteller <felix.wolfsteller@intevation.de> |
---|---|
date | Thu, 25 Jul 2013 14:58:21 +0200 |
parents | 1ca0688dddc7 |
children | 0a70a320bfca |
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; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
39 protected int[][] epoch; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
40 protected String unit; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
41 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
42 public SedimentLoadCalculation() { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
43 } |
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 public CalculationResult calculate(SedimentLoadAccess access) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
46 logger.info("SedimentLoadCalculation.calculate"); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
47 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
48 String river = access.getRiver(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
49 String yearEpoch = access.getYearEpoch(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
50 String unit = access.getUnit(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
51 int[] period = null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
52 int[][] epoch = null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
53 double kmUp = access.getUpperKM(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
54 double kmLow = access.getLowerKM(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
55 if (yearEpoch.equals("year")) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
56 period = access.getPeriod(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
57 epoch = null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
58 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
59 else if (yearEpoch.equals("epoch") || yearEpoch.equals("off_epoch")) { |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
60 epoch = access.getEpochs(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
61 period = null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
62 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
63 else { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
64 addProblem("minfo.missing.year_epoch"); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
65 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
66 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
67 if (river == null) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
68 // TODO: i18n |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
69 addProblem("minfo.missing.river"); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
70 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
71 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
72 if (period == null && epoch == null) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
73 addProblem("minfo.missing.time"); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
74 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
75 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
76 if (!hasProblems()) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
77 this.river = river; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
78 this.yearEpoch = yearEpoch; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
79 this.unit = unit; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
80 this.period = period; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
81 this.epoch = epoch; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
82 this.kmUp = kmUp; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
83 this.kmLow = kmLow; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
84 return internalCalculate(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
85 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
86 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
87 return new CalculationResult(); |
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 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
90 private CalculationResult internalCalculate() { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
91 logger.debug("internalCalulate; mode:" + yearEpoch); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
92 if (yearEpoch.equals("year")) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
93 List<SedimentLoadResult> results = |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
94 new ArrayList<SedimentLoadResult>(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
95 for (int i = period[0]; i <= period[1]; i++) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
96 SedimentLoadResult res = calculateYear(i); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
97 results.add(res); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
98 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
99 return new CalculationResult( |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
100 results.toArray(new SedimentLoadResult[results.size()]), this); |
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 else if (yearEpoch.equals("epoch")) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
103 List<SedimentLoadResult> results = |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
104 new ArrayList<SedimentLoadResult>(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
105 for (int i = 0; i < epoch.length; i++) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
106 SedimentLoadResult res = calculateEpoch(i); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
107 results.add(res); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
108 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
109 return new CalculationResult( |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
110 results.toArray(new SedimentLoadResult[results.size()]), this); |
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 else if (yearEpoch.equals("off_epoch")) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
113 List<SedimentLoadResult> results = |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
114 new ArrayList<SedimentLoadResult>(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
115 for (int i = 0; i < epoch.length; i++) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
116 SedimentLoadResult res = calculateOffEpoch(i); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
117 results.add(res); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
118 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
119 return new CalculationResult( |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
120 results.toArray(new SedimentLoadResult[results.size()]), this); |
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 null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
123 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
124 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
125 private SedimentLoadResult calculateEpoch(int i) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
126 List<SedimentLoad> epochLoads = new ArrayList<SedimentLoad>(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
127 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
|
128 epochLoads.add(SedimentLoadFactory.getLoadWithData( |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
129 this.river, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
130 this.yearEpoch, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
131 this.kmLow, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
132 this.kmUp, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
133 j, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
134 j)); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
135 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
136 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
137 SedimentLoad resLoad = new SedimentLoad(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
138 TDoubleArrayList kms = new TDoubleArrayList(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
139 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
140 for (SedimentLoad load : epochLoads) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
141 for (double km : load.getKms()) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
142 if (!kms.contains(km)) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
143 kms.add(km); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
144 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
145 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
146 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
147 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
148 for (int j = 0; j < kms.size(); j++) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
149 int cSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
150 int fmSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
151 int sSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
152 int ssSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
153 int ssbSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
154 int sseSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
155 double km = kms.get(j); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
156 for (SedimentLoad load : epochLoads) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
157 SedimentLoadFraction f = load.getFraction(km); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
158 if (f.getCoarse() > 0d) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
159 double c = resLoad.getFraction(km).getCoarse(); |
6374
48e92ff57f23
SedimentLoad*: Set range to fraction/fractions values.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5994
diff
changeset
|
160 resLoad.setCoarse(km, c + f.getCoarse(), f.getCoarseRange()); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
161 cSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
162 } |
6376
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
163 if (f.getFineMiddle() > 0d) { |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
164 double fm = resLoad.getFraction(km).getFineMiddle(); |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
165 resLoad.setFineMiddle(km, fm + f.getFineMiddle(), f.getFineMiddleRange()); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
166 fmSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
167 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
168 if (f.getSand() > 0d) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
169 double s = resLoad.getFraction(km).getSand(); |
6374
48e92ff57f23
SedimentLoad*: Set range to fraction/fractions values.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5994
diff
changeset
|
170 resLoad.setSand(km, s + f.getSand(), f.getSandRange()); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
171 sSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
172 } |
6376
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
173 if (f.getSuspSand() > 0d) { |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
174 double s = resLoad.getFraction(km).getSuspSand(); |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
175 resLoad.setSuspSand(km, s + f.getSuspSand(), f.getSuspSandRange()); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
176 ssSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
177 } |
6376
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
178 if (f.getSuspSandBed() > 0d) { |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
179 double s = resLoad.getFraction(km).getSuspSandBed(); |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
180 resLoad.setSuspSandBed(km, s + f.getSuspSandBed(), f.getSuspSandBedRange()); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
181 ssbSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
182 } |
6376
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
183 if (f.getSuspSediment() > 0d) { |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
184 double s = resLoad.getFraction(km).getSuspSediment(); |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
185 resLoad.setSuspSediment(km, s + f.getSuspSediment(), f.getSuspSedimentRange()); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
186 sseSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
187 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
188 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
189 SedimentLoadFraction fr = resLoad.getFraction(km); |
6543
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
190 // Prevent divisions by zero, the fraction defaults to 0d. |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
191 if (cSum != 0) { |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
192 resLoad.setCoarse(km, fr.getCoarse()/cSum, fr.getCoarseRange()); |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
193 } |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
194 if (fmSum != 0) { |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
195 resLoad.setFineMiddle(km, fr.getFineMiddle()/fmSum, |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
196 fr.getFineMiddleRange()); |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
197 } |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
198 if (sSum != 0) { |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
199 resLoad.setSand(km, fr.getSand()/sSum, fr.getSandRange()); |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
200 } |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
201 if (ssSum != 0) { |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
202 resLoad.setSuspSand(km, fr.getSuspSand()/ssSum, |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
203 fr.getSuspSandRange()); |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
204 } |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
205 if (ssbSum != 0) { |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
206 resLoad.setSuspSandBed(km, fr.getSuspSandBed()/ssbSum, |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
207 fr.getSuspSandBedRange()); |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
208 } |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
209 if (sseSum != 0) { |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
210 resLoad.setSuspSediment(km, fr.getSuspSediment()/sseSum, fr.getSuspSedimentRange()); |
1ca0688dddc7
issue1235: Prevent division by zeros leading to NaNs.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6500
diff
changeset
|
211 } |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
212 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
213 resLoad.setDescription(""); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
214 resLoad.setEpoch(true); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
215 |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
216 SedimentLoadResult result; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
217 SedimentLoad sl = calculateTotalLoad(resLoad, this.epoch[i][0]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
218 if (this.unit.equals("m3_per_a")) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
219 SedimentLoad slu = calculateUnit(sl, this.epoch[i][0]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
220 result = new SedimentLoadResult( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
221 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
222 this.epoch[i][1], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
223 slu); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
224 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
225 else { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
226 result = new SedimentLoadResult( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
227 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
228 this.epoch[i][1], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
229 sl); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
230 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
231 |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
232 return result; |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
233 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
234 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
235 private SedimentLoadResult calculateOffEpoch(int i) { |
5661
b6297a67823e
SedimentLoadFactory: Renamed getLoadwithData to getLoadWithData, update callers.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4543
diff
changeset
|
236 SedimentLoad load = SedimentLoadFactory.getLoadWithData( |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
237 this.river, |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
238 this.yearEpoch, |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
239 kmLow, |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
240 kmUp, |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
241 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
242 this.epoch[i][1]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
243 SedimentLoadResult result; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
244 SedimentLoad sl = calculateTotalLoad(load, this.epoch[i][0]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
245 if (unit.equals("m3_per_a")) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
246 SedimentLoad slu = calculateUnit(sl, epoch[i][0]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
247 result = new SedimentLoadResult( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
248 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
249 this.epoch[i][1], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
250 slu); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
251 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
252 else { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
253 result = new SedimentLoadResult( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
254 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
255 this.epoch[i][1], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
256 sl); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
257 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
258 |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
259 return result; |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
260 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
261 |
5788 | 262 /** Fetch loads for a single year, calculate total and |
263 * return the result containing both. */ | |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
264 private SedimentLoadResult calculateYear(int y) { |
5661
b6297a67823e
SedimentLoadFactory: Renamed getLoadwithData to getLoadWithData, update callers.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4543
diff
changeset
|
265 SedimentLoad load = SedimentLoadFactory.getLoadWithData( |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
266 this.river, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
267 this.yearEpoch, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
268 this.kmLow, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
269 this.kmUp, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
270 y, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
271 y); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
272 |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
273 SedimentLoadResult result; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
274 SedimentLoad sl = calculateTotalLoad(load, y); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
275 if (unit.equals("m3_per_a")) { |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
276 SedimentLoad slu = calculateUnit(sl, y); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
277 result = new SedimentLoadResult(y, 0, slu); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
278 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
279 else { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
280 result = new SedimentLoadResult(y, 0, sl); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
281 } |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
282 return result; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
283 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
284 |
6374
48e92ff57f23
SedimentLoad*: Set range to fraction/fractions values.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5994
diff
changeset
|
285 /** Add up the loads of a year. */ |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
286 private SedimentLoad calculateTotalLoad(SedimentLoad load, int year) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
287 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
|
288 boolean problemThisYear = false; |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
289 if (!load.hasCoarse()) { |
5742
a4a894b15c35
Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5685
diff
changeset
|
290 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
|
291 problemThisYear = true; |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
292 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
293 if (!load.hasFineMiddle()) { |
5742
a4a894b15c35
Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5685
diff
changeset
|
294 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
|
295 problemThisYear = true; |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
296 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
297 if (!load.hasSand()) { |
5742
a4a894b15c35
Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5685
diff
changeset
|
298 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
|
299 problemThisYear = true; |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
300 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
301 if (!load.hasSuspSand()) { |
5742
a4a894b15c35
Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5685
diff
changeset
|
302 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
|
303 problemThisYear = true; |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
304 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
305 if (!load.hasSuspSediment()) { |
5742
a4a894b15c35
Minor part of issue1235: dont show decimals of years.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5685
diff
changeset
|
306 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
|
307 problemThisYear = true; |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
308 } |
5685
756df79274e1
SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible).
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5661
diff
changeset
|
309 if (problemThisYear) { |
756df79274e1
SedimentLoadCalculation: Partly fix 1250 (calculate more than one total if possible).
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5661
diff
changeset
|
310 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
|
311 return load; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
312 } |
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
|
313 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
|
314 } |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
315 |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
316 |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
317 /** Returns true if all fraction values except SuspSediment are unset. */ |
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
|
318 private boolean hasOnlySuspValues(SedimentLoadFraction fraction) { |
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
|
319 return (fraction.getSuspSediment() != 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
|
320 fraction.getCoarse() == 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
|
321 fraction.getFineMiddle() == 0d && |
6411
a21f7a3e4be5
Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6405
diff
changeset
|
322 fraction.getSand() == 0d && |
a21f7a3e4be5
Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6405
diff
changeset
|
323 fraction.getSuspSand() == 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
|
324 } |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
325 |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
326 |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
327 /** Returns true if all fraction values except SuspSediment are set. */ |
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
|
328 private boolean hasButSuspValues(SedimentLoadFraction fraction) { |
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
|
329 return (fraction.getSuspSediment() == 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
|
330 fraction.getCoarse() != 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
|
331 fraction.getFineMiddle() != 0d && |
6411
a21f7a3e4be5
Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6405
diff
changeset
|
332 fraction.getSand() != 0d && |
a21f7a3e4be5
Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6405
diff
changeset
|
333 fraction.getSuspSand() != 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
|
334 } |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
335 |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
336 |
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 /** Returns true if all fraction needed for total calculation 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
|
338 private boolean complete(SedimentLoadFraction fraction) { |
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 return (fraction.getCoarse() != 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
|
340 fraction.getFineMiddle() != 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
|
341 fraction.getSand() != 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
|
342 fraction.getSuspSand() != 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
|
343 fraction.getSuspSediment() != 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
|
344 } |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
345 |
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 |
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 /** |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
348 * Set total values in 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
|
349 * Therefore, run over the kms and find ranges where either all |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
350 * 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
|
351 * Merge these ranges and add (maybe new) respective fractions to |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
352 * 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
|
353 * @param load SedimentLoad to add total values (and ranges) to. |
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 * @return input param 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 */ |
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 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
|
357 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
|
358 |
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 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
|
360 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
|
361 |
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
|
362 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
|
363 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
|
364 |
6498
686bbfeacfe5
issue1235: sort km when running over them for total sediment load.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6411
diff
changeset
|
365 TreeSet<Double> kms = new TreeSet<Double>(load.getKms()); |
686bbfeacfe5
issue1235: sort km when running over them for total sediment load.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6411
diff
changeset
|
366 |
686bbfeacfe5
issue1235: sort km when running over them for total sediment load.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6411
diff
changeset
|
367 for (double km: kms) { |
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
|
368 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
|
369 SedimentLoadFraction fraction = load.getFraction(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
|
370 if (complete(fraction)) { |
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 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
|
372 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
|
373 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
|
374 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
|
375 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
|
376 // 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
|
377 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
|
378 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
|
379 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
|
380 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
|
381 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
|
382 } |
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 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
|
384 // 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
|
385 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
|
386 < 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
|
387 // 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
|
388 // 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
|
389 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
|
390 lastSuspRange.setStart(fraction.getCoarseRange().getEnd()); |
6411
a21f7a3e4be5
Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6405
diff
changeset
|
391 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
|
392 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
|
393 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
|
394 } |
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 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
|
396 // 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
|
397 // 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
|
398 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
|
399 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
|
400 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
|
401 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
|
402 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
|
403 } |
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 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
405 } |
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
|
406 else if (hasOnlySuspValues(fraction) && 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
|
407 // 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
|
408 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
|
409 // 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
|
410 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
|
411 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
|
412 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
|
413 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
|
414 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
|
415 |
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 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
|
417 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
|
418 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
|
419 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
|
420 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
|
421 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
|
422 lastSuspRange.setStart(lastOtherRange.getEnd()); |
6411
a21f7a3e4be5
Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6405
diff
changeset
|
423 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
|
424 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
|
425 } |
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 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
|
427 // 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
|
428 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
|
429 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
|
430 } |
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 if (Math.abs(suspSedimentRange.getEnd() - lastOtherRange.getEnd()) < 0.1d) { |
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 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
|
433 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
|
434 } |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
435 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
|
436 } |
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 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
|
438 lastSuspRange = suspSedimentRange; |
6411
a21f7a3e4be5
Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6405
diff
changeset
|
439 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
|
440 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
|
441 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
442 } |
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
|
443 else if (hasButSuspValues(fraction) && 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
|
444 // 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
|
445 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
|
446 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
|
447 fraction.getSand() + |
6411
a21f7a3e4be5
Revert to rev 6404 for issue1235
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6405
diff
changeset
|
448 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
|
449 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
|
450 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
|
451 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
|
452 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
|
453 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
|
454 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
|
455 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
|
456 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
|
457 // 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
|
458 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
|
459 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
|
460 } |
df1a3826c691
SedimentLoadCalculation: Calculate and set Total Load values ind a brutal range merging manner.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6377
diff
changeset
|
461 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
|
462 // 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
|
463 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
|
464 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
|
465 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
|
466 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
|
467 } |
6499
857f4d51d859
issue1235: Avoid NPE.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6498
diff
changeset
|
468 if (lastSuspRange != null && 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
|
469 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
|
470 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
|
471 } |
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 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
|
473 } |
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 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
|
475 // 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
|
476 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
|
477 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
|
478 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
|
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 } |
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 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
|
482 // 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
|
483 // Stay as we are. |
6500
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
484 if (hasButSuspValues(fraction)) { |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
485 double total = fraction.getCoarse() + |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
486 fraction.getFineMiddle() + |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
487 fraction.getSand() + |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
488 fraction.getSuspSand(); |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
489 lastOtherRange = fraction.getCoarseRange(); |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
490 lastOtherValue = total; |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
491 lastSuspRange = null; |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
492 } |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
493 else if (hasOnlySuspValues(fraction)) { |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
494 lastSuspRange = fraction.getSuspSedimentRange(); |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
495 lastSuspValue = fraction.getSuspSediment(); |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
496 lastOtherRange = null; |
f3b3742fbdad
SedimentLoadCalculation: Handle boring cases in total sed load calculation.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6499
diff
changeset
|
497 } |
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
|
498 } |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
499 } |
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
|
500 return fairLoad; |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
501 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
502 |
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
|
503 |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
504 private SedimentLoad calculateUnit(SedimentLoad load, int year) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
505 SedimentDensity density = |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
506 SedimentDensityFactory.getSedimentDensity(river, kmLow, kmUp, year); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
507 for (double km: load.getKms()) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
508 double dens = density.getDensity(km, year); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
509 SedimentLoadFraction fraction = load.getFraction(km); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
510 double coarse = fraction.getCoarse(); |
6376
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
511 double fineMiddle = fraction.getFineMiddle(); |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
512 double sand = fraction.getSand(); |
6376
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
513 double suspSand = fraction.getSuspSand(); |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
514 double bedSand = fraction.getSuspSandBed(); |
8d6f23243d1d
SedimentLoad*: Switch from snake_case to camelCase.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
6375
diff
changeset
|
515 double sediment = fraction.getSuspSediment(); |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
516 double total = fraction.getTotal(); |
6374
48e92ff57f23
SedimentLoad*: Set range to fraction/fractions values.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
5994
diff
changeset
|
517 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
|
518 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
|
519 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
|
520 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
|
521 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
|
522 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
|
523 load.setTotal(km, (total * dens), fraction.getTotalRange()); |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
524 } |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
525 return load; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
526 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
527 } |
4483
cc6323401643
Cosmetics: Removed some WSP, obsolete imports. Documentation added.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4373
diff
changeset
|
528 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf-8 : |