SUBROUTINE getcur(curve) $$INCLUDE 'C$INC:BCS.INC' $$INCLUDE 'YBOS$LIBRARY:ERRCOD.INC' REAL curve ! curvature of muon INTEGER blocat ! YBOS named bank finder INTEGER status ! error return status C................................. curve = 0.0 status = blocat(iw,'MUON',1,ind,indmuo) IF(status .NE. yesucc) RETURN ! or do something else sensible IF(rw(indmuo+7) .NE. 0) curve = .45/rw(indmuo+7) RETURN END
The routine might have worked, most of the time, but when the bank doesn't exist, the index to the data is set to 0. The routine then reads something from the YBOS header instead of actual physics data.