sorti25.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sorti25 (n, isky, fskyi, nfskyi)

Function/Subroutine Documentation

sorti25()

subroutine sorti25	(	integer	n,
		integer, dimension(*)	isky,
			fskyi,
		integer	nfskyi )

Definition at line 28 of file sorti25.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER N, NFSKYI, ISKY(*)
      my_real fskyi(2*nfskyi,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, L, IR, IIA
      my_real rra1, rra2, rra3, rra4, rra5, rra6, rra7, rra8, rra9, rra10
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
C   test sur n=0
      IF (n==0) RETURN
      l=n/2+1
      ir=n
      IF(nfskyi==4)THEN
10     CONTINUE
        IF(l>1)THEN
          l=l-1
          iia =isky(l)
          rra1=fskyi(1,l)
          rra2=fskyi(2,l)
          rra3=fskyi(3,l)
          rra4=fskyi(4,l)
          rra5=fskyi(5,l)
          rra6=fskyi(6,l)
          rra7=fskyi(7,l)
          rra8=fskyi(8,l)
        ELSE
          iia =isky(ir)
          rra1=fskyi(1,ir)
          rra2=fskyi(2,ir)
          rra3=fskyi(3,ir)
          rra4=fskyi(4,ir)
          rra5=fskyi(5,ir)
          rra6=fskyi(6,ir)
          rra7=fskyi(7,ir)
          rra8=fskyi(8,ir)
          isky(ir)=isky(1)
          fskyi(1,ir)=fskyi(1,1)
          fskyi(2,ir)=fskyi(2,1)
          fskyi(3,ir)=fskyi(3,1)
          fskyi(4,ir)=fskyi(4,1)
          fskyi(5,ir)=fskyi(5,1)
          fskyi(6,ir)=fskyi(6,1)
          fskyi(7,ir)=fskyi(7,1)
          fskyi(8,ir)=fskyi(8,1)
          ir=ir-1
C         IF(IR==1)THEN
          IF(ir<=1)THEN
            isky(1)=iia
            fskyi(1,1)=rra1
            fskyi(2,1)=rra2
            fskyi(3,1)=rra3
            fskyi(4,1)=rra4
            fskyi(5,1)=rra5
            fskyi(6,1)=rra6
            fskyi(7,1)=rra7
            fskyi(8,1)=rra8
            RETURN
          ENDIF
        ENDIF
        i=l
        j=l+l
20      IF(j<=ir)THEN
          IF(j<ir)THEN
            IF(isky(j)<isky(j+1))j=j+1
          ENDIF
          IF(iia<isky(j))THEN
            isky(i)=isky(j)
            fskyi(1,i)=fskyi(1,j)
            fskyi(2,i)=fskyi(2,j)
            fskyi(3,i)=fskyi(3,j)
            fskyi(4,i)=fskyi(4,j)
            fskyi(5,i)=fskyi(5,j)
            fskyi(6,i)=fskyi(6,j)
            fskyi(7,i)=fskyi(7,j)
            fskyi(8,i)=fskyi(8,j)
            i=j
            j=j+j
          ELSE
            j=ir+1
          ENDIF
        GO TO 20
        ENDIF
        isky(i)=iia
        fskyi(1,i)=rra1
        fskyi(2,i)=rra2
        fskyi(3,i)=rra3
        fskyi(4,i)=rra4
        fskyi(5,i)=rra5
        fskyi(6,i)=rra6
        fskyi(7,i)=rra7
        fskyi(8,i)=rra8
       GO TO 10
      ELSE
100    CONTINUE
        IF(l>1)THEN
          l=l-1
          iia =isky(l)
          rra1=fskyi(1,l)
          rra2=fskyi(2,l)
          rra3=fskyi(3,l)
          rra4=fskyi(4,l)
          rra5=fskyi(5,l)
          rra6=fskyi(6,l)
          rra7=fskyi(7,l)
          rra8=fskyi(8,l)
          rra9=fskyi(9,l)
          rra10=fskyi(10,l)
        ELSE
          iia =isky(ir)
          rra1=fskyi(1,ir)
          rra2=fskyi(2,ir)
          rra3=fskyi(3,ir)
          rra4=fskyi(4,ir)
          rra5=fskyi(5,ir)
          rra6=fskyi(6,ir)
          rra7=fskyi(7,ir)
          rra8=fskyi(8,ir)
          rra9=fskyi(9,ir)
          rra10=fskyi(10,ir)
          isky(ir)=isky(1)
          fskyi(1,ir)=fskyi(1,1)
          fskyi(2,ir)=fskyi(2,1)
          fskyi(3,ir)=fskyi(3,1)
          fskyi(4,ir)=fskyi(4,1)
          fskyi(5,ir)=fskyi(5,1)
          fskyi(6,ir)=fskyi(6,1)
          fskyi(7,ir)=fskyi(7,1)
          fskyi(8,ir)=fskyi(8,1)
          fskyi(9,ir)=fskyi(9,1)
          fskyi(10,ir)=fskyi(10,1)
          ir=ir-1
C         IF(IR==1)THEN
          IF(ir<=1)THEN
            isky(1)=iia
            fskyi(1,1)=rra1
            fskyi(2,1)=rra2
            fskyi(3,1)=rra3
            fskyi(4,1)=rra4
            fskyi(5,1)=rra5
            fskyi(6,1)=rra6
            fskyi(7,1)=rra7
            fskyi(8,1)=rra8
            fskyi(9,1)=rra9
            fskyi(10,1)=rra10
            RETURN
          ENDIF
        ENDIF
        i=l
        j=l+l
200     IF(j<=ir)THEN
          IF(j<ir)THEN
            IF(isky(j)<isky(j+1))j=j+1
          ENDIF
          IF(iia<isky(j))THEN
            isky(i)=isky(j)
            fskyi(1,i)=fskyi(1,j)
            fskyi(2,i)=fskyi(2,j)
            fskyi(3,i)=fskyi(3,j)
            fskyi(4,i)=fskyi(4,j)
            fskyi(5,i)=fskyi(5,j)
            fskyi(6,i)=fskyi(6,j)
            fskyi(7,i)=fskyi(7,j)
            fskyi(8,i)=fskyi(8,j)
            fskyi(9,i)=fskyi(9,j)
            fskyi(10,i)=fskyi(10,j)
            i=j
            j=j+j
          ELSE
            j=ir+1
          ENDIF
        GO TO 200
        ENDIF
        isky(i)=iia
        fskyi(1,i)=rra1
        fskyi(2,i)=rra2
        fskyi(3,i)=rra3
        fskyi(4,i)=rra4
        fskyi(5,i)=rra5
        fskyi(6,i)=rra6
        fskyi(7,i)=rra7
        fskyi(8,i)=rra8
        fskyi(9,i)=rra9
        fskyi(10,i)=rra10
       GO TO 100
      ENDIF
C
      RETURN