added a function ESRI_ASCIIGrid to retrieve the cell by coordinates

impls/square-cell/vn-neighbours/1-levels/simple/ESRI_ASCIIGrid.hpp

@@ -932,9 +932,36 @@ namespace CA {
 
 
 
-		}
+		}// end of copyData
+		//--------------------------------------------------------------------------------------------------------------
+
+		//--------------------------------------------------------------------------------------------------------------
+		inline T getCellByCords(double XCord, double YCord)
+		{
+
+
+			int XCol = int(XCord - AsciiGridGeneral<T>::xllcorner) / AsciiGridGeneral<T>::cellsize;
+
+			// note , cords go from bottom left and rows,cols go grom top left, therefore need to flip this...
+
+			//int YRow = int(YCord - AsciiGridGeneral<T>::yllcorner) / AsciiGridGeneral<T>::cellsize;
+
+			int YRow = AsciiGridGeneral<T>::nrows - (int(YCord - AsciiGridGeneral<T>::yllcorner) / AsciiGridGeneral<T>::cellsize);
+
+			T returnResult = AsciiGridGeneral<T>::data[XCol + (AsciiGridGeneral<T>::ncols * YRow)];
+
+
+			std::cout << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!" << std::endl;
+			std::cout << "Local Cords ["<< XCol<<","<<YRow<<"]: " << returnResult<< std::endl;
+			std::cout << "!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!" << std::endl;
+
+			return returnResult;
+
+
+		}// end of getCellByCords
+		//--------------------------------------------------------------------------------------------------------------
+
 
-		
 
 
 	};// end of class