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