epsi
/
msspec_python3
3
0
Fork 0
Code Issues 6 Pull Requests Releases Wiki Activity

Fix lmax bug 

There was a limit in the lmax value due to too large numbers
in spherical harmonics. This was due to large arguments in
exponential function in EXPFAC_MOD and EXPFAC2_MOD. We changed
the kind from REAL to REAL*8
This commit is contained in:
Sylvain Tricot 2025-04-03 21:56:15 +02:00
parent 9f7c9f97f3
commit 9167d613d6
2 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/msspec/spec/fortran/common_sub/read_data.f
View File

@ -83,8 +83,9 @@ C
REAL TEXTE1(10),TEXTE2(10),TEXTE3(10)
REAL TEXTE4(10),TEXTE5(10),TEXTE6(10)
REAL TEXTE6B(10),TEXTE7(10)
REAL THFWD(NATP_M),THBWD(NATP_M),GLG(0:N_GAUNT),NJ(0:N_GAUNT)
REAL THFWD(NATP_M),THBWD(NATP_M)
REAL ALPHAR,BETAR,RACC
REAL*8 GLG(0:N_GAUNT),NJ(0:N_GAUNT)
C
C
C

4
src/msspec/spec/fortran/memalloc/modules.f
View File

@ -781,7 +781,7 @@ C=======================================================================
C=======================================================================
MODULE EXPFAC_MOD
IMPLICIT NONE
REAL, ALLOCATABLE, DIMENSION(:,:) :: EXPF
REAL*8, ALLOCATABLE, DIMENSION(:,:) :: EXPF
CONTAINS
SUBROUTINE ALLOC_EXPFAC()
USE DIM_MOD
@ -837,7 +837,7 @@ C=======================================================================
C=======================================================================
MODULE EXPFAC2_MOD
IMPLICIT NONE
REAL, ALLOCATABLE, DIMENSION(:,:) :: EXPF2
REAL*8, ALLOCATABLE, DIMENSION(:,:) :: EXPF2
CONTAINS
SUBROUTINE ALLOC_EXPFAC2()
USE DIM_MOD