CTITLESAPK02 - SUBROUTINE USED BY SDPICK 00010000 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** CA 00020000 CA DESIGNER P. C. LUH 00030000 CA AUTHOR P. C. LUH 00040000 CA LANGUAGE FORTRAN 00050000 CA SYSTEM IBM ONLY 00060000 CA DATE MM/DD/YY 00070000 C REVISED 03/06/86 LBL OLD NAME = DMED. 00080000 C 00090000 C 00100000 C SUBROUTINE DMED (N,B,D1,D2) 00110000 SUBROUTINE SAPK02 (N,B,D1,D2) 00120000 C 00130000 DIMENSION B(2) 00140002 C 00150000 D1=B(1) 00160000 IF(N-2) 10,15,20 00170000 10 D2=B(1) 00180000 GO TO 70 00190000 C 00200000 15 D2=B(2) 00210000 GO TO 70 00220000 C 00230000 20 M2=N/2 00240000 NP=(N+1)/2 00250000 IS=MOD(N,2)+1 00260000 NNP=N 00270000 C 00280000 DO 40 I=1,M2 00290000 X=B(I) 00300000 YK=B(I) 00310000 KX=I 00320000 KN=I 00330000 K=I+1 00340000 C 00350000 DO 30 L=K,NNP 00360000 IF(B(L).LT.X) GO TO 25 00370000 X=B(L) 00380000 KX=L 00390000 25 IF(YK.LE.B(L)) GO TO 30 00400000 YK=B(L) 00410000 KN=L 00420000 30 CONTINUE 00430000 C 00440000 Q=B(I) 00450000 B(I)=YK 00460000 B(KN)=Q 00470000 IF(KX.EQ.I) KX=KN 00480000 Q=B(NNP) 00490000 B(NNP)=X 00500000 B(KX)=Q 00510000 40 NNP=NNP-1 00520000 C 00530000 D1=B(NP) 00540000 GO TO (50,60),IS 00550000 50 D2=B(NP+1) 00560000 GO TO 70 00570000 60 D2=B(NP) 00580000 C 00590000 70 RETURN 00600000 C 00610000 END 00620000