stat_n_bcs.F File Reference

#include "implicit_f.inc"
#include "units_c.inc"
#include "task_c.inc"
#include "scr14_c.inc"
#include "spmd_c.inc"
#include "com01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	stat_n_bcs (icode, numnod, itab, itabg, leng, nodglob, iskew, nodtag)

Function/Subroutine Documentation

stat_n_bcs()

subroutine stat_n_bcs	(	integer, dimension(*)	icode,
		integer	numnod,
		integer, dimension(*)	itab,
		integer, dimension(*)	itabg,
		integer	leng,
		integer, dimension(*)	nodglob,
		integer, dimension(*)	iskew,
		integer, dimension(*)	nodtag )

Definition at line 33 of file stat_n_bcs.F.

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
.AND.             IF(ITRA/=0IROT/=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
.AND.             IF(ITRA/=0IROT/=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