DOUAR.f90

   deallocate (surface(is)%icon)
end do

deallocate (surface)
deallocate (surface0)
deallocate (mat)

if (params%isobc) then
  deallocate (bcdef%zisodisp)
  deallocate (bcdef%zisoinc)
endif
deallocate (params%materialn)
deallocate (params%surface_for_mat_props)

if (params%debug.gt.1) then
  call heap_final (threadinfo)
  close (threadinfo%Logunit)
  close (threadinfo%mem_heap_unit)
endif

call mpi_finalize (ierr)

end