Add exception handling when writing to the netCDF. This hopefully improves debugging somewhat.
Add option to control the fill value in the netCDF (now defaults to no fill value specified, rather than -999.0).
Fix the exception handling to report the shape of the arrays which do not match. Also remove the finally whose function I misunderstood.
Add the global attributes to the attributes dict if it has been requested.
Make the code mildly more pythonic.
Fix the order of the variables when reporting an IndexError.
Add note about a partial success in removing the eval approach to extracting the data.
Remove the math import.
Update the code to be compatible with Python 3 (this appears to have had no impact on the analysis).