Difference between revisions of "UCVM cvmsi tapering"

From SCECpedia
Jump to navigationJump to search
 
(11 intermediate revisions by the same user not shown)
Line 1: Line 1:
−
== UCVM tapering ==
+
== UCVM tapering using ucvm_query ==
  
−
ucvm_query flags that affects tapering:
+
Flags that affects tapering:
  
 
Adding taper interpolation to a model.
 
Adding taper interpolation to a model.
Line 11: Line 11:
 
          
 
          
 
   -L vs,vp,density   
 
   -L vs,vp,density   
−
     example, -L 200,-1,-1
+
     example, resetting vs floor to 200
−
          (only floor the vs to 200)
+
              -L 200,1700,1700
 +
            or -L 200,-1,-1 (needs to pull latest GitHub commit)
  
−
Z interpolation range. default: zmin =0, zmax = 350
+
Z interpolation range. default: zmin =0m, zmax = 350m
  
 
   -z zmin,zmax   
 
   -z zmin,zmax   
 
     example, -z 0,200
 
     example, -z 0,200
 +
 +
== UCVM tapering using C-api ==
 +
 +
int ucvm_add_model_list(const char *mlist)
 +
  ucvm_add_modellist("cvmsi,elygtl:taper")
 +
 +
int ucvm_setfloor(double *fvals)
 +
  ucvm_setfloor({500.0, 1700.0, 1700.0})
 +
 +
int ucvm_setparam(ucvm_model_param_t param, ...)
 +
  ucvm_setparam(UCVM_MODEL_PARAM_IFUNC_ZRANGE, 0.0, 150.0)
 +
 +
== Using ucvm_query calls ==
 +
 +
 
 +
===> ucvm_query -l 34.0458,-114.4995,0 -f /var/www/html/UCVM_web_sfcvm_mei/web/model/UCVM_TARGET/conf/ucvm.conf -m cvmsi
 +
Using Geo Depth coordinates as default mode.
 +
-114.4995    34.0458      0.000    236.002    710.100      cvmsi  5238.240  3073.977  2739.354      none      0.000      0.000      0.000      crust  5238.240  3073.977  2739.354
 +
 
 +
 
 +
===> ucvm_query -l 34.0458,-114.4995,0 -f /var/www/html/UCVM_web_sfcvm_mei/web/model/UCVM_TARGET/conf/ucvm.conf -m cvmsi,elygtl:taper -L 0,1700,1700
 +
Using Geo Depth coordinates as default mode.
 +
-114.4995    34.0458      0.000    236.002    710.100      cvmsi  5283.988  3108.776  2739.354    elygtl      0.000    710.100      0.000      taper  1700.000    355.050  1700.000
 +
 
 +
 
 +
===> ucvm_query -l 34.0458,-114.4995,0 -f /var/www/html/UCVM_web_sfcvm_mei/web/model/UCVM_TARGET/conf/ucvm.conf -m cvmsi,elygtl:taper -L 500,1700,1700
 +
Using Geo Depth coordinates as default mode.
 +
-114.4995    34.0458      0.000    236.002    710.100      cvmsi  5283.988  3108.776  2739.354    elygtl      0.000    710.100      0.000      taper  1700.000    500.000  1700.000
  
 
== Plots from web viewer ==
 
== Plots from web viewer ==
  
 
{|
 
{|
−
| [[FILE:cvmsi_region.png|thumb|300px|cvmsi region]]  
+
| [[FILE:taper_cvmsi_region.png|thumb|300px|cvmsi region]]  
 
|}
 
|}
  
Line 28: Line 57:
  
 
{|
 
{|
−
| [[FILE:cvmsi_vs_base.png|thumb|300px|cvmsi vs]]  
+
| [[FILE:taper_cvmsi_h_none.png|thumb|300px|cvmsi vs none]]  
−
| [[FILE:cvmsi_taper_vs_500.png|thumb|300px|cvmsi taper vs floor 500m]]
+
| [[FILE:taper_cvmsi_h_500.png|thumb|300px|cvmsi taper vs floor 500m]]  
−
| [[FILE:cvmsi_taper_vs_0.png|thumb|300px|cvmsi taper vs floor 0m]]  
 
 
|}
 
|}
 +
  
 
{|
 
{|
−
| [[FILE:cvmsi_taper_diff_0_500.png|thumb|300px|cvmsi taper diff plot]]   
+
| [[FILE:taper_cvmsi_diff_none_500.png|thumb|300px|cvmsi taper diff plot]]   
 
|}
 
|}
  
  
 
<pre>
 
<pre>
−
PROJ_LIB=/usr/local/share/anaconda2/share/proj PATH=/usr/local/share/anaconda2/bin:/usr/local/share/anaconda2/condabin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin PYTHONPATH=../model/UCVM_TARGET/ucvm_plotting plot_cross_section.py -s 0 -h 200 -d vs -c cvmsi -a sd -o ../result/UCVM_88246892_c.png -n ../model/UCVM_TARGET/conf/ucv
+
PROJ_LIB=/usr/local/share/anaconda2/share/proj PATH=/usr/local/share/anaconda2/bin:/usr/local/share/anaconda2/condab
−
m.conf -i ../model/UCVM_TARGET -v 10 -b 33.9548,-118.949 -u 34.391,-117.6306 -e 2000
+
in:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin PYTHONPATH=../model/UCVM_TARGET/ucvm_plotting plot_horizontal_s
 +
lice.py -d vs -c cvmsi -s 0.008 -a sd -o ../result/UCVM_90507555_h.png -n ../model/UCVM_TARGET/conf/uc
 +
vm.conf -i ../model/UCVM_TARGET -b 32.8542,-116.0376 -u 34.6807,-113.8623 -A 1.0,4.0
  
−
PROJ_LIB=/usr/local/share/anaconda2/share/proj PATH=/usr/local/share/anaconda2/bin:/usr/local/share/anaconda2/condabin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin PYTHONPATH=../model/UCVM_TARGET/ucvm_plotting plot_cross_section.py -s 0 -h 200 -d vs -c cvmsi,elygtl:taper -a sd -o ../result/UCVM_88246892_c.png -n ../model/UCVM_TARGET/conf/ucvm.conf -i ../model/UCVM_TARGET -v 10 -b 33.9548,-118.949 -u 34.391,-117.6306 -e 2000 -L 500,1700,1700
 
  
−
ROJ_LIB=/usr/local/share/anaconda2/share/proj PATH=/usr/local/share/anaconda2/bin:/usr/local/share/anaconda2/condabin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin PYTHONPATH=../model/UCVM_TARGET/ucvm_plotting plot_cross_section.py -s 0 -h 200 -d vs -c cvmsi,elygtl:taper -a sd -o ../result/UCVM_88246892_c.png -n ../model/UCVM_TARGET/conf/ucvm.conf -i ../model/UCVM_TARGET -v 10 -b 33.9548,-118.949 -u 34.391,-117.6306 -e 2000 -L 0,1700,1700
+
PROJ_LIB=/usr/local/share/anaconda2/share/proj PATH=/usr/local/share/anaconda2/bin:/usr/local/share/anaconda2/condab
 +
in:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin PYTHONPATH=../model/UCVM_TARGET/ucvm_plotting plot_horizontal_s
 +
lice.py -d vs -c cvmsi,elygtl:taper -s 0.008 -a sd -o ../result/UCVM_90507555_h.png -n ../model/UCVM_TARGET/conf/uc
 +
vm.conf -i ../model/UCVM_TARGET -b 32.8542,-116.0376 -u 34.6807,-113.8623 -e 500 -L 500,1700,1700 -A 1.0,4.0
 
</pre>
 
</pre>
  

Latest revision as of 23:44, 31 January 2024

UCVM tapering using ucvm_query

Flags that affects tapering:

Adding taper interpolation to a model. 

 -m MODEL,elygtl:taper
    example: -m cvmsi,elygtl:taper

Taper Interpolation floors. default: vs=500m, vp=1700m, density=1700m, 

 -L vs,vp,density  
    example, resetting vs floor to 200
              -L 200,1700,1700
           or -L 200,-1,-1 (needs to pull latest GitHub commit)

Z interpolation range. default: zmin =0m, zmax = 350m 

 -z zmin,zmax  
    example, -z 0,200

UCVM tapering using C-api

int ucvm_add_model_list(const char *mlist) 

 ucvm_add_modellist("cvmsi,elygtl:taper")

int ucvm_setfloor(double *fvals) 

  ucvm_setfloor({500.0, 1700.0, 1700.0})

int ucvm_setparam(ucvm_model_param_t param, ...) 

  ucvm_setparam(UCVM_MODEL_PARAM_IFUNC_ZRANGE, 0.0, 150.0)

Using ucvm_query calls

===> ucvm_query -l 34.0458,-114.4995,0 -f /var/www/html/UCVM_web_sfcvm_mei/web/model/UCVM_TARGET/conf/ucvm.conf -m cvmsi
Using Geo Depth coordinates as default mode.
-114.4995    34.0458      0.000    236.002    710.100      cvmsi   5238.240   3073.977   2739.354       none      0.000      0.000      0.000      crust   5238.240   3073.977   2739.354
 
 
===> ucvm_query -l 34.0458,-114.4995,0 -f /var/www/html/UCVM_web_sfcvm_mei/web/model/UCVM_TARGET/conf/ucvm.conf -m cvmsi,elygtl:taper -L 0,1700,1700
Using Geo Depth coordinates as default mode.
-114.4995    34.0458      0.000    236.002    710.100      cvmsi   5283.988   3108.776   2739.354     elygtl      0.000    710.100      0.000      taper   1700.000    355.050   1700.000
 
 
===> ucvm_query -l 34.0458,-114.4995,0 -f /var/www/html/UCVM_web_sfcvm_mei/web/model/UCVM_TARGET/conf/ucvm.conf -m cvmsi,elygtl:taper -L 500,1700,1700
Using Geo Depth coordinates as default mode.
-114.4995    34.0458      0.000    236.002    710.100      cvmsi   5283.988   3108.776   2739.354     elygtl      0.000    710.100      0.000      taper   1700.000    500.000   1700.000

Plots from web viewer

cvmsi region


cvmsi vs none
cvmsi taper vs floor 500m


cvmsi taper diff plot


PROJ_LIB=/usr/local/share/anaconda2/share/proj PATH=/usr/local/share/anaconda2/bin:/usr/local/share/anaconda2/condab
in:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin PYTHONPATH=../model/UCVM_TARGET/ucvm_plotting plot_horizontal_s
lice.py  -d vs -c cvmsi -s 0.008 -a sd -o ../result/UCVM_90507555_h.png -n ../model/UCVM_TARGET/conf/uc
vm.conf -i ../model/UCVM_TARGET  -b 32.8542,-116.0376 -u 34.6807,-113.8623 -A 1.0,4.0


PROJ_LIB=/usr/local/share/anaconda2/share/proj PATH=/usr/local/share/anaconda2/bin:/usr/local/share/anaconda2/condab
in:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin PYTHONPATH=../model/UCVM_TARGET/ucvm_plotting plot_horizontal_s
lice.py  -d vs -c cvmsi,elygtl:taper -s 0.008 -a sd -o ../result/UCVM_90507555_h.png -n ../model/UCVM_TARGET/conf/uc
vm.conf -i ../model/UCVM_TARGET  -b 32.8542,-116.0376 -u 34.6807,-113.8623 -e 500 -L 500,1700,1700 -A 1.0,4.0

See Also

Retrieved from "https://strike.scec.org/scecwiki/index.php?title=UCVM_cvmsi_tapering&oldid=28234"