Begin port the spatialDecompose.m codes to Fortran

Fortran loops are much much faster than Octave.