c
c User subroutine VFRIC_COEF to define isotropic friction coefficient
c
      subroutine vfric_coef (
c Write only - 
     *   fCoef, fCoefDeriv, 
c Read only - 
     *   nBlock, nProps, nTemp, nFields, 
     *   jFlags, rData, 
     *   surfInt, surfSec, surfMain, 
     *   props, slipData, pressure, 
     *   tempAvg, fieldAvg )
c     
c Write-only variable fCoefDeriv is not used in Explicit
c
c Surface names surfSec and surfMain are not available in 
c general contact (set to blank).

c slipData is sized as 6 for each kBlock
c Non zero quantities:
c    slipData(k,1) = slipRate     
c    slipData(k,4) = accumulated equivalent slip      
c
      include 'vaba_param.inc'
      dimension fCoef(nBlock), 
     *   fCoefDeriv(nBlock,3),
     *   props(nProps),
     *   slipData(nBlock,6),
     *   pressure(nBlock),
     *   tempAvg(nBlock), 
     *   fieldAvg(nBlock,nFields)
c
      parameter( iKStep   = 1,
     *           iKInc    = 2,
     *           nFlags   = 2 )
      parameter( iTimStep = 1,
     *           iTimGlb  = 2,
     *           iDTimCur = 3,
     *           nData    = 3 )
      dimension jFlags(nFlags), rData(nData)
      character*80 surfInt, surfSec, surfMain 
c
      us = props(1)
      uk = props(2)
      dc = props(3)
      du = us - uk
      
      do k=1, nBlock
         fCoef(k) = uk + du*exp( -dc*slipData(k,1) )
      enddo
c
      return
      end