Rename from a_r_1d_rlm_r to ar_1d_rlm

hiromatsui · hiromatsui · commit 737e17e0b680 · 2024-12-20T23:59:11.000-08:00
diff --git a/src/Fortran_libraries/MHD_src/sph_MHD/sum_coriolis_terms_rlm.f90 b/src/Fortran_libraries/MHD_src/sph_MHD/sum_coriolis_terms_rlm.f90
@@ -127,7 +127,7 @@ subroutine select_sum_coriolis_rlm_10                             &
       if(sph_rlm%istep_rlm(1) .eq. 1) then
         call sum_coriolis_rlm_10                                        &
      &     (b_trns%base, sph_rlm%nnod_rlm, sph_rlm%nidx_rlm,            &
-     &      sph_rlm%idx_gl_1d_rlm_j, sph_rlm%a_r_1d_rlm_r,              &
+     &      sph_rlm%idx_gl_1d_rlm_j, sph_rlm%ar_1d_rlm,                 &
      &      omega_sph%ws_rlm, coef_cor, gt_cor%jgi_rlm, gt_cor%jei_rlm, &
      &      gt_cor%sw_rlm, gt_cor%sr_rlm, gt_cor%tr_rlm,                &
      &      gt_cor%sh_rlm, gt_cor%hh_rlm, gt_cor%th_rlm,                &
@@ -136,7 +136,7 @@ subroutine select_sum_coriolis_rlm_10                             &
       else
         call sum_coriolis_lrm_10                                        &
      &     (b_trns%base, sph_rlm%nnod_rlm, sph_rlm%nidx_rlm,            &
-     &      sph_rlm%idx_gl_1d_rlm_j, sph_rlm%a_r_1d_rlm_r,              &
+     &      sph_rlm%idx_gl_1d_rlm_j, sph_rlm%ar_1d_rlm,                 &
      &      omega_sph%ws_rlm, coef_cor, gt_cor%jgi_rlm, gt_cor%jei_rlm, &
      &      gt_cor%sw_rlm, gt_cor%sr_rlm, gt_cor%tr_rlm,                &
      &      gt_cor%sh_rlm, gt_cor%hh_rlm, gt_cor%th_rlm,                &
diff --git a/src/Fortran_libraries/MHD_src/sph_MHD/sum_rot_coriolis_terms_rlm.f90 b/src/Fortran_libraries/MHD_src/sph_MHD/sum_rot_coriolis_terms_rlm.f90
@@ -122,13 +122,13 @@ subroutine select_sum_rot_coriolis_rlm_10(b_trns, sph_rlm, leg,   &
 !
       if(sph_rlm%istep_rlm(1) .eq. 1) then
         call sum_rot_coriolis_rlm_10(b_trns%base,                       &
-     &    sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%a_r_1d_rlm_r,     &
+     &    sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%ar_1d_rlm,        &
      &    leg%g_sph_rlm, omega_sph%ws_rlm, coef_cor,                    &
      &    gt_cor%jgi_rlm, gt_cor%jei_rlm, gt_cor%sw_rlm, gt_cor%tw_rlm, &
      &    NB, n_WR, irev_sr_rlm, WR, d_rotp_cor_rlm, d_rott_cor_rlm)
       else
         call sum_rot_coriolis_lrm_10(b_trns%base,                       &
-     &    sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%a_r_1d_rlm_r,     &
+     &    sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%ar_1d_rlm,        &
      &    leg%g_sph_rlm, omega_sph%ws_rlm, coef_cor,                    &
      &    gt_cor%jgi_rlm, gt_cor%jei_rlm, gt_cor%sw_rlm, gt_cor%tw_rlm, &
      &    NB, n_WR, irev_sr_rlm, WR, d_rotp_cor_rlm, d_rott_cor_rlm)
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/leg_b_trans_sym_matmul_big.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/leg_b_trans_sym_matmul_big.f90
@@ -117,7 +117,7 @@ subroutine leg_bwd_trans_sym_matmul_big                           &
           WK_l_bsm%time_omp(ip,0) = MPI_WTIME()
           call set_sp_rlm_vec_sym_matmul_big                            &
      &       (sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%istep_rlm,    &
-     &        sph_rlm%idx_gl_1d_rlm_j, sph_rlm%a_r_1d_rlm_r, g_sph_rlm, &
+     &        sph_rlm%idx_gl_1d_rlm_j, sph_rlm%ar_1d_rlm, g_sph_rlm,    &
      &        kst(ip), nkr(ip), jst(ip), n_jk_e(ip), n_jk_o(ip),        &
      &        ncomp, nvector, comm_rlm%irev_sr, n_WR, WR,               &
      &        WK_l_bsm%pol_e(1,ip), WK_l_bsm%tor_e(1,ip),               &
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/leg_bwd_trans_on_the_fly.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/leg_bwd_trans_on_the_fly.f90
@@ -108,13 +108,12 @@ subroutine legendre_b_trans_on_the_fly                            &
 !
         if(iflag_SDT_time) call start_elapsed_time(ist_elapsed_SDT+9)
           call set_sp_rlm_sym_mat_rin                                   &
-     &       (sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%istep_rlm,    &
-     &        sph_rlm%idx_gl_1d_rlm_j, sph_rlm%a_r_1d_rlm_r,            &
-     &        leg%g_sph_rlm, &
-     &        jst, WK_l_otf%n_jk_e(mp_rlm),  WK_l_otf%n_jk_o(mp_rlm),   &
-     &        ncomp, nvector, nscalar, comm_rlm%irev_sr, n_WR, WR,      &
-     &        WK_l_otf%Smat(1)%pol_e(1), WK_l_otf%Smat(1)%tor_e(1),     &
-     &        WK_l_otf%Smat(1)%pol_o(1), WK_l_otf%Smat(1)%tor_o(1) )
+     &      (sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%istep_rlm,     &
+     &       sph_rlm%idx_gl_1d_rlm_j, sph_rlm%ar_1d_rlm, leg%g_sph_rlm, &
+     &       jst, WK_l_otf%n_jk_e(mp_rlm),  WK_l_otf%n_jk_o(mp_rlm),    &
+     &       ncomp, nvector, nscalar, comm_rlm%irev_sr, n_WR, WR,       &
+     &       WK_l_otf%Smat(1)%pol_e(1), WK_l_otf%Smat(1)%tor_e(1),      &
+     &       WK_l_otf%Smat(1)%pol_o(1), WK_l_otf%Smat(1)%tor_o(1) )
         if(iflag_SDT_time) call end_elapsed_time(ist_elapsed_SDT+9)
 !
 !
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/leg_bwd_trans_sym_mat_jt.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/leg_bwd_trans_sym_mat_jt.f90
@@ -102,7 +102,7 @@ subroutine legendre_b_trans_sym_mat_jt                            &
         if(iflag_SDT_time) call start_elapsed_time(ist_elapsed_SDT+9)
         call set_sp_rlm_sym_mat_rin                                     &
      &     (sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%istep_rlm,      &
-     &      sph_rlm%idx_gl_1d_rlm_j, sph_rlm%a_r_1d_rlm_r, g_sph_rlm,   &
+     &      sph_rlm%idx_gl_1d_rlm_j, sph_rlm%ar_1d_rlm, g_sph_rlm,      &
      &      jst, WK_l_tsp%n_jk_e(mp_rlm),  WK_l_tsp%n_jk_o(mp_rlm),     &
      &      ncomp, nvector, nscalar, comm_rlm%irev_sr, n_WR, WR,        &
      &      WK_l_tsp%Smat(1)%pol_e(1), WK_l_tsp%Smat(1)%tor_e(1),       &
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/leg_bwd_trans_sym_mat_tj.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/leg_bwd_trans_sym_mat_tj.f90
@@ -102,7 +102,7 @@ subroutine legendre_b_trans_sym_mat_tj                            &
         if(iflag_SDT_time) call start_elapsed_time(ist_elapsed_SDT+9)
         call set_sp_rlm_sym_mat_rout                                    &
      &     (sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%istep_rlm,      &
-     &      sph_rlm%idx_gl_1d_rlm_j, sph_rlm%a_r_1d_rlm_r, g_sph_rlm,   &
+     &      sph_rlm%idx_gl_1d_rlm_j, sph_rlm%ar_1d_rlm, g_sph_rlm,      &
      &      jst, WK_l_tsp%n_jk_e(mp_rlm),  WK_l_tsp%n_jk_o(mp_rlm),     &
      &      ncomp, nvector, nscalar, comm_rlm%irev_sr, n_WR, WR,        &
      &      WK_l_tsp%Smat(1)%pol_e(1), WK_l_tsp%Smat(1)%tor_e(1),       &
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/legendre_bwd_sym_matmul.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/legendre_bwd_sym_matmul.f90
@@ -128,7 +128,7 @@ subroutine leg_b_trans_vec_sym_matmul                             &
 !          st_elapsed = MPI_WTIME()
           call set_sp_rlm_vector_sym_matmul                             &
      &       (sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%istep_rlm,    &
-     &        sph_rlm%idx_gl_1d_rlm_j, sph_rlm%a_r_1d_rlm_r, g_sph_rlm, &
+     &        sph_rlm%idx_gl_1d_rlm_j, sph_rlm%ar_1d_rlm, g_sph_rlm,    &
      &        kst(ip), nkr(ip), jst(ip), n_jk_e(ip), n_jk_o(ip),        &
      &        ncomp, comm_rlm%irev_sr, n_WR, WR,                        &
      &        WK_l_sml%pol_e(1,ip), WK_l_sml%dpoldt_e(1,ip),            &
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/legendre_bwd_trans_sym_spin.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/legendre_bwd_trans_sym_spin.f90
@@ -100,7 +100,7 @@ subroutine leg_bwd_trans_vector_sym_spin(ncomp, nvector,          &
         kst = sph_rtm%istack_rtm_kr_smp(ip-1) + 1
         ked = sph_rtm%istack_rtm_kr_smp(ip  )
         do k_rlm = kst, ked
-          a1r_1d_rlm_r = sph_rlm%a_r_1d_rlm_r(k_rlm)
+          a1r_1d_rlm_r = sph_rlm%ar_1d_rlm(k_rlm)
           a2r_1d_rlm_r = a1r_1d_rlm_r**2
           do nd = 1, nvector
             do lp = 1, idx_trns%nblock_l_rtm
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/legendre_bwd_trans_symmetry.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/legendre_bwd_trans_symmetry.f90
@@ -112,7 +112,7 @@ subroutine leg_bwd_trans_vector_sym_org(ncomp, nvector,           &
             je = 1 + jst
             jo = 1 + jst + (nj_rlm+1) / 2
             do k_rlm = kst, ked
-              a1r_1d_rlm_r = sph_rlm%a_r_1d_rlm_r(k_rlm)
+              a1r_1d_rlm_r = sph_rlm%ar_1d_rlm(k_rlm)
               a2r_1d_rlm_r = a1r_1d_rlm_r**2
               do ll = 1, nl_rtm
                 lp_rtm =  ll + lst
@@ -185,7 +185,7 @@ subroutine leg_bwd_trans_vector_sym_org(ncomp, nvector,           &
             je = 1 + jst
             jo = 1 + jst + (nj_rlm+1) / 2
             do k_rlm = kst, ked
-              a1r_1d_rlm_r = sph_rlm%a_r_1d_rlm_r(k_rlm)
+              a1r_1d_rlm_r = sph_rlm%ar_1d_rlm(k_rlm)
               a2r_1d_rlm_r = a1r_1d_rlm_r**2
 !
               ip_rtpm = 1 + (lp_rtm-1) * sph_rtm%istep_rtm(2)           &
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/legendre_bwd_trans_testloop.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/legendre_bwd_trans_testloop.f90
@@ -107,13 +107,12 @@ subroutine legendre_b_trans_vector_test                           &
 !
         if(iflag_SDT_time) call start_elapsed_time(ist_elapsed_SDT+9)
           call set_sp_rlm_sym_mat_rin                                   &
-     &       (sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%istep_rlm,    &
-     &        sph_rlm%idx_gl_1d_rlm_j, sph_rlm%a_r_1d_rlm_r,            &
-     &        leg%g_sph_rlm, &
-     &        jst, WK_l_tst%n_jk_e(mp_rlm), WK_l_tst%n_jk_o(mp_rlm),    &
-     &        ncomp, nvector, nscalar, comm_rlm%irev_sr, n_WR, WR,      &
-     &        WK_l_tst%Smat(1)%pol_e(1), WK_l_tst%Smat(1)%tor_e(1),     &
-     &        WK_l_tst%Smat(1)%pol_o(1), WK_l_tst%Smat(1)%tor_o(1) )
+     &      (sph_rlm%nnod_rlm, sph_rlm%nidx_rlm, sph_rlm%istep_rlm,     &
+     &       sph_rlm%idx_gl_1d_rlm_j, sph_rlm%ar_1d_rlm, leg%g_sph_rlm, &
+     &       jst, WK_l_tst%n_jk_e(mp_rlm), WK_l_tst%n_jk_o(mp_rlm),     &
+     &       ncomp, nvector, nscalar, comm_rlm%irev_sr, n_WR, WR,       &
+     &       WK_l_tst%Smat(1)%pol_e(1), WK_l_tst%Smat(1)%tor_e(1),      &
+     &       WK_l_tst%Smat(1)%pol_o(1), WK_l_tst%Smat(1)%tor_o(1) )
         if(iflag_SDT_time) call end_elapsed_time(ist_elapsed_SDT+9)
 !
 !
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/pole_sph_transform.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/pole_sph_transform.f90
@@ -138,7 +138,7 @@ subroutine pole_backward_transforms(ncomp, nvector, nscalar,      &
      &   sph_rlm%nnod_rlm, sph_rlm%nidx_rlm(1), sph_rlm%nidx_rlm(2),    &
      &   sph_rtp%nnod_pole, sph_rlm%istep_rlm,                          &
      &   sph_rtp%nidx_global_rtp, sph_rlm%idx_gl_1d_rlm_r,              &
-     &   sph_rlm%a_r_1d_rlm_r, comm_rlm%irev_sr,                        &
+     &   sph_rlm%ar_1d_rlm, comm_rlm%irev_sr,                           &
      &   leg%g_sph_rlm, leg%P_pole_rtm, leg%dPdt_pole_rtm,              &
      &    n_WR, WR, v_pl_local)
       call schmidt_b_trans_pole_scalar(ncomp, nvector, nscalar,         &
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/set_sp_rlm_sym_mat_tsmp.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/set_sp_rlm_sym_mat_tsmp.f90
@@ -12,11 +12,41 @@
 !!      subroutine set_sp_rlm_sym_mat_rin(nnod_rlm, nidx_rlm,           &
 !!     &          istep_rlm, idx_gl_1d_rlm_j, a_r_1d_rlm_r, g_sph_rlm,  &
 !!     &          jst, n_jk_e, n_jk_o, ncomp_recv, nvector, nscalar,    &
-!!     &          irev_sr_rlm, n_WR, WR,  pol_e, tor_e, pol_o, tor_o)
+!!     &          irev_sr_rlm, n_WR, WR, pol_e, tor_e, pol_o, tor_o)
+!!        integer(kind = kint), intent(in) :: nnod_rlm
+!!        integer(kind = kint), intent(in) :: nidx_rlm(2)
+!!        integer(kind = kint), intent(in) :: istep_rlm(2)
+!!        integer(kind = kint), intent(in)                              &
+!!     &                      :: idx_gl_1d_rlm_j(nidx_rlm(2),3)
+!!        real(kind = kreal), intent(in) :: a_r_1d_rlm_r(nidx_rlm(1))
+!!        real(kind = kreal), intent(in) :: g_sph_rlm(nidx_rlm(2),17)
+!!        integer(kind = kint), intent(in) :: jst, n_jk_e, n_jk_o
+!!        integer(kind = kint), intent(in) :: ncomp_recv
+!!        integer(kind = kint), intent(in) :: n_WR
+!!        integer(kind = kint), intent(in) :: irev_sr_rlm(nnod_rlm)
+!!        real (kind=kreal), intent(in):: WR(n_WR)
+!!        integer(kind = kint), intent(in) :: nvector, nscalar
+!!        real(kind = kreal), intent(inout)                             &
+!!     &           :: pol_e(3*nvector+nscalar,nidx_rlm(1),n_jk_e)
+!!        real(kind = kreal), intent(inout)                             &
+!!     &           :: tor_e(2*nvector,nidx_rlm(1),n_jk_e)
+!!        real(kind = kreal), intent(inout)                             &
+!!     &           :: pol_o(3*nvector+nscalar,nidx_rlm(1),n_jk_o)
+!!        real(kind = kreal), intent(inout)                             &
+!!     &           :: tor_o(2*nvector,nidx_rlm(1),n_jk_o)
 !!      subroutine set_sp_rlm_sym_mat_rout(nnod_rlm, nidx_rlm,          &
 !!     &          istep_rlm, idx_gl_1d_rlm_j, a_r_1d_rlm_r, g_sph_rlm,  &
 !!     &          jst, n_jk_e, n_jk_o, ncomp_recv, nvector, nscalar,    &
-!!     &          irev_sr_rlm, n_WR, WR,  pol_e, tor_e, pol_o, tor_o)
+!!     &          irev_sr_rlm, n_WR, WR, pol_e, tor_e, pol_o, tor_o)
+!!        integer(kind = kint), intent(in) :: nvector, nscalar
+!!        real(kind = kreal), intent(inout)                             &
+!!     &           :: pol_e(n_jk_e,nidx_rlm(1),3*nvector+nscalar)
+!!        real(kind = kreal), intent(inout)                             &
+!!     &           :: tor_e(n_jk_e,nidx_rlm(1),2*nvector)
+!!        real(kind = kreal), intent(inout)                             &
+!!     &           :: pol_o(n_jk_o,nidx_rlm(1),3*nvector+nscalar)
+!!        real(kind = kreal), intent(inout)                             &
+!!     &           :: tor_o(n_jk_o,nidx_rlm(1),2*nvector)
 !!@endverbatim
 !!
       module set_sp_rlm_sym_mat_tsmp
@@ -44,7 +74,8 @@ subroutine set_sp_rlm_sym_mat_rin(nnod_rlm, nidx_rlm,             &
       integer(kind = kint), intent(in) :: nnod_rlm
       integer(kind = kint), intent(in) :: nidx_rlm(2)
       integer(kind = kint), intent(in) :: istep_rlm(2)
-      integer(kind = kint), intent(in) :: idx_gl_1d_rlm_j(nidx_rlm(2),3)
+      integer(kind = kint), intent(in)                                  &
+     &                      :: idx_gl_1d_rlm_j(nidx_rlm(2),3)
       real(kind = kreal), intent(in) :: a_r_1d_rlm_r(nidx_rlm(1))
       real(kind = kreal), intent(in) :: g_sph_rlm(nidx_rlm(2),17)
 !
@@ -161,7 +192,8 @@ subroutine set_sp_rlm_sym_mat_rout(nnod_rlm, nidx_rlm,            &
       integer(kind = kint), intent(in) :: nnod_rlm
       integer(kind = kint), intent(in) :: nidx_rlm(2)
       integer(kind = kint), intent(in) :: istep_rlm(2)
-      integer(kind = kint), intent(in) :: idx_gl_1d_rlm_j(nidx_rlm(2),3)
+      integer(kind = kint), intent(in)                                  &
+     &                      :: idx_gl_1d_rlm_j(nidx_rlm(2),3)
       real(kind = kreal), intent(in) :: a_r_1d_rlm_r(nidx_rlm(1))
       real(kind = kreal), intent(in) :: g_sph_rlm(nidx_rlm(2),17)
 !
diff --git a/src/Fortran_libraries/SERIAL_src/SPH_SPECTR_src/t_spheric_rlm_data.f90 b/src/Fortran_libraries/SERIAL_src/SPH_SPECTR_src/t_spheric_rlm_data.f90
@@ -80,7 +80,7 @@ module t_spheric_rlm_data
 !>        1d radius data for @f$ f(r,l,m) @f$
         real(kind = kreal), allocatable :: radius_1d_rlm_r(:)
 !>        1 / radius_1d_rlm_r
-        real(kind = kreal), allocatable :: a_r_1d_rlm_r(:)
+        real(kind = kreal), allocatable :: ar_1d_rlm(:)
       end type sph_rlm_grid
 !
 ! -----------------------------------------------------------------------
@@ -112,15 +112,15 @@ subroutine alloc_sph_1d_index_rlm(sph_rlm)
       num = sph_rlm%nidx_rlm(1)
       allocate(sph_rlm%idx_gl_1d_rlm_r(num))
       allocate(sph_rlm%radius_1d_rlm_r(num))
-      allocate(sph_rlm%a_r_1d_rlm_r(num))
+      allocate(sph_rlm%ar_1d_rlm(num))
       num = sph_rlm%nidx_rlm(2)
       allocate(sph_rlm%idx_gl_1d_rlm_j(num,3))
 !
       if(sph_rlm%nidx_rlm(2) .gt. 0) sph_rlm%idx_gl_1d_rlm_j = 0
       if(sph_rlm%nidx_rlm(1) .gt. 0) then
         sph_rlm%idx_gl_1d_rlm_r = 0
         sph_rlm%radius_1d_rlm_r = 0.0d0
-        sph_rlm%a_r_1d_rlm_r    = 0.0d0
+        sph_rlm%ar_1d_rlm =       0.0d0
       end if
 !
       end subroutine alloc_sph_1d_index_rlm
@@ -164,8 +164,7 @@ subroutine dealloc_sph_1d_index_rlm(sph_rlm)
       type(sph_rlm_grid), intent(inout) :: sph_rlm
 !
 !
-      deallocate(sph_rlm%radius_1d_rlm_r)
-      deallocate(sph_rlm%a_r_1d_rlm_r   )
+      deallocate(sph_rlm%radius_1d_rlm_r, sph_rlm%ar_1d_rlm)
       deallocate(sph_rlm%idx_gl_1d_rlm_r)
       deallocate(sph_rlm%idx_gl_1d_rlm_j)
 !
@@ -194,7 +193,7 @@ subroutine set_sph_one_over_radius_rlm(sph_rlm)
 !
 !$omp parallel do private(i)
       do i = 1, sph_rlm%nidx_rlm(1)
-        sph_rlm%a_r_1d_rlm_r(i) = one / sph_rlm%radius_1d_rlm_r(i)
+        sph_rlm%ar_1d_rlm(i) = one / sph_rlm%radius_1d_rlm_r(i)
       end do
 !$omp end parallel do
 !
