accelepinch.F File Reference

#include "implicit_f.inc"
#include "vectorize.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	accelepinch (apinch, ms, mspinch, stifpinch, nodft, nodlt, dtnod, dtfac)

Function/Subroutine Documentation

accelepinch()

subroutine accelepinch	(		apinch,
			ms,
			mspinch,
			stifpinch,
		integer	nodft,
		integer	nodlt,
			dtnod,
			dtfac )

Definition at line 30 of file accelepinch.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE pinchtype_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NODFT, NODLT  
      my_real apinch(3,*),ms(*),mspinch(*),stifpinch(*)
      my_real dtnod,dtfac
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N
      my_real rtmp,mtemppinch,factor
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------    
#include      "vectorize.inc"
      DO n=nodft,nodlt
        IF(ms(n)>zero) THEN
          factor = (dtnod**2) * half / (dtfac**2)
          mspinch(n) = max(mspinch(n), stifpinch(n)*factor) 
          rtmp = one / mspinch(n)
C       
          apinch(1,n) = apinch(1,n) * rtmp
          apinch(2,n) = apinch(2,n) * rtmp
          apinch(3,n) = apinch(3,n) * rtmp
        ELSE
          apinch(1,n) = zero
          apinch(2,n) = zero
          apinch(3,n) = zero
        ENDIF
        stifpinch(n) = zero          
      ENDDO
      RETURN