diff --git a/refine_surface.f90 b/refine_surface.f90
index 727144f8389e532997303e0f6e729e01fd01b6bd..084512ff1567630c4edcae3ff43eddc9e7575bdf 100644
--- a/refine_surface.f90
+++ b/refine_surface.f90
@@ -110,9 +110,12 @@ nnmax=12 ! nnmax is maximum number of neighbours in triangulation (should be aro
 
       
 allocate (ed(nedge),stat=threadinfo%err)                       ; call heap (threadinfo,'ed','refine_surface',size(ed),'int',+1)
-allocate (nedgepernode(nsurfacen),stat=threadinfo%err)         ; call heap (threadinfo,'nedgepernode','refine_surface',size(nedgepernode),'int',+1)
-allocate (nodenodenumber(nnmax,nsurfacen),stat=threadinfo%err) ; call heap (threadinfo,'nodenodenumber','refine_surface',size(nodenodenumber),'int',+1)
-allocate (nodeedgenumber(nnmax,nsurfacen),stat=threadinfo%err) ; call heap (threadinfo,'nodeedgenumber','refine_surface',size(nodeedgenumber),'int',+1)
+allocate (nedgepernode(nsurfacen),stat=threadinfo%err)
+call heap (threadinfo,'nedgepernode','refine_surface',size(nedgepernode),'int',+1)
+allocate (nodenodenumber(nnmax,nsurfacen),stat=threadinfo%err)
+call heap (threadinfo,'nodenodenumber','refine_surface',size(nodenodenumber),'int',+1)
+allocate (nodeedgenumber(nnmax,nsurfacen),stat=threadinfo%err)
+call heap (threadinfo,'nodeedgenumber','refine_surface',size(nodeedgenumber),'int',+1)
 
 ed%t1=0
 ed%t2=0