dive4elements/river: artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/SedimentLoadDataCalculation.java annotate

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/SedimentLoadDataCalculation.java @ 8084:640fffb0ad60

Removed same package import.

author	Sascha L. Teichmann <teichmann@intevation.de>
date	Mon, 11 Aug 2014 16:45:42 +0200
parents	9ecd6267323b
children	bbad52b073a4

rev	line source
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	1 /* Copyright (C) 2014 by Bundesanstalt für Gewässerkunde
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	2 * Software engineering by Intevation GmbH
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	3 *
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	4 * This file is Free Software under the GNU AGPL (>=v3)
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	5 * and comes with ABSOLUTELY NO WARRANTY! Check out the
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	6 * documentation coming with Dive4Elements River for details.
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	7 */
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	8 package org.dive4elements.river.artifacts.model.minfo;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	9
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	10 import java.util.ArrayList;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	11 import java.util.List;
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	12 import java.util.Map;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	13 import java.util.Set;
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	14 import java.util.TreeMap;
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	15 import java.util.TreeSet;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	16
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	17 import org.dive4elements.river.artifacts.access.SedimentLoadAccess;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	18 import org.dive4elements.river.artifacts.model.Calculation;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	19 import org.dive4elements.river.artifacts.model.CalculationResult;
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	20 import org.dive4elements.river.artifacts.model.RiverFactory;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	21 import org.apache.log4j.Logger;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	22 import org.dive4elements.river.artifacts.model.minfo.SedimentLoadData.Value;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	23 import org.dive4elements.river.artifacts.model.minfo.SedimentLoadData.Station;
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	24 import org.dive4elements.river.artifacts.model.minfo.SedimentLoadDataValueFilter.And;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	25 import org.dive4elements.river.artifacts.model.minfo.SedimentLoadDataValueFilter.IsEpoch;
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	26 import org.dive4elements.river.artifacts.model.minfo.SedimentLoadDataValueFilter.IsOfficial;
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	27 import org.dive4elements.river.artifacts.model.minfo.SedimentLoadDataValueFilter.Not;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	28 import org.dive4elements.river.artifacts.model.minfo.SedimentLoadDataValueFilter.TimeRangeIntersects;
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	29 import org.dive4elements.river.model.River;
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	30 import org.dive4elements.river.utils.DoubleUtil;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	31
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	32 public class SedimentLoadDataCalculation
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	33 extends Calculation
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	34 {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	35 private static final Logger log = Logger
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	36 .getLogger(SedimentLoadDataCalculation.class);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	37
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	38 public static final int [] TOTAL_LOAD_FLYS = {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	39 SedimentLoadData.GF_COARSE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	40 SedimentLoadData.GF_FINE_MIDDLE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	41 SedimentLoadData.GF_SAND,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	42 SedimentLoadData.GF_SUSP_SEDIMENT
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	43 };
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	44
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	45 public static final int [] BED_LOAD_FLYS = {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	46 SedimentLoadData.GF_COARSE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	47 SedimentLoadData.GF_FINE_MIDDLE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	48 SedimentLoadData.GF_SAND
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	49 };
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	50
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	51 public static final int [] BED_LOAD_SUSP_SAND_FLYS = {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	52 SedimentLoadData.GF_COARSE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	53 SedimentLoadData.GF_FINE_MIDDLE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	54 SedimentLoadData.GF_SAND,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	55 SedimentLoadData.GF_SUSP_SAND
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	56 };
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	57
8050 9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	58 public static final int [] TOTAL_LOAD_BFG = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	59 SedimentLoadData.GF_TOTAL
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	60 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	61
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	62 public static final int [] BED_LOAD_BFG = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	63 SedimentLoadData.GF_BED_LOAD
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	64 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	65
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	66 public static final int [] SUSPENDED_LOAD_BFG = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	67 SedimentLoadData.GF_SUSPENDED_LOAD
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	68 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	69
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	70 public static final int [] COARSE_FLYS = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	71 SedimentLoadData.GF_COARSE
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	72 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	73
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	74 public static final int [] FINE_MIDDLE_FLYS = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	75 SedimentLoadData.GF_FINE_MIDDLE
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	76 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	77
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	78 public static final int [] SAND_FLYS = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	79 SedimentLoadData.GF_SAND
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	80 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	81
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	82 public static final int [] SUSP_SAND_FLYS = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	83 SedimentLoadData.GF_SUSP_SAND
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	84 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	85
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	86 public static final int [] SUSP_SAND_BED_FLYS = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	87 SedimentLoadData.GF_SUSP_SAND_BED
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	88 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	89
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	90 public static final int [] SUSP_SEDIMENT_FLYS = {
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	91 SedimentLoadData.GF_SUSP_SEDIMENT
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	92 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	93
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	94 public static final class GrainFraction {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	95 private String description;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	96 private int [] grainFractions;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	97
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	98 public GrainFraction(String description, int [] grainFractions) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	99 this.description = description;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	100 this.grainFractions = grainFractions;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	101 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	102 public static final GrainFraction make(String description, int [] grainFractions) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	103 return new GrainFraction(description, grainFractions);
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	104 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	105
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	106 public String getDescription() {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	107 return description;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	108 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	109
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	110 public int [] getGrainFractions() {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	111 return grainFractions;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	112 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	113 } // class GrainFraction
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	114
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	115 public static final GrainFraction [] GRAIN_FRACTIONS = {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	116 // TODO: i18n
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	117 GrainFraction.make("minfo.total.load.flys", TOTAL_LOAD_FLYS),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	118 GrainFraction.make("minfo.bed.load.flys", BED_LOAD_FLYS),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	119 GrainFraction.make("minfo.bed.load.susp.sand.flys", BED_LOAD_SUSP_SAND_FLYS),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	120 GrainFraction.make("minfo.total.load.bfg", TOTAL_LOAD_BFG),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	121 GrainFraction.make("minfo.bed.load.bfg", BED_LOAD_BFG),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	122 GrainFraction.make("minfo.suspended.load.bfg", SUSPENDED_LOAD_BFG),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	123 GrainFraction.make("minfo.coarse.flys", COARSE_FLYS),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	124 GrainFraction.make("minfo.fine.middle.flys", FINE_MIDDLE_FLYS),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	125 GrainFraction.make("minfo.sand.flys", SAND_FLYS) ,
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	126 GrainFraction.make("minfo.susp.sand.flys", SUSP_SAND_FLYS),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	127 GrainFraction.make("minfo.susp.sand.bed.flys", SUSP_SAND_BED_FLYS),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	128 GrainFraction.make("minfo.susp.sediment.flys", SUSP_SEDIMENT_FLYS),
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	129 };
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	130
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	131 public static class Sum implements Value.Visitor {
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	132
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	133 protected int n;
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	134 protected double sum;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	135
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	136 public Sum() {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	137 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	138
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	139 public double getSum() {
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	140 return sum;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	141 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	142
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	143 public int getN() {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	144 return n;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	145 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	146
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	147 public void reset() {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	148 n = 0;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	149 sum = 0.0;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	150 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	151
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	152 @Override
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	153 public void visit(Value value) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	154 sum += value.getValue();
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	155 ++n;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	156 }
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	157 } // class Sum
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	158
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	159 private String river;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	160 private String yearEpoch;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	161 private String unit;
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	162 private int [][] epochs;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	163 private int [] years;
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	164 private double from;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	165 private double to;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	166
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	167
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	168 public SedimentLoadDataCalculation() {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	169 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	170
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	171 public CalculationResult calculate(SedimentLoadAccess access) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	172 log.info("SedimentLoadDataCalculation.calculate");
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	173
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	174 String river = access.getRiverName();
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	175 String yearEpoch = access.getYearEpoch();
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	176 String unit = access.getUnit();
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	177
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	178 int [] years = null;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	179 int [][] epochs = null;
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	180
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	181 double from = access.getLowerKM();
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	182 double to = access.getUpperKM();
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	183
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	184 if (yearEpoch.equals("year")) {
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	185 years = access.getPeriod();
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	186 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	187 else if (yearEpoch.equals("epoch") \|\| yearEpoch.equals("off_epoch")) {
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	188 epochs = access.getEpochs();
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	189 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	190 else {
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	191 addProblem("minfo.missing.year_epoch");
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	192 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	193
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	194 if (river == null) {
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	195 // TODO: i18n
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	196 addProblem("minfo.missing.river");
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	197 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	198
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	199 if (years == null && epochs == null) {
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	200 addProblem("minfo.missing.time");
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	201 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	202
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	203 if (!hasProblems()) {
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	204 this.river = river;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	205 this.yearEpoch = yearEpoch;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	206 this.unit = unit;
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	207 this.years = years;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	208 this.epochs = epochs;
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	209 this.from = from;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	210 this.to = to;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	211 return internalCalculate();
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	212 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	213
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	214 return error(null);
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	215 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	216
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	217 private CalculationResult error(String msg) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	218 if (msg != null) addProblem(msg);
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	219 return new CalculationResult(this);
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	220 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	221
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	222 private CalculationResult internalCalculate() {
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	223 if ("year".equals(yearEpoch)) return calculateYears();
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	224 if ("epoch".equals(yearEpoch)) return calculateEpochs();
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	225 if ("off_epoch".equals(yearEpoch)) return calculateOffEpochs();
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	226
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	227 // TODO: i18n
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	228 return error("minfo.sediment.load.unknown.calc.mode");
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	229 }
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	230
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	231 private CalculationResult calculateYears() {
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	232 SedimentLoadData sld =
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	233 SedimentLoadDataFactory.INSTANCE.getSedimentLoadData(river);
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	234 if (sld == null) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	235 // TODO: i18n
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	236 return error("minfo.sediment.load.no.data");
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	237 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	238
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	239 SedimentLoadDataResult sldr = new SedimentLoadDataResult();
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	240
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	241 boolean isKmUp = isKmUp();
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	242 Set<Integer> missingFractions = new TreeSet<Integer>();
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	243
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	244 Not notEpochs = new Not(IsEpoch.INSTANCE);
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	245
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	246 Sum sum = new Sum();
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	247
8066 fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	248 SedimentDensity sd = getSedimentDensity();
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	249
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	250 for (int year: years) {
8068 9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	251 Value.Filter filter = new And(notEpochs)
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	252 .add(new TimeRangeIntersects(year));
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	253
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	254 for (GrainFraction gf: GRAIN_FRACTIONS) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	255 double [][] result = sum(
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	256 sld, gf.getGrainFractions(), filter, sum, isKmUp,
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	257 missingFractions);
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	258
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	259 if (result[0].length == 0 \|\| DoubleUtil.isNaN(result[1])) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	260 // TODO: resolve i18n
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	261 addProblem("minfo.sediment.load.no.fractions",
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	262 gf.getDescription());
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	263 continue;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	264 }
8066 fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	265
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	266 transformT2M3(sd, year, result);
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	267
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	268 SedimentLoadDataResult.Fraction sldrf =
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	269 new SedimentLoadDataResult.Fraction(gf.getDescription(), result);
8066 fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	270
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	271 sldr.addFraction(sldrf);
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	272 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	273 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	274 // TODO: Generate messages for missing fractions.
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	275 return new CalculationResult(sldr, this);
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	276 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	277
8066 fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	278
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	279 private CalculationResult calculateEpochs() {
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	280 SedimentLoadData sld =
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	281 SedimentLoadDataFactory.INSTANCE.getSedimentLoadData(river);
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	282 if (sld == null) {
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	283 // TODO: i18n
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	284 return error("minfo.sediment.load.no.data");
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	285 }
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	286
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	287 SedimentLoadDataResult sldr = new SedimentLoadDataResult();
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	288
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	289 boolean isKmUp = isKmUp();
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	290 Set<Integer> missingFractions = new TreeSet<Integer>();
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	291
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	292 SedimentDensity sd = getSedimentDensity();
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	293
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	294 // They are not epochs, they are single years!
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	295 Not notEpochs = new Not(IsEpoch.INSTANCE);
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	296
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	297 for (int [] epoch: epochs) {
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	298 List<double [][]> results = new ArrayList<double [][]>();
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	299
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	300 int min = Math.min(epoch[0], epoch[1]);
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	301 int max = Math.max(epoch[0], epoch[1]);
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	302
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	303 for (int year = min; year <= max; ++year) {
8068 9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	304 Value.Filter filter = new And(notEpochs)
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	305 .add(new TimeRangeIntersects(year));
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	306
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	307 Sum sum = new Sum();
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	308
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	309 for (GrainFraction gf: GRAIN_FRACTIONS) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	310 double [][] result = sum(
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	311 sld, gf.getGrainFractions(), filter, sum, isKmUp,
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	312 missingFractions);
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	313
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	314 if (result[0].length == 0 \|\| DoubleUtil.isNaN(result[1])) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	315 // TODO: resolve i18n
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	316 addProblem("minfo.sediment.load.no.fractions",
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	317 gf.getDescription());
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	318 continue;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	319 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	320
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	321 transformT2M3(sd, year, result);
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	322 results.add(result);
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	323 }
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	324 }
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	325
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	326 double [][] result = average(results);
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	327 // TODO: Optionally transform units.
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	328 SedimentLoadDataResult.Fraction sldrf =
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	329 new SedimentLoadDataResult.Fraction("TODO: nice description", result);
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	330 sldr.addFraction(sldrf);
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	331 }
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	332 // TODO: Generate messages for missing fractions.
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	333 return new CalculationResult(sldr, this);
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	334 }
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	335
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	336 private CalculationResult calculateOffEpochs() {
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	337 SedimentLoadData sld =
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	338 SedimentLoadDataFactory.INSTANCE.getSedimentLoadData(river);
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	339 if (sld == null) {
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	340 // TODO: i18n
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	341 return error("minfo.sediment.load.no.data");
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	342 }
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	343
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	344 SedimentLoadDataResult sldr = new SedimentLoadDataResult();
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	345
8068 9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	346 SedimentDensity sd = getSedimentDensity();
9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	347
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	348 boolean isKmUp = isKmUp();
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	349 Set<Integer> missingFractions = new TreeSet<Integer>();
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	350
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	351 for (int [] epoch: epochs) {
8068 9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	352 Value.Filter filter = new And(IsOfficial.INSTANCE)
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	353 .add(new TimeRangeIntersects(epoch[0], epoch[1]));
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	354
8068 9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	355 int year = Math.min(epoch[0], epoch[1]);
9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	356
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	357 Sum sum = new Sum();
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	358
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	359 for (GrainFraction gf: GRAIN_FRACTIONS) {
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	360 double [][] result = sum(
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	361 sld, gf.getGrainFractions(), filter, sum, isKmUp,
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	362 missingFractions);
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	363
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	364 if (result[0].length == 0 \|\| DoubleUtil.isNaN(result[1])) {
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	365 // TODO: resolve i18n
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	366 addProblem("minfo.sediment.load.no.fractions",
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	367 gf.getDescription());
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	368 continue;
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	369 }
8068 9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	370 transformT2M3(sd, year, result);
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	371 SedimentLoadDataResult.Fraction sldrf =
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	372 new SedimentLoadDataResult.Fraction(gf.getDescription(), result);
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	373 sldr.addFraction(sldrf);
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	374 }
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	375 }
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	376 // TODO: Generate messages for missing fractions.
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	377 return new CalculationResult(sldr, this);
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	378 }
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	379
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	380 /** Figure out flow direction of river. */
8068 9ecd6267323b Sediment load: transform units for official epochs, too. Sascha L. Teichmann <teichmann@intevation.de> parents: 8067 diff changeset	381 private final boolean isKmUp() {
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	382 River r = RiverFactory.getRiver(river);
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	383 if (r == null) {
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	384 addProblem("minfo.missing.river");
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	385 return true;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	386 }
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	387 return r.getKmUp();
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	388 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	389
8066 fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	390 private final boolean inM3() {
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	391 return unit.equals("m3_per_a");
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	392 }
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	393
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	394 private SedimentDensity getSedimentDensity() {
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	395 return inM3()
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	396 ? SedimentDensityFactory.getSedimentDensity(river, from, to)
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	397 : null;
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	398 }
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	399
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	400 private static void transformT2M3(SedimentDensity sd, int year, double [][] data) {
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	401 if (sd == null) {
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	402 return;
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	403 }
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	404 double [] kms = data[0];
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	405 double [] values = data[1];
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	406 for (int i = 0; i < kms.length; ++i) {
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	407 if (Double.isNaN(kms[i]) \|\| Double.isNaN(kms[i])) {
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	408 continue;
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	409 }
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	410 double density = sd.getDensity(kms[i], year);
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	411 values[i] /= density;
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	412 }
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	413 }
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	414
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	415 public double[][] sum(
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	416 SedimentLoadData sld,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	417 int [] grainFractions,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	418 Value.Filter filter,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	419 Sum sum,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	420 boolean isKMUp,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	421 Set<Integer> missingFractions
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	422 ) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	423 List<Station> stations = sld.findStations(from, to);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	424
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	425 double [] values = new double[grainFractions.length];
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	426
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	427 double [][] result = new double[2][stations.size()];
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	428
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	429 for (int j = 0, S = stations.size(); j < S; ++j) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	430 Station station = stations.get(j);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	431 for (int i = 0; i < grainFractions.length; ++i) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	432 int gf = grainFractions[i];
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	433 sum.reset();
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	434 station.filterGrainFraction(gf, filter, sum);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	435 if (sum.getN() == 0) { // No values found
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	436 int msType = SedimentLoadData.measurementStationType(gf);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	437 // Station of right fraction type already? No: take previous.
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	438 if (!station.isType(msType)) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	439 Station prev = station.prevByType(msType, isKMUp);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	440 if (prev != null) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	441 prev.filterGrainFraction(gf, filter, sum);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	442 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	443 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	444 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	445
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	446 if (sum.getN() == 0) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	447 missingFractions.add(gf);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	448 values[i] = Double.NaN;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	449 } else {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	450 values[i] = sum.getSum();
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	451 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	452 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	453 result[0][j] = station.getStation();
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	454 result[1][j] = DoubleUtil.sum(values);
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	455 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	456
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	457 // TODO: Handle 'virtual' measument stations 'from' and 'to'.
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	458
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	459 return result;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	460 }
8067 6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	461
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	462 private static final class XSum {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	463 private double sum;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	464 private int n;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	465 public XSum() {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	466 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	467 public void add(double v) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	468 sum += v;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	469 ++n;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	470 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	471 public double avg() {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	472 return sum/n;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	473 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	474 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	475
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	476 private static double [][] average(List<double [][]> data) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	477
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	478 TreeMap<Double, XSum> map = new TreeMap<Double, XSum>();
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	479
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	480 for (double [][] pair: data) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	481 double [] kms = pair[0];
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	482 double [] vs = pair[1];
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	483 for (int i = 0; i < kms.length; ++i) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	484 double km = kms[i];
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	485 double v = vs[i];
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	486 if (Double.isNaN(km) \|\| Double.isNaN(v)) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	487 continue;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	488 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	489 XSum xsum = map.get(km);
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	490 if (xsum == null) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	491 map.put(km, xsum = new XSum());
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	492 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	493 xsum.add(v);
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	494 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	495 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	496
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	497 double [][] result = new double[2][map.size()];
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	498 int i = 0;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	499 for (Map.Entry<Double, XSum> entry: map.entrySet()) {
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	500 result[i][0] = entry.getKey();
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	501 result[i][1] = entry.getValue().avg();
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	502 ++i;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	503 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	504
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	505 return null;
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	506 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	507
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	508 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	509 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf-8 :

Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/artifacts/model/minfo/SedimentLoadDataCalculation.java @ 8084:640fffb0ad60