Bill wants me to make four changes to his code, see source:trunk/kirkendall/A-V_phase-field.py@1457.

He needs to make four changes.

• Update to be an expression involving . This is done in r1458.

• The second item is to check source:trunk/kirkendall/A-V_phase-field.py@1457#L71 against Bill's notes. It seems correct other than that he left out the diffusion term in his notes by mistake.

• Save the data from the simulations for later use and read the saved files and plot the results independently from the simulations. See diff:@1458:1460. This now writes a data file directory of the form data/A-V_phase-field_17_36_53Thu02Feb2012. The data files are stored as numpy zip files .npz and numbered by their time step value. plotdata.py reads and plots from the latest data directory created.

The name of an alternative path can be specified as an argument

$ python plotdata.py data/A-V_phase-field_17_33_42Thu02Feb2012

This will use the specified directory, while

$ python plotdata.py

will use the latest created data directory. plotdata.py plots the data for the four main fields using the stored data.

• Add an improved method for extracting the 0.5 level set. See r1461. This is done using scipy's interp1d. After plotting the four main fields for every stored timestep, plotdata.py will throw up a plot of the interface position.