UCVM cvmsi tapering

From SCECpedia
Revision as of 23:44, 31 January 2024 by Mei (talk | contribs) (→‎UCVM tapering using ucvm_query)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigationJump to search

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