From 39c6d04310ed0217aab119db62e129e7ca8d7736 Mon Sep 17 00:00:00 2001
From: Douglas Guptill <douglas.guptill@dal.ca>
Date: Tue, 14 Jul 2009 22:17:08 +0000
Subject: [PATCH] add octreesolve_shrink_xyz after ONE of the calls to
octree_find_node_connectivity
---
embed_surface_in_octree.f90 | 2 ++
1 file changed, 2 insertions(+)
diff --git a/embed_surface_in_octree.f90 b/embed_surface_in_octree.f90
index 3222657e..b47ba765 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)
--
GitLab