error with --onlynameddiscon and source at depth, thanks Shubh

crotwell · Aug 30, 2024 · f85a549 · f85a549
1 parent 6ac8663
commit f85a549
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/src/main/java/edu/sc/seis/TauP/TauModel.java b/src/main/java/edu/sc/seis/TauP/TauModel.java
@@ -191,6 +191,7 @@ public boolean isFluidBranch(int branchNum) {
 
     /**
      * Is the given depth a "noDisconDepth"?
+     * Usually because model was split at the source or receiver depth, or an added depth used by taup pierce.
      */
     public boolean isNoDisconDepth(double noDisconDepth) {
         for(int i = 0; i < noDisconDepths.length; i++) {

diff --git a/src/main/java/edu/sc/seis/TauP/cmdline/TauP_Find.java b/src/main/java/edu/sc/seis/TauP/cmdline/TauP_Find.java
@@ -291,7 +291,9 @@ public List<Double> getExcludedDepths(TauModel tMod) {
         List<Double> exList = new ArrayList<>(getExcludeDepth(tMod.getVelocityModel()));
         double[] branchDepths = tMod.getBranchDepths();
         for (double branchDepth : branchDepths) {
-            if ((onlyNamedDiscon && branchDepth != 0 &&  !tMod.getVelocityModel().isNamedDisconDepth(branchDepth))) {
+            if ((onlyNamedDiscon && branchDepth != 0
+                    && !tMod.isNoDisconDepth(branchDepth)
+                    &&  !tMod.getVelocityModel().isNamedDisconDepth(branchDepth))) {
                 exList.add(branchDepth);
             }
         }