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 @ 8067:6d24ba2ac964

Sediment load: refactored epoch calculation to not average over all but per year.

author	Sascha L. Teichmann <teichmann@intevation.de>
date	Thu, 31 Jul 2014 17:04:41 +0200
parents	fe5ef780f8b1
children	9ecd6267323b

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;
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.Value;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	24 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	25 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	26 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	27 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	28 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	29 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	30 import org.dive4elements.river.model.River;
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	31 import org.dive4elements.river.utils.DoubleUtil;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	32
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	33 public class SedimentLoadDataCalculation
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	34 extends Calculation
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	35 {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	36 private static final Logger log = Logger
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	37 .getLogger(SedimentLoadDataCalculation.class);
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	38
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	39 public static final int [] TOTAL_LOAD_FLYS = {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	40 SedimentLoadData.GF_COARSE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	41 SedimentLoadData.GF_FINE_MIDDLE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	42 SedimentLoadData.GF_SAND,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	43 SedimentLoadData.GF_SUSP_SEDIMENT
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
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	46 public static final int [] BED_LOAD_FLYS = {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	47 SedimentLoadData.GF_COARSE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	48 SedimentLoadData.GF_FINE_MIDDLE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	49 SedimentLoadData.GF_SAND
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
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	52 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	53 SedimentLoadData.GF_COARSE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	54 SedimentLoadData.GF_FINE_MIDDLE,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	55 SedimentLoadData.GF_SAND,
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	56 SedimentLoadData.GF_SUSP_SAND
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	57 };
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	58
8050 9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	59 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	60 SedimentLoadData.GF_TOTAL
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
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	63 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	64 SedimentLoadData.GF_BED_LOAD
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
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	67 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	68 SedimentLoadData.GF_SUSPENDED_LOAD
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
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	71 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	72 SedimentLoadData.GF_COARSE
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
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	75 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	76 SedimentLoadData.GF_FINE_MIDDLE
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
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	79 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	80 SedimentLoadData.GF_SAND
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
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	83 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	84 SedimentLoadData.GF_SUSP_SAND
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
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	87 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	88 SedimentLoadData.GF_SUSP_SAND_BED
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
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	91 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	92 SedimentLoadData.GF_SUSP_SEDIMENT
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	93 };
9e79e384aa8b Sediment load: Added more 'sums' to be calculated. Sascha L. Teichmann <teichmann@intevation.de> parents: 8049 diff changeset	94
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	95 public static final class GrainFraction {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	96 private String description;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	97 private int [] grainFractions;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	98
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	99 public GrainFraction(String description, int [] grainFractions) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	100 this.description = description;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	101 this.grainFractions = grainFractions;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	102 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	103 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	104 return new GrainFraction(description, grainFractions);
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
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	107 public String getDescription() {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	108 return description;
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
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	111 public int [] getGrainFractions() {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	112 return grainFractions;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	113 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	114 } // class GrainFraction
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	115
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	116 public static final GrainFraction [] GRAIN_FRACTIONS = {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	117 // TODO: i18n
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	118 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	119 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	120 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	121 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	122 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	123 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	124 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	125 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	126 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	127 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	128 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	129 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	130 };
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	131
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	132 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	133
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	134 protected int n;
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	135 protected double sum;
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	136
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	137 public Sum() {
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
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	140 public double getSum() {
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	141 return sum;
8049 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
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	144 public int getN() {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	145 return n;
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
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	148 public void reset() {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	149 n = 0;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	150 sum = 0.0;
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
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	153 @Override
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	154 public void visit(Value value) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	155 sum += value.getValue();
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	156 ++n;
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	157 }
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	158 } // class Sum
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	159
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	160 private String river;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	161 private String yearEpoch;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	162 private String unit;
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	163 private int [][] epochs;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	164 private int [] years;
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	165 private double from;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	166 private double to;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	167
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	168
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	169 public SedimentLoadDataCalculation() {
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
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	172 public CalculationResult calculate(SedimentLoadAccess access) {
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	173 log.info("SedimentLoadDataCalculation.calculate");
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	174
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	175 String river = access.getRiverName();
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	176 String yearEpoch = access.getYearEpoch();
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	177 String unit = access.getUnit();
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	178
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	179 int [] years = null;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	180 int [][] epochs = null;
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	181
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	182 double from = access.getLowerKM();
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	183 double to = access.getUpperKM();
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	184
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	185 if (yearEpoch.equals("year")) {
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	186 years = access.getPeriod();
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	187 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	188 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	189 epochs = access.getEpochs();
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	190 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	191 else {
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	192 addProblem("minfo.missing.year_epoch");
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
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	195 if (river == null) {
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	196 // TODO: i18n
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	197 addProblem("minfo.missing.river");
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	198 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	199
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	200 if (years == null && epochs == null) {
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	201 addProblem("minfo.missing.time");
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
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	204 if (!hasProblems()) {
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	205 this.river = river;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	206 this.yearEpoch = yearEpoch;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	207 this.unit = unit;
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	208 this.years = years;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	209 this.epochs = epochs;
8052 1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	210 this.from = from;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	211 this.to = to;
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	212 return internalCalculate();
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	213 }
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	214
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	215 return error(null);
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
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	218 private CalculationResult error(String msg) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	219 if (msg != null) addProblem(msg);
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	220 return new CalculationResult(this);
8052 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
1dae69eff79d Sediment load: take arguments from sediment access. Sascha L. Teichmann <teichmann@intevation.de> parents: 8050 diff changeset	223 private CalculationResult internalCalculate() {
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	224 if ("year".equals(yearEpoch)) return calculateYears();
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	225 if ("epoch".equals(yearEpoch)) return calculateEpochs();
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	226 if ("off_epoch".equals(yearEpoch)) return calculateOffEpochs();
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	227
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	228 // TODO: i18n
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	229 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	230 }
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	231
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	232 private CalculationResult calculateYears() {
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	233 SedimentLoadData sld =
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	234 SedimentLoadDataFactory.INSTANCE.getSedimentLoadData(river);
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	235 if (sld == null) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	236 // TODO: i18n
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	237 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	238 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	239
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	240 SedimentLoadDataResult sldr = new SedimentLoadDataResult();
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	241
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	242 boolean isKmUp = isKmUp();
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	243 Set<Integer> missingFractions = new TreeSet<Integer>();
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	244
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	245 Not notEpochs = new Not(IsEpoch.INSTANCE);
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	246
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	247 Sum sum = new Sum();
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	248
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	249 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	250
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	251 for (int year: years) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	252 Value.Filter filter = new And()
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	253 .add(notEpochs)
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	254 .add(new TimeRangeIntersects(year));
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	255
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	256 for (GrainFraction gf: GRAIN_FRACTIONS) {
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	257 double [][] result = sum(
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	258 sld, gf.getGrainFractions(), filter, sum, isKmUp,
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	259 missingFractions);
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	260
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	261 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	262 // TODO: resolve i18n
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	263 addProblem("minfo.sediment.load.no.fractions",
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	264 gf.getDescription());
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	265 continue;
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	266 }
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	267
fe5ef780f8b1 Sediment load: fetch sediment density for transforming t/a to m^3/a. Sascha L. Teichmann <teichmann@intevation.de> parents: 8062 diff changeset	268 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	269
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	270 SedimentLoadDataResult.Fraction sldrf =
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	271 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	272
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	273 sldr.addFraction(sldrf);
8055 cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	274 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	275 }
cd35b76f1ef8 Sediment load. More off year based calculations. Sascha L. Teichmann <teichmann@intevation.de> parents: 8053 diff changeset	276 // 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	277 return new CalculationResult(sldr, this);
8049 d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	278 }
d49846f05108 Sediment load: (incomplete) new sediment load calculation. Sascha L. Teichmann <teichmann@intevation.de> parents: diff changeset	279
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	280
8053 72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	281 private CalculationResult calculateEpochs() {
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	282 SedimentLoadData sld =
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	283 SedimentLoadDataFactory.INSTANCE.getSedimentLoadData(river);
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	284 if (sld == null) {
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	285 // TODO: i18n
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	286 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	287 }
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	288
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	289 SedimentLoadDataResult sldr = new SedimentLoadDataResult();
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	290
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	291 boolean isKmUp = isKmUp();
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	292 Set<Integer> missingFractions = new TreeSet<Integer>();
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	293
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	294 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	295
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	296 // 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	297 Not notEpochs = new Not(IsEpoch.INSTANCE);
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	298
8057 555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	299 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	300 List<double [][]> results = new ArrayList<double [][]>();
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	301
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	302 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	303 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	304
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 for (int year = min; year <= max; ++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 Value.Filter filter = new And()
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 .add(notEpochs)
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 .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	309
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 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	311
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 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	313 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	314 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	315 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	316
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 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	318 // 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	319 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	320 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	321 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	322 }
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	323
6d24ba2ac964 Sediment load: refactored epoch calculation to not average over all but per year. Sascha L. Teichmann <teichmann@intevation.de> parents: 8066 diff changeset	324 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	325 results.add(result);
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	326 }
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	327 }
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	328
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 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	330 // 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	331 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	332 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	333 sldr.addFraction(sldrf);
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 // 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	336 return new CalculationResult(sldr, this);
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	337 }
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	338
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	339 private CalculationResult calculateOffEpochs() {
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	340 SedimentLoadData sld =
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	341 SedimentLoadDataFactory.INSTANCE.getSedimentLoadData(river);
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	342 if (sld == null) {
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	343 // TODO: i18n
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	344 return error("minfo.sediment.load.no.data");
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	345 }
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	346
8062 8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	347 SedimentLoadDataResult sldr = new SedimentLoadDataResult();
8489565ff563 Sedimen load: Added bundle for fraction results. Sascha L. Teichmann <teichmann@intevation.de> parents: 8060 diff changeset	348
8060 25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	349 boolean isKmUp = isKmUp();
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	350 Set<Integer> missingFractions = new TreeSet<Integer>();
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	351
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	352 for (int [] epoch: epochs) {
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	353 Value.Filter filter = new And()
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	354 .add(IsOfficial.INSTANCE)
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	355 .add(new TimeRangeIntersects(epoch[0], epoch[1]));
25feef564d09 Sediment load: More of official epochs. Sascha L. Teichmann <teichmann@intevation.de> parents: 8057 diff changeset	356
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 }
555dc5a9b282 Sediment load: More on calculating epoch loads. Sascha L. Teichmann <teichmann@intevation.de> parents: 8055 diff changeset	370 // TODO: Optionally transform units.
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. */
72760ca2fc2b Sediment load: dispatch calculation modes. Sascha L. Teichmann <teichmann@intevation.de> parents: 8052 diff changeset	381 private boolean isKmUp() {
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 @ 8067:6d24ba2ac964