hord3.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    hord3            ../starter/source/output/th/hord3.F
!||--- called by ------------------------------------------------------
!||    hm_read_thgrne   ../starter/source/output/th/hm_read_thgrne.F
!||--- calls      -----------------------------------------------------
!||====================================================================
        SUBROUTINE hord3(NEL,NSEL,ISK,IASC,L)
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
        END SUBROUTINE hord3