UCVM sfcvm geomodelgrid

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

at depth 0m

at depth 500m

at depth 1000m

at depth 5000m

with sfcvm,elygtl:taper,

at depth 0m with taper

at depth 500m with taper

at depth 1000m with taper

??? with sfcvm,elygtl:taper

at depth 50m with taper

at depth 100m with taper

Vertical Slice

vertical section

with sfcvm

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,

vertical slice upto 1k

vertical slice upto 50k