SUBROUTINE DSVRGN(N,RA,RB) IMPLICIT DOUBLE PRECISION(A-H,O-Z) parameter(ldiscr=10000) common/iscrb/iscrb(ldiscr) DIMENSION RA(N),RB(N) if(n.gt.ldiscr) go to 100 DO 10 I=1,N RB(I)=RA(I) 10 CONTINUE CALL DPSORT(RB,N,iscrb,2,IER) RETURN 100 call err(n,'ldiscr','iscrblk') return END