UCVM shakeout2 velocity model
From SCECpedia
Model overview
The shakeout2 velocity model for the ShakeOut Scenario 2.0 is based on CVM-S4.26.M01(with embedded GTL layer for basin regions) and GTL Tapering for all non-basin area.
Tapering's z range is from 0m to 700m, and material property floors used are vs=200m, vp=500m and density=200m a customized z coordinates list is pre-defined
Creating the NetCDF model mesh
UCVM's ucvm2mesh is used to generate 3 binary mesh data files: vs, vp and density for ROI
In 'debug mode', 3 coordinate lists are created: one for the lats, the lons and the z depths.
With these data files and lists, a simple NetCDF file is created
ucvm2mesh configuration file:
# List of CVMs to query ucvmlist=cvmsi,elygtl:taper # (optional) zranges: -Z zmin,zmax min_zrange=0.0 max_zrange=700.0 # (optional) floors: -L vs_floor,vp_floor,density_floor vs_floor=200.0 vp_floor=500.0 density_floor=200.0 # UCVM conf file ucvmconf=YOUR_UCVM_INSTALL_PATH/conf/ucvm.conf # Gridding cell: CENTER or VERTEX gridtype=VERTEX # Spacing of cells (m) spacing=1000.0 # (optional) Spacing of z depth (m) # if not set, z_spacing = spacing #z_spacing = 500.0 # (optional) Allow varying z depth # format: 1 row 1 target depth z_file = so2_zlist # Projection proj=+proj=utm +datum=WGS84 +zone=11 #to rotate, rot=-39.9 rot=0 # Origin # x,y: latlon # z: depth (m) x0=-120.033556 y0=31.869638 z0=0.0 # Number of cells along each dim nx=516 ny=468 nz=152 # how to determine distance in meter for each dim # x = nx * spacing # y = ny * spacing # z = nz * (spacing or z_spacing), or = z_file[idx]) # Partitioning of grid among processors # relevant to chunking, should be divisible into nx,ny,nz px=2 py=2 pz=2 # Vs/Vp minimum vp_min=0 vs_min=0 # Mesh grid files meshfile=so2.media gridfile=so2.grid # Mesh Grid formats # IJK-12, IJK-20, IJK-32 : single triplet file # SORD : separate vs, vp, rho files meshtype=SORD # Location of scratch dir scratch=./scratch
ucvm2mesh z depths list: so2_zlist
latitude coordinate list: so2_xlist
longitude coordinate list: so2_ylist
ucvm2mesh command used:
