From dd9e7e30ad03b7db042a93ac37cc73a0897919bf Mon Sep 17 00:00:00 2001 From: Daniel Lopez Aldama Date: Sun, 27 Jul 2025 22:08:09 +0200 Subject: [PATCH] Fixes #369 : MATXS formatted tape for constant spectrum and production --- src/matxsr.f90 | 8 ++++---- tests/03/referenceOutput | Bin 69504 -> 69504 bytes 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/src/matxsr.f90 b/src/matxsr.f90 index bbb941f8..9a71d10d 100644 --- a/src/matxsr.f90 +++ b/src/matxsr.f90 @@ -2214,12 +2214,12 @@ subroutine shufl(nscr,iref,jg1lo,jg1hi,nout3,irec3,irite) go to 145 135 do ik=1,ng2 jg2=noutg-ig2lo-ik+2 - a(icdat-1+jg2)=b(irinp+lz+ik-1) + a(icdat-1+jg2)=b(irinp+lz+nz*(ik-1)+iz-1) enddo go to 240 140 jconst=ig jg1=ning+1-ig - a(icdat-1+noutg+jg1)=b(irinp+lz+1) + a(icdat-1+noutg+jg1)=b(irinp+lz+nz+iz-1) go to 240 !--store matrix sub-blocks @@ -2359,10 +2359,10 @@ subroutine lst1io(nin,nout,nscr,a,nb,nw,namax) nwd=6-nwh if (nwd.le.0) go to 120 if (nx.ne.0) write(hfmt,& - &'(''('',i2,''a4,'',i2,''x,1p'',i1,''e11.0)'')')& + &'(''('',i2,''a4,'',i2,''x,'',i1,''e11.0)'')')& ntw,nx,nwd if (nx.eq.0) write(hfmt,& - '(''('',i2,''a4,'',''1p'',i1,''e11.0)'')')& + &'(''('',i2,''a4,'',i1,''e11.0)'')')& ntw,nwd read(nin,hfmt) (tz(i),i=1,ntw),(a(i+nw1+ntw),i=1,nwd) do i=1,ntw diff --git a/tests/03/referenceOutput b/tests/03/referenceOutput index 755cc4416127cb3b73bbb5db2fbf80a6ee2d1a6b..01435c0963044ba9fe05d4520a02eb63f46b54fb 100644 GIT binary patch delta 917 zcmZ9KPiWI{6vs(hW$PxEwiEwQ_`xE~sx?VdvxFUN4?S3@n_d)0UDk$K)^^gw!4A?6 z-EM;8i7-Uw(1UpJ@;i${59`6*z3I)12d`d+?QSy+ZQEdqBwx>HvsfA6kSB^Ag5jZJ zgwY8U`3Pt|t~M;QO^W1t&D=1Ulr37zN{+0-cv$Hc(>)L82OgPhs;`-j>Bg??FuU$p z25XqD26Ghh20V8j4>igQ8X@mdd#!5OcD=RY;+tmmp)uee!x3+2Nq1MVL+koY=r&{a zhV58Q*Y(H6@g(hfdbHv3JWzmh>P^HcvZ2a`tTGlM9{mW@*n0Sk zNJ^z4*1yfqhChx$aZH216Z3)ur~7x{&14Cx;T%t$x^e=o{c}lPiEj&D2uH>$uO%;# z%omcdJNJ^$oOyp9-r}{N(t?1&1e6NnaPuC#`>upgC*aThXM8QUPss$#R9bvT`%GdA zD)Bo1rl(xNkBsd>eU$IcO-2L-+dtvc`II0iJYRS*1u-oFKcrQ-b@vzCC?q+y389&U z-|AEDEMkwCwh6B)?ThkdEsA^6U&^$MqDrdnIXXFdMw(T19j-pu!iosQ`}1&esU&3K ou$SdqOCOd4d$Wf&KkR)8WA^!hg72>fWk3BIQ=j_XqA$z;0G?;qqyPW_ delta 568 zcmZpe&(biTMVw0^FDt)N!N|bSOwT|;LBYt-B{kX5Y+|6gpn;*jp|QT9xq_jAm5GIw zk;!C%6G5Cssmb|yMS6OAlNos&CLiEt<>pdIPt46ki0K=#a&jr8l%zpbISCl(nH6(! zDdZ-WR20Ls@W;Vvfu}vB?LH`!SkN=06eW$v`czS&VrN$K(Z)ER(%Bj3@u( z4gj--P6SMz%j3;xxcMbd0xM8RbaR|g8zYD@`JYwx<`fZW5QlqmzsbK{FO(FX{JN+=j{1;)+xaw9UDIF=qk93g&u7CI%)J5D`m|h^3x^g^8ubWZn2y zpdS7c!IRk%yn(E