[Shake-dev] Error when using Finite Fault file

[Bruce Worden]

Hi Lijam, Bob,

When did you last update your code? For the last week or so the repository has held some new code that I haven't announced yet, but it changes the way some of this is handled. 

The basic issue is that while BA08 uses Rjb, WGRW11 wants Rrup. I can try to address these errors so the system is a little more tolerant of mismatched fault types, but the best solution is to make all of your fault models into polygons as described in the software guide. That will insure that they will work with any GMPE, IPE, or GMICE (and any combination thereof). Line faults can be extended to depth pretty easily, and if a polygonal fault model is available, depths should be as well.

Again, I strongly recommend making ALL faults closed polygons with depths for each point. See "A Note About Finite Faults" in the software guide.

Bruce

P.S. As an example, suppose you have a line fault defined as:

42.0  -121.0
43.0  -122.0

You can make this a polygon by doing something like this:

42.0  -121.0  0.0
43.0  -122.0  0.0
43.0  -122.0  20.0
42.0  -121.0  20.0
42.0  -121.0  0.0

Assuming the fault ruptured to the surface. For a vertical fault, the depth doesn't matter all that much, but if you don't have a depth, you can always use Wells & Coppersmith to estimate it.





On Mar 27, 2012, at 10:01 AM, Hagos, Lijam wrote:

> Hello Bruce,
>  
> I am getting an error message when trying to  re-run ShakeMap (for Sierra El Mayor Earthquake) with a finite fault model as an input and using WGRW11 as a GMICE. It works fine when I use Wald99. The finite fault model for this earthquake is represented by two segments – a line, and polygon.  I rerun this earthquake before with an earlier release of SM3.5 and it worked fine, except for the warning message indicating that it is defaulting to rjb.  Any idea what is happening?
>  
>  
> grind: Using fixed boundary -116.5/31.5/-114/33.3
> grind: 33.3:31.5:-114:-116.5
> grind:
>  
> grind: Magnitude (7.2) exceeds cutoff (7.0); not flagging outliers.
> Warning: dist_rrup expects a polygonal fault; will use dist_rjb, but that may not be what you want! at /space/shake35/bin/../perl/lib/Shake/Distance.pm line 433.
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: Error reading 1st record of file /space/shake35/bin/../data/14607652t/output/KCjsN1Gho6.txt
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> Warning: dist_rrup expects a polygonal fault; will use dist_rjb, but that may not be what you want! at /space/shake35/bin/../perl/lib/Shake/Distance.pm line 433.
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: Error reading 1st record of file /space/shake35/bin/../data/14607652t/output/DR4I7bvSpz.txt
> Warning: dist_rrup expects a polygonal fault; will use dist_rjb, but that may not be what you want! at /space/shake35/bin/../perl/lib/Shake/Distance.pm line 433.
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: Error reading 1st record of file /space/shake35/bin/../data/14607652t/output/zzkufjdaeH.txt
> Warning: dist_rrup expects a polygonal fault; will use dist_rjb, but that may not be what you want! at /space/shake35/bin/../perl/lib/Shake/Distance.pm line 433.
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: Error reading 1st record of file /space/shake35/bin/../data/14607652t/output/IgUJxyXgx8.txt
> Warning: dist_rrup expects a polygonal fault; will use dist_rjb, but that may not be what you want! at /space/shake35/bin/../perl/lib/Shake/Distance.pm line 433.
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: ERROR: Input file has different number of rows (1) than expected (520)
> gmtmath: Error reading 1st record of file /space/shake35/bin/../data/14607652t/output/e9dYiWqDEF.txt
> Warning: dist_rrup expects a polygonal fault; will use dist_rjb, but that may not be what you want! at /space/shake35/bin/../perl/lib/Shake/Distance.pm line 526.
> write_pfile: fault file not formatted correctly at /space/shake35/bin/../perl/lib/Shake/Fault.pm line 207.
> shake_man: Error in grind: 2560
>  
>  
>  
>  
> Best regards,
>  
> Lijam
>  
> -------------
> Lijam Z. Hagos, PhD
> Associate Seismologist,
> Strong Motion Instrumentation program
> California Geological Survey                                                                                
> 801 K st. M 13 -35
> Sacramento, CA 95814                                                                                           
> Phone (w):   (916) 327-1059
> Cell: 408-420-3880
> email: Lijam.Hagos at conservation.ca.gov
>  
>  

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geohazards.usgs.gov/pipermail/shake-dev/attachments/20120328/1e2d7be1/attachment-0001.html>