Updated grid

I changed the grid to accommodate different length systems and to have
the grid points be centered in cells rather than on cell edges.
Fixed bug

Accidentally introduced bug and fixed it in line 8 of
F1D_Initial_Condition.py. I also think that sparse matrices will be
required to solve for a realistic number of time steps.
Fixed a bug

Fixed a bug that used two 40Ar concentrations instead of 36Ar and 40Ar
and led to null enrichment result.
Back up

The is a save point before getting help from Chris Patton.
Controlling F1D_Plot.py.
added a script for reading data
Working spinel measurer

One problem remains: how to measure the actual width. The part of the
code that is suppose to be doing that seems to be doing something else.
It always reports a length of about 200. No idea what is wrong. Need
help. the pixel counter seems to work.
selectivity improvement 1

Attempted to make my spinel detection algorithm selective for the main
spinel crystal, not isolated fragments
Alex edits

Alex used a slightly different method to count the pixels and fit an
equation to the resulting data. The next step will require replacing
the dummy function with the kinetics equation from Watson and Price,
and rescaling the data and function so that the results have correct
units. Almost done!
Pre Neumann

Changing boundary conditions to Neuman B. C. This is a backup save.