SUBROUTINE calcup(vector) $$INCLUDE 'C$INC:BCS.INC' $$INCLUDE 'YBOS$LIBRARY:ERRCOD.INC' REAL vector(6) ! array to be processed INTEGER wbank ! YBOS work bank creator INTEGER wdrop ! YBOS work bank dropper INTEGER windex ! work bank index INTEGER offset ! offset halfway into work bank INTEGER status ! error return status INTEGER i,j ! DO loop indices C................................. status = wbank(iw, windex, 240) IF(status .NE. yesucc) RETURN ! defeated offset = windex + 120 CALL getbnk(windex, offset) ! find JETS and load ! eta in low half and ! phi in high half of ! work bank DO 1 i = 1, 40 DO 2 j = 1, 3 vector(j) = vector(j)+rw(windex+1+i+(j-1)*40) vector(3+j) = vector(3+j) + TAN(rw(offset+1+i+(j-1)*40)) 2 CONTINUE 1 CONTINUE status = wdrop(iw, windex, 1) RETURN END
This routine has a bug. (A hint--the computation has no physics meaning, so don't try to find any.)