Fix boundary conditions
Many changes to PML and main routines