We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent d7603bf commit c2877a7Copy full SHA for c2877a7
1 file changed
SRC/simple.f90
@@ -342,16 +342,17 @@ subroutine trace_to_cut(self, si, f, z, var_cut, cut_type, ierr)
342
end do
343
end subroutine trace_to_cut
344
345
- subroutine fract_dimension(ntr,rt0,fraction)
+ subroutine fract_dimension(ntr,rt,fraction)
346
347
integer, parameter :: iunit=1003
348
integer :: itr,ntr,ngrid,nrefine,irefine,kr,kt,nboxes
349
double precision :: fraction,rmax,rmin,tmax,tmin,hr,ht
350
- double precision, dimension(2,ntr) :: rt0, rt
+ double precision, dimension(2,ntr) :: rt!0, rt
351
logical, dimension(:,:), allocatable :: free
352
353
- rt(1,:) = rt0(1,:)*cos(rt0(2,:))
354
- rt(2,:) = rt0(1,:)*sin(rt0(2,:))
+! TODO: check if this works better on tips only
+! rt(1,:) = rt0(1,:)*cos(rt0(2,:))
355
+! rt(2,:) = rt0(1,:)*sin(rt0(2,:))
356
357
rmin=minval(rt(1,:))
358
rmax=maxval(rt(1,:))
0 commit comments