青衣巷 发表于 2021-11-20 10:33
include 'aba_param.inc'
         real :: a, b, c

能运行了，谢谢

       include 'aba_param.inc'
         real :: a, b, c

交换一下位置

[Fortran] 查看源码 复制源码

      subroutine fric_coef (
C Write only - 
     *   fCoef, fCoefDeriv, 
C Read only - 
     *   nBlock, nProps, nTemp, nFields, 
     *   jFlags, rData, 
     *   surfInt, surfSlv, surfMst, 
     *   props, slipRate, pressure, 
     *   tempAvg, fieldAvg 
     *   real a, b, c )
C
      include 'aba_param.inc'
C
      dimension fCoef(nBlock), 
     *   fCoefDeriv(nBlock,3),
     *   props(nProps),
     *   slipRate(nBlock),
     *   pressure(nBlock),
     *   tempAvg(nBlock), 
     *   fieldAvg(nBlock,nFields)
C
      parameter( iKStep   = 1,
     *           iKInc    = 2,
     *           nFlags   = 2
     *           a        = 0.1,
     *           b        = -1.2,
     *           c        = 0.3) 
C
      parameter( iTimStep = 1,
     *           iTimGlb  = 2,
     *           iDTimCur = 3,
     *           nData    = 3,) 
C
      dimension jFlags(nFlags), rData(nData)
C
      character*80 surfInt, surfSlv, surfMst 
C
         xMuk = props(1)
         xMus = props(2)
         beta = props(3)
      DO  k = 1, nBlock
         fCoef(k)=(a*exp(b*pressure(k))+c)
                   *[xMuk+(xMus-xMuk)*exp(-beta*slipRate(k))]
      END DO

      return
      end