Difference between revisions of "UCVM Review"

From SCECpedia
Jump to navigationJump to search
 
(50 intermediate revisions by 3 users not shown)
Line 1: Line 1:
This entry defines test cases used to document performance of the UCVMC v17.1 release.
+
== Basin Depth plots at -118.76, 34.0 ==
 +
*[[Depth plot Review]]
  
== Test Case ucvm2mesh and ucvm2mesh-mpi tests ==
+
== Projection Correction issues ==
* The goal of this test is to construct AWP test meshes using both ucvm2mesh and ucvm2mesh-mpi.
+
*[[Final Projection Review Result]]
#Test that shows equivalent input files for ucvm2mesh and ucvm2mesh-mpi. Also includes the job submit script for mpi job, and defines the restrictions on number of nodes.
+
*[[CCA06 Projection Review]]
#Test that shows the mesh values returned by these methods are identical.
 
#Update ucvmc Installation instructions need to be updated to describe how to create mpi version if running on usc hpc.
 
  
== Test Case for basin_query serial and mpi tests ==
+
== Result of latest Basin Query ==
* The goal of this test is to evaluate the results from basin_query and basin_query_mpi to demonstrate the produce the same results.
+
*[[Final CCA06 Crossings Result]] Generated using updated UCVM v17.1 with CCA06 with updated projection fix
#Test case that queries basin_query and basin_query_mpi showing that they return the same depth values.
+
*[[Basin Query Depth Profile and Cross Section Result]]
 +
*[[Basin Query Result]] - Generated using UCVM v17.1 updated with original CCA06
  
== Test Case for Basin Query depth Algorithm  ==
+
== CCA Model Plots and Notes ==
* Goal of this test is to document the algorithms available in ucvmc, and show example return values for each.
+
*[[Basin Query Notes]]
#Test that confirms basin_query and basin_query_mpi return the same value for the same point. Document which of the basin_query values is returned by basin_query_mpi
+
*[[CCA06 Basin Depth]]
#Test that shows what basin_query and basin_query_mpi return for points outside the defined volume of the model
+
*[http://hypocenter.usc.edu/research/ucvmc/CCA0to2k.pptx CCA Plots 0k to 2k]
#Test that shows how basin_query and basin_query_mpi support model tiling.
+
*[http://hypocenter.usc.edu/research/ucvmc/CCA0to2k_gtl.pptx CCA Plots 0k to 2k w/ GTL]
 +
*[http://hypocenter.usc.edu/research/ucvmc/CCA0to70k.pptx CCA Plots 0k to 70k]
  
== Test Case for Blue Waters and Titan ==
+
== Basin Depth Algorithm Description ==
#Install UCVMC 17.1 on Blue Water and Titans.
+
*[http://hypocenter.usc.edu/research/ucvm/17.3.0/docs/availablemodels.html#operators Basin Depth Algorithm]
#Regenerate ucvm2mesh_mpi meshes and confirm them match ones from USC
+
*[https://scec.usc.edu/scecpedia/Study_17.3_Data_Products#Velocity_Maps Basin Depth Maps - OpenSHA]
 +
*[http://scec.usc.edu/scecpedia/CCA06_Basin_Depth CCA Basin Depth]
 +
*[http://scec.usc.edu/scecpedia/CVM4.26_Basin_Depth CVM426 Basin Depth]
 +
*[https://scec.usc.edu/scecpedia/File:Summary_of_z_values.pptx Summary of Z values (GMPEs) and related verification and science tasks]
  
== Test Case for Tomography Utility Projection ==
+
== Basin_Query Previous Results ==
#Test cases showing that projection issue in CCA06 is fixed. Confirm this is in latest UCVMC
+
Binary Z files used in OpenSHA, along with some metadata, are posted here:
  
== Test Case for Basin Region Specs ==
+
*4.26-M01: https://github.com/opensha/opensha-commons/tree/master/src/resources/data/site/CVM4i26
*The goal of this test is to reproduce the z1.0 and z2.5 basin models for three California velocity models. These binary files will be given to OpenSHA for use in GMPE calculations. Results we return will describe the algorithm they use to determined the returned depth value.
+
*CCA06: https://github.com/opensha/opensha-commons/tree/master/src/resources/data/site/CCAi6
*CCA: [minLat = 33.35, minLon = -123.0, nx = 1551, ny = 1201, gridSpacing = 0.005 degrees]
+
 
*CVM-S4.26 and CVM-S4.26-M01: [minLat = 31.0, minLon = -121.0, nx = 1701, ny = 1101, gridSpacing = 0.005 degrees]
+
We plan to implements a way to read these files and plot them in python which can be included with UCVM. Plotting would be similar to existing horizontal slide plotting scripts (which are likely included in both versions of UCVM, though I have not used them), and will be a useful feature going forward. Once that tool exists, comparisons will become easier. If in the meantime help is needed extracting individual values from these files, there is an OpenSHA tool which can access them (via Java).
*Output format for the basin data sets is a binary regular mesh format with separate files for Z1.0 and Z2.5 in meters
 
*Output binary files are checked by plotting and comparing with previous versions
 
*The binary files are 4 byte float, little endian, fast XY, no rotation.  
 
*The first data point is bottom left (at minX, minY).
 
*The formulas for lat/lon at a given x/y index:
 
<pre>
 
lon = xIndex * gridSpacing + minLon
 
lat = yIndex * gridSpacing + minLat
 
</pre>
 
<pre>
 
*The formula for calculating the file position for a given xIndex/yIndex:
 
filePos = 4 * (nx * yIndex + xIndex)
 
</pre>
 
  
 
== Velocity mesh plotting scripts from CyberShake codebase ==
 
== Velocity mesh plotting scripts from CyberShake codebase ==
 +
As part of the CyberShake work, we've developed a few Python scripts for plotting horizontal cross-sections (plot_horiz_slice.py), vertical cross-sections (plot_vertical_slice.py), and profile plots (plot_profile.py).  The scripts are available in the CyberShake SVN repository, at http://source.usc.edu/svn/cybershake/import/trunk/UCVM/smoothing/ .  As input, they expect a mesh in AWP format, and also a [[CyberShake_Code_Base#Coord|model coords]] file, which contains a listing of the lat/lon values for each surface point in the mesh (this is used for decimation purposes). Sample plots produced using these scripts can be seen [[Study_16.9_Velocity_Plots|here]].
  
As part of the CyberShake work, we've developed a few Python scripts for plotting horizontal cross-sections (plot_horiz_slice.py), vertical cross-sections (plot_vertical_slice.py), and profile plots (plot_profile.py).  The scripts are available in the CyberShake SVN repository, at http://source.usc.edu/svn/cybershake/import/trunk/UCVM/smoothing/ .  As input, they expect a mesh in AWP format, and also a [[CyberShake_Code_Base#Coord|model coords]] file, which contains a listing of the lat/lon values for each surface point in the mesh (this is used for decimation purposes).
+
== Review of Vp/Vs Ratio Checks in UCVM ==
 
+
*[[UCVM Vp/Vs Ratio]]
Sample plots produced using these scripts can be seen [[Study_16.9_Velocity_Plots|here]].
 
 
 
== Introductory Presentations ==
 
*[http://hypocenter.usc.edu/research/UCVM/UCVM_Intro.zip UCVM PPT Files]
 
*[http://srl.geoscienceworld.org/content/early/2017/09/04/0220170082 UCVM SRL Paper online]
 
 
 
== UCVM-Python Tutorials ==
 
Training Day 1:
 
*[http://hypocenter.usc.edu/research/Goulet/ucvmTrainingVideos/ Session 1]
 
 
 
Additional training found on hypocenter (4 videos in the root here):
 
*[http://hypocenter.usc.edu/research/ucvm/17.3.0/ Sessions 2-5]
 
 
 
== UCVM Scripts and Utilities ==
 
*[http://hypocenter.usc.edu/research/UCVM/scripts.zip UCVM Scripts]
 
*[http://hypocenter.usc.edu/research/UCVM/pbs.zip UCVM PBS Scripts]
 
  
 
== Related Entries ==
 
== Related Entries ==
 +
*[[UCVM Introduction]]
 +
*[[Planned UCVM Test Cases]]
 
*[[UCVM]]
 
*[[UCVM]]
 
*[[USR]]
 
*[[USR]]
 
*[[CVM-S]]
 
*[[CVM-S]]
 
*[[CVM-H]]
 
*[[CVM-H]]

Latest revision as of 23:32, 29 January 2019

Basin Depth plots at -118.76, 34.0

Projection Correction issues

Result of latest Basin Query

CCA Model Plots and Notes

Basin Depth Algorithm Description

Basin_Query Previous Results

Binary Z files used in OpenSHA, along with some metadata, are posted here:

We plan to implements a way to read these files and plot them in python which can be included with UCVM. Plotting would be similar to existing horizontal slide plotting scripts (which are likely included in both versions of UCVM, though I have not used them), and will be a useful feature going forward. Once that tool exists, comparisons will become easier. If in the meantime help is needed extracting individual values from these files, there is an OpenSHA tool which can access them (via Java).

Velocity mesh plotting scripts from CyberShake codebase

As part of the CyberShake work, we've developed a few Python scripts for plotting horizontal cross-sections (plot_horiz_slice.py), vertical cross-sections (plot_vertical_slice.py), and profile plots (plot_profile.py). The scripts are available in the CyberShake SVN repository, at http://source.usc.edu/svn/cybershake/import/trunk/UCVM/smoothing/ . As input, they expect a mesh in AWP format, and also a model coords file, which contains a listing of the lat/lon values for each surface point in the mesh (this is used for decimation purposes). Sample plots produced using these scripts can be seen here.

Review of Vp/Vs Ratio Checks in UCVM

Related Entries