subroutine gapcon(ak,d,flowm,temp,predef,time,ciname,
     1  slname,msname,coords,noel,node,npred,kstep,kinc)
c
      include 'aba_param.inc'
c
      character*80 ciname,slname,msname
      dimension flowm(2),temp(2),predef(2,*),
     1          ak(5),d(2),coords(3),time(2)
c
c       define gap conductivity dependent on temperature.
        tavg=0.5*(temp(1)+temp(2))
        ak(1)=15.0+25.0*tavg
        ak(4)=12.5
        ak(5)=12.5
c
      return
      end