hord3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hord3 (nel, nsel, isk, iasc, l)

Function/Subroutine Documentation

hord3()

subroutine hord3	(	integer, dimension(*)	nel,
		integer	nsel,
		integer, dimension(*)	isk,
		integer, dimension(*)	iasc,
		integer	l )

Definition at line 29 of file hord3.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 NSEL
        INTEGER NEL(*),ISK(*),IASC(*),L
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER J2, I2, I,J, NEL1,ISK1
        INTEGER LI,LJ,K,IC1
        INTEGER :: MODE,INDICE,INDICE2
        INTEGER, DIMENSION(:), ALLOCATABLE :: ITRI,INDX
        INTEGER, DIMENSION(:), ALLOCATABLE :: IASC_SAVE,ISK_SAVE
        INTEGER, DIMENSION(70000) :: WORK
!   ---------------------------------------------
        ALLOCATE( itri(nsel) )
        ALLOCATE( indx(2*nsel) )
        ALLOCATE( iasc_save(l*nsel) )
        ALLOCATE( isk_save(nsel) ) 
 
 
        itri(1:nsel) = nel(1:nsel)
        iasc_save(1:l*nsel) = iasc(1:l*nsel)
        isk_save(1:nsel) = isk(1:nsel)
        DO i=1,nsel
            indx(i) = i
        ENDDO
        mode=0
        CALL my_orders( mode, work, itri, indx, nsel , 1)
        
        DO i=1,nsel
            nel(i) = itri(indx(i))
            isk(i) = isk_save(indx(i))
        ENDDO
        DO i=1,nsel
            indice = (i-1)*l
            indice2 = (indx(i)-1)*l
            DO j=1,l
                iasc(indice+j) = iasc_save(indice2+j)
            ENDDO
        ENDDO
 
        DEALLOCATE( itri )
        DEALLOCATE( indx )
        DEALLOCATE( iasc_save )
        DEALLOCATE( isk_save ) 
 
        RETURN