From fd2671c00d2aa6659d51cd676d9d0a01016cd248 Mon Sep 17 00:00:00 2001 From: Alain Hebert Date: Sun, 19 Oct 2025 11:39:22 +0200 Subject: #8: Improve epsilon tests and error detection in module SALT: --- Dragon/src/MUSACG.f90 | 6 ++++++ 1 file changed, 6 insertions(+) (limited to 'Dragon/src/MUSACG.f90') diff --git a/Dragon/src/MUSACG.f90 b/Dragon/src/MUSACG.f90 index 4f68a52..9ceaeb8 100644 --- a/Dragon/src/MUSACG.f90 +++ b/Dragon/src/MUSACG.f90 @@ -57,6 +57,7 @@ SUBROUTINE MUSACG(ITRACK,IFTRK,IPRINT,IMACRO,NBSLIN,RCUTOF,GG,LGINF,NBNODE_MACRO INTEGER ELEM, OK, TYPE REAL(PDB) :: X1,X2,Y1,Y2,DET1,DET2 REAL(PDB) :: DGMESHX(2),DGMESHY(2) + LOGICAL :: LTEST INTEGER, DIMENSION(NSTATE) :: I_STATE,IEDIMG CHARACTER(LEN=72) :: TEXT72 !---- @@ -313,6 +314,11 @@ SUBROUTINE MUSACG(ITRACK,IFTRK,IPRINT,IMACRO,NBSLIN,RCUTOF,GG,LGINF,NBNODE_MACRO !---- ! Unfold macro geometry (many times, if required) !---- + LTEST=.TRUE. + DO IB=1,GG_MAC%NBBCDA + IF(.NOT.LFOLD(IB)) LTEST=.FALSE. + ENDDO + IF(LTEST) CALL XABORT('MUSACG: YOU CANNOT UNFOLD ALL PERIMETERS OF A MACROCELL.') DO IB=1,GG_MAC%NBBCDA IF(LFOLD(IB)) THEN ALLOCATE(IFOLD(2*GG_MAC%NB_ELEM)) -- cgit v1.2.3