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/MPOCA2.f | |
| parent | 372969d71ef6f2bdd4b468b0e997852b847a5aa2 (diff) | |
Correct a bug in SPH decorrection
Diffstat (limited to 'Dragon/src/MPOCA2.f')
| -rw-r--r-- | Dragon/src/MPOCA2.f | 9 |
1 files changed, 5 insertions, 4 deletions
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) |