The code is no different than the standard PWEM code. Search the discussion groups for code from Guo (or Gou?). From there, rather than calculate the eigen-frequencies just around the perimeter of the Brillouin zone, you will calculate the eigen-frequencie inside the entire Brillouin zone. The isofrequency contours in this space are your EFSs.