Adding GTL

Elygtl

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

Plotting depth profile with varying parameters

./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

Plotting cross sections

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