First working version. Boundary conditions are not entirely correct. Plot file now reads n automatically.