140 likes | 299 Views
Southern California Earthquake Center SCEC UCVM-based Tools. Patrick Small Research Programmer. UCVM-based Tools. 3D Meshes: MPI 3D Mesh generator 3D Viewing: Mesh to netCDF converter 2D Viewing: Python 2D map and profile plotting scripts
E N D
Southern California Earthquake CenterSCEC UCVM-basedTools Patrick Small Research Programmer
UCVM-based Tools • 3D Meshes: MPI 3D Mesh generator • 3D Viewing: Mesh to netCDF converter • 2D Viewing: Python 2D map and profile plotting scripts • Model Differencing: Differencing tool for comparing two models • CVM Evaluation: Simulation based CVM evaluation system
3D Mesh Generation • mesh-create-MPI: • Creates AWP 3/5/8 value meshes and SORD meshes from UCVM • Minimum properties: Vp, Vs, density • 3D regular grid of specified cell size • C code using MPI/IO, UCVM, and Proj.4 • Uses 3D volume partitioning among cores • Supports any Proj.4 map projection • Mesh bounding box specified with (lon, lat) origin, spacing, grid points in each dimension, rotation angle • 435 B point M8 mesh from CVM-H on 3300+ cores in 1.5 hr on NICS Kraken
3D Visualization • mesh2netcdf: • Converts AWP/SORD formatted mesh into netCDF format • C code using netCDF library • Climate and Forecast Metadata Convention used (CF netCDF) • Output file is suitable for viewing in a 3D viewer such as IDV • IDV Download: http://www.unidata.ucar.edu/software/idv/ • IDV Tutorial: http://scec.usc.edu/scecpedia/IDV • Example netCDF mesh: http://scec.usc.edu/scecpedia/CVM-SI
3D Visualization Example • Screen capture of map view from CVM-SI i5 mesh in IDV:
2D Visualization • viz-cvm: • Produces MATLAB-friendly 2D datasets and plots from UCVM or a mesh • Python scripts that use modules numpy, pyproj, matplotlib, and basemap • Since based on UCVM, any model or combination of models may be plotted • Supported plots: • Velocity map plots from UCVM/mesh • Velocity cross-section plots from UCVM/mesh • Z maps: Minimum depth to Z value from UCVM • DEM, Vs30 map plots from UCVM
2D Visualization Example Vs map overlay of USGS Bay Area, CVM-H 11.2, and Lin-Thurber models
Model Differencing • mesh-op: • Subtracts one 3D mesh from another to create a new mesh of the differences. • Supports formats AWP 3/5/8 value and SORD • C code • Resulting difference mesh may be converted to netCDF with mesh2netcdf and viewed in IDV
CVM Evaluation • cvmtest: • System based on Python, AWP-ODC, and SCEC BBand Goodness of Fit on NICS Kraken • Three simulated events with saved observed seismograms: Chino Hills 2008, Gilroy 2002, Oakland 2007 • Can submit any CVM from UCVM to any/all of those events • SCEC BBandGoF executed on the synthetic/observed seismograms from each sim and the following results are compiled: • PSA bias ln(obs/synth) GoF with plot • Map-based PSA Fit • Map-based PGV Fit • Seismogram overlay plots for all stations
CVM-H Evaluation Plots Example fit plots for CVM-H 11.2.0 / Chino Hills
CVM-S Evaluation Plots Example fit plots for CVM-S4 / Chino Hills
CVM-SI Evaluation Plots Example fit plots for CVM-SI i5 / Chino Hills
Resources • Documentation: • http://scec.usc.edu/scecpedia/CVM_Toolkit • Software Download: • cvm2mesh/mesh2netCDF SVN: https://source.usc.edu/svn/cvm2mesh • Requires MPI and netCDF library for mesh2netCDF • viz-cvm SVN: https://source.usc.edu/svn/viz-cvm • Requires Python (numpy, pyproj, matplotlib, basemap) • cvmtest SVN: https://source.usc.edu/svn/cvmtest • Requires Python (numpy, pyproj, matplotlib, basemap), AWP-ODC, SCEC BBand