modsti.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	modsti (nodft, nodlt, stifn, viscn, ms)

Function/Subroutine Documentation

modsti()

subroutine modsti	(	integer	nodft,
		integer	nodlt,
			stifn,
			viscn,
			ms )

Definition at line 28 of file modsti.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-----------------------------------------------
C
      INTEGER NODFT,NODLT
      my_real 
     .  stifn(*),viscn(*),ms(*)
C----------------------------------------
C   L o c a l   V a r i a b l e s
C----------------------------------------
      INTEGER I,N
      my_real 
     .  c2,c21,aux
C----------------------------------------
       DO i=nodft,nodlt
        IF(ms(i)>zero)THEN
         c2      = four*viscn(i)*viscn(i)
         c21     = eight*ms(i)*stifn(i)
         aux     = sqrt(c2+c21)+2.*viscn(i)
         stifn(i)= one_over_8*aux*aux/max(ms(i),em30)
C        STIFN(I)= STIFN(I)*AUX*AUX/MAX(C21,EM30)
         stifn(i)= max(stifn(i),two*viscn(i)*viscn(i)/max(ms(i),em30))
         viscn(i)= zero
        ENDIF
       ENDDO
C
      RETURN