From 25661c95dc0bf8256004f749c43808f43737cc41 Mon Sep 17 00:00:00 2001
From: Rudolf Zeller <ru.zeller@fz-juelich.de>
Date: Sat, 30 May 2020 13:07:36 +0200
Subject: [PATCH] introduced automatic determination of the factor between
chebyshev interval lengths
---
source/common/radial_solver_Chebychev/create_newmesh.f90 | 7 ++++++-
1 file changed, 6 insertions(+), 1 deletion(-)
diff --git a/source/common/radial_solver_Chebychev/create_newmesh.f90 b/source/common/radial_solver_Chebychev/create_newmesh.f90
index 8ce134ce8..ee3487922 100644
--- a/source/common/radial_solver_Chebychev/create_newmesh.f90
+++ b/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
--
GitLab