diff --git a/src/DOUAR.f90 b/src/DOUAR.f90 index 8db8e78453355bb792068644411cd6a1dce026d0..99dc93d28721871bc3fcb330617864bdd89b8679 100644 --- a/src/DOUAR.f90 +++ b/src/DOUAR.f90 @@ -430,8 +430,10 @@ do while (istep.le.params%nstep) ! refinement/improvement of osolve !------------------------------------------------------------------------| !------------------------------------------------------------------------| - call improve_osolve (osolve,ov,params,threadinfo,istep,iter,total,step,inc, & - current_level,boxes,cube_faces,increase_current_level ) + call improve_osolve (osolve,ov,params,threadinfo,istep,iter,total,step, & + inc,current_level,boxes,cube_faces) + !call improve_osolve (osolve,ov,params,threadinfo,istep,iter,total,step,inc, & + ! current_level,boxes,cube_faces,increase_current_level ) nleaves_new_mem1=osolve%octree(2) diff --git a/src/improve_osolve.f90 b/src/improve_osolve.f90 index fb8e117303f2354dfaeee529ff2cc7bbe66c425f..72f0ae1846718801286b6d329d1860e3872be3b1 100644 --- a/src/improve_osolve.f90 +++ b/src/improve_osolve.f90 @@ -17,8 +17,11 @@ !------------------------------------------------------------------------------| subroutine improve_osolve (osolve,ov,params,threadinfo,istep,iter,total,step, & - inc,refine_level,boxes,cube_faces, & - increase_current_level) + inc,refine_level,boxes,cube_faces) + +!subroutine improve_osolve (osolve,ov,params,threadinfo,istep,iter,total,step, & +! inc,refine_level,boxes,cube_faces, & +! increase_current_level) !------------------------------------------------------------------------------| !(((((((((((((((( Purpose of the routine )))))))))))))))))))))))))))))))))))))) @@ -70,22 +73,21 @@ double precision total,step,inc integer refine_level type (box) boxes(params%nboxes) type (face), dimension(6) :: cube_faces -logical increase_current_level +!logical increase_current_level !------------------------------------------------------------------------------| !(((((((((((((((( declaration of the subroutine internal variables ))))))))))))) !------------------------------------------------------------------------------| -double precision jcb(3,3),jcbi(3,3) double precision,dimension(:),allocatable :: x,y,z,vx,vy,vz double precision,dimension(:),allocatable :: dhdx,dhdy,dhdz -double precision,dimension(:),allocatable :: crit,crittemp,xc,yc,zc -double precision :: eps,r,s,t,volume,w +double precision,dimension(:),allocatable :: crit,crittemp +double precision :: eps,r,s,t,volume double precision :: exx,eyy,ezz,exy,eyz,ezx double precision critmax,xxx,yyy,zzz,b,l double precision :: x0,y0,z0,x1,y1,z1,E2d -integer err,ierr,nproc,iproc,k,i,iint,nint,mpe -integer :: j,nn,nb_of_boxes,level_box,nx,ny,nz,level,nv,nh +integer ierr,nproc,iproc,k,i,nint,mpe +integer :: j,nn,level_box,nx,ny,nz,level,nv,nh character*72 shift !------------------------------------------------------------------------------|