Final CCA06 Crossings Result

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 Basin Query Result 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 without ELY GTL

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
Raw crossing result	data	data
Metadata result	data	data

Basin Query Algorithm Result for CCA 06 with ELY GTL

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
Raw crossing result	data	data
Metadata result	data	data

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

CCA06 no GTL, section of comparable area to Kevin's zoomed in area plot :

first crossing

first or second crossing

Scatter plot :

first crossing

Related Entries

• Basin Query Notes
• UCVM
• UCVM Review
• USR