thbcs_imp.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	thbcs_imp (nodft, nodlast, a, ar, fthreac, nodreac, cptreac, fthdtm, dt3)

Function/Subroutine Documentation

thbcs_imp()

subroutine thbcs_imp	(	integer	nodft,
		integer	nodlast,
			a,
			ar,
			fthreac,
		integer, dimension(*)	nodreac,
		integer	cptreac,
			fthdtm,
			dt3 )

Definition at line 30 of file thbcs_imp.F.

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      "com01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NODREAC(*),CPTREAC,NODFT,NODLAST
      my_real
     .    fthreac(6,*),a(3,*),ar(3,*),dt3,fthdtm(6,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, N, K, NINDX,
     .        INDX(1024),L,J,FLAG_MOM
C-----------------------------------------------
C    Time History
C-----------------------------------------------
        IF (cptreac >= 1) THEN
          flag_mom = 0
          indx = 0
C
          DO 460 i = nodft,nodlast,1024
            nindx = 0
            DO 450 n=i,min(nodlast,i+1023)
              IF (nodreac(n)>0)THEN
                nindx = nindx + 1
                indx(nindx) = n
              ENDIF
  450       CONTINUE
            CALL bcs1th_imp(nindx,indx,a,fthreac,nodreac,flag_mom,
     &                      fthdtm,dt3)
  460     CONTINUE
C
          IF(iroddl/=0) THEN
            flag_mom = 3
            DO 480 i = nodft,nodlast,1024
              nindx = 0
              DO 470 n=i,min(nodlast,i+1023)
                IF (nodreac(n)>0)THEN
                  nindx = nindx + 1
                  indx(nindx) = n
                ENDIF
  470       CONTINUE
            CALL bcs1th_imp(nindx,indx,ar,fthreac,nodreac,flag_mom,
     &                      fthdtm,dt3)
  480       CONTINUE
          ENDIF
        ENDIF
C-----------------------------------------------
      RETURN