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

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

http://dive4elements.wald.intevation.org