diff options
| author | Alain Hebert <alainhebert@Alains-MacBook-Pro.local> | 2025-11-13 18:32:05 +0100 |
|---|---|---|
| committer | Alain Hebert <alainhebert@Alains-MacBook-Pro.local> | 2025-11-13 18:32:05 +0100 |
| commit | 733dff759433a8cace9805d6127b4bfa3f71e221 (patch) | |
| tree | aad98db8f19dbd1d493c069a7d03a5faafb2229c /Dragon/src/APXCA2.f | |
| parent | 372969d71ef6f2bdd4b468b0e997852b847a5aa2 (diff) | |
Correct a bug in SPH decorrection
Diffstat (limited to 'Dragon/src/APXCA2.f')
| -rw-r--r-- | Dragon/src/APXCA2.f | 7 |
1 files changed, 5 insertions, 2 deletions
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) |