Basin Query Result - Revision history

Scottcal: /* Basin Query Algorithms in Actual Implementation */

2017-12-07T15:28:37Z

Basin Query Algorithms in Actual Implementation

Maechlin at 22:41, 6 December 2017

2017-12-06T22:41:08Z

Maechlin: /* Basin Query Algorithms in Actual Implementation */

2017-12-06T22:40:36Z

Basin Query Algorithms in Actual Implementation

Maechlin at 22:34, 6 December 2017

2017-12-06T22:34:57Z

Maechlin at 22:10, 6 December 2017

2017-12-06T22:10:23Z

Mei: /* Basin Query Algorithm Result for CCA 06 */

2017-12-04T19:17:41Z

Basin Query Algorithm Result for CCA 06

Maechlin at 18:55, 4 December 2017

2017-12-04T18:55:03Z

Mei: /* Basin Query Algorithm Result for CCA 06 */

2017-12-04T06:16:35Z

Basin Query Algorithm Result for CCA 06

Mei: /* Basin Query Algorithm Result for CCA 06 */

2017-12-01T19:35:05Z

Basin Query Algorithm Result for CCA 06

Mei: /* Basin Query Algorithm Result for CCA 06 */

2017-12-01T19:34:45Z

Basin Query Algorithm Result for CCA 06

@@ Line 9: / Line 9: @@
 *model_depth = This is the depth at which we query the velocity model. We start at current-depth=0 which is the surface.
 *model_vs = This is the Vs returned by the model under evaluation, and a model_vs always has an associated model_depth.
-*crossing = each time we query the velocity_model, we check to determine whether the model_vs is greater than the target_vs. If this is the first model_vs, or if the previous model_vs was below the target_vs, and now the model_vs is larger than the target_vs, then we consider it a crossing.
+*crossing = each time we query the velocity_model, we check to determine whether the model_vs is greater than or equal to the target_vs. If this is the first model_vs, or if the previous model_vs was below the target_vs, and now the model_vs is larger than the target_vs, then we consider it a crossing.
-#Search of Vs threshold crossings of a grid point starts by looking at the Vs of this grid point from the first depth layer to the very last  depth layer.
+#Search of Vs threshold crossings of a grid point starts by looking at the Vs of this grid point from the first depth layer to the very last depth layer.
 #If the Vs is less than 0.0, it is invalid, just skip it.
-#If the Vs is greater than 0.0, check if this is a crossing point. A crossing is found if Vs of the earlier depth layer is lower than the preset threshold and the Vs of the current  depth layer is higher than the preset threshold, ie. an upward sloping crossing.
+#If the Vs is greater than 0.0, check if this is a crossing point. A crossing is found if Vs of the earlier depth layer is lower than the preset threshold and the Vs of the current depth layer is higher than or equal to the preset threshold, ie. an upward sloping crossing.
 different types of result for all the grid points are being collected and outputted into 5 different binary data files. User may selectively collect each or all data files by supplying the matching data file names on command line.

@@ Line 6: / Line 6: @@
 We define three terms.
-*target_vs = This is either 1000m (Z1.0) or 2500m (Z2.5). We are trying to determine depth at which this is Vs is observed
+*target_vs = This is either 1000m (Z1.0) or 2500m (Z2.5). Our Basin_Query script is trying to determine depth at which this is Vs is observed
 *model_depth = This is the depth at which we query the velocity model. We start at current-depth=0 which is the surface.
 *model_vs = This is the Vs returned by the model under evaluation, and a model_vs always has an associated model_depth.

@@ Line 9: / Line 9: @@
 *model_depth = This is the depth at which we query the velocity model. We start at current-depth=0 which is the surface.
 *model_vs = This is the Vs returned by the model under evaluation, and a model_vs always has an associated model_depth.
-*crossing = each time we query the velocity_model, we check to determine whether the model_vs is greater than the target_vs. If the model_vs is larger than the target_vs, then we consider it a crossing.
+*crossing = each time we query the velocity_model, we check to determine whether the model_vs is greater than the target_vs. If this is the first model_vs, or if the previous model_vs was below the target_vs, and now the model_vs is larger than the target_vs, then we consider it a crossing.
-Search of Vs threshold crossings of a grid point starts by looking at
+#Search of Vs threshold crossings of a grid point starts by looking at the Vs of this grid point from the first depth layer to the very last  depth layer.
-the Vs of this grid point from the first depth layer to the very last
+#If the Vs is less than 0.0, it is invalid, just skip it.
-depth layer.
+#If the Vs is greater than 0.0, check if this is a crossing point. A crossing is found if Vs of the earlier depth layer is lower than the preset threshold and the Vs of the current  depth layer is higher than the preset threshold, ie. an upward sloping crossing.
 different types of result for all the grid points are being collected and outputted into 5 different binary data files. User may selectively collect each or all data files by supplying the matching data file names on command line.

@@ Line 4: / Line 4: @@
 with different depths in the increment of 20m. This set of query points are submitted to the
 backend model for their matching Vs.
 Search of Vs threshold crossings of a grid point starts by looking at

@@ Line 85: / Line 85: @@
 -o cca_z2.5_first_map.png
-== Basin Query Result ==
+== Southern California Basin Query Result ==
-These result were ran on HPC and Blue Waters
+*[[Southern California Basin Query Results]]
 == Related Entries ==

@@ Line 43: / Line 43: @@
 [http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_firstOrSecond_map.png png]
 || [http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5.firstOrSecond  data]
-[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5.firstOrSecond_map.png png]
+[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_firstOrSecond_map.png png]
 |-
 ! scope="row" | Last Crossing

← Older revision		Revision as of 18:55, 4 December 2017
Line 1:		Line 1:
−	~~== Basin Query Algorithms ==~~	+	== Basin Query Algorithms in Actual Implementation ==
−	~~Both Z1.0 and Z2.5 will use the same algorithms, they will just look for different Vs values Z1.0 = Vs = 1000m/s Z2.5 = Vs = 2500m/s, which we call the Target_Z.~~
−
−	*First Crossing - If the surface Vs is > Target_Z, then depth is 0. Otherwise, query the CVM every 20m depth, until a Vs >= Target_Z is found. Return the depth at which the Vs returned from the model is greater than or equal to the Target-Z.
−
−	*Second Crossing -
−
−	*Last Crossing -
−
−	== Basin Query Algorithms in ~~actual implmentation~~ ==

	For each grid point within the user supplied boundary, a stack of query points is generated		For each grid point within the user supplied boundary, a stack of query points is generated

@@ Line 79: / Line 79: @@
 CCA06: [minLat = 33.35, minLon = -123.0, nx = 1551, ny = 1201, gridSpacing = 0.005 degrees]
-mpirun -np 20 ./basin_query_mpi -b ${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z
+mpirun -np 20 ./basin_query_mpi
-.5_n20.first,${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.firstOrSecond
+-b ${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.first,
-,${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.last -o ${TEST_TOP_PATH}/R
+${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.firstOrSecond,
-ESULT/basin_query_mpi_cca_z2.5_n20.result,${TEST_TOP_PATH}/RESULT/basin_query_mp
+${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.last,
-i_cca_z2.5_n20.secondOnly,${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.t
+${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.secondOnly,
-hreeLast -f ucvm.conf -m cca -i 20 -v 2500 -l 33.35,-123.0 -s 0.005 -x 1551 -y 1201
+${TEST_TOP_PATH}/RESULT/basin_query_mpi_cca_z2.5_n20.threeLast
 A plot command that was used to generate one of the png file,
 (data file were renamed to simplify the batch script)
-./plot_z25_map.py -b 33.35,-123.0 -u 39.35,-115.25 -c cca -s 0.005 -f cca_z2.5.first -o cca_z2.5_first_map.png
+./plot_z25_map.py -b 33.35,-123.0 -u 39.35,-115.25 -c cca -s 0.005 -f cca_z2.5.first
 == Basin Query Result ==

@@ Line 71: / Line 71: @@
 [http://hypocenter.usc.edu/research/ucvmc_result/cca_z1_threeLast_map.png png]
 || [http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5.threeLast  data]
-[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_threeLast_map.png png]|
+[http://hypocenter.usc.edu/research/ucvmc_result/cca_z2.5_threeLast_map.png png]
 |-
 |}