Skip to content
Snippets Groups Projects
Commit e279eb21 authored by Rudolf Zeller's avatar Rudolf Zeller
Browse files

Removed inconsistency in the initial definition of the B parameter

parent 59e951e7
No related branches found
No related tags found
No related merge requests found
......@@ -160,17 +160,18 @@ module RadialMeshData_mod
double precision, intent(in) :: radius_mt
double precision, intent(in) :: a_log
double precision, intent(in) :: b_log
double precision :: b_log_local
!double precision, parameter :: A = 0.025d0
integer :: ii
self%A = a_log
self%B = b_log
!self%B = radius_mt / (exp(a_log * (imt - 1)) - 1.d0)
b_log_local = radius_mt / (exp(a_log * (imt - 1)) - 1.d0)
self%B = b_log_local
do ii = 1, imt
self%r(ii) = b_log * (exp(a_log * (ii - 1)) - 1.d0)
self%drdi(ii) = a_log * b_log * exp(a_log * (ii - 1))
self%r(ii) = b_log_local * (exp(a_log * (ii - 1)) - 1.d0)
self%drdi(ii) = a_log * b_log_local * exp(a_log * (ii - 1))
enddo ! ii
self%rmt = radius_MT
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment