Difference between revisions of "UCVM Review"

From SCECpedia
Jump to navigationJump to search
 
(54 intermediate revisions by 4 users not shown)
Line 1: Line 1:
These links are information about the current UCVM software.
+
== Basin Depth plots at -118.76, 34.0 ==
 +
*[[Depth plot Review]]
  
== Example PBS scripts ==
+
== Projection Correction issues ==
*[http://hypocenter.usc.edu/research/UCVM/scripts.zip UCVM Scripts]
+
*[[Final Projection Review Result]]
*[http://hypocenter.usc.edu/research/UCVM/pbs.zip UCVM PBS Scripts]
+
*[[CCA06 Projection Review]]
  
== Test UCVMC Install on Blue Waters and Titan ==
+
== Result of latest Basin Query ==
*Install UCVMC 17.1 on Blue Water and Titans.
+
*[[Final CCA06 Crossings Result]] Generated using updated UCVM v17.1 with CCA06 with updated projection fix
*Regenerate ucvm2mesh_mpi meshes and confirm them match ones from USC
+
*[[Basin Query Depth Profile and Cross Section Result]]
 +
*[[Basin Query Result]] - Generated using UCVM v17.1 updated with original CCA06
  
== confirm en-jui projection issue fixed ==
+
== CCA Model Plots and Notes ==
DG made a change to fix projection issue in CCA06. Confirm this is in latest UCVMC
+
*[[Basin Query Notes]]
 +
*[[CCA06 Basin Depth]]
 +
*[http://hypocenter.usc.edu/research/ucvmc/CCA0to2k.pptx CCA Plots 0k to 2k]
 +
*[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]
  
== basin_query z2500 first occurence or second occurence test ==
+
== Basin Depth Algorithm Description ==
* goal of this test is to evaluate the difference between results for the z2500 test. We will test the impact of different methods for selecting the z2500 values from a velocity model. Some algorithms return the first occurence of the specified Vs, or the second occurence, or the last occurence found. We need to document the algorithms available in ucvmc, and show example return values for each.
+
*[http://hypocenter.usc.edu/research/ucvm/17.3.0/docs/availablemodels.html#operators Basin Depth Algorithm]
 +
*[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]
  
== basin_query serial and mpi tests ==
+
== Basin_Query Previous Results ==
* The goal of this test is to evaluate the results from basin_query and basin_query_mpi to demonstrate the produce the same results.
+
Binary Z files used in OpenSHA, along with some metadata, are posted here:
** We determined that the serial code returns
 
  
== mpi ucvm2mesh tests ==
+
*4.26-M01: https://github.com/opensha/opensha-commons/tree/master/src/resources/data/site/CVM4i26
* goal of this test is to construct AWP test meshes using both ucvm2mesh and ucvm2mesh_mpi.  
+
*CCA06: https://github.com/opensha/opensha-commons/tree/master/src/resources/data/site/CCAi6
**Then we want to compare the resulting meshes and make sure they are identical
 
* ucvmc Installation instructions need to be updated to describe how to create mpi version if running on usc hpc.
 
  
== Basin region Specs ==
+
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).
CCA:
 
*minLat = 33.35
 
*minLon = -123
 
*nx = 1551
 
*ny = 1201
 
*gridSpacing = 0.005 m
 
 
 
CVM-S4.26:
 
*minLat = 31;
 
*minLon = -121;
 
*nx = 1701;
 
*ny = 1101;
 
*gridSpacing = 0.005 degrees
 
 
 
# define depth
 
When we extract these basin data sets as binary regular mesh format, we will check them by plotting and comparing with previous versions. We need separate files for Z1.0 and Z2.5 in meters.
 
 
 
The binary files are 4 byte float, little endian, fast XY, no rotation. Here are the formulas for lat/lon at a given x/y index:
 
<pre>
 
lon = xIndex * gridSpacing + minLon
 
lat = yIndex * gridSpacing + minLat
 
</pre>
 
 
 
and for calculating the file position for a given xIndex/yIndex:
 
<pre>
 
filePos = 4 * (nx * yIndex + xIndex)
 
</pre>
 
 
 
The first data point is bottom left (at minX, minY).
 
  
 
== 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]]
== 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]
 
  
 
== Related Entries ==
 
== Related Entries ==
 +
*[[UCVM Introduction]]
 +
*[[Planned UCVM Test Cases]]
 
*[[UCVM]]
 
*[[UCVM]]
 +
*[[USR]]
 +
*[[CVM-S]]
 +
*[[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