Added fix to number of particles to remove for new remove_surf_pts flag

src/refine_surface.f90

@@ -210,6 +210,8 @@ if (params%remove_surf_pts) then
   allocate(remove_list(nrem))
   irem=0
   nremp=0
+else
+  nrem=0
 endif
 
 do iedge=1,nedge
@@ -271,8 +273,6 @@ deallocate (edswap)
 
 nnodemax=surface%nsurface+nadd
 nelemmax=surface%nt+(nadd-naddp)*2+naddp
-!nnodemax=surface%nsurface+(nadd-nrem)
-!nelemmax=surface%nt+((nadd-naddp)*2+naddp)-((nrem-nremp)*2+nremp)
 nedgemax=nnodemax+nelemmax-1
 
 if (nedgemax.ne.nedgen) call stop_run ('error in counting number of edges in refine_surface$')