Skip to content

plot_vs30_map

Philip Maechling edited this page Sep 1, 2022 · 1 revision

plot_vs30_map.py

A UCVM utility script called plot_vs30_map does a few things. It discritizes a region, and then it queries the specific grid points, queries data values from 0-29m, adds them using a slowness calculation, and plots a map with the vs30 values shown.

The script has a query interface, but selected capabilities, such as an output file name, and an output image name are not available through that interface. Standard interface is a command line.

Running plot_vs30_map

In addition to the png plot, a metadata file of the run and an output data file is generated. This output data can be reused as input data to plot_vs30_map with different display parameters to create a new png plot base on the same set of input data.

 -bash-4.2$ ./plot_vs30_map.py -b 30.5,-126.0 -u 42.5,-112.5 -s 0.2 -a s -c cs173 -o vs30_cs173_map.png
Using parameters:
datafile = NoneUsing parameters:
lat1 = 30.5
lat2 = 42.5
color = s
spacing = 0.2
outfile = vs30_cs173h_map.png
cvm = cs173h
lon1 = -126.0
lon2 = -112.5
Retrieving data. Please wait...

Run using Input File

Run plot_vs30_map.py with input data set generated in the earlier example, and re-create map with discrete color scale.

-bash-4.2$ ./plot_vs30_map.py -b 30.5,-126.0 -u 42.5,-112.5 -s 0.2 -a d -c cs173 -o vs30_cs173h_d_map.png -f vs30_cs173h_data.bin

Clone this wiki locally