remove lt_work

erosion.f90

@@ -67,7 +67,7 @@ integer ntmax,nhmax,npmax,nmax,nmode,nvmax,nnpnmax,int_method
 double precision,dimension(:),allocatable::field,wpoints
 double precision,dimension(:,:),allocatable::points,centres,pointsmin,pointsmax,ppoints
 integer,dimension(:,:),allocatable::vertices,neighbour
-integer,dimension(:),allocatable::hulltriangles,lt_work,ln_work
+integer,dimension(:),allocatable::hulltriangles,ln_work
 logical*1,dimension(:),allocatable::lt_work_l,ln_work_l
 integer,dimension(:),allocatable::vis_tlist,vis_elist,add_tlist
 logical clockwise,out,nnext
@@ -282,6 +282,7 @@ allocate (hulltriangles(nhmax),stat=err)      ; if (err.ne.0) call stop_run ('Er
 allocate (vis_tlist(nvmax),stat=err)          ; if (err.ne.0) call stop_run ('Error alloc vis_tlist in erode$')
 allocate (vis_elist(nvmax),stat=err)          ; if (err.ne.0) call stop_run ('Error alloc vis_elist in erode$')
 allocate (add_tlist(nvmax),stat=err)          ; if (err.ne.0) call stop_run ('Error alloc add_tlist in erode$')
+allocate (ln_work(surface%nsurface),stat=err) ; if (err.ne.0) call stop_run ('Error alloc ln_work in erode$')
 allocate (lt_work_l(ntmax),stat=err)            ; if (err.ne.0) call stop_run ('Error alloc lt_work_l in erode$')
 allocate (ln_work_l(surface%nsurface),stat=err) ; if (err.ne.0) call stop_run ('Error alloc ln_work_l in erode$')
 
@@ -342,7 +343,6 @@ if (surface%spread_surface_points.gt.0) then
 
 allocate (pointsmin(2,surface%nsurface),stat=err) ; if (err.ne.0) call stop_run ('Error alloc pointsmin in erode$')
 allocate (pointsmax(2,surface%nsurface),stat=err) ; if (err.ne.0) call stop_run ('Error alloc pointsmax in erode$')
-allocate (ln_work(surface%nsurface),stat=err) ; if (err.ne.0) call stop_run ('Error alloc ln_work in erode$')
 
 points=0.d0
 pointsmin=2.d0
@@ -392,7 +392,6 @@ deallocate (hulltriangles)
 deallocate (vis_tlist)
 deallocate (vis_elist)
 deallocate (add_tlist)
-deallocate (lt_work)
 deallocate (lt_work_l)
 deallocate (ln_work)
 deallocate (ln_work_l)