La Habra Simulation Region
Contents
Proposed Small Simulation Region
High-F modelers have proposed to design a "small" 3D model for La Habra EQ that includes the closest 3 sites. We then run tests with this model until the behavior of the codes and agreement of the results, to move forward in a much more efficient and practical manner (as opposed to trying to run test cases with the full 3D volume).
An additional benefit is that using a small model would make testing/iterating on UCVM development much more efficient. For example, modifications for RWG mesh output and testing/implementing velocity thresholds. As a potential starting point, here are the model parameters RWG used for his latest runs:
- set MODEL_LON = -117.99503 # LON at center of model region
- set MODEL_LAT = 33.86771 # LAT at center of model region
- set MODEL_ROT = -50.1 # counter-clockwise rotation of X_rwg from default azimuth of 090
- set (39.9 azimuth) XLEN = 28.0
- set (129.9) YLEN = 28.0
- set (points down) ZLEN = 12.0
- set HH = 0.020
The map below shows this region (black rectangle) along with the epicenter and recording stations. Green triangles indicate sites where we have been comparing waveforms (FUL, DLA, BRE are inside this small region). The 4 corners of the proposed region (magenta circles) are:
- c1= -118.20819 33.85173
- c2= -118.01410 34.04473
- c3= -117.78178 33.88333
- c4= -117.97604 33.69068
The above parameters give NX=1400, NY=1400, NZ=600, and and with set NT=60000 (60 sec). This took about 4 hours to run on the Caltech cluster using 864 processors, not counting queue time.
Current UCVM Configuration File
# List of CVMs to query ucvmlist=cvmsi # UCVM conf file ucvmconf=../conf/ucvm.conf # Gridding cell centered or vertex gridtype=CENTER # Spacing of cells spacing=20.0 # Projection proj=+proj=utm +datum=WGS84 +zone=11 rot=-39.9 x0=-118.20819 y0=33.85173 z0=0.0 # Number of cells along each dim nx=1400 ny=1400 nz=600 # Partitioning of grid among processors px=25 py=25 pz=20 # Vs/Vp minimum vp_min=200 vs_min=0
