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 28 of file rdtime33.F.

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 dx(3),theta(3),e0(mvsiz),th1,th2,th3,
     .   dt(jlt-jft+1), dtc(1), dtinv,  a, mass2, in2, 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-------- Masse harmonique --------------
        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