Mercurial > dive4elements > river
view flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/Calculation5.java @ 4798:39885bdfc6fc
Added calculation of the "Umhuellende" to calculation of "W fuer ungleichwertige Abfluesse".
This is done by figuring out the WST columns that imfold the data and then do
simple "gleichwertige" calculations from the start of the interval. This
is too much because only the Qs are needed for the "Umhuellende".
author | Sascha L. Teichmann <teichmann@intevation.de> |
---|---|
date | Sun, 13 Jan 2013 16:18:28 +0100 |
parents | 5642a83420f2 |
children |
line wrap: on
line source
package de.intevation.flys.artifacts.model; import de.intevation.artifacts.CallContext; import de.intevation.artifacts.CallMeta; import de.intevation.flys.utils.Formatter; import java.text.NumberFormat; import java.util.ArrayList; import java.util.Map; import org.apache.log4j.Logger; public class Calculation5 extends Calculation { private static Logger log = Logger.getLogger(Calculation5.class); protected double startKm; protected double [] endKms; public Calculation5() { } public Calculation5(double startKm, double [] endKms) { this.startKm = startKm; this.endKms = endKms; } public CalculationResult calculate( WstValueTable wst, Map<Double, Double> kms2gaugeDatums, CallContext context ) { ArrayList<WWQQ> results = new ArrayList<WWQQ>(); int numProblems = numProblems(); CallMeta meta = context.getMeta(); for (double endKm: endKms) { double [][] wws = wst.relateWs(startKm, endKm, this); int newNumProblems = numProblems(); if (wws.length == 4) { WWQQ wwqq = new WWQQ( generateName(meta, startKm, endKm), startKm, kms2gaugeDatums.get(startKm), wws[0], wws[1], endKm, kms2gaugeDatums.get(endKm), wws[2], wws[3]); if (newNumProblems > numProblems) { wwqq.removeNaNs(); } results.add(wwqq); } numProblems = newNumProblems; } return new CalculationResult( results.toArray(new WWQQ[results.size()]), this); } protected static String generateName( CallMeta meta, double startKm, double endKm ) { NumberFormat nf = Formatter.getCalculationKm(meta); return "W(km " + nf.format(startKm) + ") ~ W(km " + nf.format(endKm) + ")"; } } // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :