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

Fix call to findpaths for order > 10

This commit is contained in:
Sylvain Tricot 2025-03-18 12:42:59 +01:00
parent 183e704f34
commit c7f1223341
1 changed files with 4 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/msspec/spec/fortran/phd_se_noso_nosp_nosym/findpaths1.f
View File

@ -327,6 +327,10 @@ C
IF((ND.LT.NDIF).OR.(IPW.EQ.0)) THEN
CALL ARCSIN(COMPL1,CTROIS1,PHIJK)
ENDIF
C PRINT *,"KTYP=",KTYP,"(X,Y,Z)=",SYM_AT(1,KTYP),
C &SYM_AT(2,KTYP),SYM_AT(3,KTYP)
C PRINT *,"JTYP=",JTYP,"(X,Y,Z)=",SYM_AT(1,JTYP),
C &SYM_AT(2,JTYP),SYM_AT(3,JTYP)
CALL EULER(THJK,PHIJK,THIJ,PHIIJ,AIJK,BIJK,CIJK,IEULER)
IF((I_CP.EQ.1).AND.(ND.EQ.NDIF)) I_ABS=2
CALL MATDIF(NO,ND-1,LF2,JTYP,KTYP,JE,I_ABS,ISPEED,ISPHER,A