ingo@696: package de.intevation.flys.artifacts.model;
ingo@696: 
sascha@2196: import de.intevation.artifactdatabase.state.Facet;
ingo@696: import de.intevation.artifacts.Artifact;
ingo@696: import de.intevation.artifacts.CallContext;
felix@2122: import de.intevation.flys.artifacts.FLYSArtifact;
christian@4045: import de.intevation.flys.artifacts.math.Linear;
sascha@697: import de.intevation.flys.artifacts.states.DefaultState.ComputeType;
ingo@696: 
sascha@2196: import org.apache.log4j.Logger;
felix@2017: 
felix@1721: /**
felix@1721:  * Facet of a Waterlevel (WQKms).
felix@1721:  */
felix@2122: public class WaterlevelFacet extends DataFacet {
ingo@696: 
ingo@696:     private static Logger logger = Logger.getLogger(WaterlevelFacet.class);
ingo@696: 
ingo@696:     public WaterlevelFacet(int index, String name, String description) {
felix@2122:         super(index, name, description, ComputeType.ADVANCE, null, null);
ingo@696:     }
ingo@696: 
sascha@722:     public WaterlevelFacet(
sascha@742:         int         index,
sascha@742:         String      name,
sascha@722:         String      description,
sascha@722:         ComputeType type,
sascha@722:         String      stateID,
sascha@722:         String      hash
sascha@722:     ) {
felix@2125:         super(index, name, description, type, hash, stateID);
felix@2122:     }
felix@2122: 
felix@2122:     public WaterlevelFacet() {
sascha@722:     }
ingo@696: 
sascha@3461:     protected WQKms [] getWQKms(CalculationResult res) {
sascha@3461:         return (WQKms [])res.getData();
sascha@3461:     }
felix@1136: 
felix@1136:     /**
felix@1136:      * Get waterlevel data.
felix@1136:      * @return a WQKms at given index.
felix@1136:      */
christian@4045:     @Override
ingo@696:     public Object getData(Artifact artifact, CallContext context) {
sascha@2196: 
sascha@2196:         if (logger.isDebugEnabled()) {
sascha@2196:             logger.debug("Get data for waterlevels at index: " + index +
sascha@2196:                 " /stateId: " + stateId);
sascha@2196:         }
ingo@696: 
christian@4045:         if (artifact == null) {
christian@4045:             logger.error("WaterlevelFacet.getData: artifact is null");
christian@4045:             return null;
christian@4045:         }
christian@4045: 
felix@3455:         FLYSArtifact winfo = (FLYSArtifact) artifact;
ingo@696: 
sascha@709:         CalculationResult res = (CalculationResult)
felix@2122:             winfo.compute(context, hash, stateId, type, false);
sascha@709: 
felix@3455:         if (res == null) {
felix@3455:             logger.error("WaterlevelFacet.getData: null result");
felix@3455:             return null;
felix@3455:         }
felix@3455: 
sascha@3461:         WQKms [] wqkms = getWQKms(res);
sascha@3597:         Object KM = context.getContextValue("currentKm");
sascha@3597:         if (KM != null) {
felix@3596:             logger.debug("interpolate at given km");
felix@3596:             // TODO handle exact match.
felix@3596: 
felix@3596:             WQKms wqkmsI = wqkms[index];
sascha@3597:             double km = (Double)KM;
felix@3596: 
felix@3596:             // TODO employ DataUtils interface to TDoubleArraList
felix@3596:             int size = wqkmsI.size();
sascha@3597:             boolean kmIncreasing = wqkmsI.getKm(0) < wqkmsI.getKm(size-1);
felix@3596:             int mod = kmIncreasing ? +1 : -1;
felix@3596:             int idx = 0;
felix@3596:             if (!kmIncreasing) {
felix@3596:                 while (idx < size && wqkmsI.getKm(idx) < km) {
felix@3596:                     idx++;
felix@3596:                 }
felix@3596:             }
felix@3596:             else {
felix@3596:                 idx = size-1;
felix@3596:                 while (idx > 0 && wqkmsI.getKm(idx) > km) {
felix@3596:                     idx--;
felix@3596:                 }
felix@3596:             }
christian@4045: 
felix@3596:             WQKms resultWQKms = new WQKms();
christian@4045:             if ((idx != -1) && (idx < size) && (idx - mod != -1) && (idx - mod < size)) {
sascha@3597:                 double inW = Linear.linear(
sascha@3597:                     km,
christian@4045:                     wqkmsI.getKm(idx), wqkmsI.getKm(idx - mod),
christian@4045:                     wqkmsI.getW(idx), wqkmsI.getW(idx - mod));
sascha@3597:                 double inQ = Linear.linear(
sascha@3597:                     km,
christian@4045:                     wqkmsI.getKm(idx), wqkmsI.getKm(idx - mod),
christian@4045:                     wqkmsI.getQ(idx), wqkmsI.getQ(idx - mod));
felix@3596:                 resultWQKms.add(inW, inQ, km);
felix@3596:             }
christian@4045: 
felix@3596:             return resultWQKms;
felix@3596:         }
felix@3455: 
sascha@3461:         return wqkms != null ? wqkms[index] : null;
ingo@696:     }
sascha@1061: 
felix@1136: 
felix@1136:     /** Copy deeply. */
felix@1136:     @Override
sascha@1061:     public Facet deepCopy() {
sascha@1061:         WaterlevelFacet copy = new WaterlevelFacet();
sascha@1061:         copy.set(this);
sascha@1061:         copy.type    = type;
felix@2125:         copy.hash    = hash;
felix@2122:         copy.stateId = stateId;
sascha@1061:         return copy;
sascha@1061:     }
ingo@696: }
ingo@696: // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :