csip-hydrotools Displaying differences for changeset

Displaying differences for changeset

display as Allow wrapping

src/java/bin/lin-amd64/HECRAS_Converter_V31.exe

src/java/bin/lin-amd64/NOAA_Converter_V31_non.exe

src/java/bin/lin-amd64/NOAA_Converter_V31_sm.exe

src/java/bin/lin-amd64/NRCC_Converter_V31_non.exe

src/java/bin/lin-amd64/NRCC_Converter_V31_sm.exe

src/java/bin/lin-amd64/TR20.exe

src/java/m/efh2/EFH2HydrologyModel.java

 @@ -1,3 +1,4 @@
++
   package m.efh2;
   import m.efh2.EFH2RainfallCoefficients.RfC;
 @@ -10,218 +11,235 @@
    */
   public class EFH2HydrologyModel {
--    /**
--     * storm type parameter, from values in EFH2RainfallCoefficients
--     */
--    private String stormType = "II";
--    /**
--     * total precipitation parameter, inches
--     */
--    private double precip = 1.0;
--    /**
--     * watershed length parameter, feet
--     */
--    private double watershedLength = 200.0;
--    /**
--     * watershed slope parameter, percent
--     */
--    private double watershedSlope = 0.1;
--    /**
--     * runoff curve number parameter, for lookups
--     */
--    private int runoffCurveNumber = 40;
--    // outputs
--    /**
--     * runoff depth output, in
--     */
--    private double runoffQ = Double.NaN;
--    /**
--     * time of concentration, hours
--     */
--    private double Tc;
--    /**
--     * unit peak discharge output, cfs/ac-in
--     */
--    private double qu;
++  /**
++   * storm type parameter, from values in EFH2RainfallCoefficients
++   */
++  private String stormType = "II";
++  /**
++   * total precipitation parameter, inches
++   */
++  private double precip = 1.0;
++  /**
++   * watershed length parameter, feet
++   */
++  private double watershedLength = 200.0;
++  /**
++   * watershed slope parameter, percent
++   */
++  private double watershedSlope = 0.1;
++  /**
++   * runoff curve number parameter, for lookups
++   */
++  private int runoffCurveNumber = 40;
++  // outputs
++  /**
++   * runoff depth output, in
++   */
++  private double runoffQ = Double.NaN;
++  /**
++   * time of concentration, hours
++   */
++  private double Tc;
++  /**
++   * unit peak discharge output, cfs/ac-in
++   */
++  private double qu;
--    public void simulate() {
--        // lookup Q based on P and CN, using Table 2.2
--        runoffQ = EFH2RunoffTable.lookupQ(precip, runoffCurveNumber);
--        computeTc();
--        computeQu();
++
++  public void simulate() {
++    // lookup Q based on P and CN, using Table 2.2
++    runoffQ = EFH2RunoffTable.lookupQ(precip, runoffCurveNumber);
++    computeTc();
++    computeQu();
++  }
++
++
++  private void computeTc() {
++    // compute Tc using Eq. 2-5
++    Tc = Math.pow(watershedLength, 0.8) * Math.pow((1000.0 / runoffCurveNumber) - 9.0, 0.7)
++        / (1140.0 * Math.sqrt(watershedSlope));
++    if (Tc < 0.1) {
++      Tc = 0.1;
++    } else if (Tc > 10) {
++      throw new IllegalArgumentException(
++          "Computed Tc > 10 hours; EFH2 model is not valid for this range.\n" + "Please estimate peak discharge using NEH-4 procedure, TR-20 program");
++    }
++  }
++
++
++  private void computeQu() {
++    // lookup rainfall coefficients for qu computation
++    double Ia = EFH2IaTable.lookupIa(runoffCurveNumber);
++    double SAPG = Ia / precip;    // Ia/P
++    List<RfC> rfcs = EFH2RainfallCoefficients.getRainfallCoeff(stormType, SAPG);
++    if (rfcs.size() < 1) {
++      //Logger logger = Logger.getLogger(this.getClass());
++      System.err.println("Failed to lookup rainfall coefficients for storm " + stormType + " and Ia/P = " + SAPG);
++      qu = Double.NaN;
++      return;
       }
--    private void computeTc() {
--        // compute Tc using Eq. 2-5
--        Tc = Math.pow(watershedLength, 0.8) * Math.pow((1000.0 / runoffCurveNumber) - 9.0, 0.7)
--                / (1140.0 * Math.sqrt(watershedSlope));
--        if (Tc < 0.1) {
--            Tc = 0.1;
--        } else if (Tc > 10) {
--            throw new IllegalArgumentException(
--                    "Computed Tc > 10 hours; EFH2 model is not valid for this range.\n" + "Please estimate peak discharge using NEH-4 procedure, TR-20 program");
--        }
++    // compute qu, per algorithm in Basic source
++    double logTc = Math.log10(Tc);
++    RfC rfc = rfcs.get(0);
++    double logLower = rfc.getRfC2() + (rfc.getRfC3() + rfc.getRfC4() * logTc) * logTc;
++    double lowerQug = Math.pow(10.0, logLower);
++    if (rfcs.size() == 1) {
++
++      // Ia/P matched one of the standard values, use just one set of coeffs
++      qu = lowerQug / 640.0;
++
++    } else {
++      // Ia/P fell between two standard values, read second RfC coeffs
++      double lowerSAPG = rfc.getSAPG();
++      rfc = rfcs.get(1);
++      double upperSAPG = rfc.getSAPG();
++
++      // compute second Qug result and interpolate
++      double logUpper = rfc.getRfC2() + (rfc.getRfC3() + rfc.getRfC4() * logTc) * logTc;
++      double upperQug = Math.pow(10.0, logUpper);
++      double Qug = lowerQug + (upperQug - lowerQug) * (SAPG - lowerSAPG) / (upperSAPG - lowerSAPG);
++      qu = Qug / 640.0;
       }
++  }
--    private void computeQu() {
--        // lookup rainfall coefficients for qu computation
--        double Ia = EFH2IaTable.lookupIa(runoffCurveNumber);
--        double SAPG = Ia / precip;    // Ia/P
--        List<RfC> rfcs = EFH2RainfallCoefficients.getRainfallCoeff(stormType, SAPG);
--        if (rfcs.size() < 1) {
--            //Logger logger = Logger.getLogger(this.getClass());
--            System.err.println("Failed to lookup rainfall coefficients for storm " + stormType + " and Ia/P = " + SAPG);
--            qu = Double.NaN;
--            return;
--        }
--        // compute qu, per algorithm in Basic source
--        double logTc = Math.log10(Tc);
--        RfC rfc = rfcs.get(0);
--        double logLower = rfc.getRfC2() + (rfc.getRfC3() + rfc.getRfC4() * logTc) * logTc;
--        double lowerQug = Math.pow(10.0, logLower);
--        if (rfcs.size() == 1) {
++  /**
++   * get storm type parameter, see enumerated values
++   */
++  public String getStormType() {
++    return stormType;
++  }
--            // Ia/P matched one of the standard values, use just one set of coeffs
--            qu = lowerQug / 640.0;
--        } else {
--            // Ia/P fell between two standard values, read second RfC coeffs
--            double lowerSAPG = rfc.getSAPG();
--            rfc = rfcs.get(1);
--            double upperSAPG = rfc.getSAPG();
++  /**
++   * set storm type parameter
++   *
++   * @param stormType
++   */
++  public void setStormType(String stormType) {
++    this.stormType = stormType;
++  }
--            // compute second Qug result and interpolate
--            double logUpper = rfc.getRfC2() + (rfc.getRfC3() + rfc.getRfC4() * logTc) * logTc;
--            double upperQug = Math.pow(10.0, logUpper);
--            double Qug = lowerQug + (upperQug - lowerQug) * (SAPG - lowerSAPG) / (upperSAPG - lowerSAPG);
--            qu = Qug / 640.0;
--        }
++
++  public List<String> getStormTypes() {
++    return EFH2RainfallCoefficients.getStormTypes();
++  }
++
++
++  /**
++   * get total precipitation parameter
++   *
++   * @return precip in inches
++   */
++  public double getPrecip() {
++    return precip;
++  }
++
++
++  /**
++   * set total precipitation parameter
++   *
++   * @param precip total storm precip in inches
++   */
++  public void setPrecip(double precip) {
++    if (precip < 1.0 || precip > 15) {
++      throw new IllegalArgumentException("Precip out of range: (1.0 .. 15.0) : " + precip);
       }
++    this.precip = Math.abs(precip);
++  }
--    /**
--     * get storm type parameter, see enumerated values
--     */
--    public String getStormType() {
--        return stormType;
++
++  /**
++   * get watershed length parameter, feet
++   */
++  public double getWatershedLength() {
++    return watershedLength;
++  }
++
++
++  /**
++   * set watershed length parameter
++   *
++   * @param watershedLength feet
++   */
++  public void setWatershedLength(double watershedLength) {
++    if (this.watershedLength < 200 || this.watershedLength > 26000) {
++      throw new IllegalArgumentException(
++          "Watershed flow length " + watershedLength
++          + " out of bounds; please use TR-55 program to compute Tc and qu");
       }
++    this.watershedLength = Math.abs(watershedLength);
++  }
--    /**
--     * set storm type parameter
--     *
--     * @param stormType
--     */
--    public void setStormType(String stormType) {
--        this.stormType = stormType;
++
++  /**
++   * get watershed slope parameter
++   *
++   * @return slope, ft/ft
++   */
++  public double getWatershedSlope() {
++    return watershedSlope;
++  }
++
++
++  /**
++   * set watershed slope parameter
++   *
++   * @param watershedSlope ft/ft
++   */
++  public void setWatershedSlope(double watershedSlope) {
++    if (Math.abs(watershedSlope) > 64.0 || Math.abs(watershedSlope) < 0.1) {
++      throw new IllegalArgumentException(
++          "Watershed slope must be in percent: 0.1% .. 64%");
       }
++    this.watershedSlope = Math.abs(watershedSlope);
++  }
--    public List<String> getStormTypes() {
--        return EFH2RainfallCoefficients.getStormTypes();
++
++  /**
++   * get runoff curve number parameter
++   *
++   * @return curve number
++   */
++  public int getRunoffCurveNumber() {
++    return runoffCurveNumber;
++  }
++
++
++  /**
++   * set runoff curve number parameter
++   *
++   * @param runoffCurveNumber curve number
++   */
++  public void setRunoffCurveNumber(int runoffCurveNumber) {
++    if (runoffCurveNumber < 40 || runoffCurveNumber > 95) {
++      throw new IllegalArgumentException("Runoff Curve Number: 40 .. 95");
       }
++    this.runoffCurveNumber = runoffCurveNumber;
++  }
--    /**
--     * get total precipitation parameter
--     *
--     * @return precip in inches
--     */
--    public double getPrecip() {
--        return precip;
--    }
--    /**
--     * set total precipitation parameter
--     *
--     * @param precip total storm precip in inches
--     */
--    public void setPrecip(double precip) {
--        if (precip < 1.0 || precip > 15) {
--            throw new IllegalArgumentException("Precip out of range: (1.0 .. 15.0) : " + precip);
--        }
--        this.precip = Math.abs(precip);
--    }
++  /**
++   * get computed runoff, in inches
++   */
++  public double getRunoffQ() {
++    return runoffQ;
++  }
--    /**
--     * get watershed length parameter, feet
--     */
--    public double getWatershedLength() {
--        return watershedLength;
--    }
--    /**
--     * set watershed length parameter
--     *
--     * @param watershedLength feet
--     */
--    public void setWatershedLength(double watershedLength) {
--        if (this.watershedLength < 200 || this.watershedLength > 26000) {
--            throw new IllegalArgumentException(
--                    "Watershed flow length " + watershedLength
--                    + " out of bounds; please use TR-55 program to compute Tc and qu");
--        }
--        this.watershedLength = Math.abs(watershedLength);
--    }
++  /**
++   * get computed Tc time of concentration, hours
++   */
++  public double getTimeOfConcentration() {
++    return Tc;
++  }
--    /**
--     * get watershed slope parameter
--     *
--     * @return slope, ft/ft
--     */
--    public double getWatershedSlope() {
--        return watershedSlope;
--    }
--    /**
--     * set watershed slope parameter
--     *
--     * @param watershedSlope ft/ft
--     */
--    public void setWatershedSlope(double watershedSlope) {
--        if (Math.abs(watershedSlope) > 64.0 || Math.abs(watershedSlope) < 0.1) {
--            throw new IllegalArgumentException(
--                    "Watershed slope must be in percent: 0.1% .. 64%");
--        }
--        this.watershedSlope = Math.abs(watershedSlope);
--    }
--
--    /**
--     * get runoff curve number parameter
--     *
--     * @return curve number
--     */
--    public int getRunoffCurveNumber() {
--        return runoffCurveNumber;
--    }
--
--    /**
--     * set runoff curve number parameter
--     *
--     * @param runoffCurveNumber curve number
--     */
--    public void setRunoffCurveNumber(int runoffCurveNumber) {
--        if (runoffCurveNumber < 40 || runoffCurveNumber > 95) {
--            throw new IllegalArgumentException("Runoff Curve Number: 40 .. 95");
--        }
--        this.runoffCurveNumber = runoffCurveNumber;
--    }
--
--    /**
--     * get computed runoff, in inches
--     */
--    public double getRunoffQ() {
--        return runoffQ;
--    }
--
--    /**
--     * get computed Tc time of concentration, hours
--     */
--    public double getTimeOfConcentration() {
--        return Tc;
--    }
--
--    /**
--     * get computed peak discharge rate, cfs/ac-in??
--     */
--    public double getUnitPeakDischarge() {
--        return qu;
--    }
++  /**
++   * get computed peak discharge rate, cfs/ac-in??
++   */
++  public double getUnitPeakDischarge() {
++    return qu;
++  }
   }

src/java/m/hydrotools/efh2/V1_0.java

@@ -19,7 +19,7 @@

static final String RUNOFF = "runoff";

static final String RUNOFFCURVENUMBER = "runoffcurvenumber";

static final String STORMTYPE = "stormtype";

- static final String TIMEOFCONCENTRATION = "timeofconcentration";

+ static final String TIMEOFCONCENTRATION = "Tc";

static final String UNITPEAKDISCHARGE = "unitpeakdischarge";

static final String WATERSHEDLENGTH = "watershedlength";

static final String WATERSHEDSLOPE = "watershedslope";

@@ -45,9 +45,9 @@

@Override

protected void postProcess() throws Exception {

- putResult(RUNOFF, model.getRunoffQ(), "Runoff depth output", "inches");

- putResult(TIMEOFCONCENTRATION, model.getTimeOfConcentration(), "Time of concentration", "hours" );

- putResult(UNITPEAKDISCHARGE, model.getUnitPeakDischarge(), "Unit Peak Discharge", "cfs/ac-in");

+ putResult(RUNOFF, model.getRunoffQ(), "Runoff depth output", "inch");

+ putResult(TIMEOFCONCENTRATION, model.getTimeOfConcentration(), "Time in hours for runoff to flow from the most hydraulically remote point to the project site", "hours" );

+ putResult(UNITPEAKDISCHARGE, model.getUnitPeakDischarge(), "Unit Peak Discharge", "cfs/acre/inch");

}

src/java/m/hydrotools/efh2/V1_0.json

@@ -3,7 +3,7 @@

"parameter": [

{

"name": "precip",

- "description": "precip",

+ "description": "The volume of excess precipitation, in inches on the watershed, not absorbed by the soil nor held on the vegetation.",

"unit": "inch",

"min": "1",

"max": "15",

@@ -11,7 +11,7 @@

{

"name": "runoffcurvenumber",

- "description": "runoffcurvenumber",

+ "description": "Weighted NRCS curve number representing the portion of precipitation that runs off the land.",

"unit": "",

"min": "40",

"max": "95",

@@ -24,15 +24,15 @@

{

"name": "watershedlength",

- "description": "watershedlength",

- "unit": "foot",

+ "description": "Length in feet along the flow path from the hydraulically most distant point to the project site.",

+ "unit": "ft",

"min": "200",

"max": "26000",

"value": 1500

{

"name": "watershedslope",

- "description": "watershedslope",

+ "description": "Average slope in percent of the all the contributing land within the watershed boundary.",

"unit": "%",

"min": "0.1",

"max": "64",

src/java/m/hydrotools/efh2/V2_0.java

 @@ -90,20 +90,20 @@
       st_name = getStringParam(STORM_NAME);
       double dr_area = getDoubleParam(DRAINAGE_AREA);
--    checkRange(dr_area, DRAINAGE_AREA, 0, 2000);
--    checkGT(dr_area, DRAINAGE_AREA, 0);
++//    checkRange(dr_area, DRAINAGE_AREA, 0, 3.125);
++//    checkGT(dr_area, DRAINAGE_AREA, 0);
       double cu_number = getDoubleParam(CURVE_NUMBER);
--    checkRange(cu_number, CURVE_NUMBER, 30, 100);
++//    checkRange(cu_number, CURVE_NUMBER, 30, 100);
       double ws_slope = getDoubleParam(WATERSHED_SLOPE);
--    checkRange(ws_slope, WATERSHED_SLOPE, 0.5, 64);
++//    checkRange(ws_slope, WATERSHED_SLOPE, 0.5, 64);
       double ws_length = getDoubleParam(WATERSHED_LENGTH);
--    checkRange(ws_length, WATERSHED_LENGTH, 200, 26000);
++//    checkRange(ws_length, WATERSHED_LENGTH, 200, 26000);
       double ra_depth = getDoubleParam(RAIN_DEPTH);
--    checkRange(ra_depth, RAIN_DEPTH, 0, 26);
++//    checkRange(ra_depth, RAIN_DEPTH, 0, 26);
       String ra_dist = getStringParam(RAIN_DIST);
 @@ -128,38 +128,44 @@
       int ret = tr20.exec();
       if (ret != 0) {
         File err = tr20.stderr();
--      if (err.exists()) {
--        new ServiceException(FileUtils.readFileToString(err));
++      if (err.exists() && err.length() > 0) {
++        throw new ServiceException(FileUtils.readFileToString(err));
         }
--      new ServiceException("Error executing TR20");
++      throw new ServiceException("Error executing TR20");
       }
++
++    File err = new File(getWorkspaceDir(), "TR20.err");
++    if (err.exists() && err.length() > 0) {
++      throw new ServiceException(FileUtils.readFileToString(err));
++    }
++
       File out = new File(getWorkspaceDir(), "TR20.out");
       if (!out.exists()) {
--      new ServiceException("Missing File : TR20.out");
++      throw new ServiceException("Missing File : TR20.out");
       }
       File dbg = new File(getWorkspaceDir(), "TR20.dbg");
       if (!dbg.exists()) {
--      new ServiceException("Missing File : TR20.dbg");
++      throw new ServiceException("Missing File : TR20.dbg");
       }
       res = parseOutputFile(st_name, out);
       if (res == null || res.length != 3) {
--      new ServiceException("Missing output.");
++      throw new ServiceException("Missing output.");
       }
       toc = parseDebugFile(dbg);
       if (toc == null) {
--      new ServiceException("Missing tocs.");
++      throw new ServiceException("Missing tocs.");
       }
     }
     @Override
     protected void postProcess() throws Exception {
--    putResult(RUNOFF_VOLUME, res[0], "inches");
--    putResult(PEAK_TIME, res[1], "h");
--    putResult(PEAK_DISCHARGE, res[2], "cfs");
--    putResult(TOC, toc, "h");
++    putResult(RUNOFF_VOLUME, res[0], "Runoff depth output", "inch");
++    putResult(PEAK_TIME, res[1], "Peak Time", "h");
++    putResult(PEAK_DISCHARGE, res[2], "Unit Peak Discharge", "cfs/acre/inch");
++    putResult(TOC, toc, "Time of Concentration", "h");
     }

src/java/m/hydrotools/efh2/V2_0.json

 @@ -7,8 +7,9 @@
     }, {
      "name": "drainage_area",
      "description": "Drainage Area",
--   "min": "1",
--   "max": "2000",
++   "unit": "miles^2",
++   "min": "0",
++   "max": "3.125",
      "value": "0.01"
     }, {
      "name": "curve_number",

src/java/m/hydrotools/tr20/V2_0.java

@@ -32,7 +32,6 @@

@Resource(file = "/bin/lin-amd64/NRCC_Converter_V31_non.exe", type = EXECUTABLE, id = "nrcc")

@Resource(file = "/bin/lin-amd64/NRCC_Converter_V31_sm.exe", type = EXECUTABLE, id = "nrcc_sm")

@Resource(file = "*.out *.hyd *.err *.dbg stdout.txt stderr.txt", type = OUTPUT)

public class V2_0 extends ModelDataService {

String conv;

@@ -79,5 +78,4 @@

throw new ServiceException(" TR20 error " + conv);

}