Test Cases 1D Models
The model/1d/1d.conf file defines the SCEC 1d model. The material properties from querying 1d model is interpolated. The material properties will vary gradually between layers, showing interpolation. This model has several layers between 0-10km.
The input file, test_latlons_1d_pts.txt is a list of lat/lon with depths from 0 to 11K at 100m steps.
./bin/ucvm_query -f ./conf/ucvm.conf -m 1d < ./tests/inputs/test_latlons_1d_pts.txt > 1d_res.txt
The depth plot of 1D model shows smoothly varying material properties between each defined layer.
./plot_depth_profile.py -s 34,-118 -b 0 -e 20000 -d vs,vp,density -v 100 -c 1d -o 1d.png
The bbp1d model uses the defaults in model/1d/bbp1d.conf model configuration file.
The default interpolation in bbp1d.conf is set to linear,
# Interpolation method. # Set interpolation=none to turn off interpolation # Set interpolation=linear to turn on interpolation interpolation=linear
Query for material properties using BBP 1d model,
./bin/ucvm_query -f ./conf/ucvm.conf -m bbp1d < ./tests/inputs/test_latlons_1d_pts.txt > bbp1d_res.txt
The depth plot of BBP 1D model also has Layers with smoothly varying material properties values,
./plot_depth_profile.py -s 34,-118 -b 0 -e 20000 -d vs,vp,density -v 100 -c bbp1d -o bbp1d.png
When the bb1d.conf interpolation flag is changed to none, The material properties will show discrete steps in returning material properties in between layers.
# Interpolation method. # Set interpolation=none to turn off interpolation # Set interpolation=linear to turn on interpolation interpolation=none
Running query with BBP 1d model,
./bin/ucvm_query -f ./conf/ucvm.conf -m bbp1d < ./tests/inputs/test_latlons_1d_pts.txt > bbp1d_interpolation_none_res.txt
Plotting the depth profile shows a distinct step between two layers,
./plot_depth_profile.py -s 34,-118 -b 0 -e 20000 -d vs,vp,density -v 100 -c bbp1d -o bbp1d_none.png
CVM-S4 model returns default interpolated 1D model values when outside of the defined CVM-S4 region.
Using Phoenix point which is outside of CVM-S4 region to verify the returned material values are the interpolated 1D values.
./bin/ucvm_query -f ./conf/ucvm.conf -m cvms < ./tests/inputs/test_latlons_cvms_bkg_pts.txt > cvms_bkg_res.txt
- Input Points Outside of CVM-S4 defined Region
- Results from CVM-S4 for Points outside of defined Region
Plotting the cross section across the region shows a distinct change in material properties,
./plot_cross_section.py -b 33.5,-118 -u 33.5,-112.1 -h 1000 -v 100 -d vs -c cvms -a s -s 0.0 -e 20000 -o cvms_1d.png
The CVM-S4.26 model does not have a background model. We define a point outside of the model region, and expect the model to return no values.
-bash-4.2$ ./bin/ucvm_query -f ./conf/ucvm.conf -m cvms5 < ./tests/inputs/test_latlons_cvms_bkg_pts.txt > cvms5_bkg_res.txt
Plotting the cross section across the region,
./plot_cross_section.py -b 33.5,-118 -u 33.5,-112.1 -h 1000 -v 100 -d vs -c cvms5 -a s -s 0.0 -e 20000 -o cvms5_1d.png
We run ucvm.query, specifying the cvm-s4.26.m01 model (cvmsi), for depths from 0-10km. We expect that model will return interpolated values from a 1d background model
We will use a point outside of the define CVM-S4 region. We will use the Phoenix point defined earlier, input to ucvm_query, and expect CVM-S4.26.M01 to return interpolated 1D values.
-bash-4.2$ ./bin/ucvm_query -f ./conf/ucvm.conf -m cvmi < ./tests/inputs/test_latlons_cvms_bkg_pts.txt > cvmsi_bkg_res.txt
Plotting the cross section across the region and now the right side of cross section plot this used to be empty is now filled with the interpolated background 1D values,
./plot_cross_section.py -b 33.5,-118 -u 33.5,-112.1 -h 1000 -v 100 -d vs -c cvmsi -a s -s 0.0 -e 20000 -o cvmsi_1d.png
