diff --git a/embed_surface_in_octree.f90 b/embed_surface_in_octree.f90
index 3222657e65ce2666a5a75f8863d1c5c0cb59c2f8..b47ba765d9e29febd8dcdef9d8e099f6c1cbac25 100644
--- a/embed_surface_in_octree.f90
+++ b/embed_surface_in_octree.f90
@@ -200,6 +200,8 @@ call heap (threadinfo,'osolve%z',   'embed_surf...',size(osolve%z),'dp',+1)
 call octree_find_node_connectivity (osolve%octree,osolve%noctree, &
                                     osolve%icon,osolve%nleaves,   &
                                     osolve%x,osolve%y,osolve%z,osolve%nnode)
+! osolve%nnode has been changed by octree_find_node_connectivity, re-size x, y, z
+call octreesolve_shrink_xyz(threadinfo, osolve, osolve%nnode)
 
 ! now that osolve%nnode is known we can allocate osolve%lsf
 allocate (osolve%lsf(osolve%nnode,osolve%nlsf),stat=threadinfo%err)