dum(1) = RDATA(jlast)
Several months ago - when I started the COVID stuff I had this line. I made a mistake and created jlast as a real, and I did not assign it an initial value, but the thing gave an answer - which is actually correct
Next loop it got assigned the correct update and looped through RDATA ok.
I added one line in another routine
lambda = Constants%lambda * ((0.1*randdata(nrhs)) + 0.99)
and it tells me the error exists.
Sometimes, there is weirdness in the land.
I realize there is no answer.