Adding GTL - SCECcode/ucvm GitHub Wiki
To add the external Ely-Jordan GTL method in ucvm_query, the user would supply elygtl:ely in the model list and may include the optional z-range for the interpolation range. If not supplied, a default of 0 to 350 is assumed.
./ucvm_query -f ../conf/ucvm.conf -m cvms5,elygtl:ely -z 0,350
with input data (cvms5 with use_gtl=off)
-118.0000 34.0000 0.000 -118.0000 34.0000 50.000 -118.0000 34.0000 100.000 -118.0000 34.0000 250.000 -118.0000 34.0000 300.000 -118.0000 34.0000 350.000 -118.0000 34.0000 400.000
result:
-118.0000 34.0000 0.000 281.668 468.400 cvms5 2330.334 1187.775 2175.202 elygtl 0.000 468.400 0.000 ely 884.192 234.200 1143.546 -118.0000 34.0000 50.000 281.668 468.400 cvms5 2330.334 1187.775 2175.202 elygtl 0.000 468.400 0.000 ely 2203.489 711.726 1990.240 -118.0000 34.0000 100.000 281.668 468.400 cvms5 2330.334 1187.775 2175.202 elygtl 0.000 468.400 0.000 ely 2393.050 874.489 2058.621 -118.0000 34.0000 250.000 281.668 468.400 cvms5 2330.334 1187.775 2175.202 elygtl 0.000 468.400 0.000 ely 2387.138 1117.435 2056.630 -118.0000 34.0000 300.000 281.668 468.400 cvms5 2330.334 1187.775 2175.202 elygtl 0.000 468.400 0.000 ely 2353.512 1158.983 2045.136 -118.0000 34.0000 350.000 281.668 468.400 cvms5 2330.334 1187.775 2175.202 none 0.000 0.000 0.000 crust 2330.334 1187.775 2175.202 -118.0000 34.0000 400.000 281.668 468.400 cvms5 2397.913 1227.419 2188.256 none 0.000 0.000 0.000 crust 2397.913 1227.419 2188.256
./plot_depth_profile.py -s 34,-118 -b 0 -e 1000 -d vs -v 50 -c cvms5 -z 0,350 -o cvms5_noGTL.png ./plot_depth_profile.py -s 34,-118 -b 0 -e 1000 -d vs -v 50 -c cvms5,elygtl:ely -z 0,350 -o cvms5_noGTL_elygtl_full.png ./plot_depth_profile.py -s 34,-118 -b 0 -e 1000 -d vs -v 50 -c cvms5,elygtl:ely -z 200,350 -o cvms5_noGTL_elygtl_part.png
Baseline CVMS5 with no internal or external GTLs,
With external Elgtl GTL and z-range from 0 to 350,
With external Elygtl GTL and z-range from 200 to 350,
Using internal GTL in cvms5 (use_gtl=on),
./plot_depth_profile.py -s 34,-118 -b 0 -e 1000 -d vs -v 50 -c cvms5 -z 0,350 -o cvms5_GTL.png
Baseline with no gtl and no Elygtl,
./plot_cross_section.py -b 34.5,-116 -u 34.5,-113 -h 1000 -v 10 -d vs -c cvms5 -z 0,3 50 -a d -s 0 -e 1000 -o cross_cvms5_nogtl.png
Enable Elygtl with range of 0-350,
./plot_cross_section.py -b 34.5,-116 -u 34.5,-113 -h 1000 -v 10 -d vs -c cvms5,elygtl :ely -z 0,350 -a d -s 0 -e 1000 -o cross_cvms5_nogtl_elygtl.png
Enable Elygtl and use Z-range of 0-800,
./plot_cross_section.py -b 34.5,-116 -u 34.5,-113 -h 1000 -v 10 -d vs -c cvms5,elygtl :ely -z 0,800 -a d -s 0 -e 1000 -o cross_cvms5_nogtl_elygtl_2.png
Turn on GTL in model/cvms5/data/config with no Elygtl,
./plot_cross_section.py -b 34.5,-116 -u 34.5,-113 -h 1000 -v 10 -d vs -c cvms5 -z 0,3 50 -a d -s 0 -e 1000 -o cross_cvms5_nogtl.png
Turn on GTL in model/cvms5/data/config with Elygtl and Z-range of 0-350,
./plot_cross_section.py -b 34.5,-116 -u 34.5,-113 -h 1000 -v 10 -d vs -c cvms5,elygtl :ely -z 0,350 -a d -s 0 -e 1000 -o cross_cvms5_gtl_elygtl.png
Turn on GTL in model/cvms5/data/config with Elygtl and Z-range of 0-800,
./plot_cross_section.py -b 34.5,-116 -u 34.5,-113 -h 1000 -v 10 -d vs -c cvms5,elygtl :ely -z 0,800 -a d -s 0 -e 1000 -o cross_cvms5_gtl_elygtl.png
