Added ImmersedField : copy the value of the field just outside the domain. Should act as an effective null flux boundary condition. Added python code to do 2d animated scatter plots.