Difference between revisions of "CTM"

From SCECpedia
Jump to navigationJump to search
Line 17: Line 17:
 
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.
 
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 Data Format ==
+
== 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.
 
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.
 
<pre>
 
<pre>
Line 28: Line 28:
 
-122.166584 29.731010 -100000.000000 1.000000 1316.400024 1.000000 -4157.000000
 
-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
 
-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
 +
</pre>
 +
 +
== 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
 +
<pre>
 +
-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
 +
</pre>
 +
 +
Same points in GFM cell file format
 +
<pre>
 +
long. lat. Z regionID temp_median heatRegionID elevation
 
-120.00  35.00  -4532.0 1.000000 1316.400024 1.000000 468.40
 
-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
 
-119.00  35.00  -4556.0 1.000000 1316.400024 1.000000 444.00

Revision as of 22:01, 24 July 2020

We are using jupyter notebook developed by Liz Hearn to calculate CTM values at points defined for the GFM.

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.

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:

Related Entries