diff --git a/src/msspec/spec/fortran/phd_se_noso_nosp_nosym/findpaths1.f b/src/msspec/spec/fortran/phd_se_noso_nosp_nosym/findpaths1.f index 47c0139..3083730 100644 --- a/src/msspec/spec/fortran/phd_se_noso_nosp_nosym/findpaths1.f +++ b/src/msspec/spec/fortran/phd_se_noso_nosp_nosym/findpaths1.f @@ -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