Removed unused variables from seismic._wavefd
Tried a different free surface bc for P-SV. Still unstable
Simplified interface to wavefd.elastic_sh. Sources now take the x, z coordinates and convert to indexes