Skip to content
Snippets Groups Projects

Issues running DisaggCalc for certain CEUS sites

    • Closed Issue created by Andrew J. Makdisi @amakdisi

    I cloned and compiled the most recent version of nshmp-haz and nshm-conus this past Friday (Jan 7), and have started testing out the calculation capabilities a bit more today beyond the documented examples. One of the things I was hoping to do was perform disaggregation calcs for multiple sites across the US in the same sites.csv file. So I ran a slightly modified version of the disagg example in etc/examples/7-disaggregation using a site file that includes some WUS (SF, Seattle, LA, Salt Lake), and CEUS (Memphis and Charleston) sites, for five IMs (PGA, SA0P1, SA0P5, SA1P0, SA2P0).

    The WUS worked perfectly, but I ran into some issues with Charleston and Memphis, ultimately crashing with the following stack:

    INFO
   INFO Sites file: C:\Users\amakdisi\Projects\nshm-project\latest-version-2\nshmp-haz\etc\examples\7-disagg-CEUS-test\sites.csv
   INFO Sites: Sites [size=6]
               Los Angeles CA              -118.250 34.050 Vs30=760.0 inferred  Z1.0=NaN Z2.5=NaN
               San Francisco CA            -122.400 37.750 Vs30=760.0 inferred  Z1.0=NaN Z2.5=NaN
               Seattle WA                  -122.300 47.600 Vs30=760.0 inferred  Z1.0=NaN Z2.5=NaN
               Salt Lake City UT           -111.900 40.750 Vs30=760.0 inferred  Z1.0=NaN Z2.5=NaN
               Charleston SC               -79.930 32.780 Vs30=760.0 inferred  Z1.0=NaN Z2.5=NaN
               ... and 1 more ...
   INFO Threads: 16
   INFO DisaggCalc: calculating ...
 SEVERE @ HazardCalc.handleError()
** Exiting **
DisaggCalc: error
 Arguments: [../../../../nshm-conus, sites.csv, config.json]
java.lang.ArrayIndexOutOfBoundsException: Index -1 out of bounds for length 50
        at gov.usgs.earthquake.nshmp.data.IntervalTable.get(IntervalTable.java:90)
        at gov.usgs.earthquake.nshmp.calc.DisaggExport.createSummaryElements(DisaggExport.java:343)
        at gov.usgs.earthquake.nshmp.calc.DisaggExport.<init>(DisaggExport.java:79)
        at gov.usgs.earthquake.nshmp.calc.HazardExport.writeDisagg(HazardExport.java:334)
        at gov.usgs.earthquake.nshmp.calc.HazardExport.write(HazardExport.java:207)
        at gov.usgs.earthquake.nshmp.DisaggCalc.calc(DisaggCalc.java:153)
        at gov.usgs.earthquake.nshmp.DisaggCalc.run(DisaggCalc.java:106)
        at gov.usgs.earthquake.nshmp.DisaggCalc.main(DisaggCalc.java:63)
DisaggCalc [TODO get version from resource]
Usage:
  java -cp nshmp-haz.jar gov.usgs.earthquake.nshmp.DisaggCalc model sites [config]
Where:
  'model' is a model directory
  'sites' is a *.csv file or *.geojson file of sites and data
     - site class and basin terms are optional
  'config' (optional) supplies a calculation configuration
For more information, see:
  https://code.usgs.gov/ghsc/nshmp/nshmp-haz/-/tree/main/docs

    For Charleston, I was able to at least obtain the hazard curves for all IMs, and the Disagg results for PGA, SA0P1, and SA0P5. Memphis didn't work for any of the IMs. I should note that I was able to run these sites through the HazardCalc module with no issues.

    Edited

    Designs

      Child items ...

      2. Activity

      • Andrew J. Makdisi changed the description

        changed the description

      • Andrew J. Makdisi added Bug label

        added Bug label

      • Powers, Peter M.
        Powers, Peter M. @pmpowers ·
        Owner

        @amakdisi There were significant revisions to the CLI applications over the past few weeks. I just ran disaggs with the latest nshmp-haz and nshm-conus models (main branches) for the set of IMTs above and at sites that includes Charleston and Memphis; 2475 return periop. I didn't encounter any errors and the disaggCalc application output both hazard and disaggregation results. Can you pull the latest and rerun?

        Edited by Powers, Peter M.
      • Andrew J. Makdisi
        Andrew J. Makdisi @amakdisi ·
        Author Developer

        I just cloned and compiled the latest version, and attempted to run the example hazard calc in etc/examples/5-complex-model/ (with and without the -Dfile.encoding="UTF-8" specification in the HazardCalc call, but am now getting an error and crash at the model parsing stage.

        Based on the relevant portion of the stacktrace below, and diving into the latest nshm-conus, it looks like there's a subfolder (R1-1300-yr) missing from UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\cluster\scaled\δ+15° (or it's trying to parse something it's not supposed to?). Not totally sure.

        ...
    tree: [2716] UT\West Cache (Junction Hills).geojson
    tree: [2715] UT\West Cache (Wellsville).geojson
    tree: [2799] UT\Wasatch
  branch:        UT\Wasatch\segmented
  branch:        UT\Wasatch\segmented\Brigham City
  branch:        UT\Wasatch\segmented\Brigham City\full-rupture
  branch:        UT\Wasatch\segmented\Brigham City\full-rupture\scaled
  branch:        UT\Wasatch\segmented\Brigham City\full-rupture\scaled\δ+15°
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Brigham City\full-rupture\scaled\--δ--
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Brigham City\full-rupture\scaled\δ-15°
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Brigham City\partial-rupture
  branch:        UT\Wasatch\segmented\Brigham City\partial-rupture\scaled
  branch:        UT\Wasatch\segmented\Brigham City\partial-rupture\scaled\δ+15°
  branch:        UT\Wasatch\segmented\Brigham City\partial-rupture\scaled\--δ--
  branch:        UT\Wasatch\segmented\Brigham City\partial-rupture\scaled\δ-15°
  branch:        UT\Wasatch\segmented\Weber
  branch:        UT\Wasatch\segmented\Weber\full-rupture
  branch:        UT\Wasatch\segmented\Weber\full-rupture\scaled
  branch:        UT\Wasatch\segmented\Weber\full-rupture\scaled\δ+15°
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Weber\full-rupture\scaled\--δ--
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Weber\full-rupture\scaled\δ-15°
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Weber\partial-rupture
  branch:        UT\Wasatch\segmented\Weber\partial-rupture\scaled
  branch:        UT\Wasatch\segmented\Weber\partial-rupture\scaled\δ+15°
  branch:        UT\Wasatch\segmented\Weber\partial-rupture\scaled\--δ--
  branch:        UT\Wasatch\segmented\Weber\partial-rupture\scaled\δ-15°
  branch:        UT\Wasatch\segmented\Salt Lake City
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geodetic
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geodetic\δ+15°
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geodetic\--δ--
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geodetic\δ-15°
                 WARNING: Multiple SINGLE MFDs with epistemic uncertainty
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\geologic
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\geologic\δ+15°
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\geologic\--δ--
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\geologic\δ-15°
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\connected
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\connected\scaled
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\connected\scaled\δ+15°
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\connected\scaled\--δ--
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\connected\scaled\δ-15°
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\cluster
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\cluster\scaled
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\cluster\scaled\δ+15°
  branch:        UT\Wasatch\segmented\Salt Lake City\full-rupture\scaled\geologic\cluster\scaled\δ+15°\R1-1300-yr
 SEVERE @ ModelLoader.load()Error loading model: C:\Users\amakdisi\Projects\nshm-project\latest-version\nshm-conus

java.nio.file.InvalidPathException: Illegal char <:> at index 34: Wasatch (Salt Lake City, cluster) : R1-1300-yr
        at java.base/sun.nio.fs.WindowsPathParser.normalize(WindowsPathParser.java:182)
        at java.base/sun.nio.fs.WindowsPathParser.parse(WindowsPathParser.java:153)
        at java.base/sun.nio.fs.WindowsPathParser.parse(WindowsPathParser.java:77)
        at java.base/sun.nio.fs.WindowsPath.parse(WindowsPath.java:92)
        at java.base/sun.nio.fs.WindowsFileSystem.getPath(WindowsFileSystem.java:229)
        at java.base/java.nio.file.Path.of(Path.java:147)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processClusterBranch(ModelLoader.java:639)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processBranch(ModelLoader.java:590)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processBranch(ModelLoader.java:579)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processBranch(ModelLoader.java:579)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processBranch(ModelLoader.java:579)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processBranch(ModelLoader.java:579)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processBranch(ModelLoader.java:579)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processBranch(ModelLoader.java:579)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.processBranch(ModelLoader.java:579)
        at gov.usgs.earthquake.nshmp.model.ModelLoader$Fault.loadSourceTree(ModelLoader.java:542)
        at gov.usgs.earthquake.nshmp.model.ModelLoader.loadSourceDirectory(ModelLoader.java:298)
        at gov.usgs.earthquake.nshmp.model.ModelFiles.lambda$readNestedDirectories$2(ModelFiles.java:125)
        at java.base/java.util.stream.ReferencePipeline$3$1.accept(ReferencePipeline.java:195)
        at java.base/java.util.stream.ReferencePipeline$2$1.accept(ReferencePipeline.java:177)
        at java.base/java.util.Iterator.forEachRemaining(Iterator.java:133)
        at java.base/java.util.Spliterators$IteratorSpliterator.forEachRemaining(Spliterators.java:1801)
        at java.base/java.util.stream.AbstractPipeline.copyInto(AbstractPipeline.java:484)
        at java.base/java.util.stream.AbstractPipeline.wrapAndCopyInto(AbstractPipeline.java:474)
        at java.base/java.util.stream.ForEachOps$ForEachOp.evaluateSequential(ForEachOps.java:150)
        at java.base/java.util.stream.ForEachOps$ForEachOp$OfRef.evaluateSequential(ForEachOps.java:173)
        at java.base/java.util.stream.AbstractPipeline.evaluate(AbstractPipeline.java:234)
        at java.base/java.util.stream.ReferencePipeline.forEach(ReferencePipeline.java:497)
        at gov.usgs.earthquake.nshmp.model.ModelFiles.readNestedDirectories(ModelFiles.java:126)
        at gov.usgs.earthquake.nshmp.model.ModelLoader.loadSourceDirectory(ModelLoader.java:307)
        at gov.usgs.earthquake.nshmp.model.ModelFiles.lambda$readNestedDirectories$2(ModelFiles.java:125)
        at java.base/java.util.stream.ReferencePipeline$3$1.accept(ReferencePipeline.java:195)
        at java.base/java.util.stream.ReferencePipeline$2$1.accept(ReferencePipeline.java:177)
        at java.base/java.util.Iterator.forEachRemaining(Iterator.java:133)
        at java.base/java.util.Spliterators$IteratorSpliterator.forEachRemaining(Spliterators.java:1801)
        at java.base/java.util.stream.AbstractPipeline.copyInto(AbstractPipeline.java:484)
        at java.base/java.util.stream.AbstractPipeline.wrapAndCopyInto(AbstractPipeline.java:474)
        at java.base/java.util.stream.ForEachOps$ForEachOp.evaluateSequential(ForEachOps.java:150)
        at java.base/java.util.stream.ForEachOps$ForEachOp$OfRef.evaluateSequential(ForEachOps.java:173)
        at java.base/java.util.stream.AbstractPipeline.evaluate(AbstractPipeline.java:234)
        at java.base/java.util.stream.ReferencePipeline.forEach(ReferencePipeline.java:497)
        at gov.usgs.earthquake.nshmp.model.ModelFiles.readNestedDirectories(ModelFiles.java:126)
        at gov.usgs.earthquake.nshmp.model.ModelLoader.loadSourceDirectory(ModelLoader.java:307)
        at gov.usgs.earthquake.nshmp.model.ModelLoader.faultSources(ModelLoader.java:328)
        at gov.usgs.earthquake.nshmp.model.ModelLoader.crustalSources(ModelLoader.java:190)
        at gov.usgs.earthquake.nshmp.model.ModelLoader.loadTectonicSetting(ModelLoader.java:166)
        at gov.usgs.earthquake.nshmp.model.ModelLoader.loadTectonicSettings(ModelLoader.java:150)
        at gov.usgs.earthquake.nshmp.model.ModelLoader.load(ModelLoader.java:128)
        at gov.usgs.earthquake.nshmp.model.HazardModel.load(HazardModel.java:103)
        at gov.usgs.earthquake.nshmp.HazardCalc.run(HazardCalc.java:102)
        at gov.usgs.earthquake.nshmp.HazardCalc.main(HazardCalc.java:75)
      • Powers, Peter M.
        Powers, Peter M. @pmpowers ·
        Owner

        See nshmp-lib#141 (closed)

      • Powers, Peter M. marked this issue as related to nshmp-lib#141 (closed)

        marked this issue as related to nshmp-lib#141 (closed)

      • Powers, Peter M. changed milestone to %NSHM Development & Support

        changed milestone to %NSHM Development & Support

      • Clayton, Brandon Scott added To Do label

        added To Do label

      • Powers, Peter M.
        Powers, Peter M. @pmpowers ·
        Owner

        Resolved in !631 (merged)
        Please use latest nshmp-haz and nshm-conus
        Latest nshm-conus is 2023 development model. Use tag 5.1.2 for 2018.

      • Powers, Peter M. closed

        closed

      • Powers, Peter M. removed To Do label

        removed To Do label

      Please register or sign in to reply