Fixed issue with material properties right above model boundary.

Author: scecdavid

.gitignore

@@ -30,3 +30,5 @@
 
 # Debug files
 *.dSYM/
+
+.DS_Store

.project

@@ -0,0 +1,17 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<projectDescription>
+	<name>CCA</name>
+	<comment></comment>
+	<projects>
+	</projects>
+	<buildSpec>
+		<buildCommand>
+			<name>org.python.pydev.PyDevBuilder</name>
+			<arguments>
+			</arguments>
+		</buildCommand>
+	</buildSpec>
+	<natures>
+		<nature>org.python.pydev.pythonNature</nature>
+	</natures>
+</projectDescription>

.pydevproject

@@ -0,0 +1,5 @@
+<?xml version="1.0" encoding="UTF-8" standalone="no"?>
+<?eclipse-pydev version="1.0"?><pydev_project>
+<pydev_property name="org.python.pydev.PYTHON_PROJECT_INTERPRETER">Default</pydev_property>
+<pydev_property name="org.python.pydev.PYTHON_PROJECT_VERSION">python 2.7</pydev_property>
+</pydev_project>

data/config

@@ -14,11 +14,11 @@ depth = 50000
 depth_interval = 500
 
 # Corners, in UTM projection.
-top_left_corner_e		= 504472.106530
+top_left_corner_e	= 504472.106530
 top_left_corner_n   	= 4050290.088321  
 bottom_left_corner_e  	= 779032.901477
 bottom_left_corner_n  	= 3699450.983449
-top_right_corner_e	  	= 905367.666914 
-top_right_corner_n		= 4366503.431060
+top_right_corner_e  	= 905367.666914 
+top_right_corner_n	= 4366503.431060
 bottom_right_corner_e 	= 1181157.928029
 bottom_right_corner_n	= 4014713.414724

src/cca.c

@@ -137,15 +137,15 @@ int cca_query(cca_point_t *points, cca_properties_t *data, int numpoints) {
 		point_utm_n = sin_rotation_angle * temp_utm_e + cos_rotation_angle * temp_utm_n;
 
 		// Which point base point does that correspond to?
-		load_x_coord = floor(point_utm_e / total_width_m * configuration->nx);
-		load_y_coord = floor(point_utm_n / total_height_m * configuration->ny);
+		load_x_coord = floor(point_utm_e / total_width_m * (configuration->nx - 1));
+		load_y_coord = floor(point_utm_n / total_height_m * (configuration->ny - 1));
 
 		// And on the Z-axis?
 		load_z_coord = (configuration->depth / configuration->depth_interval - 1) -
 					   floor(points[i].depth / configuration->depth_interval);
 
 		// Are we outside the model's X and Y boundaries?
-		if (load_x_coord > configuration->nx - 1 || load_y_coord > configuration->ny - 1 || load_x_coord < 0 || load_y_coord < 0) {
+		if (load_x_coord > configuration->nx - 2 || load_y_coord > configuration->ny - 2 || load_x_coord < 0 || load_y_coord < 0) {
 			data[i].vp = -1;
 			data[i].vs = -1;
 			data[i].rho = -1;