summaryrefslogtreecommitdiff
path: root/Dragon/src/MUSACG.f90
diff options
context:
space:
mode:
authorAlain Hebert <alainhebert@Alains-MacBook-Pro.local>2025-10-19 11:39:22 +0200
committerAlain Hebert <alainhebert@Alains-MacBook-Pro.local>2025-10-19 11:39:22 +0200
commitfd2671c00d2aa6659d51cd676d9d0a01016cd248 (patch)
tree446377010774826f3d24fa3e8f52905a14c0a7ea /Dragon/src/MUSACG.f90
parentf00eea9ff4dd094374672f5ef73939ed6a4dd9d5 (diff)
#8: Improve epsilon tests and error detection in module SALT:
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))