SUBROUTINE delphi(dphi) $$INCLUDE 'C$INC:BCS.INC' $$INCLUDE 'YBOS$LIBRARY:ERRCOD.INC' REAL dphi ! difference in phi's REAL phi1 ! phi of muon REAL phi2 ! phi of electron INTEGER blocat, status C................................. dphi = -9999. status = blocat(iw,'MUON',1,ind,indmuo) IF(status .NE. yesucc) RETURN status = blocat(iw,'ELEC',1,ind,indele) IF(status .NE. yesucc) RETURN phi1 = rw(indmuo+15) phi2 = rw(indele+5) dphi = ABS(phi1-phi2) RETURN END
This routine would have crashed, because the YBOS return codes were carefully selected to be illegal as REAL numbers. The YBOS functions return INTEGER values--make sure the variable set to receive them is also INTEGER.