I found a paper that describes this calculation, you would need to do some post-processing of the field data on your own.
http://fy.chalmers.se/~klavs/papers/Nano_2006_field-enhance.pdf
You might also have to run a simulation with a simpler geometry as a kind of baseline to get the macroscopic applied field they talk about in the paper.
Regards