Skip to content
Snippets Groups Projects
Commit 03f170aa authored by Dave Whipp's avatar Dave Whipp
Browse files

Added initial allocation of ov%octree

parent 6a92023e
No related branches found
No related tags found
No related merge requests found
...@@ -48,7 +48,7 @@ type (nest_info) nest ...@@ -48,7 +48,7 @@ type (nest_info) nest
!(((((((((((((((( declaration of the subroutine internal variables ))))))))))))) !(((((((((((((((( declaration of the subroutine internal variables )))))))))))))
!------------------------------------------------------------------------------| !------------------------------------------------------------------------------|
type (octreev) ovls type (octreev) :: ovls
integer :: i,j,k,nface,nlsf,np,kfixls,matnum,err integer :: i,j,k,nface,nlsf,np,kfixls,matnum,err
logical :: is_plas logical :: is_plas
double precision :: eps,xls,yls,zls,s,crit,e2d,evisc,dilatr,xlsf,current_time double precision :: eps,xls,yls,zls,s,crit,e2d,evisc,dilatr,xlsf,current_time
...@@ -59,6 +59,9 @@ double precision :: eps,xls,yls,zls,s,crit,e2d,evisc,dilatr,xlsf,current_time ...@@ -59,6 +59,9 @@ double precision :: eps,xls,yls,zls,s,crit,e2d,evisc,dilatr,xlsf,current_time
eps=1.d-10 eps=1.d-10
ovls%noctree=params%noctreemax ovls%noctree=params%noctreemax
allocate (ov%octree(ov%noctree),stat=err)
!if (params%debug.gt.1) call heap (threadinfo,'ov%octree','main',size(ov%octree),'int',+1)
! Read in coarse model octree; needs to be stored in params @ end of step ! Read in coarse model octree; needs to be stored in params @ end of step
open (9,file=trim(nest%lsoutfile),status='old',form='unformatted') open (9,file=trim(nest%lsoutfile),status='old',form='unformatted')
read (9) ovls%octree(3),& read (9) ovls%octree(3),&
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment