[cig-commits] [commit] Legendre_transform_w_symmetry: Add maximum numbers for blocked operation (537bef9)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/calypso

On branch  : Legendre_transform_w_symmetry
Link       : https://github.com/geodynamics/calypso/compare/0000000000000000000000000000000000000000...537bef902fcbd7167c78a1b35a7a41d022964d74

>---------------------------------------------------------------

commit 537bef902fcbd7167c78a1b35a7a41d022964d74
Author: Hiroaki Matsui <h_kemono at mac.com>
Date:   Fri Aug 29 09:12:57 2014 -0700

    Add maximum numbers for blocked operation


>---------------------------------------------------------------

537bef902fcbd7167c78a1b35a7a41d022964d74
 src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/init_sph_trans.f90   | 2 ++
 .../PARALLEL_src/SPH_SHELL_src/m_work_4_sph_trans.f90                 | 4 ++++
 2 files changed, 6 insertions(+)

diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/init_sph_trans.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/init_sph_trans.f90
index 689aa41..6bb7ca0 100644
--- a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/init_sph_trans.f90
+++ b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/init_sph_trans.f90
@@ -115,10 +115,12 @@
       mdx_p_rlm_rtm(1:nidx_rlm(2)) = 0
       mdx_n_rlm_rtm(1:nidx_rlm(2)) = 0
       lstack_rlm(0) = 0
+      maxdegree_rlm = 0
       do m = 1, nidx_rtm(3)
         mm = idx_gl_1d_rtm_m(m,2)
         mdx_rlm_rtm(mm) = m
         lstack_rlm(m) = lstack_rlm(m-1) + (l_truncation - abs(mm) + 1)
+        maxdegree_rlm = max(maxdegree_rlm,(l_truncation - abs(mm) + 1))
       end do
 !
       do m = 1, nidx_rtm(3)
diff --git a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/m_work_4_sph_trans.f90 b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/m_work_4_sph_trans.f90
index 7a6fbda..2f5cdf2 100644
--- a/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/m_work_4_sph_trans.f90
+++ b/src/Fortran_libraries/PARALLEL_src/SPH_SHELL_src/m_work_4_sph_trans.f90
@@ -84,6 +84,8 @@
 !
 !>      End address of spherical harmonics order for SMP parallelization
       integer(kind = kint), allocatable :: lstack_rlm(:)
+!>      MAximum number of spherical harmonics order for SMP parallelization
+      integer(kind = kint), allocatable :: maxdegree_rlm
 !
 !>      Data size for Legendre transform to check work area
       integer(kind = kint), private :: iflag_sph_trans = -1
@@ -129,6 +131,7 @@
       allocate(sin_theta_1d_rtp(nidx_rtp(2)))
       allocate(cot_theta_1d_rtp(nidx_rtp(2)))
 !
+      maxdegree_rlm = 0
       lstack_rlm = 0
       mdx_p_rlm_rtm = 0
       mdx_n_rlm_rtm = 0
@@ -161,6 +164,7 @@
       deallocate(sp_rj, vr_rtp)
       deallocate(sp_rlm, vr_rtm)
 !
+      maxdegree_rlm =   0
       iflag_sph_trans = 0
 !
       end subroutine deallocate_work_4_sph_trans