dive4elements/river: flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/minfo/SedimentLoadCalculation.java annotate

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/minfo/SedimentLoadCalculation.java @ 4483:cc6323401643

Cosmetics: Removed some WSP, obsolete imports. Documentation added.

author	Felix Wolfsteller <felix.wolfsteller@intevation.de>
date	Tue, 13 Nov 2012 16:52:01 +0100
parents	1fb224bb2c6b
children	a9753f717b3d

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.Iterator;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	7 import java.util.List;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	8 import java.util.Set;
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 org.apache.log4j.Logger;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	11 import org.jfree.util.Log;
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 de.intevation.flys.artifacts.access.SedimentLoadAccess;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	14 import de.intevation.flys.artifacts.model.Calculation;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	15 import de.intevation.flys.artifacts.model.CalculationResult;
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
4483 cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 4373 diff changeset	18 /** Calculate sediment load. */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	19 public class SedimentLoadCalculation
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	20 extends Calculation
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	21 {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	22
4483 cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 4373 diff changeset	23 /** Private logger. */
4373 1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	24 private static final Logger logger = Logger
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	25 .getLogger(SedimentLoadCalculation.class);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	26
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	27 protected String river;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	28 protected String yearEpoch;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	29 protected double kmUp;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	30 protected double kmLow;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	31 protected int[] period;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	32 protected int[][] epoch;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	33 protected String unit;
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 SedimentLoadCalculation() {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	36 }
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 public CalculationResult calculate(SedimentLoadAccess access) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	39 logger.info("SedimentLoadCalculation.calculate");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	40
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	41 String river = access.getRiver();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	42 String yearEpoch = access.getYearEpoch();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	43 String unit = access.getUnit();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	44 int[] period = null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	45 int[][] epoch = null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	46 double kmUp = access.getUpperKM();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	47 double kmLow = access.getLowerKM();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	48 if (yearEpoch.equals("year")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	49 period = access.getPeriod();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	50 epoch = null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	51 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	52 else if (yearEpoch.equals("epoch")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	53 epoch = access.getEpochs();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	54 period = null;
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 else {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	57 addProblem("minfo.missing.year_epoch");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	58 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	59
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	60 if (river == null) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	61 // TODO: i18n
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	62 addProblem("minfo.missing.river");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	63 }
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 if (period == null && epoch == null) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	66 addProblem("minfo.missing.time");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	67 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	68
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	69 if (!hasProblems()) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	70 this.river = river;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	71 this.yearEpoch = yearEpoch;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	72 this.unit = unit;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	73 this.period = period;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	74 this.epoch = epoch;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	75 this.kmUp = kmUp;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	76 this.kmLow = kmLow;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	77 return internalCalculate();
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 return new CalculationResult();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	81 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	82
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	83 private CalculationResult internalCalculate() {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	84 logger.debug("internalCalulate; mode:" + yearEpoch);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	85 if (yearEpoch.equals("year")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	86 List<SedimentLoadResult> results =
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	87 new ArrayList<SedimentLoadResult>();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	88 for (int i = period[0]; i <= period[1]; i++) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	89 logger.debug("calculating for year: " + i);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	90 SedimentLoadResult res = calculateYear(i);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	91 results.add(res);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	92 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	93 return new CalculationResult(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	94 results.toArray(new SedimentLoadResult[results.size()]), this);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	95 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	96 else if (yearEpoch.equals("epoch")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	97 List<SedimentLoadResult> results =
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	98 new ArrayList<SedimentLoadResult>();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	99 for (int i = 0; i < epoch.length; i++) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	100 SedimentLoadResult res = calculateEpoch(i);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	101 results.add(res);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	102 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	103 return new CalculationResult(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	104 results.toArray(new SedimentLoadResult[results.size()]), this);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	105 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	106 else if (yearEpoch.equals("off_epoch")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	107 List<SedimentLoadResult> results =
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	108 new ArrayList<SedimentLoadResult>();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	109 for (int i = 0; i < epoch.length; i++) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	110 SedimentLoadResult res = calculateOffEpoch(i);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	111 results.add(res);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	112 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	113 return new CalculationResult(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	114 results.toArray(new SedimentLoadResult[results.size()]), this);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	115 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	116 return null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	117 }
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 private SedimentLoadResult calculateEpoch(int i) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	120 List<SedimentLoad> epochLoads = new ArrayList<SedimentLoad>();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	121 for (int j = epoch[i][0]; j < epoch[i][1]; j++) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	122 epochLoads.add(SedimentLoadFactory.getLoadwithData(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	123 this.river,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	124 this.yearEpoch,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	125 this.kmLow,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	126 this.kmUp,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	127 j,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	128 j));
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
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	131 SedimentLoad resLoad = new SedimentLoad();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	132 TDoubleArrayList kms = new TDoubleArrayList();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	133
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	134 for (SedimentLoad load : epochLoads) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	135 for (double km : load.getKms()) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	136 if (!kms.contains(km)) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	137 kms.add(km);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	138 }
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 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	141
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	142 for (int j = 0; j < kms.size(); j++) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	143 int cSum = 0;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	144 int fmSum = 0;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	145 int sSum = 0;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	146 int ssSum = 0;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	147 int ssbSum = 0;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	148 int sseSum = 0;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	149 double km = kms.get(j);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	150 for (SedimentLoad load : epochLoads) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	151 SedimentLoadFraction f = load.getFraction(km);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	152 if (f.getCoarse() > 0d) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	153 double c = resLoad.getFraction(km).getCoarse();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	154 resLoad.setCoarse(km, c + f.getCoarse());
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	155 cSum++;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	156 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	157 if (f.getFine_middle() > 0d) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	158 double fm = resLoad.getFraction(km).getFine_middle();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	159 resLoad.setFineMiddle(km, fm + f.getFine_middle());
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	160 fmSum++;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	161 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	162 if (f.getSand() > 0d) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	163 double s = resLoad.getFraction(km).getSand();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	164 resLoad.setSand(km, s + f.getSand());
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	165 sSum++;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	166 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	167 if (f.getSusp_sand() > 0d) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	168 double s = resLoad.getFraction(km).getSand();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	169 resLoad.setSuspSand(km, s + f.getSusp_sand());
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	170 ssSum++;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	171 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	172 if (f.getSusp_sand_bed() > 0d) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	173 double s = resLoad.getFraction(km).getSusp_sand_bed();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	174 resLoad.setSuspSandBed(km, s + f.getSusp_sand_bed());
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	175 ssbSum++;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	176 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	177 if (f.getSusp_sediment() > 0d) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	178 double s = resLoad.getFraction(km).getSusp_sediment();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	179 resLoad.setSuspSediment(km, s + f.getSusp_sediment());
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	180 sseSum++;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	181 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	182 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	183 SedimentLoadFraction fr = resLoad.getFraction(km);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	184 resLoad.setCoarse(km, fr.getCoarse()/cSum);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	185 resLoad.setFineMiddle(km, fr.getFine_middle()/fmSum);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	186 resLoad.setSand(km, fr.getSand()/sSum);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	187 resLoad.setSuspSand(km, fr.getSusp_sand()/ssSum);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	188 resLoad.setSuspSandBed(km, fr.getSusp_sand_bed()/ssbSum);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	189 resLoad.setSuspSediment(km, fr.getSusp_sediment()/sseSum);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	190 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	191 resLoad.setDescription("");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	192 resLoad.setEpoch(true);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	193
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	194 SedimentLoad sl = calculateTotalLoad(resLoad);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	195 return new SedimentLoadResult(epoch[i][0], epoch[i][1], sl);
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
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	198 private SedimentLoadResult calculateOffEpoch(int i) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	199 return null;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	200 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	201
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	202 private SedimentLoadResult calculateYear(int y) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	203 SedimentLoad load = SedimentLoadFactory.getLoadwithData(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	204 this.river,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	205 this.yearEpoch,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	206 this.kmLow,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	207 this.kmUp,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	208 y,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	209 y);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	210
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	211 SedimentLoad sl = calculateTotalLoad(load);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	212 if (unit.equals("m3_per_a")) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	213 SedimentLoad slu = calculateUnit(sl);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	214 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	215
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	216 SedimentLoadResult result = new SedimentLoadResult(
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	217 y,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	218 0,
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	219 sl);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	220 return result;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	221 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	222
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	223 private SedimentLoad validateData(SedimentLoad load) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	224 SedimentLoad clean = new SedimentLoad();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	225 Set<Double> kms = load.getKms();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	226 for (double km : kms) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	227 SedimentLoadFraction fraction = load.getFraction(km);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	228 if (fraction.getCoarse() > 0 &&
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	229 fraction.getFine_middle() > 0 &&
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	230 fraction.getSand() > 0 &&
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	231 fraction.getSusp_sand() > 0&&
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	232 fraction.getSusp_sand_bed() > 0 &&
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	233 fraction.getSusp_sediment() > 0) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	234 clean.addKm(km, fraction);
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 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	237 return clean;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	238 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	239
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	240 private SedimentLoad calculateTotalLoad(SedimentLoad load) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	241 Log.debug("calculateTotalLoad");
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	242 SedimentLoad clean = validateData(load);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	243 if (clean.getKms().size() > 0) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	244 Iterator<Double> iter = clean.getKms().iterator();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	245 while (iter.hasNext()) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	246 double km = iter.next();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	247 SedimentLoadFraction f =clean.getFraction(km);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	248 double total = f.getCoarse() +
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	249 f.getFine_middle() +
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	250 f.getSand() +
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	251 f.getSusp_sand() +
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	252 f.getSusp_sediment();
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	253 load.setTotal(km, total);
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	254 }
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 return load;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	257 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	258
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	259 private SedimentLoad calculateUnit(SedimentLoad load) {
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	260 //TODO implement me!
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	261 return load;
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	262 }
1fb224bb2c6b Added sediment load calculation. Raimund Renkert <rrenkert@intevation.de> parents: diff changeset	263 }
4483 cc6323401643 Cosmetics: Removed some WSP, obsolete imports. Documentation added. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 4373 diff changeset	264 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf-8 :

Mercurial > dive4elements > river

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/minfo/SedimentLoadCalculation.java @ 4483:cc6323401643