Difference between revisions of "UCVM sfcvm geomodelgrid"

From SCECpedia
Jump to navigationJump to search
Line 244: Line 244:
  
 
{|
 
{|
−
| [[FILE:sfcvm_gabbro_yes_map.png|thumb|300px | sfcvm gabbro map]]
+
| [[FILE:sfcvm_gabbro_yes_map.png|thumb|500px | sfcvm gabbro map]]
−
| [[FILE:sfcvm_gabbro_yes.png|thumb|300px|sfcvm gabbro 0m]]
+
| [[FILE:sfcvm_gabbro_yes.png|thumb|400px|sfcvm gabbro 0m]]
 
|}
 
|}
  
Line 253: Line 253:
  
 
{|
 
{|
−
| [[FILE:cencal-o-map.png|thumb|300px|cencal map]]
+
| [[FILE:cencal-o-map.png|thumb|500px|cencal map]]
−
| [[FILE:cencal-o.png|thumb|300px|cencal plot]]
+
| [[FILE:cencal-o.png|thumb|400px|cencal plot]]
 
|}
 
|}

Revision as of 19:40, 27 February 2024

UCVM with SFCVM

Software stack for supporting geomodelgrid based format:

  1. curl (curl-7.87.0)
  2. hdf5 (hdf5-1.10.9)
  3. openssl (openssl-3.0.8)
  4. proj (proj-8.1.1)
  5. sqlite (sqlite-autoconf-3400100)
  6. tiff (tiff-4.5.0)

SFCVM in UCVM

geomodelgrid based SFCVM is set with 

   
    geomodelgrids_squery_setSquashing(sfcvm_geo_query_object, GEOMODELGRIDS_SQUASH_TOPOGRAPHY_BATHYMETRY);
    geomodelgrids_squery_setSquashMinElev(sfcvm_geo_query_object, SFCVM_SquashMinElev);

where SFCVM_SquashMinElev is default to -5000m (sfcvm is elevation based, depth is negative value)

NOTE: For locations that are below water, their geomodelgrids_squery_queryTopoBathyElevation() could return a NA value and so a grid height is added across the board to always have an usable value.

Plots from web viewer

sfcvm regions
sfcvm cross section
sfcvm depth 0m
sfcvm depth 100m
sfcvm depth 1000m
sfcvm depth 5000m
sfcvm depth 6000m

New command line option for model specific configuration

UCVM option: 

-P Optional model configuration(model_param:key,value)


Example: 

ucvm_query -l 37.455000,-121.941,3000 -m sfcvm -f ../../conf/ucvm.conf -P sfcvm_param:SquashMinElev,-5000.0


Allows multiple -P options, 

ucvm_query ... -P cvmh_param:USE_1D_BKG,False -P cvmh_param:USE_GTL,False


More detailed plots from SFCVM

Internal calls from web viewer, (set e to 0,100,1000,5000,6000) 

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_horizontal_slice.py  -d vs -c sfcvm -s 0.010 -a sd -o ../result/UCVM_15167424_h.png -n ../model/UCVM_TARGET/conf/ucvm.conf -i ../model/UCVM_TARGET  -b 36.5416,-123.8269 -u 38.8108,-120.6299 -e 0
sfcvm depth 0m
sfcvm depth 100m
sfcvm depth 1000m
sfcvm depth 5000m
sfcvm depth 6000m

Detail plots from Cencal

cencal depth 0m
cencal depth 100m
cencal depth 1000m
cencal depth 5000m
cencal depth 6000m
cencal cross section


High Velocity surface regions in SFCVM

Horizontal region,

horizontal region
plot_horizontal_slice.py  -d vs -c sfcvm,elygtl:taper -s 0.0005 -a sd -o ../result/UCVM_18933492_h.png -n ../model/UCVM_TARGET/conf/ucvm.conf -i ../model/UCVM_TARGET  -L 500,1700,1700 -z 0,700 -b 37.6021,-122.1711 -u 37.7922,-122.0021 -e 500 -A 0.0,4.0

with sfcvm (only detailed region)

at depth 0m
at depth 500m
at depth 1000m
at depth 5000m

with sfcvm,elygtl:taper (only detailed region)

at depth 0m with taper
at depth 500m with taper
at depth 1000m with taper

with sfcvm,elygtl:taper -A 0,2 (only detailed region)

at depth 0m with taper
at depth 50m with taper
at depth 100m with taper

SFCVM(detailed + regional regions) with Gabbro modification

Upper San Leandro Gabbro region Map,

regions

Rectangle horizontal box,

no Gabbro
Gabbro

Lower Cross section,

no Gabbro
Gabbro

Diagonal Cross section,

no Gabbro
Gabbro

Depth Profile at the blue marker,

no Gabbro
Gabbro


Comparing depth profile plot (red=no gabbro, blue=gabbro)

vs comparison
zoomed in vs comparison
vp comparison


with Gabbro modification (full detailed +regional regions)

at depth 0m with Gabbro
at depth 500m with Gabbro
at depth 1000m with Gabbro
at depth 5000m with Gabbro
at depth 7500m with Gabbro

Vertical Slice

vertical section

with sfcvm (detailed region) 

plot_cross_section.py -s 0 -h 48 -d vs -c sfcvm,elygtl:taper -a sd -o ../result/UCVM_23820778_c.png -n ../model/UCVM_TARGET/conf/ucvm.conf -i ../model/UCVM_TARGET -v 25 -L 500,1700,1700 -z 0,700 -b 37.7026,-122.1728 -u 37.701,-121.9984 -e 5000 -A 0.0,4.0
vertical slice upto 1k
vertical slice upto 50k

with sfcvm,elygtl:taper (detailed region)

vertical slice upto 1k
vertical slice upto 50k

More Study

SFCVM (detailed+regional)

sfcvm coverage


With Water-Step1-Step2 modifications 

 surf = geomodelgrids_squery_queryTopoBathyElevation
 surf < 0:
    call geomodelgrids_squery_query with elevation= N
    if vs<0 || vp<0
       call geomodelgrids_squery_query with elevation= N-25 (first step down grid height)
       if vs < 0 || vp < 0
          call geomodelgrids_squery_query with elevation= N-125 (second step down grid height)

PROJ_LIB=../model/UCVM_TARGET/lib/proj/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_horizontal_slice.py  -d vs -c sfcvm
 -s 0.010 -a sd -o ../result/UCVM_93552953_h.png -n ../model/UCVM_TARGET/conf/ucvm.con
f -i ../model/UCVM_TARGET  -b 34.5541,-126.9141 -u 41.4695,-117.9492 -e 0 -A 0,4
sfcvm gabbro map
sfcvm gabbro 0m


As baseline, original CENCAL

cencal map
cencal plot
Retrieved from "https://strike.scec.org/scecwiki/index.php?title=UCVM_sfcvm_geomodelgrid&oldid=28301"