CTM
We are using jupyter notebook developed by Liz Hearn to calculate CTM values at points defined for the GFM.
Contents
GFM Points
AP has created a discretized version of the GFM. The coverage region for the GFM is similar to the CVM-H GTL coverage region. Corner points from the GFM grid are listed below.
The resolution of the GFM 3D voxel data set is 10km x 10km x 1km (depth).
Each level of the GFM has 8256 points, and the depth goes from -100km to +4km (8256 x 105 = 866880) points.
Corner Points from GFM file:
- -122.681200, 29.709900 -123.205500, 37.335600 -112.481200, 37.412300 -112.864200, 29.767300
A white-space delimited file that defines the GFM region types for each of these points is posted here. This is the updated file that include elevation for each point.
- GFM Cell Points and Data (Warning 13MB zip file download. 65Mb expanded as text file)
The cell location in this file use the CVM-H convention of elevation relative to sea level in meter, with positive values above sea level, and negative values depth below sea level. The cells extend from -100km to 4km.
Example GFM File Data Format
This shows the format of the GFM cell data, before it is updated with CTM heatRegionID and CTM temps. The elevation data is from the GFM topography surface. The last 9 lines are CTM test points, with known results, used to verify that the script used to populate the heatRegionID and Temp are producing the correct results for the GFM points.
long. lat. Z regionID temp_median heatRegionID elevation -122.681222 29.709863 -100000.000000 1.000000 1316.400024 1.000000 -4253.000000 -122.578339 29.714251 -100000.000000 1.000000 1316.400024 1.000000 -4312.000000 -122.475433 29.718560 -100000.000000 1.000000 1316.400024 1.000000 -4310.586426 -122.372505 29.722790 -100000.000000 1.000000 1316.400024 1.000000 -4233.651367 -122.269555 29.726940 -100000.000000 1.000000 1316.400024 1.000000 -4059.687012 -122.166584 29.731010 -100000.000000 1.000000 1316.400024 1.000000 -4157.000000 -122.063592 29.735001 -100000.000000 1.000000 1316.400024 1.000000 -4212.185059 -120.00 35.00 -4532.0 1.000000 1316.400024 1.000000 468.40 -119.00 35.00 -4556.0 1.000000 1316.400024 1.000000 444.00 -118.00 35.00 -4290.0 1.000000 1316.400024 1.000000 710.00 -120.00 34.00 -5820.0 1.000000 1316.400024 1.000000 180.00 -119.00 34.00 -5820.0 1.000000 1316.400024 1.000000 180.00 -118.00 34.00 -5532.0 1.000000 1316.400024 1.000000 468.00 -120.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00 -119.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00 -118.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00
CTM Test Points in GFM format
The CTM notebook defines a series of test points. The query points are given as depth below surface. We've converted these test points into the GFM file format, that expresses them as depth below earth surface. We used the elevation information to express the test points in GFM format. We'll use this GFM test file, to confirm our CTM query script is returning the same CTM results, for both smoothed and unsmoothed modes, as the CTm notebook.
CTM QueryPoint.txt format
-120.00 35.00 5.0 -119.00 35.00 5.0 -118.00 35.00 5.0 -120.00 34.00 6.0 -119.00 34.00 6.0 -118.00 34.00 6.0 -120.00 33.00 8.0 -119.00 33.00 8.0 -118.00 33.00 8.0
Same points in GFM cell file format
long. lat. Z regionID temp_median heatRegionID elevation -120.00 35.00 -4532.0 1.000000 1316.400024 1.000000 468.40 -119.00 35.00 -4556.0 1.000000 1316.400024 1.000000 444.00 -118.00 35.00 -4290.0 1.000000 1316.400024 1.000000 710.00 -120.00 34.00 -5820.0 1.000000 1316.400024 1.000000 180.00 -119.00 34.00 -5820.0 1.000000 1316.400024 1.000000 180.00 -118.00 34.00 -5532.0 1.000000 1316.400024 1.000000 468.00 -120.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00 -119.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00 -118.00 33.00 -7820.0 1.000000 1316.400024 1.000000 180.00
KML file showing model coverage regions
Heatflow Region IDs and Names
ID's for the 14 original CTM heat flow regions, from Thatcher and Chapman (2020): 1 WBR (Western Basin and Range) 2 SN (Sierra Nevada) 3 CCR (Central Coast Ranges) 4 PR (Peninsular Ranges) 5 ST (Outer Salton Trough, aka oST) 6 WTR (Western Transverse Ranges) 7 LA (Los Angeles) 8 SG (San Gabriels) 9 MD (Mojave Desert) 10 GV (Great Valley) 11 ETR (Eastern Transverse Ranges) 12 iST (Inner Salton Trough, also called STi5 in some files) 13 iCB (Inner Continental Borderland) 14 VB (Ventura Basin) (the following were added to expand CTM region to cover the Geologic Framework region) 15 Ocean 16 CCR extension 17 GV extension 18 SN extension 19 WBR extension 20 ST extension 21 PR extension
Google Pages
Google Colab notebook: CTM_tools.ipynb
files used by the query code:
