diff --git a/build_surface_octree.f90 b/build_surface_octree.f90 index bb99d1d95f85ed3163f31613b846388795019b20..8f1084a7e75b635dff2e24f3e9d526e62464c1e5 100644 --- a/build_surface_octree.f90 +++ b/build_surface_octree.f90 @@ -136,7 +136,8 @@ elseif (surface%criterion.eq.2) then ! anglemax=surface%anglemaxoctree/(0.5d0**surface%leveloct) angle=acos(cosd) anglemax=surface%anglemaxoctree - level=max(level,params%leveluniform_oct+ nint(min(angle/surface%anglemaxoctree,1.d0)*(params%levelmax_oct-params%leveluniform_oct))) + level=max(level,params%leveluniform_oct+ & + nint(min(angle/surface%anglemaxoctree,1.d0)*(params%levelmax_oct-params%leveluniform_oct))) ! level=max(level,params%leveluniform_oct+ nint(min(angle/surface%anglemaxoctree,1.d0)*(surface%leveloct-params%leveluniform_oct))) ! level=max(level,params%leveluniform_oct+ nint(min(angle/surface%anglemax,1.d0)*(surface%leveloct-params%leveluniform_oct))) enddo