summaryrefslogtreecommitdiff
path: root/Dragon/src/MUSACG.f90
diff options
context:
space:
mode:
authorHEBERT Alain <alain.hebert@polymtl.ca>2025-10-19 07:29:26 -0400
committerHEBERT Alain <alain.hebert@polymtl.ca>2025-10-19 07:29:26 -0400
commit46a6e3604f7ea10bd35d032bc282c6e031da0779 (patch)
tree446377010774826f3d24fa3e8f52905a14c0a7ea /Dragon/src/MUSACG.f90
parentf00eea9ff4dd094374672f5ef73939ed6a4dd9d5 (diff)
parentfd2671c00d2aa6659d51cd676d9d0a01016cd248 (diff)
Merge branch '8-improve-epsilon-tests-and-error-detection-in-module-salt' into 'main'
#8: Improve epsilon tests and error detection in module SALT: See merge request dragon/5.1!14
Diffstat (limited to 'Dragon/src/MUSACG.f90')
-rw-r--r--Dragon/src/MUSACG.f906
1 files changed, 6 insertions, 0 deletions
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))