subroutine gapcon(ak,d,flowm,temp,predef,time,ciname, 1 secname,mainname,coords,noel,node,npred,kstep,kinc) c include 'aba_param.inc' c character*80 ciname,secname,mainname 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