stat_n_bcs.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/.
!||====================================================================
!||    stat_n_bcs        ../engine/source/output/sta/stat_n_bcs.F
!||--- called by ------------------------------------------------------
!||    genstat           ../engine/source/output/sta/genstat.F
!||--- calls      -----------------------------------------------------
!||    spmd_istat_gath   ../engine/source/mpi/output/spmd_stat.f
!||    strs_txt50        ../engine/source/output/sta/sta_txt.F
!||--- uses       -----------------------------------------------------
!||    my_alloc_mod      ../common_source/tools/memory/my_alloc.F90
!||====================================================================
      SUBROUTINE stat_n_bcs(ICODE,NUMNOD,ITAB,ITABG,LENG,
     .   NODGLOB,ISKEW,NODTAG)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE my_alloc_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "units_c.inc"
#include      "task_c.inc"
#include      "scr14_c.inc"
#include      "spmd_c.inc"
#include      "com01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ICODE(*),ITAB(*), NUMNOD,ITABG(*),LENG,
     .        nodglob(*),iskew(*),nodtag(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,N, IC, IC1, IC2, IROT, ITRA, ISK
      INTEGER,DIMENSION(:),ALLOCATABLE::ICODEG
      INTEGER,DIMENSION(:),ALLOCATABLE::ISKEWG
      INTEGER,DIMENSION(:),ALLOCATABLE::NODTAGLOB
      CHARACTER TRA*3,ROT*3
      CHARACTER*100 LINE
C-----------------------------------------------
      CALL my_alloc(icodeg,leng)
      CALL my_alloc(iskewg,leng)
      CALL my_alloc(nodtaglob,leng)
C
      IF (nspmd > 1) THEN
        CALL spmd_istat_gath(icode,nodglob,icodeg)
        CALL spmd_istat_gath(iskew,nodglob,iskewg)
        CALL spmd_istat_gath(nodtag,nodglob,nodtaglob)
        IF (ispmd/=0) RETURN
      ENDIF
 
 
      IF (izipstrs == 0) THEN
        WRITE(iugeo,'(A)')'/NBCS/1   '
        WRITE(iugeo,'(A)')' INITIAL BOUNDARY CONDITIONS'
        WRITE(iugeo,'(2A)')'#   TRAROT      SKEW      NODE'
      ELSE
        WRITE(line,'(A)')'/NBCS/1'
               CALL strs_txt50(line,100)
        WRITE(line,'(A)')' INITIAL BOUNDARY CONDITIONS'
               CALL strs_txt50(line,100)
        WRITE(line,'(A)')
     .      '#   TRAROT      SKEW      NODE'
               CALL strs_txt50(line,100)
      ENDIF
 
 
      IF (nspmd == 1) THEN
        DO n=1,numnod
           IF(nodtag(n) /= 0) THEN
             ic=icode(n)       
             ic1=ic/512
             itra = 1
             irot = 1
             tra = '000'
             rot = '000'
             IF(ic1 == 7) THEN
                tra = '111'
             ELSEIF(ic1 == 6) THEN
                tra = '110'
             ELSEIF(ic1 == 5) THEN
                tra = '101'
             ELSEIF(ic1 == 4) THEN
                tra = '100'
             ELSEIF(ic1 == 3) THEN
                tra = '011'
             ELSEIF(ic1 == 2) THEN
                tra = '010'
             ELSEIF(ic1 == 1) THEN
                tra = '001'
             ELSE
                itra = 0    
             ENDIF
             IF(iroddl>0)THEN 
                ic2=(ic-512*ic1)/64
                IF(ic2 == 7) THEN
                   rot = '111'
                ELSEIF(ic2 == 6) THEN
                   rot = '110'
                ELSEIF(ic2 == 5) THEN
                   rot = '101'
                ELSEIF(ic2 == 4) THEN
                   rot = '100'
                ELSEIF(ic2 == 3) THEN
                   rot = '011'
                ELSEIF(ic2 == 2) THEN
                   rot = '010'
                ELSEIF(ic2 == 1) THEN
                   rot = '001'
                ELSE
                   irot = 0        
                ENDIF
             ELSE
                irot = 0        
             ENDIF
C
             IF(itra/=0.AND.irot/=0)THEN
                isk = iskew(n) -1
                IF (izipstrs == 0) THEN
                    WRITE(iugeo,'(3X,A3,1X,A3,I10,I10)')
     .                tra,rot,isk,itab(n)
                ELSE
                    WRITE(line,'(3X,A3,1X,A3,I10,I10)')
     .                tra,rot,isk,itab(n)
                     CALL strs_txt50(line,100)
                ENDIF
             ENDIF
           ENDIF
        END DO
      ELSE
        DO n=1,numnodg
           IF(nodtaglob(n) /= 0) THEN
             ic=icodeg(n)       
             ic1=ic/512
             itra = 1
             irot = 1
             tra = '000'
             rot = '000'
             IF(ic1 == 7) THEN
                tra = '111'
             ELSEIF(ic1 == 6) THEN
                tra = '110'
             ELSEIF(ic1 == 5) THEN
                tra = '101'
             ELSEIF(ic1 == 4) THEN
                tra = '100'
             ELSEIF(ic1 == 3) THEN
                tra = '011'
             ELSEIF(ic1 == 2) THEN
                tra = '010'
             ELSEIF(ic1 == 1) THEN
                tra = '001'
             ELSE
                itra = 0    
             ENDIF
             IF(iroddl>0)THEN 
                ic2=(ic-512*ic1)/64
                IF(ic2 == 7) THEN
                   rot = '111'
                ELSEIF(ic2 == 6) THEN
                   rot = '110'
                ELSEIF(ic2 == 5) THEN
                   rot = '101'
                ELSEIF(ic2 == 4) THEN
                   rot = '100'
                ELSEIF(ic2 == 3) THEN
                   rot = '011'
                ELSEIF(ic2 == 2) THEN
                   rot = '010'
                ELSEIF(ic2 == 1) THEN
                   rot = '001'
                ELSE
                   irot = 0        
                ENDIF
             ELSE
                irot = 0        
             ENDIF
C
             IF(itra/=0.AND.irot/=0)THEN
               isk = iskewg(n) -1
               IF (izipstrs == 0) THEN
                 WRITE(iugeo,'(3X,A3,1X,A3,I10,I10)')
     .           tra,rot,isk,itabg(n)
               ELSE
                 WRITE(line,'(3X,A3,1X,A3,I10,I10)')
     .           tra,rot,isk,itabg(n)
                 CALL strs_txt50(line,100)
               END IF
           END IF
          ENDIF
        END DO
      ENDIF
 
C
      DEALLOCATE(icodeg)
      DEALLOCATE(iskewg)
      DEALLOCATE(nodtaglob)
 
      RETURN
      END