rdtime33.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "scr02_c.inc"
#include "sms_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	rdtime33 (jft, jlt, dt2t, neltst, ityptst, ixr, eint, sti, stir, off, xkm, xkr, xcm, xcr, umas, uiner, fx, fy, fz, xmom, ymom, zmom, rot1, rot2, msrt, dmelrt, nuvar, uvar, jntyp, jsms)

Function/Subroutine Documentation

rdtime33()

subroutine rdtime33	(	integer	jft,
		integer	jlt,
			dt2t,
		integer	neltst,
		integer	ityptst,
		integer, dimension(nixr,*)	ixr,
			eint,
			sti,
			stir,
			off,
			xkm,
			xkr,
			xcm,
			xcr,
			umas,
			uiner,
			fx,
			fy,
			fz,
			xmom,
			ymom,
			zmom,
			rot1,
			rot2,
			msrt,
			dmelrt,
		integer	nuvar,
			uvar,
		integer	jntyp,
		integer, intent(in)	jsms )

Definition at line 30 of file rdtime33.F.

      use element_mod , only : nixr
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr02_c.inc"
#include      "sms_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: JSMS
      INTEGER JFT, JLT, IXR(NIXR,*), NELTST, ITYPTST,
     .   NUVAR,JNTYP
      my_real dt2t, sti(3,*), stir(3,*), off(*), eint(*) ,  
     .   xkm(*),xkr(*),xcm(*),xcr(*),umas(*),uiner(*),
     .   fx(*), fy(*), fz(*), xmom(*), ymom(*),zmom(*),
     .   rot1(*), rot2(*), msrt(*), dmelrt(*),
     .   uvar(nuvar,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real 
     .   dt(jlt-jft+1), dta, dtb, mx2,
     .   ms, in   
C=======================================================================
C     nodal time step
C--------------------------------------------
      IF((idtmins/=2).AND.(jntyp==33)) nodadt = 1
C
      DO i=jft,jlt
        sti(1,i) = xkm(i) 
        stir(1,i) = xkr(i)
        sti(2,i) = sti(1,i) 
        stir(2,i) = stir(1,i)
 
C-------- Harmonic mass --------------
        ms = (uvar(34,i)*uvar(35,i))/max(em20,uvar(34,i)+uvar(35,i))
        in = (uvar(36,i)*uvar(37,i))/max(em20,uvar(36,i)+uvar(37,i))
 
C-------- Node 1--------------
          IF (ms>em15)
     .    sti(1,i) = ((xcm(i)+sqrt(xcm(i)**2+xkm(i)*ms))**2)/ms
          IF (in>em15)
     .    stir(1,i) = ((xcr(i)+sqrt(xcr(i)**2+xkr(i)*in))**2)/in
     
C-------- Node 2--------------
        sti(2,i) = sti(1,i) 
        stir(2,i) = stir(1,i)
           
      ENDDO
C
      IF(idtmins==2.AND.jsms/=0)THEN
C
        dta=dtmins/dtfacs
        dtb=dta*dta
        DO i=jft,jlt
         IF(off(i)<=zero) cycle
          xkm(i)   = max(em15,xkm(i))
          dmelrt(i)=max(dmelrt(i),
     .      xcm(i)*dta+half*xkm(i)*dtb-half*msrt(i))
C
C         MX2 = 2*(Mn+2*DeltaM)
          mx2 =msrt(i)+two*dmelrt(i)
          dt(i)=dtfacs*
     .      mx2 /max(em15,sqrt(xcm(i)*xcm(i)+mx2*xkm(i))+xcm(i))
          dt(i)=dtmins
        ENDDO
C
        DO i=jft,jlt
         IF(off(i)<=zero) cycle
         IF(dt(i)<dt2t) THEN
          dt2t=dt(i)
          neltst =ixr(nixr,i)
          ityptst=6
         ENDIF
        ENDDO
C
      ENDIF
C-----------------------------------------------
      RETURN