add octreesolve_shrink_xyz after ONE of the calls to octree_find_node_connectivity

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)