From 733dff759433a8cace9805d6127b4bfa3f71e221 Mon Sep 17 00:00:00 2001 From: Alain Hebert Date: Thu, 13 Nov 2025 18:32:05 +0100 Subject: Correct a bug in SPH decorrection --- Dragon/src/APXCA2.f | 7 +++++-- Dragon/src/MPOCA2.f | 9 +++++---- 2 files changed, 10 insertions(+), 6 deletions(-) (limited to 'Dragon') diff --git a/Dragon/src/APXCA2.f b/Dragon/src/APXCA2.f index e541cc8..97aa69e 100644 --- a/Dragon/src/APXCA2.f +++ b/Dragon/src/APXCA2.f @@ -129,7 +129,7 @@ CALL LCMGET(IPEDIT,'VOLUME',VOLMIX) JPEDIT=LCMGID(IPEDIT,'GROUP') LSPH=.FALSE. - ALLOCATE(SPH(NMIL,NG)) + ALLOCATE(SPH(NMIL+NALBP,NG)) DO 80 IGR=1,NG KPEDIT=LCMGIL(JPEDIT,IGR) CALL LCMLEN(KPEDIT,'NSPH',ILONG,ITYLCM) @@ -139,8 +139,11 @@ DO 70 IMIL=1,NMIL SPH(IMIL,IGR)=1.0/WORK1(IMIL) 70 CONTINUE + DO 75 IALB=1,NALBP + SPH(NMIL+IALB,IGR)=1.0 + 75 CONTINUE ELSE - SPH(:NMIL,IGR)=1.0 + SPH(:NMIL+NALBP,IGR)=1.0 ENDIF 80 CONTINUE CALL LCMSIX(IPEDIT,' ',2) diff --git a/Dragon/src/MPOCA2.f b/Dragon/src/MPOCA2.f index 3fd0927..70a24a8 100644 --- a/Dragon/src/MPOCA2.f +++ b/Dragon/src/MPOCA2.f @@ -246,7 +246,7 @@ CALL LCMSIX(IPEDIT,'MACROLIB',1) JPEDIT=LCMGID(IPEDIT,'GROUP') LSPH=.FALSE. - ALLOCATE(SPH(NMIL,NG)) + ALLOCATE(SPH(NMIL+NALBP,NG)) DO 120 IGR=1,NG KPEDIT=LCMGIL(JPEDIT,IGR) CALL LCMLEN(KPEDIT,'NSPH',ILONG,ITYLCM) @@ -256,10 +256,11 @@ DO 100 IMIL=1,NMIL SPH(IMIL,IGR)=1.0/WORK1(IMIL) 100 CONTINUE + DO 110 IALB=1,NALBP + SPH(NMIL+IALB,IGR)=1.0 + 110 CONTINUE ELSE - DO 110 IMIL=1,NMIL - SPH(IMIL,IGR)=1.0 - 110 CONTINUE + SPH(:NMIL+NALBP,IGR)=1.0 ENDIF 120 CONTINUE CALL LCMSIX(IPEDIT,' ',2) -- cgit v1.2.3