C SUBROUTINE SORT C SORTS COL NK+1 AND NK+2 BY NK+1 C ONLY FIRST NT ELEMENTS SUBROUTINE SORT(NT) COMMON /DATA/XL COMMON IAN,KOL1,KOL2,NK,N DIMENSION XL(1) C SOME DOLOOP OTIMIZATIONS KX=NK*N KY=(NK+1)*N DO 100 I=1,NT-1 DO 50 J=I+1,NT X1=XL(KX+J) X2=XL(KX+I) IF (X1.LE.X2)GOTO 50 XL(KX+I)=X1 XL(KX+J)=X2 X1=XL(KY+I) XL(KY+I)=XL(KY+J) XL(KY+J)=X1 50 CONTINUE 100 CONTINUE RETURN END