View source for Registering CS173 into UCVM

This page details the steps involved in delivering a velocity mesh to John Shaw and Andreas Plesch, so that the Harvard group's Central California basins can be integrated into a velocity model, which can be fed back into UCVM and used for CyberShake.

== Test Points ==
To evaluate the registration of CS17.3 into UCVM we are querying a set of point using Scott's original CS17.3 and the recently updated UCVMC version. The following test points are sites from the CS17.3 study. We will query each point at 5 depths 0km, 100m, 1km, 10km, and 40km

<pre>
<name>Osito Audit: Castaic Lake Dam (OSI)</name>
<coordinates>-118.723500,34.614500</coordinates>

<name>Santa Barbara (SBR)</name>
<coordinates>-119.698189,34.420831</coordinates>

<name>s1356 (s1356)</name>
<coordinates>-120.582740,34.968060</coordinates>

<name>Calstate Bakersfield (BAK)</name>	
<coordinates>-119.104450,35.344440</coordinates>

<name>Lime Mtn. Quarry (LMQ)</name>
<coordinates>-120.990669,35.668835</coordinates>

<name>Parkfield (PARK)</name>
<coordinates>-120.432800,35.899700</coordinates>

<name>Huron (HUR)</name>
<coordinates>-120.102917,36.202731</coordinates>

<name>Lemoore (LEM)</name>
<coordinates>-119.782911,36.300783</coordinates>

</pre>

== Test points and depths list ==
<pre>
-118.723500 34.614500 0.0
-118.723500 34.614500 100.0
-118.723500 34.614500 1000.0
-118.723500 34.614500 10000.0
-118.723500 34.614500 40000.0
-119.698189 34.420831 0.0
-119.698189 34.420831 100.0
-119.698189 34.420831 1000.0
-119.698189 34.420831 10000.0
-119.698189 34.420831 40000.0
-120.582740 34.968060 0.0
-120.582740 34.968060 100.0
-120.582740 34.968060 1000.0
-120.582740 34.968060 10000.0
-120.582740 34.968060 40000.0
-119.104450 35.344440 0.0
-119.104450 35.344440 100.0
-119.104450 35.344440 1000.0
-119.104450 35.344440 10000.0
-119.104450 35.344440 40000.0
-120.990669 35.668835 0.0
-120.990669 35.668835 100.0
-120.990669 35.668835 1000.0
-120.990669 35.668835 10000.0
-120.990669 35.668835 40000.0
-120.432800 35.899700 0.0
-120.432800 35.899700 100.0
-120.432800 35.899700 1000.0
-120.432800 35.899700 10000.0
-120.432800 35.899700 40000.0
-120.102917 36.202731 0.0
-120.102917 36.202731 100.0
-120.102917 36.202731 1000.0
-120.102917 36.202731 10000.0
-120.102917 36.202731 40000.0
-119.782911 36.300783 0.0
-119.782911 36.300783 100.0
-119.782911 36.300783 1000.0
-119.782911 36.300783 10000.0
-119.782911 36.300783 40000.0
</pre>



ucvm_query result from UCVM-C 18.5 alpha

<pre>
../bin/ucvm_query -f ../conf/ucvm.conf -m cs173 << EOF &> cs173-it.log
-118.723500 34.614500 0.0
...
EOF
</pre>

<pre>
WARNING: Could not load model into memory. Reading the model from the
hard disk may result in slow performance.
Using Geo Depth coordinates as default mode.
 -118.7235    34.6145      0.000    713.454    390.000      cs173   2047.660   1089.404   2139.522       none      0.000      0.000      0.000      crust   2047.660   1089.404   2139.522
 -118.7235    34.6145    100.000    713.454    390.000      cs173   2734.936   1518.460   2232.448       none      0.000      0.000      0.000      crust   2734.936   1518.460   2232.448
 -118.7235    34.6145   1000.000    713.454    390.000      cs173   5460.268   3237.855   2605.311       none      0.000      0.000      0.000      crust   5460.268   3237.855   2605.311
 -118.7235    34.6145  10000.000    713.454    390.000      cs173   6561.615   3871.144   2890.995       none      0.000      0.000      0.000      crust   6561.615   3871.144   2890.995
 -118.7235    34.6145  40000.000    713.454    390.000      cs173   7545.088   4324.382   3141.390       none      0.000      0.000      0.000      crust   7545.088   4324.382   3141.390
 -119.6982    34.4208      0.000     15.086    387.000      cs173   4670.389   2585.928   2430.952       none      0.000      0.000      0.000      crust   4670.389   2585.928   2430.952
 -119.6982    34.4208    100.000     15.086    387.000      cs173   4668.188   2583.502   2430.503       none      0.000      0.000      0.000      crust   4668.188   2583.502   2430.503
 -119.6982    34.4208   1000.000     15.086    387.000      cs173   4625.225   2526.679   2420.277       none      0.000      0.000      0.000      crust   4625.225   2526.679   2420.277
 -119.6982    34.4208  10000.000     15.086    387.000      cs173   5915.015   3283.019   2621.902       none      0.000      0.000      0.000      crust   5915.015   3283.019   2621.902
 -119.6982    34.4208  40000.000     15.086    387.000      cs173   7619.421   4327.448   3143.100       none      0.000      0.000      0.000      crust   7619.421   4327.448   3143.100
 -120.5827    34.9681      0.000     18.443    280.000      cs173   4767.372   2397.728   2399.038       none      0.000      0.000      0.000      crust   4767.372   2397.728   2399.038
 -120.5827    34.9681    100.000     18.443    280.000      cs173   4771.665   2403.488   2399.941       none      0.000      0.000      0.000      crust   4771.665   2403.488   2399.941
 -120.5827    34.9681   1000.000     18.443    280.000      cs173   4833.515   2478.506   2412.065       none      0.000      0.000      0.000      crust   4833.515   2478.506   2412.065
 -120.5827    34.9681  10000.000     18.443    280.000      cs173   5818.455   3107.401   2560.882       none      0.000      0.000      0.000      crust   5818.455   3107.401   2560.882
 -120.5827    34.9681  40000.000     18.443    280.000      cs173   7883.487   4559.606   3267.956       none      0.000      0.000      0.000      crust   7883.487   4559.606   3267.956
 -119.1044    35.3444      0.000    114.259    280.000      cs173   3372.060   1940.204   2335.459       none      0.000      0.000      0.000      crust   3372.060   1940.204   2335.459
 -119.1044    35.3444    100.000    114.259    280.000      cs173   3378.782   1943.201   2335.872       none      0.000      0.000      0.000      crust   3378.782   1943.201   2335.872
 -119.1044    35.3444   1000.000    114.259    280.000      cs173   3507.057   1995.574   2343.018       none      0.000      0.000      0.000      crust   3507.057   1995.574   2343.018
 -119.1044    35.3444  10000.000    114.259    280.000      cs173   6903.631   4181.306   3061.306       none      0.000      0.000      0.000      crust   6903.631   4181.306   3061.306
 -119.1044    35.3444  40000.000    114.259    280.000      cs173   7822.489   4575.821   3276.215       none      0.000      0.000      0.000      crust   7822.489   4575.821   3276.215
 -120.9907    35.6688      0.000    602.291    782.000      cs173   4317.442   2337.356   2389.908       none      0.000      0.000      0.000      crust   4317.442   2337.356   2389.908
 -120.9907    35.6688    100.000    602.291    782.000      cs173   4324.663   2342.934   2390.741       none      0.000      0.000      0.000      crust   4324.663   2342.934   2390.741
 -120.9907    35.6688   1000.000    602.291    782.000      cs173   4501.126   2466.457   2410.151       none      0.000      0.000      0.000      crust   4501.126   2466.457   2410.151
 -120.9907    35.6688  10000.000    602.291    782.000      cs173   5709.418   3078.566   2551.708       none      0.000      0.000      0.000      crust   5709.418   3078.566   2551.708
 -120.9907    35.6688  40000.000    602.291    782.000      cs173   7870.905   4547.270   3261.620       none      0.000      0.000      0.000      crust   7870.905   4547.270   3261.620
 -120.4328    35.8997      0.000    467.078    349.000      cs173   4388.868   2370.534   2394.871       none      0.000      0.000      0.000      crust   4388.868   2370.534   2394.871
 -120.4328    35.8997    100.000    467.078    349.000      cs173   4390.334   2371.916   2395.081       none      0.000      0.000      0.000      crust   4390.334   2371.916   2395.081
 -120.4328    35.8997   1000.000    467.078    349.000      cs173   4371.586   2367.767   2394.453       none      0.000      0.000      0.000      crust   4371.586   2367.767   2394.453
 -120.4328    35.8997  10000.000    467.078    349.000      cs173   5795.831   3301.412   2628.905       none      0.000      0.000      0.000      crust   5795.831   3301.412   2628.905
 -120.4328    35.8997  40000.000    467.078    349.000      cs173   7586.191   4296.916   3126.092       none      0.000      0.000      0.000      crust   7586.191   4296.916   3126.092
 -120.1029    36.2027      0.000    113.530    280.000      cs173   2851.014   1495.107   2259.639       none      0.000      0.000      0.000      crust   2851.014   1495.107   2259.639
 -120.1029    36.2027    100.000    113.530    280.000      cs173   2859.429   1498.931   2260.468       none      0.000      0.000      0.000      crust   2859.429   1498.931   2260.468
 -120.1029    36.2027   1000.000    113.530    280.000      cs173   3050.430   1580.739   2277.315       none      0.000      0.000      0.000      crust   3050.430   1580.739   2277.315
 -120.1029    36.2027  10000.000    113.530    280.000      cs173   5612.277   3107.839   2561.259       none      0.000      0.000      0.000      crust   5612.277   3107.839   2561.259
 -120.1029    36.2027  40000.000    113.530    280.000      cs173   7696.131   4393.316   3179.477       none      0.000      0.000      0.000      crust   7696.131   4393.316   3179.477
 -119.7829    36.3008      0.000     69.541    280.000      cs173   3031.683   1489.972   2258.508       none      0.000      0.000      0.000      crust   3031.683   1489.972   2258.508
 -119.7829    36.3008    100.000     69.541    280.000      cs173   3040.815   1494.782   2259.555       none      0.000      0.000      0.000      crust   3040.815   1494.782   2259.555
 -119.7829    36.3008   1000.000     69.541    280.000      cs173   3253.721   1604.038   2281.744       none      0.000      0.000      0.000      crust   3253.721   1604.038   2281.744
 -119.7829    36.3008  10000.000     69.541    280.000      cs173   5852.581   3303.596   2629.943       none      0.000      0.000      0.000      crust   5852.581   3303.596   2629.943
 -119.7829    36.3008  40000.000     69.541    280.000      cs173   7657.497   4359.726   3160.989       none      0.000      0.000      0.000      crust   7657.497   4359.726   3160.989
</pre>

== Another set of Test Points ==
<pre>
-118.72254 34.61444 0
-118.72254 34.61444 175
-118.72254 34.61444 1050
-118.72254 34.61444 9975
-118.72254 34.61444 40075
-119.69846 34.42187 0
-119.69846 34.42187 175
-119.69846 34.42187 1050
-119.69846 34.42187 9975
-119.69846 34.42187 40075
-120.58229 34.96759 0
-120.58229 34.96759 175
-120.58229 34.96759 1050
-120.58229 34.96759 9975
-120.58229 34.96759 40075
-119.10517 35.34494 0
-119.10517 35.34494 175
-119.10517 35.34494 1050
-119.10517 35.34494 9975
-119.10517 35.34494 40075
-120.99114 35.66917 0
-120.99114 35.66917 175
-120.99114 35.66917 1050
-120.99114 35.66917 9975
-120.99114 35.66917 40075
-120.43380 35.89934 0
-120.43380 35.89934 175
-120.43380 35.89934 1050
-120.43380 35.89934 9975
-120.43380 35.89934 40075
-120.50437 35.81076 0
-120.50437 35.81076 175
-120.50437 35.81076 1050
-120.50437 35.81076 9975
-120.50437 35.81076 40075
-119.78270 36.30157 0
-119.78270 36.30157 175
-119.78270 36.30157 1050
-119.78270 36.30157 9975
-119.78270 36.30157 40075
</pre>


ucvm_query result

<pre>
-118.7225    34.6144      0.000    720.976    390.000      cs173   2043.192   1085.995   2138.201       none      0.000      0.000      0.000      crust   2043.192   1085.995   2138.201
 -118.7225    34.6144    175.000    720.976    390.000      cs173   3248.173   1837.993   2301.267       none      0.000      0.000      0.000      crust   3248.173   1837.993   2301.267
 -118.7225    34.6144   1050.000    720.976    390.000      cs173   5459.535   3237.474   2605.180       none      0.000      0.000      0.000      crust   5459.535   3237.474   2605.180
 -118.7225    34.6144   9975.000    720.976    390.000      cs173   6556.426   3867.625   2889.135       none      0.000      0.000      0.000      crust   6556.426   3867.625   2889.135
 -118.7225    34.6144  40075.000    720.976    390.000      cs173   7551.173   4330.115   3144.578       none      0.000      0.000      0.000      crust   7551.173   4330.115   3144.578
 -119.6985    34.4219      0.000     16.241    387.000      cs173   4656.098   2577.254   2429.342       none      0.000      0.000      0.000      crust   4656.098   2577.254   2429.342
 -119.6985    34.4219    175.000     16.241    387.000      cs173   4652.338   2573.052   2428.572       none      0.000      0.000      0.000      crust   4652.338   2573.052   2428.572
 -119.6985    34.4219   1050.000     16.241    387.000      cs173   4609.664   2514.169   2418.086       none      0.000      0.000      0.000      crust   4609.664   2514.169   2418.086
 -119.6985    34.4219   9975.000     16.241    387.000      cs173   5913.173   3281.206   2621.216       none      0.000      0.000      0.000      crust   5913.173   3281.206   2621.216
 -119.6985    34.4219  40075.000     16.241    387.000      cs173   7622.714   4330.343   3144.709       none      0.000      0.000      0.000      crust   7622.714   4330.343   3144.709
 -120.5823    34.9676      0.000     18.570    280.000      cs173   4765.085   2395.066   2398.626       none      0.000      0.000      0.000      crust   4765.085   2395.066   2398.626
 -120.5823    34.9676    175.000     18.570    280.000      cs173   4772.822   2405.407   2400.245       none      0.000      0.000      0.000      crust   4772.822   2405.407   2400.245
 -120.5823    34.9676   1050.000     18.570    280.000      cs173   4837.920   2483.522   2412.900       none      0.000      0.000      0.000      crust   4837.920   2483.522   2412.900
 -120.5823    34.9676   9975.000     18.570    280.000      cs173   5814.321   3105.312   2560.201       none      0.000      0.000      0.000      crust   5814.321   3105.312   2560.201
 -120.5823    34.9676  40075.000     18.570    280.000      cs173   7883.104   4559.139   3267.717       none      0.000      0.000      0.000      crust   7883.104   4559.139   3267.717
 -119.1052    35.3449      0.000    114.266    280.000      cs173   3371.219   1940.716   2335.530       none      0.000      0.000      0.000      crust   3371.219   1940.716   2335.530
 -119.1052    35.3449    175.000    114.266    280.000      cs173   3382.928   1945.926   2336.248       none      0.000      0.000      0.000      crust   3382.928   1945.926   2336.248
-119.1052    35.3449   1050.000    114.266    280.000      cs173   3516.983   2000.467   2343.681       none      0.000      0.000      0.000      crust   3516.983   2000.467   2343.681
 -119.1052    35.3449   9975.000    114.266    280.000      cs173   6906.483   4185.805   3063.830       none      0.000      0.000      0.000      crust   6906.483   4185.805   3063.830
 -119.1052    35.3449  40075.000    114.266    280.000      cs173   7823.070   4576.244   3276.429       none      0.000      0.000      0.000      crust   7823.070   4576.244   3276.429
 -120.9911    35.6692      0.000    612.945    782.000      cs173   4317.940   2337.993   2390.002       none      0.000      0.000      0.000      crust   4317.940   2337.993   2390.002
 -120.9911    35.6692    175.000    612.945    782.000      cs173   4330.583   2347.756   2391.460       none      0.000      0.000      0.000      crust   4330.583   2347.756   2391.460
 -120.9911    35.6692   1050.000    612.945    782.000      cs173   4516.919   2477.496   2411.950       none      0.000      0.000      0.000      crust   4516.919   2477.496   2411.950
 -120.9911    35.6692   9975.000    612.945    782.000      cs173   5711.902   3079.826   2552.101       none      0.000      0.000      0.000      crust   5711.902   3079.826   2552.101
 -120.9911    35.6692  40075.000    612.945    782.000      cs173   7870.165   4546.529   3261.238       none      0.000      0.000      0.000      crust   7870.165   4546.529   3261.238
 -120.4338    35.8993      0.000    468.744    349.000      cs173   4401.550   2376.688   2395.809       none      0.000      0.000      0.000      crust   4401.550   2376.688   2395.809
 -120.4338    35.8993    175.000    468.744    349.000      cs173   4403.934   2378.959   2396.156       none      0.000      0.000      0.000      crust   4403.934   2378.959   2396.156
 -120.4338    35.8993   1050.000    468.744    349.000      cs173   4379.906   2371.956   2395.088       none      0.000      0.000      0.000      crust   4379.906   2371.956   2395.088
 -120.4338    35.8993   9975.000    468.744    349.000      cs173   5789.239   3302.620   2629.347       none      0.000      0.000      0.000      crust   5789.239   3302.620   2629.347
 -120.4338    35.8993  40075.000    468.744    349.000      cs173   7586.570   4297.113   3126.202       none      0.000      0.000      0.000      crust   7586.570   4297.113   3126.202
 -120.5044    35.8108      0.000    553.646    455.000      cs173   5124.013   2859.918   2490.301       none      0.000      0.000      0.000      crust   5124.013   2859.918   2490.301
 -120.5044    35.8108    175.000    553.646    455.000      cs173   5115.474   2852.077   2488.361       none      0.000      0.000      0.000      crust   5115.474   2852.077   2488.361
 -120.5044    35.8108   1050.000    553.646    455.000      cs173   5039.287   2789.464   2473.323       none      0.000      0.000      0.000      crust   5039.287   2789.464   2473.323
 -120.5044    35.8108   9975.000    553.646    455.000      cs173   5692.855   3347.414   2646.662       none      0.000      0.000      0.000      crust   5692.855   3347.414   2646.662
-120.5044    35.8108  40075.000    553.646    455.000      cs173   7612.118   4317.427   3137.526       none      0.000      0.000      0.000      crust   7612.118   4317.427   3137.526
 -119.7827    36.3016      0.000     69.588    280.000      cs173   3032.590   1490.514   2258.628       none      0.000      0.000      0.000      crust   3032.590   1490.514   2258.628
 -119.7827    36.3016    175.000     69.588    280.000      cs173   3048.569   1498.927   2260.458       none      0.000      0.000      0.000      crust   3048.569   1498.927   2260.458
 -119.7827    36.3016   1050.000     69.588    280.000      cs173   3273.004   1613.803   2283.640       none      0.000      0.000      0.000      crust   3273.004   1613.803   2283.640
 -119.7827    36.3016   9975.000     69.588    280.000      cs173   5843.598   3296.038   2626.953       none      0.000      0.000      0.000      crust   5843.598   3296.038   2626.953
 -119.7827    36.3016  40075.000     69.588    280.000      cs173   7657.531   4359.616   3160.929       none      0.000      0.000      0.000      crust   7657.531   4359.616   3160.929
</pre>



== Comparison Plots CS173 ==

=== Raw mesh plots ===

Below are cross-section plots produced by using the CyberShake mesh plotting tools on the raw cs173 velocity mesh files.  The color scale is basemap.cm.GMT_seis, which was originally defined in GMT and then ported to basemap.

==== Horizontal cross-sections ====

The surface cross-section includes lines indicating where the vertical cross-section plots are.

{|
| [[File:cs173_horiz_slice_0z.png|thumb|400px|z=0, depth=0 km, with vertical cross-section locations indicated]]
| [[File:cs173_horiz_slice_6z.png|thumb|400px|z=6, depth=1.05 km]]
| [[File:cs173_horiz_slice_57z.png|thumb|400px|z=57, depth=9.975 km]]
| [[File:cs173_horiz_slice_114z.png|thumb|400px|z=114, depth=19.95 km]]
|}

==== Vertical cross-sections ====

Parallel to Y-axis

X=1000: (32.550700, -115.409380) to (38.446980, -126.695380)

X=2000: (33.824270, -114.301790) to (39.813370, -125.716460)

{|
| [[File:cs173_vert_slice_1000x.png|thumb|400px|X=1000 (blue line on horizontal plot)]]
| [[File:cs173_vert_slice_2000x.png|thumb|400px|X=2000 (orange line on horizontal plot)]]
|}

Parallel to X-axis

Y=1740: (32.824810, -119.141480) to (36.959180, -115.655030)

Y=3480: (34.320700, -121.888010) to (38.531550, -118.487980)

Y=5220: (35.744050, -124.726900) to (40.029510, -121.442660)

{|
| [[File:cs173_vert_slice_1740y.png|thumb|400px|Y=1740 (green line on horizontal plot)]]
| [[File:cs173_vert_slice_3480y.png|thumb|400px|Y=3480 (red line on horizontal plot)]]
| [[File:cs173_vert_slice_5220y.png|thumb|400px|Y=5220 (purple line on horizontal plot)]]
|}

=== UCVM plots ===

{|
| [[File:scott_cs173_hort_0.png|thumb|400px|z=0, depth=0 km]]
| [[File:scott_cs173_hort_1050.png|thumb|400px|z=6, depth=1.05 km]]
|}
{|
| [[File:scott_cs173_hort_9975.png|thumb|400px|z=57, depth=9.975 km]]
| [[File:scott_cs173_hort_19950.png|thumb|400px|z=114, depth=19.95 km]]
|}
{|
| [[File:scott_cs173_hort_50400.png|thumb|400px|z=287, depth=50.4 km]]
|}

Parallel to Y-axis

{|
| [[File:scott_cs173_cross_x1.png|thumb|400px|X=1000]]
| [[File:scott_cs173_cross_x2.png|thumb|400px|X=2000]]
|}

Parallel to X-axis

{|
| [[File:scott_cs173_cross_y1.png|thumb|400px|Y=1740]]
| [[File:scott_cs173_cross_y2.png|thumb|400px|Y=3480]]
|}
{|
| [[File:scott_cs173_cross_y3.png|thumb|400px|Y=5220]]
|}

== Example Comparison Data Points ==
* CS173
* UCVM CS173

Vs at top 6 layers(counting down at steps of 175m) at lon/lats,(36.02718,-119.446), that translates to grid point (1600, 3000) on the mesh

    x 1600 y 3000 z 287  >>  vs 1776.336060
    x 1600 y 3000 z 286  >>  vs 1789.887817
    x 1600 y 3000 z 285  >>  vs 1803.439453
    x 1600 y 3000 z 284  >>  vs 1818.368774
    x 1600 y 3000 z 283  >>  vs 1841.563721
    x 1600 y 3000 z 282  >>  vs 1864.758545

== Velocity mesh region ==

The region we will extract is 560 km x 1218 km x 50.4 km (as measured in the WGS84 coordinate system), with 55 degree rotation counter-clockwise.  The mesh will have 175m grid point spacing.

== Mesh extraction process ==

To deliver the mesh, we will follow the following procedure:

#CyberShake processing tools will be used to construct a mesh, with the location of each point defined as km offsets from the center of the mesh at the surface.
## The mesh will evenly discretize the region in x,y,z, but due to the rotation, this is not North,South. The resolution will be 175m as measured in the UTM coordinate system.
##A number of CyberShake codes, originally from RWG, are involved in this step. These codes use his great circle projection.  These codes figure out the outer bounds of the mesh given the location of a site and fault surfaces. They return a lat/lon of the mesh center. 
## The codes used are in version control here: [http://source.usc.edu/svn/cybershake/import/trunk/PreCVM].
##Then we use the Proj 4 library to convert this center point to UTM Zone 11, which is the southern California UTM zone. While some of the mesh region extends north into UTM zone 10, we will use Zone 11 since the majority of our region of interest is in that zone.  We'll do the remaining next steps in UTM.  
#The location of each mesh point will be redefined as a UTM easting and northing, by translating the center of the mesh from lat/lon to UTM coordinates in Zone 11 and then applying the offsets.
##The UTM mesh points will be written to a file, in fast y, x, z format, starting with all the mesh points are the surface, first changing y values, then x values, then moving down one the Z axis, and listing the next layer of points.  These points will be in the same order as the points in the velocity mesh.
##This file in internal, and will not be part of the data sent to the USR group.
#The pyproj package will be used to convert the UTM coordinates at each point into a latitude and longitude value, to 5 decimal places.
##The proj4 coordinate system to coordinate system (cs2cs) conversion string used will be proj='utm', zone=11, ellps='WGS84'.
#This list of points will be used to query UCVM v15.10 through the C API, using CyberShake's [[CyberShake_Code_Base#UCVM|UCVM wrapper]], same as used in CyberShake Study 17.3, and create a velocity mesh.  For a given point, we will first query the CCA-06 model, then CVM-S4.26, then USGS Bay Area, and finally a background 1D model.  Additionally, we will query the USGS topo model in UCVM 17.3 to obtain the surface elevation at each point and write it to a file.
#Smoothing will be applied to the velocity mesh for all points within 10 km of an interface, by averaging the neighbors for 10km in the N, S, E, and W directions.  We will use the [[CyberShake_Code_Base#Smoothing|CyberShake smoothing code]].  Sample plots illustrating the effect of this smoothing are available [[Study_16.9_Velocity_Plots|here]].
#Scott will deliver three meshes (one each with Vp, Vs, and rho) in fast y, x, z binary format with 4-byte floats, a list of the elevation of each surface point, a list of the coordinates of each surface point in lat/lon and UTM zone 11, and MD5 sums for all products.

== Mesh region details ==

The actual mesh we extracted has its center at (36.44287, -120.22745), or 210707 easting, 4037914 northing in UTM Zone 11.

The actual corner points are:
 S: (31.26643, -116.48562)
 W: (37.08416, -127.65648)
 E: (35.33518, -112.92896)
 N: (41.45284, -124.51032)

Using proj4 project, we project these points into UTM Zone 10 (since most of the mesh is in northern California:

Lat/Lon WGS84 to UTM Zone 10
 S: (31.26643, -116.48562) : 1120786.790041  3477495.533344
 W: (37.08416, -127.65648) : 86016.577605    4114365.092550
 E: (35.33518, -112.92896)  : 1416797.228564  3957099.088864
 N: (41.45284, -124.51032) : 373851.268121 4590129.392433


Lat/Lon WGS84 to UTM Zone 11
 S: (31.26643, -116.48562) : 548969.079292 3459243.769232
 W: (37.08416, -127.65648) : -448610.671359 4157755.975881
 E: (35.33518, -112.92896) : 870071.761707 3917825.705794
 N: (41.45284, -124.51032) :  -127508.133290 4616338.305196

== Content of the config file ==

  # UTM Zone
  utm_zone = 11
  #
  # Model directory
  model_dir = cs173
  #
  # Number of cells in x, y, and z.
  nx = 3200
  ny = 6960
  nz = 288
  #
  # Depth, in meters
  depth = 50400
  depth_interval = 175
  # 
  # Four corners
  bottom_left_corner_e = 548969.079292
  bottom_left_corner_n = 3459243.769232
  top_left_corner_e    = -448610.671359
  top_left_corner_n    = 4157755.975881
  bottom_right_corner_e = 870071.761707
  bottom_right_corner_n = 3917825.705794
  top_right_corner_e   = -127508.133290
  top_right_corner_n   = 4616338.305196
  #
  # Seek method, fast-X or fast-Y
  seek_axis = fast-Y
  seek_direction = top-down

== Mesh integration process ==
# Returned mesh will be registered into UCVM
#We will develop evaluation plan to determine whether returned points are registered properly.

== Model Coverage Regions ==
[[File:Velocity_Mesh_with_models.png|thumb|left|600px|The mesh is depicted in white; the region covered by CCA-06 is in blue, by CVM-S4.26 in red, and by USGS Bay Area in green.]]

== Related Entries ==

*[[USR]]
*[[CVM-S]]
*[[CVM-H]]
*[[UCVM]]
*[[Proj4]]