dive4elements/river: artifacts/src/main/java/org/dive4elements/river/exports/fixings/FixATWriter.java annotate

annotate artifacts/src/main/java/org/dive4elements/river/exports/fixings/FixATWriter.java @ 9801:1d7a72a50183 3.2.x tip

Assume Compose V2, consistently

author	Tom Gottfried <tom@intevation.de>
date	Thu, 23 Nov 2023 10:14:13 +0100
parents	0a5239a1e46e
children

rev	line source
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	1 /* Copyright (C) 2011, 2012, 2013 by Bundesanstalt für Gewässerkunde
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	2 * Software engineering by Intevation GmbH
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	3 *
5994 af13ceeba52a Removed trailing whitespace. Sascha L. Teichmann <teichmann@intevation.de> parents: 5863 diff changeset	4 * This file is Free Software under the GNU AGPL (>=v3)
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	5 * and comes with ABSOLUTELY NO WARRANTY! Check out the
5994 af13ceeba52a Removed trailing whitespace. Sascha L. Teichmann <teichmann@intevation.de> parents: 5863 diff changeset	6 * documentation coming with Dive4Elements River for details.
5863 4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	7 */
4897a58c8746 River artifacts: Added new copyright headers. Sascha L. Teichmann <teichmann@intevation.de> parents: 5838 diff changeset	8
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	9 package org.dive4elements.river.exports.fixings;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	11 import org.dive4elements.artifacts.CallMeta;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	13 import org.dive4elements.river.artifacts.math.fitting.Function;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	15 import org.dive4elements.river.artifacts.model.Parameters;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	17 import org.dive4elements.river.artifacts.resources.Resources;
bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	18
bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	19 import org.dive4elements.river.exports.ATWriter;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20
6325 ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	21 import org.dive4elements.river.model.Gauge;
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	22 import org.dive4elements.river.model.River;
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	23
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	24 import java.io.IOException;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25 import java.io.PrintWriter;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26 import java.io.Writer;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 import java.util.Locale;
7687 f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	29 import java.util.Arrays;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30
9726 0a5239a1e46e Upgrade to Log4j 2 Tom Gottfried <tom@intevation.de> parents: 8724 diff changeset	31 import org.apache.logging.log4j.Logger;
0a5239a1e46e Upgrade to Log4j 2 Tom Gottfried <tom@intevation.de> parents: 8724 diff changeset	32 import org.apache.logging.log4j.LogManager;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33
3905 8bbb9e173297 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 3643 diff changeset	34 /** Export Fixation Analysis Results to AT. */
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35 public class FixATWriter
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36 {
8202 e4606eae8ea5 sed src/*/.java 's/logger/log/g' Sascha L. Teichmann <teichmann@intevation.de> parents: 7687 diff changeset	37 /** Private log. */
9726 0a5239a1e46e Upgrade to Log4j 2 Tom Gottfried <tom@intevation.de> parents: 8724 diff changeset	38 private static Logger log = LogManager.getLogger(FixATWriter.class);
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40 public static final String I18N_HEADER_KEY =
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 "fix.export.at.header";
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42
6325 ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	43 public static final String I18N_GAUGE_HEADER_KEY =
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	44 "fix.export.at.gauge.header";
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	45
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46 public static final String I18N_HEADER_DEFAULT =
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47 "Exported fixings discharge curve for {0} {0}-km: {1}";
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48
6325 ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	49 public static final String I18N_GAUGE_HEADER_DEFAULT =
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	50 "Exported fixings discharge curve for {0}, gauge: {1} datum[{3}] = {2}";
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	51
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 public static final String [] Q_MAX_COLUMN = new String [] { "max_q" };
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	54 private static final int MAX_ITERATIONS = 10000;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	55 private static final double EPSILON = 1e-8;
3643 171db4d5d3cb FixA: AT export: prevent some numerical problems with steep functions around zero. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3642 diff changeset	56 private static final double MIN_Q = 1e-4;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58 protected Function function;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 protected Parameters parameters;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61 public FixATWriter() {
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62 }
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 public FixATWriter(Function function, Parameters parameters) {
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 this.function = function;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66 this.parameters = parameters;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 }
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 public void write(
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 Writer writer,
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 CallMeta meta,
6325 ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	72 River river,
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 double km
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74 )
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 throws IOException {
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 PrintWriter out = new PrintWriter(writer);
6325 ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	77
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	78 int subtractPNP = 0;
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	79 // Special case handling for at's at gauges
8724 47199406994a (issue1801) Determine gauge at a station always with same tolerance. Tom Gottfried <tom@intevation.de> parents: 8202 diff changeset	80 Gauge gauge = river.determineGaugeAtStation(km);
47199406994a (issue1801) Determine gauge at a station always with same tolerance. Tom Gottfried <tom@intevation.de> parents: 8202 diff changeset	81 if (gauge != null) {
6325 ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	82 printGaugeHeader(out, meta, river, gauge);
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	83 subtractPNP = (int)Math.round(gauge.getDatum().doubleValue() * 100);
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	84 } else {
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	85 printHeader(out, meta, river.getName(), km);
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	86 }
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 double [] coeffs = parameters.interpolate(
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 "km", km, function.getParameterNames());
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 double [] qMax = parameters.interpolate(
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 "km", km, Q_MAX_COLUMN);
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 if (coeffs == null \|\| qMax == null) {
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 log.debug("No data found at km " + km + ".");
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96 return;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 }
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	99 org.dive4elements.river.artifacts.math.Function funcInst =
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	100 function.instantiate(coeffs);
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101
3642 83c0735092a3 FixA: Increase Q max about 5% if storing as AT file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3641 diff changeset	102 // Increase Q max about 5%.
83c0735092a3 FixA: Increase Q max about 5% if storing as AT file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3641 diff changeset	103 qMax[0] += Math.abs(qMax[0])*0.05;
83c0735092a3 FixA: Increase Q max about 5% if storing as AT file. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3641 diff changeset	104
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 double wMax = funcInst.value(qMax[0]);
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	107 if (Double.isNaN(wMax) \|\| wMax < 0d) {
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	108 log.debug("function '" + function.getName() +
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	109 "' eval failed at " + wMax);
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	110 return;
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	111 }
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 Function inverse = function.getInverse();
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	115 org.dive4elements.river.artifacts.math.Function invInst =
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	116 inverse.instantiate(coeffs);
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	118 double wMin = minW(invInst, wMax, qMax[0]);
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119
6637 65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155) Andre Heinecke <aheinecke@intevation.de> parents: 6325 diff changeset	120 double wMinCM = wMin * 100d - subtractPNP;
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	121 double wMaxCM = wMax * 100d;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	123 if ((wMinCM - (int)wMinCM) > 0d) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	124 wMinCM = (int)wMinCM + 1d;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	125 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	126
7687 f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	127 int wRow = ((int)wMinCM / 10) * 10;
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	128
6637 65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155) Andre Heinecke <aheinecke@intevation.de> parents: 6325 diff changeset	129 double w = (wMinCM + subtractPNP) / 100.0;
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	130
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	131 int wcm = ((int)wMinCM) % 10;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	132
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	133 if (log.isDebugEnabled()) {
7687 f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	134 log.debug("km: " + km);
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	135 log.debug("wMinCM: " + wMinCM);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	136 log.debug("wMaxCM: " + wMaxCM);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	137 log.debug("wcm: " + wcm);
6637 65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155) Andre Heinecke <aheinecke@intevation.de> parents: 6325 diff changeset	138 log.debug("subtractPNP: " + subtractPNP);
7687 f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	139 log.debug("coeffs: " + Arrays.toString(coeffs));
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	140 log.debug("function description: " + inverse.getDescription());
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	141 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	142
6637 65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155) Andre Heinecke <aheinecke@intevation.de> parents: 6325 diff changeset	143 out.printf(Locale.US, "%8d", wRow);
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	144
6209 1642e08d14d2 Fix empty column logic in FixAT export Andre Heinecke <aheinecke@intevation.de> parents: 6208 diff changeset	145 for (int i = 0; i < wcm; i++) {
1642e08d14d2 Fix empty column logic in FixAT export Andre Heinecke <aheinecke@intevation.de> parents: 6208 diff changeset	146 out.print(ATWriter.EMPTY);
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 }
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	149 for (;;) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	150 while (wcm++ < 10) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	151 if (w > wMax) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	152 break;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	153 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	154 double q = invInst.value(w);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	155 if (Double.isNaN(w)) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	156 out.print(ATWriter.EMPTY);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	157 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	158 else {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	159 ATWriter.printQ(out, q);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	160 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	161 w += 0.01d;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	162 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	163 out.println();
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	164 if (w > wMax) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	165 break;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	166 }
6637 65a5deaee769 Subtract PNP before calculating wcm from wMinCM (issue1155) Andre Heinecke <aheinecke@intevation.de> parents: 6325 diff changeset	167 out.printf(Locale.US, "%8d", (wRow += 10));
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	168 wcm = 0;
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 }
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171 out.flush();
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 }
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 protected void printHeader(
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 PrintWriter out,
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 CallMeta meta,
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177 String river,
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178 double km
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 ) {
6208 715e010d40de Add * in the first line of AT Exports. Andre Heinecke <aheinecke@intevation.de> parents: 5994 diff changeset	180 out.println("*" + Resources.format(
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181 meta,
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 I18N_HEADER_KEY,
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 I18N_HEADER_DEFAULT,
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 river, km));
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 }
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	186
6325 ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	187 protected void printGaugeHeader(
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	188 PrintWriter out,
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	189 CallMeta meta,
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	190 River river,
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	191 Gauge gauge
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	192 ) {
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	193 out.println("*" + Resources.format(
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	194 meta,
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	195 I18N_GAUGE_HEADER_KEY,
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	196 I18N_GAUGE_HEADER_DEFAULT,
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	197 new Object[] { river.getName(), gauge.getName(),
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	198 gauge.getDatum(), river.getWstUnit().getName() }));
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	199 }
ee701de9c875 Add special case handling for Fixing AT export at a gauge Andre Heinecke <aheinecke@intevation.de> parents: 6209 diff changeset	200
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	201 private static double minW(
5831 bd047b71ab37 Repaired internal references Sascha L. Teichmann <teichmann@intevation.de> parents: 3905 diff changeset	202 org.dive4elements.river.artifacts.math.Function function,
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	203 double maxW,
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	204 double maxQ
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	205 ) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	206 double stepWidth = 10d;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	207
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	208 double lastW = maxW;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	209 double lastQ = maxQ;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	210
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	211 for (int i = 0; i < MAX_ITERATIONS; ++i) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	212 double w = lastW - stepWidth;
7687 f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	213
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	214 if (w <= 0) {
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	215 return 0;
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	216 }
f4b3861712a4 (issue1666) Calculate wRow after rounding the first value. Andre Heinecke <aheinecke@intevation.de> parents: 6637 diff changeset	217
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	218 double q = function.value(w);
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	219
3643 171db4d5d3cb FixA: AT export: prevent some numerical problems with steep functions around zero. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3642 diff changeset	220 if (Double.isNaN(q) \|\| q > lastQ \|\| q < MIN_Q) {
3641 1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	221 if (stepWidth < EPSILON) {
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	222 break;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	223 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	224 stepWidth *= 0.5d;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	225 continue;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	226 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	227
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	228 lastW = w;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	229 lastQ = q;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	230 }
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	231
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	232 return lastW;
1dc904370a64 Fixed issue687: Rewrote AT export in FixA Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3220 diff changeset	233 }
3220 0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	234 }
0c8a6145098b FixA: Added AT writer for fitted functions. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	235 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :

Mercurial > dive4elements > river

annotate artifacts/src/main/java/org/dive4elements/river/exports/fixings/FixATWriter.java @ 9801:1d7a72a50183 3.2.x tip