HighF v14.12 Data Comparison

From SCECpedia
Revision as of 16:42, 17 May 2016 by Maechlin (talk | contribs)
Jump to navigationJump to search

In this page we show a comparison between seismograms generated by Hercules and AWP for some of the stations in the La Habra station list.

Contents

Simulation Results - 17 May 2016

AWP simulation results compared to Center for Engineering Strong Motion Data (CESMD) processed V2 data for La Habra Earthquake for large simulation region. Matching synthetics and observations for 217 sites.

Simulation Results - 01 May 2016

Data from the small region against observations.

Simulation Results - 21 March 2016

Same as 3 March but using corrected AWP data. Updated results using the small domain region. There are two sets of AWP results, one using RWG's rules for clamping Vp/Vs and another one using Bill's rules, which can be described by conditions where Vp/Vs is limited between 1.45 and 3. VpMin is clamped at 1500 m/s and VsMin is clamped at 500 m/s. RWG's rules, on the other hand can be described with conditions that limits Vp/VsMin to 1.45 and VsMin to 500 m/s. Hercules data was calculated with the following simulation parameters: Max Frequency: 4 Hz, Minimum Vs: 500 m/s, Minimum Vp:1500 m/s, Q = 100 Vs, Domain: 28*28*14 km. RWG data was calculated with the small model region. These use Vsmin = 500 m/s, Vpmin = 1500 m/s.

First set of data includes RWG, Hercules, and AWP-RWG:

Simulation Results - 3 March 2016

Don't use these plots, use data from 21 March above instead.

Updated results using the small domain region. There are two sets of AWP results, one using RWG's rules for clamping Vp/Vs and another one using Bill's rules, which can be described by conditions where Vp/Vs is limited between 1.45 and 3. VpMin is clamped at 1500 m/s and VsMin is clamped at 500 m/s. RWG's rules, on the other hand can be described with conditions that limits Vp/VsMin to 1.45 and VsMin to 500 m/s. Hercules data was calculated with the following simulation parameters: Max Frequency: 4 Hz, Minimum Vs: 500 m/s, Minimum Vp:1500 m/s, Q = 100 Vs, Domain: 28*28*14 km. RWG data was calculated with the small model region. These use Vsmin = 500 m/s, Vpmin = 1500 m/s.

First set of data includes RWG, Hercules, and AWP-RWG:

Second set of plots also includes AWP-WHS:

Simulation Results - February 22 - Using small domain

The following results correspond to runs done with Hercules comparing previous results obtained with the large simulation domain against new results obtained with the smaller simulation domain. The former are shown in red, whereas the latter are shown for two conditions in blue and green. The results shown in blue are equivalent to those in red with the only difference that the red results were obtained with the original large 180 km x 135 km domain, but the blue were obtained with the smaller 28 km x 28 km domain. In both these cases the minimum Vs is equal to 500 m/s but there is no cap on the minimum Vp. Then the results shown in green correspond to a run done with the smaller domain but now applying a constrain to a minimum Vp equal to 1500 m/s.

Proposed Small Simulation Region

Simulation Results - 5 January 2016

Same data as 14 December 2015, but updated RWG data for 3 closest stations computed using a 3D velocity model generated directly from "ucvm_query".

Simulation Results - 15 December 2015

Same data as 14 December 2015, but adding observation data.

Simulation Results - 14 December 2015

Same data from 7 December 2015, updated RWG data for the three closest stations only.

Simulation Results - 7 December 2015

Same data from 20 November 2015. Here are the steps used to process the simulated time series:

  • Read all seismograms for a given station, find out largest dt, it becomes new_dt
  • Filter at 10Hz back and forward, 4 poles, using filtfilt
  • Interpolate using 1D linear to new_dt (in this case, AWP's 0.025 dt)
  • Trim longest seismograms to match the shortest one
  • Taper at the last 2 seconds: hanning with 161 points, then take 81 last values and apply this to last 81 points of each time series
  • Pad last 5 seconds with zeroes
  • Filter at 4Hz 4 poles back and forward with filtfilt
  • Plot time series and FFT using 1.25 filter frequency
  • Repeat last 2 steps with filter frequency set to 0.5, 1, 2, and 3 Hz

Simulation Results - 20 November 2015

AWP and Hercules data from 5 November, but updated RWG results to correct issue in the mesh conversion code. Added results for far-distance stations as RWG data now goes to 100s.

Simulation Results - 05 November 2015

Results from RWG, updated SDSU results by increasing the number of points in the sponge zones for the Cerjan boundary conditions from 20 points (previous runs) to 90 points.

Simulation Results - 14 October 2015

La Habra earthquake with Hercules at 14 ppw, 4Hz, Min Vs = 500m/s, CVM-S4.26, 3D, with Q

Simulation Results - 2 October 2015

La Habra earthquake with constant Q based on CVM-S4.26 w Min Vs 500m/s. Hercules results down with 10ppw (same as 28 September) but reference frequency is 0.0 Hz (it helps to decrease the phase dispersion).

Simulation Results - 28 September 2015

La Habra earthquake with constant Q based on CVM-S4.26 w Min Vs 500m/s. Hercules results down with 10ppw

Simulation Results - 31 August 2015

La Habra-scale earthquake w/ constant Q in 1D Velocity. Adding Hercules results to previous AWP-ODC and RWG results from August 13. Filtered plots were generated using SAC's "lowpass" command, with 4 poles and 4Hz as the corner frequency.

Simulation Results - 18 August 2015

Updated SDSU results, comparing La Habra synthetics for AWP, with Q(f) in the 1D model. Here, Q(f)=Q0 f^0.8, where Q0 are the values that we have been using for the constant-Q case. Filtered plots were generated using SAC's "lowpass" command, with 4 poles and 4Hz as the corner frequency.

Simulation Results - 13 August 2015

La Habra-scale earthquake w/ constant Q in 1D Velocity. Updated results for AWP-ODC and RWG. Filtered plots were generated using SAC's "lowpass" command, with 4 poles and 4Hz as the corner frequency.

Simulation Results - 3 August 2015

La Habra-scale earthquake w/ constant Q in 1D Velocity

Simulation Results - 20 July 2015 - Half Space

Updated results from all groups. AWP-ODC results shifted right by 0.025s, as instructed by KO. Hercules team provided two sets of results, with 10 and 20 points per wavelength. These are plotted separately as two sets:

Simulation Results - 20 July 2015

Updated results:

  • AWP/ODC: With corrected station positions and a fix for a small time delay. Please note that the orientations are relative to the rotated domain; i.e. '000' means N40dE, and 090 N130dE
  • RWG: Short, mid and far sites computed using emod3d. This is the same as what I sent earlier (for eight of the 9 sites), except that this package includes the response for station CI_Q0049.
  • Hercules: New set of results obtained with Hercules for La Habra with the 1D profile. This set was obtained in a reduced domain and only for the near and mid-range stations, using 20 points per wavelength. I did some plots myself and in my opinion this is the best we are ever going to do, and I think it is pretty good. I still want to see the comparisons done by you, independently. I recommend you do the following, plot this against AWP and Graves results, without filtering. And then plot the same, filtering at 4 Hz. No need to do other lower filters as you did before, 4 Hz will be enough.

Simulation Results - 29 June 2015

Results from AWP, Hercules, RWG, and Hercules-Refined. The Hercules refined results were done with 10 p/w without filtering the source, lossless, refined mesh for upper ~60 meters.

Simulation Results - 25 June 2015

Here's an update on the preliminary modeling I've done for La Habra. I took the 1D velocity model and ran an FK calculation for the 3 nearby sites. I also convolved the FK response with the same slip-rate function used for the FD calculation. The attached plots compare the FD and FK waveforms for these 3 sites. I did not apply any additional filtering to the waveforms, so they should contain frequencies up to 5 Hz (based on my understanding of the freq. content of the slip-rate function). The agreement between FD and FK is very good. This is the level of agreement that I think we can obtain as long as the model inputs are defined consistently. (Note: the late arriving blips in the FD are boundary reflections.)

On this later point, I used a velocity structure for the FK calculation that is consistent with the model representation used for my FD implementation. This means that the layer boundaries are at depths of z=0.01 km + J*0.02 km for some integer J. In other words, the shallowest layer has its bottom at 10 meters depth, and then all subsequent layers are at depths of 30, 50, 70, ... meters. I've attached a couple of plots (different resolutions) showing Vp (heavy line), Vs (thin line) and Den (dashed line) for the upper portion of the model I'm using. While this representation of the 1D structure is valid, I'm not sure if this structure is strictly the same as that "seen" by AWP and Hercules due to the way they discretize the medium.

Results without attenuation for AWP, Hercules, and RWG.

Simulation Results - 16 June 2015

Simulation Results - May 18th 2015

Results without attenuation for both AWP and Hercules. Filtered plots were generated using SAC's "lowpass" command, with 4 poles and the corresponding corner frequency.

Data from the two methods are available here:

Simulation Results - May 13th 2015

AWP results are the same as in May 11th, Hercules results were updated with Reference Frequency = 0.5Hz, source without filtering.

Data from the two methods are available here:

Simulation Results - May 11th 2015

Hercules data uses Reference Frequency of 1Hz.

Background Information

Station Selection

We picked the following stations for this comparison:

  • Short distance:
    • FUL - 5.6 km
    • BRE - 13.5 km
    • DLA - 17.4 km
  • Mid distance
    • 24945 - 40.2 km
    • GR2 - 40.5 km
    • BHP - 40.5 km
  • Far distance
    • Q0049 - 80.1 km
    • AGO - 81.1 km
    • PTD - 81.4 km

You can see the above stations on a map by downloading the kml file and opening with Google Earth.

Hercules Data Files

The Hercules results are located on Blue Waters at ‘/u/sciteam/khoshnev/scratch/N_Lahabra_sim_4Hz/ ‘. The stations are in '../outputfiles/stations/' folder. X and Y components are according to the simulation domain (rotation angle 39.9). For Z component simply multiply it by -1 to get the same positive direction with data. There are 356 stations, which are same as the SCEC wiki page but are ordered according to the epicentral distance (smaller to larger, please see the station list below).

Hercules File format

Each file has 10 columns as follows:

  • C1 : Time steps
  • C2,C3,C4 : Displacement ( Directions : | , --- , . ) (m)
  • C5,C6,C7 : Velocity ( Directions : | , --- , . )(m/s)
  • C8,C9,C10 : Acceleration ( Directions : | , --- , . )(m/s2)

Hercules Station List

The stations below are sorted by epicentral distance:

station.0       33.909  -117.931
station.1       33.931  -117.956
station.2       33.9454 -117.9237
station.3       33.933  -117.896
station.4       33.889  -117.887
station.5       33.8717 -117.9225
station.6       33.947  -117.991
station.7       33.866  -117.959
station.8       33.946  -118.001
station.9       33.905  -118.011
station.10      33.989  -117.89
station.11      34      -117.935
station.12      33.84   -117.957
station.13      34.0095 -117.9508
station.14      33.853  -117.858
station.15      33.889  -118.047
station.16      34.001  -118.0133
station.17      33.927  -117.803
station.18      33.928  -118.058
station.19      34.021  -117.866
station.20      34.036  -117.951
station.21      33.853  -117.818
station.22      33.962  -118.066
station.23      33.8078 -117.9812
station.24      34.043  -117.912
station.25      33.819  -117.85
station.26      34.002  -117.814
station.27      33.899  -117.777
station.28      34.032  -118.019
station.29      33.789  -117.921
station.30      33.791  -117.968
station.31      33.785  -117.89
station.32      33.939  -118.099
station.33      33.837  -118.072
station.34      34.001  -118.078
station.35      34.065  -117.893
station.36      33.8285 -117.7894
station.37      34.073  -117.938
station.38      33.7882 -117.844
station.39      34.067  -117.985
station.40      33.767  -117.918
station.41      33.8482 -118.0962
station.42      33.883  -118.118
station.43      33.92   -118.127
station.44      33.846  -118.107
station.45      33.777  -117.834
station.46      33.757  -117.956
station.47      33.863  -117.742
station.48      33.768  -117.85
station.49      33.786  -117.809
station.50      33.991  -118.121
station.51      34.096  -117.906
station.52      33.749  -117.875
station.53      34.075  -118.042
station.54      34.0508 -118.0808
station.55      33.818  -118.108
station.56      33.798  -118.0884
station.57      33.798  -118.0884
station.58      33.742  -117.894
station.59      33.772  -117.803
station.60      34.1047 -117.9796
station.61      33.969  -118.149
station.62      33.748  -117.844
station.63      33.99   -117.718
station.64      33.9753 -118.1492
station.65      34.091  -117.809
station.66      34.092  -117.81
station.67      33.924  -118.168
station.68      34.0918 -117.8071
station.69      34.056  -117.749
station.70      34.0497 -118.12
station.71      34.114  -117.851
station.72      34.017  -118.147
station.73      33.727  -118.001
station.74      34.041  -118.13
station.75      33.8811 -118.1757
station.76      33.873  -118.175
station.77      33.757  -118.085
station.78      34.132  -117.911
station.79      33.834  -118.161
station.80      33.81   -118.15
station.81      33.723  -118.026
station.82      34.109  -118.058
station.83      34.084  -118.101
station.84      33.818  -118.158
station.85      33.999  -118.178
station.86      33.7043 -117.8858
station.87      33.9988 -117.6804
station.88      34.141  -117.967
station.89      33.777  -118.133
station.90      33.723  -118.054
station.91      33.926  -117.658
station.92      33.7    -117.88
station.93      33.696  -117.953
station.94      34.103  -118.1
station.95      34.096  -117.745
station.96      34.147  -117.999
station.97      33.898  -118.209
station.98      33.759  -118.131
station.99      33.84   -118.195
station.100     34.057  -117.698
station.101     34.037  -118.178
station.102     33.931  -118.214
station.103     34.071  -118.151
station.104     33.976  -118.209
station.105     33.6845 -117.943
station.106     34.125  -117.777
station.107     34.053  -118.172
station.108     34.05   -118.176
station.109     34.097  -118.131
station.110     34.149  -118.03
station.111     33.785  -118.173
station.112     34.066  -118.169
station.113     33.727  -117.748
station.114     33.676  -117.963
station.115     33.876  -117.635
station.116     33.763  -118.168
station.117     34.121  -117.738
station.118     34.031  -118.208
station.119     33.924  -118.239
station.120     33.827  -118.217
station.121     33.7903 -118.1988
station.122     33.663  -117.902
station.123     33.876  -118.24
station.124     34.108  -118.154
station.125     34.076  -117.671
station.126     34.062  -118.201
station.127     34.062  -118.201
station.128     33.802  -118.215
station.129     34.057  -118.206
station.130     33.662  -117.997
station.131     34.027  -118.225
station.132     34.1368 -118.122
station.133     34.058  -118.209
station.134     34.1504 -118.1016
station.135     34.058  -118.21
station.136     34.1369 -118.126
station.137     34.1369 -118.126
station.138     33.769  -118.197
station.139     34.026  -118.23
station.140     34.1368 -118.1283
station.141     34.137  -118.128
station.142     33.9487 -118.2554
station.143     34.1401 -118.1275
station.144     33.929  -118.261
station.145     34.112  -117.681
station.146     33.867  -118.262
station.147     33.876  -117.593
station.148     34.075  -118.218
station.149     33.678  -117.754
station.150     33.678  -117.753
station.151     34.013  -118.257
station.152     34.055  -118.236
station.153     34.155  -118.133
station.154     34.032  -117.608
station.155     33.634  -117.903
station.156     34.041  -118.248
station.157     34.052  -118.242
station.158     34.033  -117.607
station.159     33.664  -117.7686
station.160     34.149  -118.151
station.161     33.66   -117.775
station.162     33.7953 -117.6091
station.163     34.058  -118.246
station.164     33.955  -117.576
station.165     34.059  -118.246
station.166     34.059  -118.247
station.167     33.624  -117.93
station.168     33.624  -117.929
station.169     34.056  -118.252
station.170     34.056  -118.252
station.171     34.058  -118.252
station.172     34.037  -118.264
station.173     34.051  -118.26
station.174     34.048  -118.262
station.175     33.969  -118.291
station.176     34.166  -118.15
station.177     33.904  -117.561
station.178     34.188  -118.115
station.179     34.004  -118.287
station.180     33.939  -118.3
station.181     34.043  -118.272
station.182     33.811  -118.277
station.183     33.8757 -117.5611
station.184     34.0192 -118.2863
station.185     34.025  -118.285
station.186     34.1207 -117.6381
station.187     33.843  -117.565
station.188     34.021  -118.29
station.189     34.056  -118.274
station.190     34.0046 -117.5616
station.191     33.605  -117.89
station.192     34.14   -118.211
station.193     33.909  -118.315
station.194     34.2236 -118.0583
station.195     34.224  -118.058
station.196     34.145  -117.651
station.197     34.08   -118.268
station.198     33.99   -118.312
station.199     33.747  -118.26
station.200     33.69   -117.655
station.201     34.1714 -118.1852
station.202     33.5974 -117.8791
station.203     33.922  -118.327
station.204     33.856  -117.541
station.205     33.938  -118.327
station.206     33.856  -117.54
station.207     33.732  -118.258
station.208     33.749  -118.272
station.209     34.077  -118.287
station.210     33.75   -118.278
station.211     33.8689 -118.3314
station.212     34.021  -118.325
station.213     34.1    -118.288
station.214     34.104  -117.574
station.215     34.134  -117.598
station.216     33.884  -117.514
station.217     34.063  -118.316
station.218     33.993  -118.343
station.219     34.084  -118.308
station.220     33.976  -118.353
station.221     34.014  -117.511
station.222     34.163  -118.255
station.223     34.1182 -118.3002
station.224     33.9905 -118.3617
station.225     33.919  -117.488
station.226     33.869  -118.368
station.227     33.7126 -118.2938
station.228     33.605  -117.705
station.229     34.008  -118.362
station.230     33.886  -118.373
station.231     33.975  -117.487
station.232     33.961  -118.379
station.233     33.945  -118.381
station.234     34.069  -118.346
station.235     33.963  -118.38
station.236     34.0003 -118.3779
station.237     34.208  -118.233
station.238     34.123  -117.54
station.239     33.7403 -118.3352
station.240     33.921  -118.391
station.241     34.037  -118.376
station.242     33.77   -117.501
station.243     34.2388 -117.6585
station.244     34.036  -118.378
station.245     34.051  -118.384
station.246     33.849  -117.455
station.247     33.951  -117.445
station.248     33.917  -117.443
station.249     33.9405 -118.4207
station.250     33.643  -117.569
station.251     34.0762 -118.3959
station.252     33.528  -117.768
station.253     34.072  -118.399
station.254     34.3334 -118.0259
station.255     34.115  -118.38
station.256     33.9627 -118.437
station.257     33.5526 -117.6617
station.258     33.7435 -118.4041
station.259     34.1764 -118.3597
station.260     34.0216 -118.4463
station.261     34.125  -117.454
station.262     34.0996 -117.4388
station.263     33.951  -117.388
station.264     34.2001 -118.3604
station.265     33.731  -117.428
station.266     33.979  -117.374
station.267     34.069  -117.398
station.268     33.9935 -117.3755
station.269     34.2535 -118.3338
station.270     34.1062 -118.455
station.271     34.1062 -118.455
station.272     34.172  -117.428
station.273     33.979  -117.344
station.274     34.193  -117.427
station.275     34.323  -117.565
station.276     33.9731 -117.3272
station.277     33.9733 -117.3267
station.278     34.382  -117.6782
station.279     34.3148 -117.5464
station.280     33.895  -117.322
station.281     34.103  -117.349
station.282     34.033  -117.312
station.283     34.277  -118.397
station.284     34.188  -118.477
station.285     34.4489 -117.9554
station.286     34.0928 -117.3169
station.287     34.2712 -117.4245
station.288     34.206  -117.367
station.289     34.135  -117.322
station.290     34.066  -117.293
station.291     34.163  -117.333
station.292     34.066  -117.28
station.293     34.249  -118.479
station.294     34.096  -117.287
station.295     33.968  -117.252
station.296     34.1054 -117.2829
station.297     34.439  -117.647
station.298     34.0838 -118.599
station.299     34.282  -118.4792
station.300     34.282  -118.4792
station.301     34.4837 -118.1178
station.302     34.165  -117.289
station.303     34.466  -118.204
station.304     33.699  -117.266
station.305     34.426  -117.544
station.306     34.133  -117.253
station.307     34.3053 -118.4881
station.308     34.197  -117.285
station.309     34.5215 -117.991
station.310     34.2287 -118.5583
station.311     33.8616 -117.2053
station.312     33.912  -117.196
station.313     34.3108 -118.4968
station.314     34.0593 -118.6461
station.315     34.222  -118.571
station.316     34.066  -117.215
station.317     34.166  -117.251
station.318     34.1213 -117.219
station.319     34.54   -118.06
station.320     33.921  -117.173
station.321     33.831  -117.181
station.322     34.145  -117.219
station.323     33.3706 -117.5574
station.324     34.5298 -118.1445
station.325     34.307  -117.313
station.326     34.3656 -117.3668
station.327     34.1717 -118.6497
station.328     34.026  -117.155
station.329     34.2324 -117.2348
station.330     34.3956 -118.4984
station.331     34.0355 -117.1124
station.332     33.748  -117.127
station.333     34.1065 -117.0982
station.334     34.197  -117.126
station.335     34.1465 -118.767
station.336     34.0044 -118.8069
station.337     34.442  -118.5821
station.338     34.046  -117.035
station.339     34.614  -118.2893
station.340     34.617  -118.285
station.341     34.2709 -118.7441
station.342     34.6884 -117.8242
station.343     34.6306 -118.305
station.344     33.9214 -116.9806
station.345     33.93   -116.973
station.346     33.96   -116.953
station.347     34.6871 -118.2995
station.348     34.049  -116.9403
station.349     34.7031 -118.4257
station.350     34.281  -118.9047
station.351     34.6822 -118.574
station.352     34.6145 -118.7235
station.353     34.1082 -119.0659
station.354     34.8788 -118.2286
station.355     34.9906 -118.3416

Hercules Station Mapping

Since the Hercules station list was sorted by epicentral distance, we need to map the desired station names into the Hercules stations.

  • Short distance
    • FUL --> Hercules station.5
    • BRE --> Hercules station.23
    • DLA --> Hercules station.41
  • Mid distance
    • 24945 --> Hercules station.222
    • GR2 --> Hercules station.223
    • BHP --> Hercules station.224
  • Large distance
    • Q0049 --> Hercules station.334
    • AGO --> Hercules station.335
    • PTD --> Hercules station.336

AWP Data Files

The AWP results are located at '/u/sciteam/wsavran/scratch/la_habra/1d’ on Blue Waters. You will find 3 files named seis_[x,y,z]0075000, which correspond to the x, y, and z components respectively. These components are aligned with the simulation domain, so they will need to be rotated into the E-W and N-S directions if desired.

The AWP files are 4-byte little-endian floating point numbers stored fast time, then station. The time-series are 3000 points long with a dt of 0.025s resulting in 75.0 seconds of simulation time. We have extracted seismograms at the 356 stations listed in the order they are present on the SCEC wiki page:

The files labeled seis_[x,y,z]0075000 are the aggregated seismograms, meaning they are already extracted from the simulation output. They are 4 byte little-endian floating point numbers with fast-time. There are 356 stations with 3,000 time points each. The ordering of the stations is the same as the station list provided at the SCEC High-F wiki page. They are the [x,y,z] components of ground motion aligned with the simulation domain.

AWP File Format

The 3 files above were processed by William, and split into 356 separate files, one for each station. These files were names <station_id>-AWP-ODC.dat, where station_id is the id listed for each station in the SCEC station wiki page linked above. The processed files are velocity seismograms in ASCII format, using the following format:

  • C1 (Time)
  • C2 (X_Velocity m/s)
  • C3 (Y_Velocity m/s)
  • C4 (Z_Velocity m/s)

They are in a raw form, meaning they have not been rotated or filtered as this should be done carefully to be the exact same between both simulations.

Generating Comparison Plots

In order to generate comparison plots using SAC, the following steps are needed.

Processing AWP files

AWP files are already in a 4-column timestamp, data, data, data format, but need to be converted so that X and Y components match the order in the BBP files. We can use the awp2bbp tool for that:

  • awp2bbp.py CI_FUL-AWP-ODC.dat CI_FUL-AWP-ODC.bbp
  • awp2bbp.py CI_BRE-AWP-ODC.dat CI_BRE-AWP-ODC.bbp
  • awp2bbp.py CI_DLA-AWP-ODC.dat CI_DLA-AWP-ODC.bbp

Then, we can use the bbp2sac utility to convert these files into 3 sac files (X/Y/Z), using the following commands:

  • bbp2sac.py CI_FUL-AWP-ODC.bbp
  • bbp2sac.py CI_DLA-AWP-ODC.bbp
  • bbp2sac.py CI_BRE-AWP-ODC.bbp

Processing Hercules files

Hercules files first need to be converted into a 4-colum timestamp, data, data, data format. We use a new "hercules2bbp.py" script that reads each line and selects the columns containing the velocity data:

  • hercules2bbp.py station.5 CI_FUL-Hercules-vel.bbp
  • hercules2bbp.py station.23 CI_BRE-Hercules-vel.bbp
  • hercules2bbp.py station.41 CI_DLA-Hercules-vel.bbp

After this step is completed, we can use the bbp2sac.py tool to convert the bbp files into sac format:

  • bbp2sac.py CI_FUL-Hercules-vel.bbp
  • bbp2sac.py CI_BRE-Hercules-vel.bbp
  • bbp2sac.py CI_DLA-Hercules-vel.bbp

Plotting using SAC

$ sac
SAC>  r ./CI_FUL-*.000.sac
SAC> xlim 0 8
SAC> ylim -0.25 0.30
SAC> xlabel "Time (s)"
SAC> ylabel "m/s"
SAC> color r
SAC> grid on
SAC> title "STATION FUL - Y component"
SAC> p1
SAC> saveimg FUL-Y.pdf
SAC> r ./CI_FUL-*.090.sac
SAC> title "STATION FUL - X component"
SAC> p1
SAC> saveimg FUL-X.pdf
SAC> r ./CI_FUL-*.ver.sac
SAC> title "STATION FUL - Z component"
SAC> p1
SAC> saveimg FUL-Z.pdf

Plotting FFT using SAC

$ sac
SAC> r ./CI_FUL-*.000.sac
SAC> color r increment list red blue
SAC> xlim 0.1 25
SAC> xlabel "Frequency"
SAC> ylabel "FAS"
SAC> loglog
SAC> fft
SAC> title "Station FUL - Y component"
SAC> p2
SAC> saveimg FUL-Y-FFT.pdf
SAC> r ./CI_FUL-*.090.sac
SAC> fft
SAC> title "Station FUL - X component"
SAC> p2
SAC> saveimg FUL-X-FFT.pdf
SAC> r ./CI_FUL-*.ver.sac
SAC> fft
SAC> title "Station FUL - Z component"
SAC> p2
SAC> saveimg FUL-Z-FFT.pdf

Basic Simulation Parameters

Hercules

  • Computer used: Blue Waters
  • Max freq: 4Hz
  • Min Vs: 500 m/s
  • Number of cores used: 300 * 32 = 9600
  • Time to complete simulation: ( 11 hours 5 minutes 8.9 seconds )
  • Estimated Cost (Node hours or SUs): 300*11 = 3300 (Node hours)
  • simulated time: 100 s
  • Total elements : 6178209792
  • Number of time steps : 40000
  • Compute/Communicate Ratio : 13.84
  • Minimum size of elements : 10.9863 m
  • Maximum size of elements : 87.8906 m

AWP / ODP

  • Computer used: Blue Waters
  • Max freq: 4Hz
  • Min Vs: 500 m/s
  • Number of cores used: 97,920
  • Time to complete simulation: ( 10 hours 30 minutes )
  • Estimated Cost (Node hours or SUs): 32,000 SU
  • simulated time: s
  • Total elements : 97,920,000,000
  • Number of time steps : 100,000
  • Compute/Communicate Ratio : n/a (run stopped short, did not write i/o or compute statistics.
  • Spatial Discretization: 20 m

SORD

  • Computer used: Blue Waters
  • Max freq: ~8 Hz
  • Min Vs: 860 m/s
  • Number of cores used: 16,384
  • Machine time: 2.6 Hrs
  • Node hours: 1024*2.6=2662
  • Simulation time: 12 sec
  • Total elements: 14.35 billion
  • Number of time steps: 9600
  • Average element size: ~16 m

BW Storage Paths

  • /home/scec-01/wsavran/qs010_qp2qs_1d_smth/
  • /home/scec-00/scottcal/N_Lahabra_sim_4Hz_May13/outputfiles/stations/

Related Entries