Difference between revisions of "Final CCA06 Crossings Result"

From SCECpedia
Jump to navigationJump to search
Line 1: Line 1:
 +
As the result of this crossing review study, it is apparent that it is not possible to
 +
have a single crossing algorithm that is applicable to all the velocity models and
 +
a reasonable crossing algorithm requires domain scientists' post processing evaluation.
 +
Therefore, we are retaining the 5 algorithms that were produced from the *[[ ]] and
 +
also create a new output file that contains raw crossing information that can be used
 +
in later time.
 +
 +
== raw crossing result file ==
  
 
An option was added in basin_query_mpi(-o) to create an ascii file of top 3 crossings plus the last crossing per grid point within ROI region. This result file can be post processed to determine the best
 
An option was added in basin_query_mpi(-o) to create an ascii file of top 3 crossings plus the last crossing per grid point within ROI region. This result file can be post processed to determine the best

Revision as of 21:27, 19 March 2018

As the result of this crossing review study, it is apparent that it is not possible to have a single crossing algorithm that is applicable to all the velocity models and a reasonable crossing algorithm requires domain scientists' post processing evaluation. Therefore, we are retaining the 5 algorithms that were produced from the *[[ ]] and also create a new output file that contains raw crossing information that can be used in later time.

raw crossing result file

An option was added in basin_query_mpi(-o) to create an ascii file of top 3 crossings plus the last crossing per grid point within ROI region. This result file can be post processed to determine the best strategy to use in selecting the suitable crossing to use for the model.

The format of the content:

lon lat count crossing1 crossing2 crossing3 crossing4

examples,

 *no crossing
     lon  lat  0  -1  -1  -1  -1
 *1 crossing (c1 = first crossing = last crossing)
     lon  lat  1  c1  -1  -1  -1
 *2 crossings (c1 = first, c2 = 2nd = last)
     lon  lat  2  c1  c2  -1  -1
 *3 crossings (c1 = first, c2 = 2nd, c3 = third = last)
     lon  lat  3  c1  c2  c3  -1
 *4 crossings (c1 = first, c2 = 2nd, c3 = third, c4 = fourth = last)
     lon  lat  4  c1  c2  c3  c4
 *5 crossings (c1 = first, c2 = 2nd, c3 = third, c5 = fifth = last)
     lon  lat  5  c1  c2  c3  c5
 *6 crossings (c1 = first, c2 = 2nd, c3 = third, c6 = sixth = last)
     lon  lat  6  c1  c2  c3  c6


Following is the final set of CCA06 crossing result with updated projection configuration from CVM_Projection_Issue_Discussion.

The target crossing to be used for CCA06 is the 2nd row, the firstOrSecond set.


Basin Query Algorithm Result for CCA 06

Algorithm Threshold = 1000m/s Threshold = 2500m/s
First Crossing data

png

data

png

First Or Second Crossing data

png

data

png

Last Crossing data

png

data

png|

Second Only Crossing data

png

data

png

Last but

not First or Second Crossing

data

png

data

png

The basin query command that was used to generate the result:

 mpirun -np 20 ./basin_query_mpi 
 -b ${TEST_TOP_PATH}/RESULT/cca_z2.5.first,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.firstOrSecond,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.last,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.secondOnly,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.threeLast 
 -o ${TEST_TOP_PATH}/RESULT/cca_z2.5.result,
 ${TEST_TOP_PATH}/RESULT/cca_z2.5.json 
 -f ucvm.conf -m cca -i 20 -v 2500 -l 33.35,-123.0 -s 0.005 -x 1551 -y 1201

and,

 mpirun -np 20 ./basin_query_mpi 
 -b ${TEST_TOP_PATH}/RESULT/cca_z1.0.first,
 ${TEST_TOP_PATH}/RESULT/cca_z1.0.firstOrSecond,
 ${TEST_TOP_PATH}/RESULT/cca_z1.0.last,
 ${TEST_TOP_PATH}/RESULT/cca_z1.0.secondOnly,
 ${TEST_TOP_PATH}/RESULT/cca_z1.0.threeLast 
 -o ${TEST_TOP_PATH}/RESULT/cca_z1.0.result 
 -f ucvm.conf -m cca -i 20 -v 1000 -l 33.35,-123.0 -s 0.005 -x 1551 -y 1201 


and the plotting commands that were used are:

 ./plot_z25_map.py -f cca_z2.5.firstOrSecond 
 -o cca_z2.5_firstOrSecond_map.png 
 -e "(first or second crossing)" -b 33.349998,-123.000000 
 -u 39.355000,-115.245003 -c cca -s 0.005000 -x 1551 -y 1201

and,

 ./plot_z10_map.py -f cca_z1.firstOrSecond
 -o cca_z1_firstOrSecond_map.png -e "(first or second crossing)"
 -b 33.349998,-123.000000 -u 39.355000,-115.245003 -c cca 
 -s 0.005000 -x 1551 -y 1201


Related Entries