Mercurial > dive4elements > river
annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/minfo/SedimentLoadCalculation.java @ 4630:63368dcc3f94
Moved MovingAverage into math package
author | Sascha L. Teichmann <teichmann@intevation.de> |
---|---|
date | Mon, 03 Dec 2012 18:44:48 +0100 |
parents | e357716cf346 |
children | b6297a67823e |
rev | line source |
---|---|
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
1 package de.intevation.flys.artifacts.model.minfo; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
2 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
3 import gnu.trove.TDoubleArrayList; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
4 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
5 import java.util.ArrayList; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
6 import java.util.List; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
7 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
8 import org.apache.log4j.Logger; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
9 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
10 import de.intevation.flys.artifacts.access.SedimentLoadAccess; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
11 import de.intevation.flys.artifacts.model.Calculation; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
12 import de.intevation.flys.artifacts.model.CalculationResult; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
13 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
14 |
4483
cc6323401643
Cosmetics: Removed some WSP, obsolete imports. Documentation added.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4373
diff
changeset
|
15 /** Calculate sediment load. */ |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
16 public class SedimentLoadCalculation |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
17 extends Calculation |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
18 { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
19 |
4483
cc6323401643
Cosmetics: Removed some WSP, obsolete imports. Documentation added.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4373
diff
changeset
|
20 /** Private logger. */ |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
21 private static final Logger logger = Logger |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
22 .getLogger(SedimentLoadCalculation.class); |
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 protected String river; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
25 protected String yearEpoch; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
26 protected double kmUp; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
27 protected double kmLow; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
28 protected int[] period; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
29 protected int[][] epoch; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
30 protected String unit; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
31 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
32 public SedimentLoadCalculation() { |
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 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
35 public CalculationResult calculate(SedimentLoadAccess access) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
36 logger.info("SedimentLoadCalculation.calculate"); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
37 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
38 String river = access.getRiver(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
39 String yearEpoch = access.getYearEpoch(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
40 String unit = access.getUnit(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
41 int[] period = null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
42 int[][] epoch = null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
43 double kmUp = access.getUpperKM(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
44 double kmLow = access.getLowerKM(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
45 if (yearEpoch.equals("year")) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
46 period = access.getPeriod(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
47 epoch = null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
48 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
49 else if (yearEpoch.equals("epoch") || yearEpoch.equals("off_epoch")) { |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
50 epoch = access.getEpochs(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
51 period = null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
52 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
53 else { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
54 addProblem("minfo.missing.year_epoch"); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
55 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
56 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
57 if (river == null) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
58 // TODO: i18n |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
59 addProblem("minfo.missing.river"); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
60 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
61 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
62 if (period == null && epoch == null) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
63 addProblem("minfo.missing.time"); |
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 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
66 if (!hasProblems()) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
67 this.river = river; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
68 this.yearEpoch = yearEpoch; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
69 this.unit = unit; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
70 this.period = period; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
71 this.epoch = epoch; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
72 this.kmUp = kmUp; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
73 this.kmLow = kmLow; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
74 return internalCalculate(); |
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 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
77 return new CalculationResult(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
78 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
79 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
80 private CalculationResult internalCalculate() { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
81 logger.debug("internalCalulate; mode:" + yearEpoch); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
82 if (yearEpoch.equals("year")) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
83 List<SedimentLoadResult> results = |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
84 new ArrayList<SedimentLoadResult>(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
85 for (int i = period[0]; i <= period[1]; i++) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
86 SedimentLoadResult res = calculateYear(i); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
87 results.add(res); |
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 results.toArray(new SedimentLoadResult[results.size()]), this); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
91 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
92 else if (yearEpoch.equals("epoch")) { |
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 = 0; i < epoch.length; i++) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
96 SedimentLoadResult res = calculateEpoch(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("off_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 = calculateOffEpoch(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 return null; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
113 } |
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 private SedimentLoadResult calculateEpoch(int i) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
116 List<SedimentLoad> epochLoads = new ArrayList<SedimentLoad>(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
117 for (int j = epoch[i][0]; j < epoch[i][1]; j++) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
118 epochLoads.add(SedimentLoadFactory.getLoadwithData( |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
119 this.river, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
120 this.yearEpoch, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
121 this.kmLow, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
122 this.kmUp, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
123 j, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
124 j)); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
125 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
126 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
127 SedimentLoad resLoad = new SedimentLoad(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
128 TDoubleArrayList kms = new TDoubleArrayList(); |
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 for (SedimentLoad load : epochLoads) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
131 for (double km : load.getKms()) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
132 if (!kms.contains(km)) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
133 kms.add(km); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
134 } |
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 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
138 for (int j = 0; j < kms.size(); j++) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
139 int cSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
140 int fmSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
141 int sSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
142 int ssSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
143 int ssbSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
144 int sseSum = 0; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
145 double km = kms.get(j); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
146 for (SedimentLoad load : epochLoads) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
147 SedimentLoadFraction f = load.getFraction(km); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
148 if (f.getCoarse() > 0d) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
149 double c = resLoad.getFraction(km).getCoarse(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
150 resLoad.setCoarse(km, c + f.getCoarse()); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
151 cSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
152 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
153 if (f.getFine_middle() > 0d) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
154 double fm = resLoad.getFraction(km).getFine_middle(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
155 resLoad.setFineMiddle(km, fm + f.getFine_middle()); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
156 fmSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
157 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
158 if (f.getSand() > 0d) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
159 double s = resLoad.getFraction(km).getSand(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
160 resLoad.setSand(km, s + f.getSand()); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
161 sSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
162 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
163 if (f.getSusp_sand() > 0d) { |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
164 double s = resLoad.getFraction(km).getSusp_sand(); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
165 resLoad.setSuspSand(km, s + f.getSusp_sand()); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
166 ssSum++; |
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.getSusp_sand_bed() > 0d) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
169 double s = resLoad.getFraction(km).getSusp_sand_bed(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
170 resLoad.setSuspSandBed(km, s + f.getSusp_sand_bed()); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
171 ssbSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
172 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
173 if (f.getSusp_sediment() > 0d) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
174 double s = resLoad.getFraction(km).getSusp_sediment(); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
175 resLoad.setSuspSediment(km, s + f.getSusp_sediment()); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
176 sseSum++; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
177 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
178 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
179 SedimentLoadFraction fr = resLoad.getFraction(km); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
180 resLoad.setCoarse(km, fr.getCoarse()/cSum); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
181 resLoad.setFineMiddle(km, fr.getFine_middle()/fmSum); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
182 resLoad.setSand(km, fr.getSand()/sSum); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
183 resLoad.setSuspSand(km, fr.getSusp_sand()/ssSum); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
184 resLoad.setSuspSandBed(km, fr.getSusp_sand_bed()/ssbSum); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
185 resLoad.setSuspSediment(km, fr.getSusp_sediment()/sseSum); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
186 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
187 resLoad.setDescription(""); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
188 resLoad.setEpoch(true); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
189 |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
190 SedimentLoadResult result; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
191 SedimentLoad sl = calculateTotalLoad(resLoad, this.epoch[i][0]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
192 if (this.unit.equals("m3_per_a")) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
193 SedimentLoad slu = calculateUnit(sl, this.epoch[i][0]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
194 result = new SedimentLoadResult( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
195 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
196 this.epoch[i][1], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
197 slu); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
198 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
199 else { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
200 result = new SedimentLoadResult( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
201 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
202 this.epoch[i][1], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
203 sl); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
204 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
205 |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
206 return result; |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
207 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
208 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
209 private SedimentLoadResult calculateOffEpoch(int i) { |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
210 SedimentLoad load = SedimentLoadFactory.getLoadwithData( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
211 this.river, |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
212 this.yearEpoch, |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
213 kmLow, |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
214 kmUp, |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
215 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
216 this.epoch[i][1]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
217 SedimentLoadResult result; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
218 SedimentLoad sl = calculateTotalLoad(load, this.epoch[i][0]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
219 if (unit.equals("m3_per_a")) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
220 SedimentLoad slu = calculateUnit(sl, epoch[i][0]); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
221 result = new SedimentLoadResult( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
222 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
223 this.epoch[i][1], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
224 slu); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
225 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
226 else { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
227 result = new SedimentLoadResult( |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
228 this.epoch[i][0], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
229 this.epoch[i][1], |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
230 sl); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
231 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
232 |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
233 return result; |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
234 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
235 |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
236 private SedimentLoadResult calculateYear(int y) { |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
237 SedimentLoad load = SedimentLoadFactory.getLoadwithData( |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
238 this.river, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
239 this.yearEpoch, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
240 this.kmLow, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
241 this.kmUp, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
242 y, |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
243 y); |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
244 |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
245 SedimentLoadResult result; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
246 SedimentLoad sl = calculateTotalLoad(load, y); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
247 if (unit.equals("m3_per_a")) { |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
248 SedimentLoad slu = calculateUnit(sl, y); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
249 result = new SedimentLoadResult(y, 0, slu); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
250 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
251 else { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
252 result = new SedimentLoadResult(y, 0, sl); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
253 } |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
254 return result; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
255 } |
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 private SedimentLoad calculateTotalLoad(SedimentLoad load, int year) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
258 logger.debug("calculateTotalLoad"); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
259 if (!load.hasCoarse()) { |
4543
e357716cf346
Removed some debug logs and changed report strings.
Raimund Renkert <rrenkert@intevation.de>
parents:
4522
diff
changeset
|
260 addProblem(year, "missing.fraction.coarse"); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
261 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
262 if (!load.hasFineMiddle()) { |
4543
e357716cf346
Removed some debug logs and changed report strings.
Raimund Renkert <rrenkert@intevation.de>
parents:
4522
diff
changeset
|
263 addProblem(year, "missing.fraction.fine_middle"); |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
264 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
265 if (!load.hasSand()) { |
4543
e357716cf346
Removed some debug logs and changed report strings.
Raimund Renkert <rrenkert@intevation.de>
parents:
4522
diff
changeset
|
266 addProblem(year, "missing.fraction.sand"); |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
267 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
268 if (!load.hasSuspSand()) { |
4543
e357716cf346
Removed some debug logs and changed report strings.
Raimund Renkert <rrenkert@intevation.de>
parents:
4522
diff
changeset
|
269 addProblem(year, "missing.fraction.susp_sand"); |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
270 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
271 if (!load.hasSuspSediment()) { |
4543
e357716cf346
Removed some debug logs and changed report strings.
Raimund Renkert <rrenkert@intevation.de>
parents:
4522
diff
changeset
|
272 addProblem(year, "missing.fraction.susp_sediment"); |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
273 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
274 if (hasProblems()) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
275 return load; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
276 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
277 for(double km : load.getKms()) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
278 SedimentLoadFraction fraction = load.getFraction(km); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
279 double total = 0d; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
280 if ((fraction.getCoarse() <= 0d && load.hasCoarse())){ |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
281 addProblem(km, "missing.data.coarse"); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
282 continue; |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
283 } |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
284 if (fraction.getFine_middle() <= 0d && load.hasFineMiddle()) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
285 addProblem(km, "missing.data.fine_middle"); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
286 continue; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
287 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
288 if (fraction.getSand() <= 0d && load.hasSand()) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
289 addProblem(km, "missing data.sand"); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
290 continue; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
291 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
292 if (fraction.getSusp_sand() <= 0d && load.hasSuspSand()) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
293 addProblem(km, "missing.data.susp_sand"); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
294 continue; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
295 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
296 if (fraction.getSusp_sediment() <= 0d && load.hasSuspSediment()) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
297 addProblem(km, "missing.data.susp_sediment"); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
298 continue; |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
299 } |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
300 total += fraction.getCoarse() + |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
301 fraction.getFine_middle() + |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
302 fraction.getSand() + |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
303 fraction.getSusp_sand() + |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
304 fraction.getSusp_sediment(); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
305 load.setTotal(km, total); |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
306 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
307 return load; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
308 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
309 |
4522
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
310 private SedimentLoad calculateUnit(SedimentLoad load, int year) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
311 SedimentDensity density = |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
312 SedimentDensityFactory.getSedimentDensity(river, kmLow, kmUp, year); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
313 for (double km: load.getKms()) { |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
314 double dens = density.getDensity(km, year); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
315 SedimentLoadFraction fraction = load.getFraction(km); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
316 double coarse = fraction.getCoarse(); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
317 double fineMiddle = fraction.getFine_middle(); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
318 double sand = fraction.getSand(); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
319 double suspSand = fraction.getSusp_sand(); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
320 double bedSand = fraction.getSusp_sand_bed(); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
321 double sediment = fraction.getSusp_sediment(); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
322 double total = fraction.getTotal(); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
323 load.setCoarse(km, (coarse * dens)); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
324 load.setFineMiddle(km, (fineMiddle * dens)); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
325 load.setSand(km, (sand * dens)); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
326 load.setSuspSand(km, (suspSand * dens)); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
327 load.setSuspSandBed(km, (bedSand * dens)); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
328 load.setSuspSediment(km, (sediment * dens)); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
329 load.setTotal(km, (total * dens)); |
a9753f717b3d
Completed and fixed sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
4483
diff
changeset
|
330 } |
4373
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
331 return load; |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
332 } |
1fb224bb2c6b
Added sediment load calculation.
Raimund Renkert <rrenkert@intevation.de>
parents:
diff
changeset
|
333 } |
4483
cc6323401643
Cosmetics: Removed some WSP, obsolete imports. Documentation added.
Felix Wolfsteller <felix.wolfsteller@intevation.de>
parents:
4373
diff
changeset
|
334 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf-8 : |