introduced automatic determination of the factor between chebyshev interval lengths

source/common/radial_solver_Chebychev/create_newmesh.f90

@@ -72,7 +72,7 @@ contains
 
     ! .. Local variables
 !   real (kind=dp), parameter :: fac=2e0_dp
-    real (kind=dp), parameter :: fac=1.2e0_dp
+    real (kind=dp) :: fac
     integer :: npan_inst, i_stat, i_all
     integer :: i1, ir2, ip, icell
     integer :: imin, imax, iminnew, imaxnew, lm1
@@ -105,6 +105,7 @@ contains
       ! log panel
       rmin = rmesh(2, i1)
       rmax = r_log
+      fac = (rmax/rmin)**(1.d0/npan_log)
       rval = 0e0_dp
       ishift = 0
       if (r_log>rmesh(irmin(i1),i1)) then
@@ -162,6 +163,10 @@ contains
 
       call chebmesh(npan_tot(i1), ncheb, rpan_intervall(0:,i1), rnew(1,i1))
 
+!     do ir2 = 1,50
+!     write(6,*) ir2, rnew(ir2,i1)
+!     end do
+
       ! do interpolation only when optional arguments are given
       if (present(ntcell)) then
         ! interpolate shape function THETAS to new shape function THETASNEW