Skip to content
This repository has been archived by the owner on Feb 8, 2022. It is now read-only.

Test Cases Tiling

meihuisu edited this page May 29, 2019 · 25 revisions

Model Tiling

UCVM support tiling multiple CVMs together. Tiling specified the query order. The ucvm_query interface will try to retrieve material properties from each specified CVM, in the order given and returns the material properties from the first model that defines properties for the specified location.

Model tiling at a Lemore site

A point at Lemoore, CA is a location where three models overlap, CVM-S4, CVM-H, and Cencal. A tiled CVMs query at Lemoore returns different material properties based on the CVM tiling order. The material properties returned are from the first model in the list that has valid material properties.

Test Points in test_latlons_lemoore_pts.txt,

-119.7785 36.3121 0.0
-119.7785 36.3121 100.0
-119.7785 36.3121 200.0
-119.7785 36.3121 500.0
-119.7785 36.3121 1000.0
-119.7785 36.3121 5000.0
-119.7785 36.3121 10000.0
-119.7785 36.3121 25000.0
-119.7785 36.3121 50000.0
-119.7785 36.3121 100000.0

Order : CVM-S4, CVM-H, Cencal

./bin/ucvm_query -f ./conf/ucvm.conf -m cvms,cvmh,cencal < ../test/inputs/test_latlons_lemoore_pts.txt > cvms_cvmh_cencal.txt

Results

-119.7785    36.3121      0.000     70.983    280.000       cvms   5508.030   3196.060   2734.718       none      0.000      0.000 
     0.000      crust   5508.030   3196.060   2734.718
 -119.7785    36.3121    100.000     70.983    280.000       cvms   5508.030   3196.060   2734.718       none      0.000      0.000 
     0.000      crust   5508.030   3196.060   2734.718
 -119.7785    36.3121    200.000     70.983    280.000       cvms   5508.030   3196.060   2734.718       none      0.000      0.000 
     0.000      crust   5508.030   3196.060   2734.718
 -119.7785    36.3121    500.000     70.983    280.000       cvms   5508.030   3196.060   2734.718       none      0.000      0.000 
     0.000      crust   5508.030   3196.060   2734.718
 -119.7785    36.3121   1000.000     70.983    280.000       cvms   5508.061   3196.078   2734.723       none      0.000      0.000 
     0.000      crust   5508.061   3196.078   2734.723
 -119.7785    36.3121   5000.000     70.983    280.000       cvms   6013.406   3491.435   2814.517       none      0.000      0.000 
     0.000      crust   6013.406   3491.435   2814.517
 -119.7785    36.3121  10000.000     70.983    280.000       cvms   6561.158   3855.761   2901.007       none      0.000      0.000 
     0.000      crust   6561.158   3855.761   2901.007
 -119.7785    36.3121  25000.000     70.983    280.000       cvms   6760.734   4024.399   2932.520       none      0.000      0.000 
     0.000      crust   6760.734   4024.399   2932.520
 -119.7785    36.3121  50000.000     70.983    280.000       cvms   7920.000   4572.614   3115.568       none      0.000      0.000 
     0.000      crust   7920.000   4572.614   3115.568
 -119.7785    36.3121 100000.000     70.983    280.000       cvms   8116.667   4686.160   3146.622       none      0.000      0.000 
     0.000      crust   8116.667   4686.160   3146.622

Order : CVM-H, Cencal, CVM-S4

-bash-4.2$ ./bin/ucvm_query -f ./conf/ucvm.conf -m cvmh,cencal,cvms < lemore_test_pts.txt > cvmh_cencal_cvms.txt

Results

 -119.7785    36.3121      0.000     70.983    280.000       cvmh   5258.205   2938.976   2576.084       none      0.000      0.000 
     0.000      crust   5258.205   2938.976   2576.084
 -119.7785    36.3121    100.000     70.983    280.000       cvmh   5266.341   2943.600   2577.442       none      0.000      0.000 
     0.000      crust   5266.341   2943.600   2577.442
 -119.7785    36.3121    200.000     70.983    280.000       cvmh   5275.458   2948.787   2578.968       none      0.000      0.000 
     0.000      crust   5275.458   2948.787   2578.968
 -119.7785    36.3121    500.000     70.983    280.000       cvmh   5325.177   2976.568   2587.370       none      0.000      0.000 
     0.000      crust   5325.177   2976.568   2587.370
 -119.7785    36.3121   1000.000     70.983    280.000       cvmh   5395.763   3016.100   2599.540       none      0.000      0.000 
     0.000      crust   5395.763   3016.100   2599.540
 -119.7785    36.3121   5000.000     70.983    280.000       cvmh   5923.339   3312.400   2700.424       none      0.000      0.000 
     0.000      crust   5923.339   3312.400   2700.424
 -119.7785    36.3121  10000.000     70.983    280.000       cvmh   6118.245   3427.836   2742.522       none      0.000      0.000 
     0.000      crust   6118.245   3427.836   2742.522
 -119.7785    36.3121  25000.000     70.983    280.000       cvmh   6584.390   3699.164   2854.530       none      0.000      0.000 
     0.000      crust   6584.390   3699.164   2854.530
 -119.7785    36.3121  50000.000     70.983    280.000       cvmh   7688.477   4332.974   3183.304       none      0.000      0.000 
     0.000      crust   7688.477   4332.974   3183.304
 -119.7785    36.3121 100000.000     70.983    280.000       cvmh   7910.747   4483.328   3259.537       none      0.000      0.000 
     0.000      crust   7910.747   4483.328   3259.537

Order : Cencal, CVM-S4, CVM-H

./bin/ucvm_query -f ./conf/ucvm.conf -m cencal,cvms,cvmh < lemore_test_pts.txt > cencal_cvms_cvmh.txt

Results

 -119.7785    36.3121      0.000     70.983    280.000     cencal   1820.000    510.000   2020.000       none      0.000      0.000 
     0.000      crust   1820.000    510.000   2020.000
 -119.7785    36.3121    100.000     70.983    280.000     cencal   1940.000    570.000   2060.000       none      0.000      0.000 
     0.000      crust   1940.000    570.000   2060.000
 -119.7785    36.3121    200.000     70.983    280.000     cencal   2060.000    650.000   2090.000       none      0.000      0.000 
     0.000      crust   2060.000    650.000   2090.000
 -119.7785    36.3121    500.000     70.983    280.000     cencal   2420.000    920.000   2170.000       none      0.000      0.000 
     0.000      crust   2420.000    920.000   2170.000
 -119.7785    36.3121   1000.000     70.983    280.000     cencal   2860.000   1290.000   2260.000       none      0.000      0.000 
     0.000      crust   2860.000   1290.000   2260.000
 -119.7785    36.3121   5000.000     70.983    280.000     cencal   5900.000   3500.000   2880.000       none      0.000      0.000 
     0.000      crust   5900.000   3500.000   2880.000
 -119.7785    36.3121  10000.000     70.983    280.000     cencal   5900.000   3500.000   2880.000       none      0.000      0.000 
     0.000      crust   5900.000   3500.000   2880.000
 -119.7785    36.3121  25000.000     70.983    280.000     cencal   6740.000   3880.000   3040.000       none      0.000      0.000 
     0.000      crust   6740.000   3880.000   3040.000
 -119.7785    36.3121  50000.000     70.983    280.000       cvms   7920.000   4572.614   3115.568       none      0.000      0.000 
     0.000      crust   7920.000   4572.614   3115.568
 -119.7785    36.3121 100000.000     70.983    280.000       cvms   8116.667   4686.160   3146.622       none      0.000      0.000 
     0.000      crust   8116.667   4686.160   3146.622

Cencal's depth profile at Lemore site,

./plot_depth_profile.py -s 36.3121,-119.7785 -b 0  -e 60000 -d vs -v 1000 -c cencal -o cencal_depth.png

CVM-S4's depth profile,

./plot_depth_profile.py -s 36.3121,-119.7785 -b 0  -e 60000 -d vs -v 1000 -c cvms -o cvms_depth.png

CVM-H's depth profile,

./plot_depth_profile.py -s 36.3121,-119.7785 -b 0  -e 60000 -d vs -v 1000 -c cvmh -o cvmh_depth.png

Cencal's depth boundary ends at 45km while CVM-S4's and CVM-H extend beyond 45km. The tiled CVMs extends and combined the overlapping models in next two depth profiles.

Depth profile of tiling Cencal, CVM-S4, CVM-H,

./plot_depth_profile.py -s 36.3121,-119.7785 -b 0  -e 60000 -d vs -v 1000 -c cencal,cvms,cvmh -o cencal_cvms_cvmh_depth.png

Depth profile of tiling Cencal, CVM-H, CVM-S4,

./plot_depth_profile.py -s 36.3121,-119.7785 -b 0  -e 60000 -d vs -v 1000 -c cencal,cvms,cvmh -o cencal_cvmh_cvms_depth.png

Add a 1D background model to CVM-S4.26 using tiling

CVM-S4.26 has no background model. It is undefined outside a specific region.

./plot_cross_section.py -b 34.5,-116 -u 34.5,-112 -h 1000 -v 100 -d vs -c cvms5 -a d -s 0.0 -e 15000 -o cvms5.png

Tiling with 1d model fills in material properties for the pointed that is not defined in CVM-S4.26.

./plot_cross_section.py -b 34.5,-116 -u 34.5,-112 -h 1000 -v 100 -d vs -c cvms5,1d -a d -s 0.0 -e 15000 -o cvms5_1d.png 

Using query to show the same effect,

Test points

-120.0 37 0 
-120.0 37.5 0 
-120.0 38 0 

No Tiling

./bin/ucvm_query -f ./conf/ucvm.conf -m cvms5

Result,

 -120.0000    37.0000      0.000     92.219    280.000      cvms5    734.408    140.000   1460.649       none      0.000      0.000      0.000      crust    734.408    140.000   1460.649
 -120.0000    37.5000      0.000    816.159    660.136       none      0.000      0.000      0.000       none      0.000      0.000      0.000       none      0.000      0.000      0.000
 -120.0000    38.0000      0.000   1662.498    515.000       none      0.000      0.000      0.000       none      0.000      0.000      0.000       none      0.000      0.000      0.000

Tiling CVM-S4.26 with 1d

./bin/ucvm_query -f ./conf/ucvm.conf -m cvms5,1d

Result,

 -120.0000    37.0000      0.000     92.219    280.000      cvms5    734.408    140.000   1460.649       none      0.000      0.000      0.000      crust    734.408    140.000   1460.649
 -120.0000    37.5000      0.000    816.159    660.136         1d   5000.000   2886.751   2654.500       none      0.000      0.000      0.000      crust   5000.000   2886.751   2654.500
 -120.0000    38.0000      0.000   1662.498    515.000         1d   5000.000   2886.751   2654.500       none      0.000      0.000      0.000      crust   5000.000   2886.751   2654.500

Model order impacts tiled query results

The material properties returned changes depending on the models tile order in a query.

Test point

-120.0 35.0 100

Tiling order: cca,cvms5

./bin/ucvm_query -f ./conf/ucvm.conf -m cca,cvms5

Result

 -120.0000    35.0000    100.000   1368.254    515.000        cca   2499.033   1078.751   2135.382       none      0.000      0.000      0.000      crust   2499.033   1078.751   2135.382

Tiling order: cvms5,cca

./bin/ucvm_query -f ./conf/ucvm.conf -m cvms5,cca

Result

 -120.0000    35.0000    100.000   1368.254    515.000      cvms5   2470.229   1040.747   2120.054       none      0.000      0.000      0.000      crust   2470.229   1040.747   2120.054
Clone this wiki locally