La Habra Simulation Region

From SCECpedia
Jump to navigationJump to search

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

Simulation Region Map

  • Base lahabra-win1.png

Related Entries